diff --git a/vignettes/tutorial.Rmd b/vignettes/tutorial.Rmd
index 593b25bda082ebeacb694eb877a2f662d313dac7..cddcec2d09fc505ff8b42c352f152d3c857ba4ed 100644
--- a/vignettes/tutorial.Rmd
+++ b/vignettes/tutorial.Rmd
@@ -44,8 +44,7 @@ pub_metadata <- read.table(file = pub_metadata_fn,  header = TRUE)
 
 To estimate plausible parameter distributions, we will fit a generalized linear  model (GLM) to the publicly available data. This GLM will allow us to estimate plausible coefficients betaG, betaE, and betaGxE for 6000 genes. The `tidy_pub_fit` object contains the estimate obtain by HTRfit for each coefficient and each genes.
 
-```{ r fit_pub_data, eval=FALSE, message=FALSE, warning=FALSE }
-
+```{r fit_pub_data, eval=FALSE, message=FALSE, warning=FALSE, include=TRUE}
 ## -- prepare data &  fit a model 
 data2fit = prepareData2fit(countMatrix = pub_data, 
                            metadata = pub_metadata, 
@@ -58,7 +57,6 @@ l_tmb <- fitModelParallel(formula = kij ~ genotype + environment + genotype:envi
                           n.cores = 1)
 ## -- extract results
 tidy_pub_fit <- tidy_results(l_tmb)
-
 ```
 
 ### Estimating parameter distributions