diff --git a/vignettes/rnaseq_notOnly.Rmd b/vignettes/rnaseq_notOnly.Rmd
index a1bc5d957f1cfdc898ae7bab5f185b8cc1a5fea3..8d01b7696d29a7ee575c83c7e71c8cefcaaf5733 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. 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://www.rdocumentation.org/packages/stats/versions/3.6.2/topics/family).
 
 
 ```{r example-fitModelParallel_nonRNA, warning = FALSE, message = FALSE}
@@ -45,7 +45,6 @@ l_tmb_iris <- fitModelParallel(
                 formula =  Sepal.Width ~ Sepal.Length ,
                 data = iris,
                 group_by = "Species",
-                family = gaussian(),
                 n.cores = 1)
 
 tidy_results(l_tmb_iris)