diff --git a/vignettes/rnaseq_notOnly.Rmd b/vignettes/rnaseq_notOnly.Rmd
index 8d01b7696d29a7ee575c83c7e71c8cefcaaf5733..e750ff6fa8601e57a5edbb0673a2cea1fbaa23e7 100644
--- a/vignettes/rnaseq_notOnly.Rmd
+++ b/vignettes/rnaseq_notOnly.Rmd
@@ -37,7 +37,7 @@ data("iris")
 plot(density(unlist(iris$Sepal.Width)))
 ```
 
-Given this observation, we opt to fit a model for each species using a Gaussian family model (default). Details about model family [here](https://www.rdocumentation.org/packages/stats/versions/3.6.2/topics/family).
+Given this observation, we opt to fit a model for each species using a Gaussian family model (default). Details about model family [here](https://search.r-project.org/CRAN/refmans/glmmTMB/html/nbinom2.html).
 
 
 ```{r example-fitModelParallel_nonRNA, warning = FALSE, message = FALSE}
@@ -45,7 +45,10 @@ l_tmb_iris <- fitModelParallel(
                 formula =  Sepal.Width ~ Sepal.Length ,
                 data = iris,
                 group_by = "Species",
+                family = stats::gaussian(),
                 n.cores = 1)
+```
 
+```{r example-print } 
 tidy_results(l_tmb_iris)
 ```