diff --git a/dev/manuscript_differential_transcript_abundance.Rmd b/dev/manuscript_differential_transcript_abundance.Rmd
index c5038614..48d6e8ed 100644
--- a/dev/manuscript_differential_transcript_abundance.Rmd
+++ b/dev/manuscript_differential_transcript_abundance.Rmd
@@ -156,7 +156,7 @@ tt_mds_adj_mds %>%
 ```
 
 ```{r}
-tt_test <- tt_adj %>% test_differential_abundance(~ condition + type)
+tt_test <- tt_adj %>% test_differential_abundance(~ condition + type, method = "edgeR_quasi_likelihood")
 
 # MA plot
 tt_test %>%
diff --git a/vignettes/introduction.Rmd b/vignettes/introduction.Rmd
index 55d6e1fb..68d32b6e 100755
--- a/vignettes/introduction.Rmd
+++ b/vignettes/introduction.Rmd
@@ -439,7 +439,7 @@ TidyTranscriptomics
 ```{r de, message=FALSE, warning=FALSE, results='hide'}
 se_mini.de =
 	se_mini |>
-	test_differential_abundance( ~ condition, action="get")
+	test_differential_abundance( ~ condition, method = "edgeR_quasi_likelihood", action="get")
 se_mini.de
 ```
 </div>
@@ -468,6 +468,7 @@ se_mini.de =
 	identify_abundant(factor_of_interest = condition) |>
 	test_differential_abundance(
 		~ 0 + condition,                  
+		method = "edgeR_quasi_likelihood",
 		.contrasts = c( "conditionTRUE - conditionFALSE"),
 		action="get"
 	)
@@ -681,6 +682,7 @@ se_mini.norm.SNN |>
 	mutate(factor_of_interest = `cluster_SNN` == 3) |>
 	test_differential_abundance(
     ~ factor_of_interest,
+    method = "edgeR_quasi_likelihood",
     action="get"
    )
 ```