@@ -147,53 +147,6 @@ airway = airway |>
147147 test_differential_expression(~ dex , method = " limma_voom_sample_weights" , prefix = " voom_weights__" )
148148```
149149
150- ## Warning: The `.abundance` argument of `test_differential_abundance()` is deprecated as
151- ## of tidybulk 2.0.0.
152- ## ℹ Please use the `abundance` argument instead.
153- ## ℹ The deprecated feature was likely used in the tidybulk package.
154- ## Please report the issue at <https://github.com/stemangiola/tidybulk/issues>.
155- ## This warning is displayed once every 8 hours.
156- ## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
157- ## generated.
158-
159- ## =====================================
160- ## tidybulk says: All testing methods use raw counts, irrespective of if scale_abundance
161- ## or adjust_abundance have been calculated. Therefore, it is essential to add covariates
162- ## such as batch effects (if applicable) in the formula.
163- ## =====================================
164- ## tidybulk says: The design column names are "(Intercept), dexuntrt"
165- ##
166- ## tidybulk says: to access the DE object do `metadata(.)$tidybulk$edgeR_quasi_likelihood_object`
167- ## tidybulk says: to access the raw results (fitted GLM) do `metadata(.)$tidybulk$edgeR_quasi_likelihood_fit`
168- ## tidybulk says: The design column names are "(Intercept), dexuntrt"
169- ##
170- ## tidybulk says: to access the DE object do `metadata(.)$tidybulk$edger_robust_likelihood_ratio_object`
171- ## tidybulk says: to access the raw results (fitted GLM) do `metadata(.)$tidybulk$edger_robust_likelihood_ratio_fit`
172- ## estimating size factors
173- ##
174- ## estimating dispersions
175- ##
176- ## gene-wise dispersion estimates
177- ##
178- ## mean-dispersion relationship
179- ##
180- ## final dispersion estimates
181- ##
182- ## fitting model and testing
183- ##
184- ## tidybulk says: to access the DE object do `metadata(.)$tidybulk$DESeq2_object`
185- ## tidybulk says: to access the raw results (fitted GLM) do `metadata(.)$tidybulk$DESeq2_fit`
186- ## tidybulk says: The design column names are "(Intercept), dexuntrt"
187- ##
188- ## tidybulk says: to access the DE object do `metadata(.)$tidybulk$limma_voom_object`
189- ## tidybulk says: to access the raw results (fitted GLM) do `metadata(.)$tidybulk$limma_voom_fit`
190- ## tidybulk says: The design column names are "(Intercept), dexuntrt"
191- ##
192- ## tidybulk says: to access the DE object do `metadata(.)$tidybulk$limma_voom_sample_weights_object`
193- ## tidybulk says: to access the raw results (fitted GLM) do `metadata(.)$tidybulk$limma_voom_sample_weights_fit`
194- ## This message is displayed once per session.
195-
196-
197150
198151### Histograms of p-values across methods
199152
@@ -275,34 +228,6 @@ airway |>
275228 labs(title = " Pairplot of p-values across methods" )
276229```
277230
278- ## Warning in ggally_statistic(data = data, mapping = mapping, na.rm = na.rm, :
279- ## Removed 30 rows containing missing values
280-
281- ## Warning in ggally_statistic(data = data, mapping = mapping, na.rm = na.rm, :
282- ## Removed 30 rows containing missing values
283- ## Warning in ggally_statistic(data = data, mapping = mapping, na.rm = na.rm, :
284- ## Removed 30 rows containing missing values
285- ## Warning in ggally_statistic(data = data, mapping = mapping, na.rm = na.rm, :
286- ## Removed 30 rows containing missing values
287-
288- ## Warning in scale_y_continuous(trans = tidybulk::log10_reverse_trans()):
289- ## log10_reverse transformation introduced infinite values.
290-
291- ## Warning: Removed 30 rows containing missing values or values outside the scale range
292- ## (`geom_point()`).
293- ## Removed 30 rows containing missing values or values outside the scale range
294- ## (`geom_point()`).
295- ## Removed 30 rows containing missing values or values outside the scale range
296- ## (`geom_point()`).
297- ## Removed 30 rows containing missing values or values outside the scale range
298- ## (`geom_point()`).
299-
300- ## Warning: Removed 30 rows containing non-finite outside the scale range
301- ## (`stat_density()`).
302-
303- ## Warning in scale_y_continuous(trans = tidybulk::log10_reverse_trans()):
304- ## log10_reverse transformation introduced infinite values.
305-
306231![ ] ( README_files/figure-gfm/differential-expression-pvalue-pairplot-1.png ) <!-- -->
307232
308233### Pairplot of effect sizes across methods (GGpairs)
0 commit comments