ModuleCmd_Load.c(213):ERROR:105: Unable to locate a modulefile for 'load' ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ── ✔ ggplot2 3.3.5 ✔ purrr 0.3.4 ✔ tibble 3.1.5 ✔ dplyr 1.0.7 ✔ tidyr 1.1.4 ✔ stringr 1.4.0 ✔ readr 2.0.2 ✔ forcats 0.5.1 ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ── ✖ dplyr::filter() masks stats::filter() ✖ dplyr::lag() masks stats::lag() Loading required package: S4Vectors Loading required package: stats4 Loading required package: BiocGenerics Attaching package: ‘BiocGenerics’ The following objects are masked from ‘package:dplyr’: combine, intersect, setdiff, union The following objects are masked from ‘package:stats’: IQR, mad, sd, var, xtabs The following objects are masked from ‘package:base’: anyDuplicated, 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, sort, table, tapply, union, unique, unsplit, which.max, which.min Attaching package: ‘S4Vectors’ The following objects are masked from ‘package:dplyr’: first, rename The following object is masked from ‘package:tidyr’: expand The following objects are masked from ‘package:base’: expand.grid, I, unname Loading required package: IRanges Attaching package: ‘IRanges’ The following objects are masked from ‘package:dplyr’: collapse, desc, slice The following object is masked from ‘package:purrr’: reduce Loading required package: GenomicRanges Loading required package: GenomeInfoDb 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 converting counts to integer mode seqname source feature start end score strand frame 1 ##gff-version 3 NA NA 2 1 FIG CDS 1 1368 . + 1 3 1 FIG CDS 1379 2479 . + 2 4 1 FIG CDS 2485 3567 . + 1 5 1 FIG CDS 3631 4776 . - 1 6 1 FIG CDS 4813 5616 . - 1 att 1 2 ID=fig|729.2200.peg.1;Name=Chromosomal replication initiator protein DnaA 3 ID=fig|729.2200.peg.2;Name=DNA polymerase III beta subunit (EC 2.7.7.7);Ontology_term=KEGG_ENZYME:2.7.7.7 4 ID=fig|729.2200.peg.3;Name=DNA recombination and repair protein RecF 5 ID=fig|729.2200.peg.4;Name=N-acetylglucosamine-6-phosphate deacetylase (EC 3.5.1.25);Ontology_term=KEGG_ENZYME:3.5.1.25 6 ID=fig|729.2200.peg.5;Name=Glucosamine-6-phosphate deaminase (EC 3.5.99.6);Ontology_term=KEGG_ENZYME:3.5.99.6 seqname source feature start end score strand frame 2 1 FIG CDS 138 1231 . + 1 3 1 FIG CDS 1489 2369 . + 2 4 1 FIG CDS 2593 3459 . + 1 5 1 FIG CDS 3746 4662 . - 1 6 1 FIG CDS 4893 5536 . - 1 7 1 FIG CDS 5795 6497 . - 1 att 2 ID=fig|729.2200.peg.1;Name=Chromosomal replication initiator protein DnaA 3 ID=fig|729.2200.peg.2;Name=DNA polymerase III beta subunit (EC 2.7.7.7);Ontology_term=KEGG_ENZYME:2.7.7.7 4 ID=fig|729.2200.peg.3;Name=DNA recombination and repair protein RecF 5 ID=fig|729.2200.peg.4;Name=N-acetylglucosamine-6-phosphate deacetylase (EC 3.5.1.25);Ontology_term=KEGG_ENZYME:3.5.1.25 6 ID=fig|729.2200.peg.5;Name=Glucosamine-6-phosphate deaminase (EC 3.5.99.6);Ontology_term=KEGG_ENZYME:3.5.99.6 7 ID=fig|729.2200.peg.6;Name=N-acetylneuraminate lyase (EC 4.1.3.3);Ontology_term=KEGG_ENZYME:4.1.3.3 [1] "Generating pseudo-datasets" Joining, by = "Pos" [1] 1 Joining, by = "Pos" [1] 2 Joining, by = "Pos" [1] 3 Joining, by = "Pos" [1] 4 Joining, by = "Pos" [1] 5 Joining, by = "Pos" [1] 6 Joining, by = "Pos" [1] 7 Joining, by = "Pos" [1] 8 Joining, by = "Pos" [1] 9 Joining, by = "Pos" [1] 10 Joining, by = "Pos" [1] 11 Joining, by = "Pos" [1] 12 Joining, by = "Pos" [1] 13 Joining, by = "Pos" [1] 14 Joining, by = "Pos" [1] 15 Joining, by = "Pos" [1] 16 Joining, by = "Pos" [1] 17 Joining, by = "Pos" [1] 18 Joining, by = "Pos" [1] 19 Joining, by = "Pos" [1] 20 Joining, by = "Pos" [1] 21 Joining, by = "Pos" [1] 22 Joining, by = "Pos" [1] 23 Joining, by = "Pos" [1] 24 Joining, by = "Pos" [1] 25 Joining, by = "Pos" [1] 26 Joining, by = "Pos" [1] 27 Joining, by = "Pos" [1] 28 Joining, by = "Pos" [1] 29 Joining, by = "Pos" [1] 30 Joining, by = "Pos" [1] 31 Joining, by = "Pos" [1] 32 Joining, by = "Pos" [1] 33 Joining, by = "Pos" [1] 34 Joining, by = "Pos" [1] 35 Joining, by = "Pos" [1] 36 Joining, by = "Pos" [1] 37 Joining, by = "Pos" [1] 38 Joining, by = "Pos" [1] 39 Joining, by = "Pos" [1] 40 Joining, by = "Pos" [1] 41 Joining, by = "Pos" [1] 42 Joining, by = "Pos" [1] 43 Joining, by = "Pos" [1] 44 Joining, by = "Pos" [1] 45 Joining, by = "Pos" [1] 46 Joining, by = "Pos" [1] 47 Joining, by = "Pos" [1] 48 Joining, by = "Pos" [1] 49 Joining, by = "Pos" [1] 50 Joining, by = "Pos" [1] 51 Joining, by = "Pos" [1] 52 Joining, by = "Pos" [1] 53 Joining, by = "Pos" [1] 54 Joining, by = "Pos" [1] 55 Joining, by = "Pos" [1] 56 Joining, by = "Pos" [1] 57 Joining, by = "Pos" [1] 58 Joining, by = "Pos" [1] 59 Joining, by = "Pos" [1] 60 Joining, by = "Pos" [1] 61 Joining, by = "Pos" [1] 62 Joining, by = "Pos" [1] 63 Joining, by = "Pos" [1] 64 Joining, by = "Pos" [1] 65 Joining, by = "Pos" [1] 66 Joining, by = "Pos" [1] 67 Joining, by = "Pos" [1] 68 Joining, by = "Pos" [1] 69 Joining, by = "Pos" [1] 70 Joining, by = "Pos" [1] 71 Joining, by = "Pos" [1] 72 Joining, by = "Pos" [1] 73 Joining, by = "Pos" [1] 74 Joining, by = "Pos" [1] 75 Joining, by = "Pos" [1] 76 Joining, by = "Pos" [1] 77 Joining, by = "Pos" [1] 78 Joining, by = "Pos" [1] 79 Joining, by = "Pos" [1] 80 Joining, by = "Pos" [1] 81 Joining, by = "Pos" [1] 82 Joining, by = "Pos" [1] 83 Joining, by = "Pos" [1] 84 Joining, by = "Pos" [1] 85 Joining, by = "Pos" [1] 86 Joining, by = "Pos" [1] 87 Joining, by = "Pos" [1] 88 Joining, by = "Pos" [1] 89 Joining, by = "Pos" [1] 90 Joining, by = "Pos" [1] 91 Joining, by = "Pos" [1] 92 Joining, by = "Pos" [1] 93 Joining, by = "Pos" [1] 94 Joining, by = "Pos" [1] 95 Joining, by = "Pos" [1] 96 Joining, by = "Pos" [1] 97 Joining, by = "Pos" [1] 98 Joining, by = "Pos" [1] 99 Joining, by = "Pos" [1] 100 Joining, by = "Pos" [1] "Binning read counts by gene boundaries" converting counts to integer mode Warning message: In DESeqDataSet(se, design = design, ignoreRank) : some variables in design formula are characters, converting to factors estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing -- replacing outliers and refitting for 9 genes -- DESeq argument 'minReplicatesForReplace' = 7 -- original counts are preserved in counts(dds) estimating dispersions fitting model and testing log2 fold change (MLE): condition el1 vs Expected Wald test p-value: condition el1 vs Expected DataFrame with 6 rows and 6 columns baseMean log2FoldChange lfcSE stat pvalue padj 2 3422.79 -11.683926 0.936120 -12.481223 9.45240e-36 1.81290e-34 3 3442.73 -11.692542 0.926982 -12.613554 1.77790e-36 3.76508e-35 4 2756.74 -2.658549 0.519282 -5.119667 3.06076e-07 1.02852e-06 5 3333.18 -1.480146 0.522502 -2.832806 4.61414e-03 9.55249e-03 6 3043.79 -0.796511 0.585638 -1.360074 1.73806e-01 2.36875e-01 7 3104.43 -0.149578 0.481163 -0.310867 7.55901e-01 8.13524e-01 Package 'mclust' version 5.4.7 Type 'citation("mclust")' for citing this R package in publications. Attaching package: ‘mclust’ The following object is masked from ‘package:purrr’: map ---------------------------------------------------- Gaussian finite mixture model fitted by EM algorithm ---------------------------------------------------- Mclust V (univariate, unequal variance) model with 2 components: log-likelihood n df BIC ICL -5123.418 2033 5 -10284.92 -10468.88 Clustering table: 1 2 641 1392 Mixing probabilities: 1 2 0.3315414 0.6684586 Means: 1 2 -7.8372946 0.5209397 Variances: 1 2 13.029941 1.568428 ------------------------------------------------------- Density estimation via Gaussian finite mixture modeling ------------------------------------------------------- Mclust V (univariate, unequal variance) model with 2 components: log-likelihood n df BIC ICL -5123.421 2033 5 -10284.93 -10468.96 null device 1 [1] "Reduced" "Reduced" "Reduced" "Unchanged" "Unchanged" "Unchanged" [7] "Reduced" "Unchanged" "Unchanged" "Unchanged" converting counts to integer mode Warning message: In DESeqDataSet(se, design = design, ignoreRank) : some variables in design formula are characters, converting to factors estimating size factors estimating dispersions gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates fitting model and testing -- replacing outliers and refitting for 9 genes -- DESeq argument 'minReplicatesForReplace' = 7 -- original counts are preserved in counts(dds) estimating dispersions fitting model and testing using 'apeglm' for LFC shrinkage. If used in published research, please cite: Zhu, A., Ibrahim, J.G., Love, M.I. (2018) Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics. https://doi.org/10.1093/bioinformatics/bty895 log2 fold change (MAP): condition Expected vs el1 Wald test p-value: condition Expected vs el1 DataFrame with 6 rows and 5 columns baseMean log2FoldChange lfcSE pvalue padj 2 3422.79 11.527401 0.913260 9.45240e-36 1.81290e-34 3 3442.73 11.550559 0.905255 1.77790e-36 3.76508e-35 4 2756.74 2.495622 0.557320 3.06076e-07 1.02852e-06 5 3333.18 1.303268 0.548541 4.61414e-03 9.55249e-03 6 3043.79 0.621907 0.562536 1.73806e-01 2.36875e-01 7 3104.43 0.138064 0.437156 7.55901e-01 8.13524e-01 ---------------------------------------------------- Gaussian finite mixture model fitted by EM algorithm ---------------------------------------------------- Mclust V (univariate, unequal variance) model with 2 components: log-likelihood n df BIC ICL -4882.578 2033 5 -9803.243 -9940.514 Clustering table: 1 2 1414 619 Mixing probabilities: 1 2 0.6829948 0.3170052 Means: 1 2 -0.4209497 7.8608404 Variances: 1 2 1.157739 12.145192 ------------------------------------------------------- Density estimation via Gaussian finite mixture modeling ------------------------------------------------------- Mclust V (univariate, unequal variance) model with 2 components: log-likelihood n df BIC ICL -4882.543 2033 5 -9803.173 -9939.386 null device 1 [1] "Reduced" "Reduced" "Reduced" "Unchanged" "Unchanged" "Unchanged" [7] "Reduced" "Unchanged" "Unchanged" "Unchanged"