__STDOUT LOG__ Job runs on fgcz-h-166 at /scratch/PostSamsa2Analysis_after--over--before_2024-12-13--10-50-21_temp2136341 Starting EzAppPostSamsa2Analysis PostSamsa2Analysis PostSamsa2Analysis_after--over--before_2024-12-13--10-50-21_temp2136341 2024-12-13 10:50:27 mkdir orgresults mkdir funcresults cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/06_006_a.RefSeq.annotated_org.txt orgresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/06_006_b.RefSeq.annotated_org.txt orgresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/07_001_a.RefSeq.annotated_org.txt orgresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/07_001_b.RefSeq.annotated_org.txt orgresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/15_006_a.RefSeq.annotated_org.txt orgresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/15_006_b.RefSeq.annotated_org.txt orgresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/17_004_a.RefSeq.annotated_org.txt orgresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/17_004_b.RefSeq.annotated_org.txt orgresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/06_006_a.RefSeq.annotated_func.txt funcresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/06_006_b.RefSeq.annotated_func.txt funcresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/07_001_a.RefSeq.annotated_func.txt funcresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/07_001_b.RefSeq.annotated_func.txt funcresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/15_006_a.RefSeq.annotated_func.txt funcresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/15_006_b.RefSeq.annotated_func.txt funcresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/17_004_a.RefSeq.annotated_func.txt funcresults/ cp /srv/GT/analysis/course_sushi/public/gstore/projects/p8002/Samsa2_2023-12-12--09-56-50/17_004_b.RefSeq.annotated_func.txt funcresults/ [1] "06_006_a.RefSeq.annotated_org.txt" "06_006_b.RefSeq.annotated_org.txt" [3] "07_001_a.RefSeq.annotated_org.txt" "07_001_b.RefSeq.annotated_org.txt" [5] "15_006_a.RefSeq.annotated_org.txt" "15_006_b.RefSeq.annotated_org.txt" [7] "17_004_a.RefSeq.annotated_org.txt" "17_004_b.RefSeq.annotated_org.txt" [[1]] [1] "p8002/Samsa2_2023-12-12--09-56-50/06_006_a.RefSeq.annotated.txt" [[2]] [1] "p8002/Samsa2_2023-12-12--09-56-50/07_001_a.RefSeq.annotated.txt" [[3]] [1] "p8002/Samsa2_2023-12-12--09-56-50/15_006_a.RefSeq.annotated.txt" [[4]] [1] "p8002/Samsa2_2023-12-12--09-56-50/17_004_a.RefSeq.annotated.txt" [1] "06_006_a.RefSeq.annotated.txt" "07_001_a.RefSeq.annotated.txt" [3] "15_006_a.RefSeq.annotated.txt" "17_004_a.RefSeq.annotated.txt" [1] "before" [1] "Medication" [1] "06_006_a" "07_001_a" "15_006_a" "17_004_a" [1] "06_006_a.RefSeq.annotated_org.txt" "07_001_a.RefSeq.annotated_org.txt" [3] "15_006_a.RefSeq.annotated_org.txt" "17_004_a.RefSeq.annotated_org.txt" [1] "06_006_a.RefSeq.annotated_func.txt" "06_006_b.RefSeq.annotated_func.txt" [3] "07_001_a.RefSeq.annotated_func.txt" "07_001_b.RefSeq.annotated_func.txt" [5] "15_006_a.RefSeq.annotated_func.txt" "15_006_b.RefSeq.annotated_func.txt" [7] "17_004_a.RefSeq.annotated_func.txt" "17_004_b.RefSeq.annotated_func.txt" [[1]] [1] "p8002/Samsa2_2023-12-12--09-56-50/06_006_a.RefSeq.annotated.txt" [[2]] [1] "p8002/Samsa2_2023-12-12--09-56-50/07_001_a.RefSeq.annotated.txt" [[3]] [1] "p8002/Samsa2_2023-12-12--09-56-50/15_006_a.RefSeq.annotated.txt" [[4]] [1] "p8002/Samsa2_2023-12-12--09-56-50/17_004_a.RefSeq.annotated.txt" [1] "06_006_a.RefSeq.annotated.txt" "07_001_a.RefSeq.annotated.txt" [3] "15_006_a.RefSeq.annotated.txt" "17_004_a.RefSeq.annotated.txt" [1] "before" [1] "Medication" [1] "06_006_a" "07_001_a" "15_006_a" "17_004_a" [1] "06_006_a.RefSeq.annotated_func.txt" "07_001_a.RefSeq.annotated_func.txt" [3] "15_006_a.RefSeq.annotated_func.txt" "17_004_a.RefSeq.annotated_func.txt" rename s/_/-/g */*.txt rename s/control-/control_/g */*.txt rename s/experimental-/experimental_/g */*.txt bash /usr/local/ngseq/src/samsa2/R_scripts/run_all_Rscripts.sh orgresults/ [1] "USAGE: $ diversity_stats.R -I working_directory/ " Working directory is orgresults/ [1] "USAGE: $ diversity_graphs.R -I working_directory/ -O save.filename" Working directory is orgresults/ Saving diversity graphs as diversity_graph.rds Success! Saving diversity graphs as diversity_graph.rds now. [1] "USAGE: $ get_raw_counts_table.R -I working_directory/ -O save.filename" Working directory is orgresults/ [1] "WARNING: no raw counts file specified, skipping this info for DESeq analysis." Now saving raw counts table as raw_counts_table.tab [1] "USAGE: $ get_normalized_counts_table.R -I working_directory/ -O save.filename" Working directory is orgresults/ [1] "WARNING: no raw counts file specified, skipping this info for DESeq analysis." Success! Saving normalized counts table as normalized_counts_table.tab [1] "USAGE: $ run_DESeq_stats.R -I working_directory/ -O save.filename" Working directory is orgresults/ [1] "WARNING: no raw counts file specified, skipping this info for DESeq analysis." Success! Saving results file as DESeq_results.tab [1] "USAGE: $ run_DESeq_stats.R -I working_directory/ -O save.filename -L level (1,2,3,4)" Working directory is orgresults/ Saving results as PCA_plot.tab.rds Calculating DESeq results for hierarchy level Saving PCA plot as PCA_plot.tab.rds now. [1] "USAGE: $ make_DESeq_heatmap -I working_directory/ -O save.filename -N number of row entries (50)" Working directory is orgresults/ Saving results as DESeq_heatmap.rds Loading packages, may take a second... [1] "./control_07-001-a.RefSeq.annotated-org.txt" [1] "./control_15-006-a.RefSeq.annotated-org.txt" [1] "./control_17-004-a.RefSeq.annotated-org.txt" [1] "./experimental_07-001-b.RefSeq.annotated-org.txt" [1] "./experimental_15-006-b.RefSeq.annotated-org.txt" [1] "./experimental_17-004-b.RefSeq.annotated-org.txt" Calculating DESeq results, may take a second... Saving heatmap matrix as DESeq_heatmap.rds now. [1] "USAGE: $ make_volcano_plot.R -D working_directory/ -O save.filename -G title -F fc_cutoff -P p_cutoff -C TRUE" Working directory is orgresults/ Saving results as volcano_plot.rds Organism.Name baseMean controlMean experimentalMean log2FoldChange 1 Gardnerella_vaginalis 170843.80 165204.68 176482.93 0.09527971 2 Gardnerella 61614.91 73868.78 49361.04 -0.58158296 3 Atopobium_vaginae 27831.54 18682.14 36980.94 0.98511870 4 Bifidobacterium_longum 27702.25 31810.91 23593.60 -0.43111832 5 Sneathia_sanguinegens 25395.52 23690.11 27100.93 0.19408269 6 Gardnerella_sp._30-4 20279.91 24548.17 16011.65 -0.61648232 lfcSE stat pvalue padj 1 0.7766250 0.1226843 0.9023571 0.9993366 2 1.0506848 -0.5535275 0.5799022 0.9993366 3 0.9849079 1.0002140 0.3172070 0.9993366 4 0.8061202 -0.5348065 0.5927837 0.9993366 5 1.3982497 0.1388040 0.8896050 0.9993366 6 1.1899430 -0.5180772 0.6044044 0.9993366 'data.frame': 5549 obs. of 3 variables: $ Organism.Name : chr "Gardnerella_vaginalis" "Gardnerella" "Atopobium_vaginae" "Bifidobacterium_longum" ... $ log2FoldChange: num 0.0953 -0.5816 0.9851 -0.4311 0.1941 ... $ padj : num 0.999 0.999 0.999 0.999 0.999 ... Organism.Name log2FoldChange padj Length:5549 Min. :-6.7533 Min. :0.9993 Class :character 1st Qu.:-0.6601 1st Qu.:0.9993 Mode :character Median :-0.1664 Median :0.9993 Mean :-0.1043 Mean :0.9993 3rd Qu.: 0.7610 3rd Qu.:0.9993 Max. : 4.9681 Max. :1.0000 NA's :51 Saving volcano plot as volcano_plot.rds now. [1] "USAGE: $ make_combined_graphs.R -D working_directory/ -O save.filename -N org_number -G title" Working directory is orgresults/ Saving results as combined_graphs_ctrl_vs_exp.rds Saving PCA plot as combined_graphs_ctrl_vs_exp.rds now. [1] "USAGE: $ make_combined_graphs.R -D working_directory/ -O save.filename -N org_number -G title" Working directory is orgresults/ Saving results as combined_graph.rds Saving PCA plot as combined_graph.rds now. bash /usr/local/ngseq/src/samsa2/R_scripts/run_all_Rscripts.sh funcresults/ [1] "USAGE: $ diversity_stats.R -I working_directory/ " Working directory is funcresults/ [1] "USAGE: $ diversity_graphs.R -I working_directory/ -O save.filename" Working directory is funcresults/ Saving diversity graphs as diversity_graph.rds Success! Saving diversity graphs as diversity_graph.rds now. [1] "USAGE: $ get_raw_counts_table.R -I working_directory/ -O save.filename" Working directory is funcresults/ [1] "WARNING: no raw counts file specified, skipping this info for DESeq analysis." Now saving raw counts table as raw_counts_table.tab [1] "USAGE: $ get_normalized_counts_table.R -I working_directory/ -O save.filename" Working directory is funcresults/ [1] "WARNING: no raw counts file specified, skipping this info for DESeq analysis." Success! Saving normalized counts table as normalized_counts_table.tab [1] "USAGE: $ run_DESeq_stats.R -I working_directory/ -O save.filename" Working directory is funcresults/ [1] "WARNING: no raw counts file specified, skipping this info for DESeq analysis." Success! Saving results file as DESeq_results.tab [1] "USAGE: $ run_DESeq_stats.R -I working_directory/ -O save.filename -L level (1,2,3,4)" Working directory is funcresults/ Saving results as PCA_plot.tab.rds Calculating DESeq results for hierarchy level Saving PCA plot as PCA_plot.tab.rds now. [1] "USAGE: $ make_DESeq_heatmap -I working_directory/ -O save.filename -N number of row entries (50)" Working directory is funcresults/ Saving results as DESeq_heatmap.rds Loading packages, may take a second... [1] "./control_07-001-a.RefSeq.annotated-func.txt" [1] "./control_15-006-a.RefSeq.annotated-func.txt" [1] "./control_17-004-a.RefSeq.annotated-func.txt" [1] "./experimental_07-001-b.RefSeq.annotated-func.txt" [1] "./experimental_15-006-b.RefSeq.annotated-func.txt" [1] "./experimental_17-004-b.RefSeq.annotated-func.txt" Calculating DESeq results, may take a second... Saving heatmap matrix as DESeq_heatmap.rds now. [1] "USAGE: $ make_volcano_plot.R -D working_directory/ -O save.filename -G title -F fc_cutoff -P p_cutoff -C TRUE" Working directory is funcresults/ Saving results as volcano_plot.rds Organism.Name baseMean controlMean 1 hypothetical_protein 90023.364 91335.202 2 ABC_transporter_substrate-binding_protein 9055.276 9090.082 3 hypothetical_protein,_partial 6405.709 5904.871 4 peptidase 5339.821 5304.153 5 ABC_transporter_ATP-binding_protein 5308.037 5180.366 6 MULTISPECIES:_hypothetical_protein 4643.327 5086.042 experimentalMean log2FoldChange lfcSE stat pvalue padj 1 88711.526 -0.04204911 0.3314697 -0.1268566 0.8990539 0.9998296 2 9020.469 -0.01108259 0.4967186 -0.0223116 0.9821994 0.9998296 3 6906.548 0.22605579 0.4784065 0.4725183 0.6365569 0.9998296 4 5375.489 0.01926109 0.6233563 0.0308990 0.9753501 0.9998296 5 5435.708 0.06942312 0.3385138 0.2050821 0.8375080 0.9998296 6 4200.613 -0.27596087 0.5728664 -0.4817194 0.6300053 0.9998296 'data.frame': 4747 obs. of 3 variables: $ Organism.Name : chr "hypothetical_protein" "ABC_transporter_substrate-binding_protein" "hypothetical_protein,_partial" "peptidase" ... $ log2FoldChange: num -0.042 -0.0111 0.2261 0.0193 0.0694 ... $ padj : num 1 1 1 1 1 ... Organism.Name log2FoldChange padj Length:4747 Min. :-22.97997 Min. :0.5912 Class :character 1st Qu.: -0.54550 1st Qu.:0.9998 Mode :character Median : 0.04611 Median :0.9998 Mean : 0.11510 Mean :0.9997 3rd Qu.: 0.71415 3rd Qu.:0.9998 Max. : 6.08286 Max. :0.9999 NA's :1 NA's :106 Saving volcano plot as volcano_plot.rds now. [1] "USAGE: $ make_combined_graphs.R -D working_directory/ -O save.filename -N org_number -G title" Working directory is funcresults/ Saving results as combined_graphs_ctrl_vs_exp.rds Saving PCA plot as combined_graphs_ctrl_vs_exp.rds now. [1] "USAGE: $ make_combined_graphs.R -D working_directory/ -O save.filename -N org_number -G title" Working directory is funcresults/ Saving results as combined_graph.rds Saving PCA plot as combined_graph.rds now. ezRun tag: 81a203b2a83152a315eaa7463292c4b4d9afaa20 ezRun github link: https://github.com/uzh/ezRun/tree/81a203b2a83152a315eaa7463292c4b4d9afaa20 R version 4.4.0 (2024-04-24) Platform: x86_64-pc-linux-gnu Running under: Debian GNU/Linux 12 (bookworm) Matrix products: default BLAS: /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.11.0 LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.11.0 locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=en_US.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C time zone: Europe/Zurich tzcode source: system (glibc) attached base packages: [1] stats4 stats graphics grDevices utils datasets methods [8] base other attached packages: [1] htmltools_0.5.8.1 DT_0.33 [3] pheatmap_1.0.12 rmarkdown_2.27 [5] vegan_2.6-6.1 lattice_0.22-6 [7] permute_0.9-7 scales_1.3.0 [9] reshape2_1.4.4 plyr_1.8.9 [11] optparse_1.7.5 knitr_1.47 [13] gridExtra_2.3 genefilter_1.86.0 [15] RColorBrewer_1.1-3 EnhancedVolcano_1.22.0 [17] ggrepel_0.9.5 DESeq2_1.44.0 [19] SummarizedExperiment_1.34.0 Biobase_2.64.0 [21] MatrixGenerics_1.16.0 matrixStats_1.3.0 [23] ezRun_3.19.1 lubridate_1.9.3 [25] forcats_1.0.0 stringr_1.5.1 [27] dplyr_1.1.4 purrr_1.0.2 [29] readr_2.1.5 tidyr_1.3.1 [31] tibble_3.2.1 ggplot2_3.5.1 [33] tidyverse_2.0.0 GenomicRanges_1.56.1 [35] Biostrings_2.72.1 GenomeInfoDb_1.40.1 [37] XVector_0.44.0 IRanges_2.38.0 [39] S4Vectors_0.42.0 BiocGenerics_0.50.0 [41] data.table_1.15.4 loaded via a namespace (and not attached): [1] DBI_1.2.3 rlang_1.1.4 magrittr_2.0.3 [4] compiler_4.4.0 RSQLite_2.3.7 mgcv_1.9-1 [7] png_0.1-8 vctrs_0.6.5 pkgconfig_2.0.3 [10] crayon_1.5.2 fastmap_1.2.0 labeling_0.4.3 [13] utf8_1.2.4 tzdb_0.4.0 UCSC.utils_1.0.0 [16] bit_4.0.5 xfun_0.45 zlibbioc_1.50.0 [19] cachem_1.1.0 jsonlite_1.8.8 blob_1.2.4 [22] highr_0.11 DelayedArray_0.30.1 BiocParallel_1.38.0 [25] parallel_4.4.0 cluster_2.1.6 R6_2.5.1 [28] bslib_0.7.0 stringi_1.8.4 jquerylib_0.1.4 [31] Rcpp_1.0.12 Matrix_1.7-0 splines_4.4.0 [34] timechange_0.3.0 tidyselect_1.2.1 yaml_2.3.8 [37] abind_1.4-5 codetools_0.2-20 withr_3.0.0 [40] KEGGREST_1.44.0 evaluate_0.24.0 survival_3.7-0 [43] getopt_1.20.4 pillar_1.9.0 plotly_4.10.4 [46] generics_0.1.3 hms_1.1.3 munsell_0.5.1 [49] xtable_1.8-4 glue_1.7.0 lazyeval_0.2.2 [52] tools_4.4.0 annotate_1.82.0 locfit_1.5-9.9 [55] XML_3.99-0.16.1 grid_4.4.0 crosstalk_1.2.1 [58] AnnotationDbi_1.66.0 colorspace_2.1-1 nlme_3.1-165 [61] GenomeInfoDbData_1.2.12 cli_3.6.2 fansi_1.0.6 [64] viridisLite_0.4.2 S4Arrays_1.4.1 gtable_0.3.5 [67] sass_0.4.9 digest_0.6.35 SparseArray_1.4.8 [70] farver_2.1.2 htmlwidgets_1.6.4 memoise_2.0.1 [73] lifecycle_1.0.4 httr_1.4.7 bit64_4.0.5 [76] MASS_7.3-61 Finished EzAppPostSamsa2Analysis PostSamsa2Analysis PostSamsa2Analysis_after--over--before_2024-12-13--10-50-21_temp2136341 2024-12-13 10:53:12 [1] "/scratch/PostSamsa2Analysis_after--over--before_2024-12-13--10-50-21_temp2136341/PostSamsa2Analysis/00index.html" __SCRIPT END__ __STDERR LOG__ Loading required package: data.table Loading required package: Biostrings Loading required package: BiocGenerics Attaching package: ‘BiocGenerics’ The following objects are masked from ‘package:stats’: IQR, mad, sd, var, xtabs The following objects are masked from ‘package:base’: anyDuplicated, aperm, append, as.data.frame, basename, cbind, colnames, dirname, do.call, duplicated, eval, evalq, Filter, Find, get, grep, grepl, intersect, is.unsorted, lapply, Map, mapply, match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank, rbind, Reduce, rownames, sapply, setdiff, table, tapply, union, unique, unsplit, which.max, which.min Loading required package: S4Vectors Loading required package: stats4 Attaching package: ‘S4Vectors’ The following objects are masked from ‘package:data.table’: first, second The following object is masked from ‘package:utils’: findMatches The following objects are masked from ‘package:base’: expand.grid, I, unname Loading required package: IRanges Attaching package: ‘IRanges’ The following object is masked from ‘package:data.table’: shift Loading required package: XVector Loading required package: GenomeInfoDb Attaching package: ‘Biostrings’ The following object is masked from ‘package:base’: strsplit Loading required package: GenomicRanges Loading required package: tidyverse ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ── ✔ dplyr 1.1.4 ✔ readr 2.1.5 ✔ forcats 1.0.0 ✔ stringr 1.5.1 ✔ ggplot2 3.5.1 ✔ tibble 3.2.1 ✔ lubridate 1.9.3 ✔ tidyr 1.3.1 ✔ purrr 1.0.2 ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ── ✖ lubridate::%within%() masks IRanges::%within%() ✖ dplyr::between() masks data.table::between() ✖ dplyr::collapse() masks Biostrings::collapse(), IRanges::collapse() ✖ dplyr::combine() masks BiocGenerics::combine() ✖ purrr::compact() masks XVector::compact() ✖ dplyr::desc() masks IRanges::desc() ✖ tidyr::expand() masks S4Vectors::expand() ✖ dplyr::filter() masks stats::filter() ✖ dplyr::first() masks S4Vectors::first(), data.table::first() ✖ lubridate::hour() masks data.table::hour() ✖ lubridate::isoweek() masks data.table::isoweek() ✖ dplyr::lag() masks stats::lag() ✖ dplyr::last() masks data.table::last() ✖ lubridate::mday() masks data.table::mday() ✖ lubridate::minute() masks data.table::minute() ✖ lubridate::month() masks data.table::month() ✖ ggplot2::Position() masks BiocGenerics::Position(), base::Position() ✖ lubridate::quarter() masks data.table::quarter() ✖ purrr::reduce() masks GenomicRanges::reduce(), IRanges::reduce() ✖ dplyr::rename() masks S4Vectors::rename() ✖ lubridate::second() masks S4Vectors::second(), data.table::second() ✖ lubridate::second<-() masks S4Vectors::second<-() ✖ dplyr::slice() masks XVector::slice(), IRanges::slice() ✖ purrr::transpose() masks data.table::transpose() ✖ lubridate::wday() masks data.table::wday() ✖ lubridate::week() masks data.table::week() ✖ lubridate::yday() masks data.table::yday() ✖ lubridate::year() masks data.table::year() ℹ Use the conflicted package () to force all conflicts to become errors unknown param: partition unknown param: sushi_app unknown param: isLastJob Loading required package: SummarizedExperiment Loading required package: MatrixGenerics Loading required package: matrixStats Attaching package: ‘matrixStats’ The following object is masked from ‘package:dplyr’: count Attaching package: ‘MatrixGenerics’ The following objects are masked from ‘package:matrixStats’: colAlls, colAnyNAs, colAnys, colAvgsPerRowSet, colCollapse, colCounts, colCummaxs, colCummins, colCumprods, colCumsums, colDiffs, colIQRDiffs, colIQRs, colLogSumExps, colMadDiffs, colMads, colMaxs, colMeans2, colMedians, colMins, colOrderStats, colProds, colQuantiles, colRanges, colRanks, colSdDiffs, colSds, colSums2, colTabulates, colVarDiffs, colVars, colWeightedMads, colWeightedMeans, colWeightedMedians, colWeightedSds, colWeightedVars, rowAlls, rowAnyNAs, rowAnys, rowAvgsPerColSet, rowCollapse, rowCounts, rowCummaxs, rowCummins, rowCumprods, rowCumsums, rowDiffs, rowIQRDiffs, rowIQRs, rowLogSumExps, rowMadDiffs, rowMads, rowMaxs, rowMeans2, rowMedians, rowMins, rowOrderStats, rowProds, rowQuantiles, rowRanges, rowRanks, rowSdDiffs, rowSds, rowSums2, rowTabulates, rowVarDiffs, rowVars, rowWeightedMads, rowWeightedMeans, rowWeightedMedians, rowWeightedSds, rowWeightedVars Loading required package: Biobase Welcome to Bioconductor Vignettes contain introductory material; view with 'browseVignettes()'. To cite Bioconductor, see 'citation("Biobase")', and for packages 'citation("pkgname")'. Attaching package: ‘Biobase’ The following object is masked from ‘package:MatrixGenerics’: rowMedians The following objects are masked from ‘package:matrixStats’: anyMissing, rowMedians Loading required package: ggrepel Attaching package: ‘genefilter’ The following objects are masked from ‘package:MatrixGenerics’: rowSds, rowVars The following objects are masked from ‘package:matrixStats’: rowSds, rowVars The following object is masked from ‘package:readr’: spec Attaching package: ‘gridExtra’ The following object is masked from ‘package:Biobase’: combine The following object is masked from ‘package:dplyr’: combine The following object is masked from ‘package:BiocGenerics’: combine ------------------------------------------------------------------------------ You have loaded plyr after dplyr - this is likely to cause problems. If you need functions from both plyr and dplyr, please load plyr first, then dplyr: library(plyr); library(dplyr) ------------------------------------------------------------------------------ Attaching package: ‘plyr’ The following object is masked from ‘package:matrixStats’: count The following objects are masked from ‘package:dplyr’: arrange, count, desc, failwith, id, mutate, rename, summarise, summarize The following object is masked from ‘package:purrr’: compact The following object is masked from ‘package:XVector’: compact The following object is masked from ‘package:IRanges’: desc The following object is masked from ‘package:S4Vectors’: rename Attaching package: ‘reshape2’ The following object is masked from ‘package:tidyr’: smiths The following objects are masked from ‘package:data.table’: dcast, melt Attaching package: ‘scales’ The following object is masked from ‘package:purrr’: discard The following object is masked from ‘package:readr’: col_factor Loading required package: permute Loading required package: lattice This is vegan 2.6-6.1 Loading required package: rmarkdown converting counts to integer mode estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing converting counts to integer mode estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing converting counts to integer mode estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing using ntop=500 top features by variance converting counts to integer mode estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing Warning message: Removed 1 row containing missing values or values outside the scale range (`geom_bar()`). converting counts to integer mode estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing converting counts to integer mode estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing converting counts to integer mode estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing using ntop=500 top features by variance converting counts to integer mode estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing Warning messages: 1: In scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : EOF within quoted string 2: In scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : number of items read is not a multiple of the number of columns Warning message: Removed 1 row containing missing values or values outside the scale range (`geom_bar()`).