Maaslin2 1.13.0
I was running the ready-made example in the function Maaslin2
help page:
input_data <- system.file(
'extdata','HMP2_taxonomy.tsv', package="Maaslin2")
input_metadata <-system.file(
'extdata','HMP2_metadata.tsv', package="Maaslin2")
fit_data <- Maaslin2(
input_data, input_metadata,'demo_output', transform = "AST",
fixed_effects = c('diagnosis', 'dysbiosisnonIBD','dysbiosisUC','dysbiosisCD', 'antibiotics', 'age'),
random_effects = c('site', 'subject'),
normalization = 'NONE',
reference = 'diagnosis,nonIBD',
standardize = FALSE)
This leads to the following error:
....
2023-04-28 10:41:54.896843 INFO::Writing heatmap of significant results to file: demo_output/heatmap.pdf
Error in xtfrm.data.frame(x) : cannot xtfrm data frames
In addition: Warning messages:
1: Model failed to converge with 1 negative eigenvalue: -5.6e+00
2: Model failed to converge with 1 negative eigenvalue: -1.1e+01
3: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 0.00291214 (tol = 0.002, component 1)
4: Model failed to converge with 1 negative eigenvalue: -2.1e+02
5: Model failed to converge with 1 negative eigenvalue: -2.2e+02
Information my R session:
> sessionInfo()
R version 4.3.0 (2023-04-21)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 22.04.2 LTS
Matrix products: default
BLAS: /home/xxx/bin/R-4.3.0/lib/libRblas.so
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.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/Mariehamn
tzcode source: system (glibc)
attached base packages:
[1] stats4 stats graphics grDevices utils datasets methods
[8] base
other attached packages:
[1] doRNG_1.8.6 rngtools_1.5.2
[3] foreach_1.5.2 ANCOMBC_2.2.0
[5] lubridate_1.9.2 forcats_1.0.0
[7] stringr_1.5.0 dplyr_1.1.2
[9] purrr_1.0.1 readr_2.1.4
[11] tidyr_1.3.0 tibble_3.2.1
[13] ggplot2_3.4.2 tidyverse_2.0.0
[15] knitr_1.42 MicrobiomeStat_1.1
[17] Maaslin2_1.13.0 ALDEx2_1.32.0
[19] zCompositions_1.4.0-1 truncnorm_1.0-9
[21] NADA_1.6-1.1 survival_3.5-5
[23] MASS_7.3-59 tidySummarizedExperiment_1.10.0
[25] patchwork_1.1.2.9000 mia_1.8.0
[27] MultiAssayExperiment_1.26.0 TreeSummarizedExperiment_2.8.0
[29] Biostrings_2.68.0 XVector_0.40.0
[31] SingleCellExperiment_1.22.0 SummarizedExperiment_1.30.0
[33] Biobase_2.60.0 GenomicRanges_1.52.0
[35] GenomeInfoDb_1.36.0 IRanges_2.34.0
[37] S4Vectors_0.38.0 BiocGenerics_0.46.0
[39] MatrixGenerics_1.12.0 matrixStats_0.63.0
[41] BiocStyle_2.28.0 rebook_1.9.0
loaded via a namespace (and not attached):
[1] bitops_1.0-7 DirichletMultinomial_1.42.0
[3] doParallel_1.0.17 httr_1.4.5
[5] numDeriv_2016.8-1.1 backports_1.4.1
[7] tools_4.3.0 utf8_1.2.3
[9] R6_2.5.1 vegan_2.6-4
[11] lazyeval_0.2.2 mgcv_1.8-42
[13] rhdf5filters_1.12.0 permute_0.9-7
[15] withr_2.5.0 gridExtra_2.3
[17] cli_3.6.1.9000 logging_0.10-108
[19] biglm_0.9-2.1 sandwich_3.0-2
[21] mvtnorm_1.1-3 robustbase_0.95-1
[23] pbapply_1.7-0 proxy_0.4-27
[25] yulab.utils_0.0.6 foreign_0.8-84
[27] scater_1.28.0 decontam_1.20.0
[29] readxl_1.4.2 rstudioapi_0.14
[31] RSQLite_2.3.1 generics_0.1.3
[33] Matrix_1.5-4 biomformat_1.28.0
[35] ggbeeswarm_0.7.1 fansi_1.0.4
[37] DescTools_0.99.48 DECIPHER_2.28.0
[39] lifecycle_1.0.3 multcomp_1.4-23
[41] yaml_2.3.7 rhdf5_2.44.0
[43] grid_4.3.0 blob_1.2.4
[45] crayon_1.5.2 dir.expiry_1.8.0
[47] lattice_0.21-8 beachmat_2.16.0
[49] CodeDepends_0.6.5 pillar_1.9.0
[51] optparse_1.7.3 statip_0.2.3
[53] boot_1.3-28.1 gld_2.6.6
[55] estimability_1.4.1 codetools_0.2-19
[57] glue_1.6.2 data.table_1.14.8
[59] Rdpack_2.4 vctrs_0.6.2
[61] treeio_1.24.0 cellranger_1.1.0
[63] gtable_0.3.3 cachem_1.0.7
[65] xfun_0.39 rbibutils_2.2.13
[67] Rfast_2.0.7 coda_0.19-4
[69] pcaPP_2.0-3 modeest_2.4.0
[71] timeDate_4022.108 iterators_1.0.14
[73] statmod_1.5.0 gmp_0.7-1
[75] TH.data_1.1-2 ellipsis_0.3.2
[77] nlme_3.1-162 phyloseq_1.44.0
[79] bit64_4.0.5 filelock_1.0.2
[81] fBasics_4022.94 irlba_2.3.5.1
[83] vipor_0.4.5 rpart_4.1.19
[85] colorspace_2.1-0 DBI_1.1.3
[87] Hmisc_5.0-1 nnet_7.3-18
[89] ade4_1.7-22 Exact_3.2
[91] tidyselect_1.2.0 emmeans_1.8.5
[93] timeSeries_4021.105 bit_4.0.5
[95] compiler_4.3.0 graph_1.78.0
[97] htmlTable_2.4.1 BiocNeighbors_1.18.0
[99] expm_0.999-7 DelayedArray_0.25.0
[101] plotly_4.10.1 checkmate_2.2.0
[103] scales_1.2.1 DEoptimR_1.0-12
[105] spatial_7.3-16 digest_0.6.31
[107] minqa_1.2.5 rmarkdown_2.21.3
[109] base64enc_0.1-3 htmltools_0.5.5
[111] pkgconfig_2.0.3 lme4_1.1-33
[113] sparseMatrixStats_1.12.0 lpsymphony_1.28.0
[115] stabledist_0.7-1 fastmap_1.1.1
[117] rlang_1.1.0 htmlwidgets_1.6.2
[119] DelayedMatrixStats_1.22.0 energy_1.7-11
[121] zoo_1.8-12 jsonlite_1.8.4
[123] BiocParallel_1.34.0 BiocSingular_1.16.0
[125] RCurl_1.98-1.12 magrittr_2.0.3
[127] Formula_1.2-5 scuttle_1.10.0
[129] GenomeInfoDbData_1.2.10 Rhdf5lib_1.22.0
[131] munsell_0.5.0 Rcpp_1.0.10
[133] ape_5.7-1 viridis_0.6.2
[135] RcppZiggurat_0.1.6 CVXR_1.0-11
[137] stringi_1.7.12 rootSolve_1.8.2.3
[139] stable_1.1.6 zlibbioc_1.46.0
[141] plyr_1.8.8 parallel_4.3.0
[143] ggrepel_0.9.3 lmom_2.9
[145] splines_4.3.0 hash_2.2.6.2
[147] multtest_2.56.0 hms_1.1.3
[149] igraph_1.4.2 reshape2_1.4.4
[151] ScaledMatrix_1.7.1 rmutil_1.1.10
[153] XML_3.99-0.14 evaluate_0.20
[155] BiocManager_1.30.20 nloptr_2.0.3
[157] tzdb_0.3.0 getopt_1.20.3
[159] clue_0.3-64 rsvd_1.0.5
[161] xtable_1.8-4 Rmpfr_0.9-2
[163] e1071_1.7-13 tidytree_0.4.2
[165] viridisLite_0.4.1 class_7.3-21
[167] gsl_2.1-8 lmerTest_3.1-3
[169] memoise_2.0.1 beeswarm_0.4.0
[171] cluster_2.1.4 timechange_0.2.0