diff --git a/404.html b/404.html index dece77f..4b51e36 100644 --- a/404.html +++ b/404.html @@ -32,7 +32,7 @@
diff --git a/LICENSE-text.html b/LICENSE-text.html index 69ca8b3..0d3b856 100644 --- a/LICENSE-text.html +++ b/LICENSE-text.html @@ -17,7 +17,7 @@ diff --git a/LICENSE.html b/LICENSE.html index 832e1d0..fddb811 100644 --- a/LICENSE.html +++ b/LICENSE.html @@ -17,7 +17,7 @@ diff --git a/articles/HMP_2012_16S_gingival_V13.html b/articles/HMP_2012_16S_gingival_V13.html index 98c58f9..7968c90 100644 --- a/articles/HMP_2012_16S_gingival_V13.html +++ b/articles/HMP_2012_16S_gingival_V13.html @@ -32,7 +32,7 @@ @@ -333,7 +333,7 @@## user system elapsed
-## 48.049 0.508 48.291
+## 48.440 0.520 48.687
## Running method 1: DA_DESeq2.1 - 2024-12-04 23:53:06.44456
-## Running method 2: DA_DESeq2.1 - 2024-12-04 23:53:10.536615
-## Running method 3: DA_limma.1 - 2024-12-04 23:53:15.657624
-## Running method 4: DA_limma.1 - 2024-12-04 23:53:15.68115
-## Running method 5: DA_metagenomeSeq.1 - 2024-12-04 23:53:15.697864
-## Running method 6: DA_ALDEx2.1 - 2024-12-04 23:53:16.169794
-## Running method 7: DA_MAST.1 - 2024-12-04 23:53:17.964839
-## Running method 8: DA_Seurat.1 - 2024-12-04 23:53:19.185622
+## Running method 1: DA_DESeq2.1 - 2024-12-05 05:53:37.353985
+## Running method 2: DA_DESeq2.1 - 2024-12-05 05:53:41.482595
+## Running method 3: DA_limma.1 - 2024-12-05 05:53:46.682774
+## Running method 4: DA_limma.1 - 2024-12-05 05:53:46.706542
+## Running method 5: DA_metagenomeSeq.1 - 2024-12-05 05:53:46.723061
+## Running method 6: DA_ALDEx2.1 - 2024-12-05 05:53:47.202794
+## Running method 7: DA_MAST.1 - 2024-12-05 05:53:49.012416
+## Running method 8: DA_Seurat.1 - 2024-12-05 05:53:50.24576
## For a (much!) faster implementation of the Wilcoxon Rank Sum Test,
## (default method for FindMarkers) please install the presto package
## --------------------------------------------
@@ -246,23 +246,23 @@ Differential abundance analysis## After installation of presto, Seurat will automatically use the more
## efficient implementation (no further action necessary).
## This message will be shown once per session
-## Running method 9: ancombc.1 - 2024-12-04 23:53:19.740719
-## Running method 10: wilcox.3 - 2024-12-04 23:53:19.789119
-## Running method 11: wilcox.4 - 2024-12-04 23:53:19.995617
-## Running method 12: ZINQ.9 - 2024-12-04 23:53:20.125108
-## Running method 13: ZINQ.10 - 2024-12-04 23:53:21.568033
-## Running method 14: lefse.12 - 2024-12-04 23:53:22.718059
+## Running method 9: ancombc.1 - 2024-12-05 05:53:50.791175
+## Running method 10: wilcox.3 - 2024-12-05 05:53:50.839201
+## Running method 11: wilcox.4 - 2024-12-05 05:53:51.048368
+## Running method 12: ZINQ.9 - 2024-12-05 05:53:51.178807
+## Running method 13: ZINQ.10 - 2024-12-05 05:53:52.633139
+## Running method 14: lefse.12 - 2024-12-05 05:53:53.787877
## Running lefse original
## The outcome variable is specified as 'body_subsite' and the reference category is 'subgingival_plaque'.
## See `?factor` or `?relevel` to change the reference category.
-## Running method 15: lefse.13 - 2024-12-04 23:53:23.143877
+## Running method 15: lefse.13 - 2024-12-05 05:53:54.218306
## Running lefse original
## The outcome variable is specified as 'body_subsite' and the reference category is 'subgingival_plaque'.
## See `?factor` or `?relevel` to change the reference category.
tim
## user system elapsed
-## 16.917 0.092 16.995
+## 17.082 0.100 17.154
## user system elapsed
-## 8.849 0.127 8.993
+## 8.866 0.128 9.010
tss_fun <- function(x) (x + 1) / sum((x + 1))
+# tss_fun <- function(x) (x + 1) / sum((x + 1))
# tss_fun <- function(x) log((x + 1) / sum((x + 1)))
-# tss_fun <- function(x) (x) / sum((x))
gmn_fun <- function(x) (x + 1) / exp(mean(log((x + 1))))
# gnm_fun <- function(x) (x + 1) / prod((x + 1)^(1 / length(x)))
# gmn_fun <- function(x) log((x + 1) / exp(mean(log((x + 1)))))
-assay(tse, "TSS") <- apply(assay(tse, 'counts'), 2, tss_fun)
+assay(tse, "TSS") <- apply(assay(tse, 'counts'), 2, tss_fun)
assay(tse, "GMN") <- apply(assay(tse, 'counts'), 2, gmn_fun)
Get coefficient of variation for all features:
gmn_mat <- assay(tse, "GMN")
@@ -336,7 +336,7 @@ Ranking
bind_rows(.id = "feature")
})
#> user system elapsed
-#> 82.329 0.004 82.335
Add rankings
tss <- tss_cv |>
@@ -348,6 +348,7 @@ Ranking
ranking = min_rank(cv)
)
+
tss$mean <- apply(tss_mat, 1, mean)
gmn$mean <- apply(gmn_mat, 1, mean)
@@ -381,10 +382,10 @@ Ranking
spikein = factor(spikein, levels = c("Spike-in", "Non-spike-in"))
) |>
ggplot(aes(x = tss_ranking, y = gmn_ranking)) +
- geom_abline(slope = 1, intercept = 0, linetype = 2, linewidth = 0.3) +
geom_point(
aes(color = spikein, shape = spikein, size = spikein)
) +
+ geom_abline(slope = 1, intercept = 0, linetype = 2, linewidth = 0.3) +
geom_label_repel(
data = filter(dat, feature_label != ""),
mapping = aes(label = feature_label),
@@ -392,6 +393,8 @@ Ranking
)+
scale_fill_viridis_c(option = "C") +
scale_size_manual(values = c(3, 1)) +
+ scale_color_manual(values = c("gray20", "gray70")) +
+ scale_shape_manual(values = c(20, 4)) +
scale_x_continuous(labels = \(x) format(x, big.mark = ",")) +
scale_y_continuous(labels = \(x) format(x, big.mark = ",")) +
labs(
@@ -406,9 +409,9 @@ Ranking
p_b <- dat |>
ggplot(aes(x = tss_ranking, y = gmn_ranking)) +
geom_hex(
- color = "white"
+ color = "black", alpha = 0.9
) +
- geom_abline(slope = 1, intercept = 0, linetype = 2, linewidth = 0.3) +
+ geom_abline(slope = 1, intercept = 0, linetype = 2, linewidth = 0.3, color = "red") +
# geom_label_repel(
# data = filter(dat, feature_label != ""),
# mapping = aes(label = feature_label),
@@ -432,16 +435,17 @@ Ranking
plotlist = list(p_a, p_b),
labels = c("A", "B"),
nrow = 1
-)
-plts_1
+plts_1
hexFun <- function(myDat, xvar, yvar) {
xvar <- enquo(xvar)
yvar <- enquo(yvar)
myDat |>
ggplot(aes(!!xvar, !!yvar)) +
- geom_hex(color = "white") +
+ geom_hex(color = "black", alpha = 0.9) +
geom_label_repel(
data = filter(myDat, feature_label != ""),
mapping = aes(label = feature_label),
@@ -487,7 +491,7 @@ Ranking
ncol = 2, nrow = 2
)
plts
+diff --git a/authors.html b/authors.html index 2cd3aa4..e9a1c0f 100644 --- a/authors.html +++ b/authors.html @@ -17,7 +17,7 @@sessioninfo::session_info() #> ─ Session info ─────────────────────────────────────────────────────────────── #> setting value @@ -511,7 +515,7 @@
Session info#> collate en_US.UTF-8 #> ctype en_US.UTF-8 #> tz Etc/UTC -#> date 2024-12-04 +#> date 2024-12-05 #> pandoc 3.5 @ /usr/bin/ (via rmarkdown) #> #> ─ Packages ─────────────────────────────────────────────────────────────────── @@ -559,7 +563,7 @@
Session info#> fs 1.6.5 2024-10-30 [1] RSPM (R 4.4.0) #> generics 0.1.3 2022-07-05 [1] RSPM (R 4.4.0) #> GenomeInfoDb * 1.42.1 2024-11-28 [1] Bioconductor 3.20 (R 4.4.2) -#> GenomeInfoDbData 1.2.13 2024-12-04 [1] Bioconductor +#> GenomeInfoDbData 1.2.13 2024-12-05 [1] Bioconductor #> GenomicRanges * 1.58.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2) #> ggplot2 * 3.5.1 2024-04-23 [1] RSPM (R 4.4.0) #> ggpubr * 0.6.0 2023-02-10 [1] RSPM (R 4.4.0) @@ -587,7 +591,7 @@
Session info#> matrixStats * 1.4.1 2024-09-08 [1] RSPM (R 4.4.0) #> memoise 2.0.1 2021-11-26 [1] RSPM (R 4.4.0) #> MicrobiomeBenchmarkData * 1.8.0 2024-10-31 [1] Bioconductor 3.20 (R 4.4.2) -#> MicrobiomeBenchmarkDataAnalyses * 0.99.11 2024-12-04 [1] local +#> MicrobiomeBenchmarkDataAnalyses * 0.99.12 2024-12-05 [1] local #> modeltools 0.2-23 2020-03-05 [1] RSPM (R 4.4.0) #> multcomp 1.4-26 2024-07-18 [1] RSPM (R 4.4.0) #> munsell 0.5.1 2024-04-01 [1] RSPM (R 4.4.0) diff --git a/articles/Stammler_2016_16S_spikein_files/figure-html/unnamed-chunk-15-1.png b/articles/Stammler_2016_16S_spikein_files/figure-html/unnamed-chunk-15-1.png index 727d987..9549473 100644 Binary files a/articles/Stammler_2016_16S_spikein_files/figure-html/unnamed-chunk-15-1.png and b/articles/Stammler_2016_16S_spikein_files/figure-html/unnamed-chunk-15-1.png differ diff --git a/articles/Stammler_2016_16S_spikein_files/figure-html/unnamed-chunk-17-1.png b/articles/Stammler_2016_16S_spikein_files/figure-html/unnamed-chunk-17-1.png new file mode 100644 index 0000000..4347448 Binary files /dev/null and b/articles/Stammler_2016_16S_spikein_files/figure-html/unnamed-chunk-17-1.png differ diff --git a/articles/index.html b/articles/index.html index c8ab9cb..7ea349b 100644 --- a/articles/index.html +++ b/articles/index.html @@ -17,7 +17,7 @@
Gamboa-Tuz S (2024). MicrobiomeBenchmarkDataAnalyses: Benchmarking analyses with the MicrobiomeBenchmarkData package. -R package version 0.99.11, https://github.com/waldronlab/MicrobiomeBenchmarkDataAnalyses. +R package version 0.99.12, https://github.com/waldronlab/MicrobiomeBenchmarkDataAnalyses.
@Manual{, title = {MicrobiomeBenchmarkDataAnalyses: Benchmarking analyses with the MicrobiomeBenchmarkData package}, author = {Samuel David Gamboa-Tuz}, year = {2024}, - note = {R package version 0.99.11}, + note = {R package version 0.99.12}, url = {https://github.com/waldronlab/MicrobiomeBenchmarkDataAnalyses}, }diff --git a/index.html b/index.html index f09bef1..c1ce858 100644 --- a/index.html +++ b/index.html @@ -33,7 +33,7 @@