Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

returnSamples() error when setting 'counts = TRUE' for 'relative_abundance' dataType #306

Open
Guan06 opened this issue Sep 28, 2024 · 1 comment
Open

returnSamples() error when setting 'counts = TRUE' for 'relative_abundance' dataType #306

Guan06 opened this issue Sep 28, 2024 · 1 comment
Assignees
@schifferl
Labels
Bug Report

Comments

@Guan06
Copy link

Guan06 commented Sep 28, 2024

Describe the bug
I am tring to get the count table of relative abundance of some disease related samples, all good when loading with counts = FALSE but errors occur when setting counts to TRUE.

To Reproduce
I have version curatedMetagenomicData_3.8.0

Blow is the scripts that got errors:

study <- "IaniroG_2022"

  samples <- filter(sampleMetadata, 
                 study_name == study) |>
    filter(body_site=="stool") |>
    select(where(~ !all(is.na(.x))))
  
  print(study)
  filter_lst <- c("acute_diarrhoea", "CDI", "FMT", "IBD", "control")
  samples <- samples[(samples$study_condition %in% filter_lst) &
                       (samples$country != "CHN"), ]
  ra_2022 <- returnSamples(samples, counts = TRUE,
                      dataType = "relative_abundance", rownames = "long")

Expected behavior
Return of a count table.

Additional context
Error message:
Error in validObject(.Object) :
invalid class “SummarizedExperiment” object:
nb of cols in 'assay' (0) must equal nb of rows in 'colData' (165)
@Guan06 Guan06 changed the title returnSamples() error when setting 'counts = TRUE' for 'relative_abundance' aataType returnSamples() error when setting 'counts = TRUE' for 'relative_abundance' dataType Sep 28, 2024
@lwaldron
Copy link
Member

Thank you for the report, and I'm terribly sorry for the delay. @schifferl I can reproduce this bug, it seems to be an invalid replacement of the assay in a SummarizedExperiment. I added to the example that it doesn't happen with the ?returnSamples example, maybe because it merges multiples datasets instead of using just one.

library(curatedMetagenomicData)
library(dplyr)

study <- "IaniroG_2022"

samples <- dplyr::filter(sampleMetadata, 
                  study_name == study) |>
  filter(body_site=="stool") |>
  select(where(~ !all(is.na(.x))))

filter_lst <- c("acute_diarrhoea", "CDI", "FMT", "IBD", "control")

samples <- samples[(samples$study_condition %in% filter_lst) &
                     (samples$country != "CHN"), ]

returnSamples(samples, counts = TRUE, dataType = "relative_abundance")
#> 
#> $`2022-10-19.IaniroG_2022.relative_abundance`
#> dropping rows without rowTree matches:
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Atopobiaceae|g__Olsenella|s__Olsenella_profusa
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Collinsella|s__Collinsella_stercoris
#>   k__Bacteria|p__Firmicutes|c__Bacilli|o__Bacillales|f__Bacillales_unclassified|g__Gemella|s__Gemella_bergeri
#>   k__Bacteria|p__Firmicutes|c__Bacilli|o__Lactobacillales|f__Carnobacteriaceae|g__Granulicatella|s__Granulicatella_elegans
#>   k__Bacteria|p__Firmicutes|c__Clostridia|o__Clostridiales|f__Ruminococcaceae|g__Ruminococcus|s__Ruminococcus_champanellensis
#>   k__Bacteria|p__Firmicutes|c__Erysipelotrichia|o__Erysipelotrichales|f__Erysipelotrichaceae|g__Bulleidia|s__Bulleidia_extructa
#>   k__Bacteria|p__Proteobacteria|c__Betaproteobacteria|o__Burkholderiales|f__Sutterellaceae|g__Sutterella|s__Sutterella_parvirubra
#>   k__Bacteria|p__Synergistetes|c__Synergistia|o__Synergistales|f__Synergistaceae|g__Cloacibacillus|s__Cloacibacillus_evryensis
#> Error in validObject(.Object): invalid class "SummarizedExperiment" object: 
#>     nb of cols in 'assay' (0) must equal nb of rows in 'colData' (165)
returnSamples(samples, counts = FALSE, dataType = "relative_abundance")
#> 
#> $`2022-10-19.IaniroG_2022.relative_abundance`
#> dropping rows without rowTree matches:
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Atopobiaceae|g__Olsenella|s__Olsenella_profusa
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Collinsella|s__Collinsella_stercoris
#>   k__Bacteria|p__Firmicutes|c__Bacilli|o__Bacillales|f__Bacillales_unclassified|g__Gemella|s__Gemella_bergeri
#>   k__Bacteria|p__Firmicutes|c__Bacilli|o__Lactobacillales|f__Carnobacteriaceae|g__Granulicatella|s__Granulicatella_elegans
#>   k__Bacteria|p__Firmicutes|c__Clostridia|o__Clostridiales|f__Ruminococcaceae|g__Ruminococcus|s__Ruminococcus_champanellensis
#>   k__Bacteria|p__Firmicutes|c__Erysipelotrichia|o__Erysipelotrichales|f__Erysipelotrichaceae|g__Bulleidia|s__Bulleidia_extructa
#>   k__Bacteria|p__Proteobacteria|c__Betaproteobacteria|o__Burkholderiales|f__Sutterellaceae|g__Sutterella|s__Sutterella_parvirubra
#>   k__Bacteria|p__Synergistetes|c__Synergistia|o__Synergistales|f__Synergistaceae|g__Cloacibacillus|s__Cloacibacillus_evryensis
#> class: TreeSummarizedExperiment 
#> dim: 637 135 
#> metadata(0):
#> assays(1): relative_abundance
#> rownames(637):
#>   k__Bacteria|p__Firmicutes|c__Clostridia|o__Clostridiales|f__Clostridiaceae|g__Clostridium|s__Clostridium_sp_CAG_413
#>   k__Bacteria|p__Firmicutes|c__Clostridia|o__Clostridiales|f__Ruminococcaceae|g__Ruminococcus|s__Ruminococcus_bromii
#>   ...
#>   k__Bacteria|p__Firmicutes|c__Negativicutes|o__Selenomonadales|f__Selenomonadaceae|g__Centipeda|s__Centipeda_periodontii
#>   k__Bacteria|p__Bacteroidetes|c__Bacteroidia|o__Bacteroidales|f__Prevotellaceae|g__Prevotella|s__Prevotella_buccalis
#> rowData names(7): superkingdom phylum ... genus species
#> colnames(135): Sample_S25223_FMT_01_t7 Sample_S25235_FMT_09_t7 ...
#>   R5_290108_t7 R6_271006_t30
#> colData names(17): study_name subject_id ... subcohort fmt_id
#> reducedDimNames(0):
#> mainExpName: NULL
#> altExpNames(0):
#> rowLinks: a LinkDataFrame (637 rows)
#> rowTree: 1 phylo tree(s) (10430 leaves)
#> colLinks: NULL
#> colTree: NULL

# However the dataset from examples in ?returnSamples doesn't have the same problem:

sampleMetadata |>
  dplyr::filter(age >= 18) |>
  dplyr::filter(!base::is.na(alcohol)) |>
  dplyr::filter(body_site == "stool") |>
  dplyr::select(where(~ !base::all(base::is.na(.x)))) |>
  returnSamples(dataType="relative_abundance", counts = TRUE)
#> 
#> $`2021-10-14.KaurK_2020.relative_abundance`
#> dropping rows without rowTree matches:
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Atopobiaceae|g__Olsenella|s__Olsenella_profusa
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Collinsella|s__Collinsella_stercoris
#>   k__Bacteria|p__Firmicutes|c__Clostridia|o__Clostridiales|f__Ruminococcaceae|g__Ruminococcus|s__Ruminococcus_champanellensis
#>   k__Bacteria|p__Proteobacteria|c__Betaproteobacteria|o__Burkholderiales|f__Sutterellaceae|g__Sutterella|s__Sutterella_parvirubra
#> $`2021-03-31.KeohaneDM_2020.relative_abundance`
#> dropping rows without rowTree matches:
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Collinsella|s__Collinsella_stercoris
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Enorma|s__[Collinsella]_massiliensis
#>   k__Bacteria|p__Firmicutes|c__Clostridia|o__Clostridiales|f__Ruminococcaceae|g__Ruminococcus|s__Ruminococcus_champanellensis
#>   k__Bacteria|p__Proteobacteria|c__Betaproteobacteria|o__Burkholderiales|f__Sutterellaceae|g__Sutterella|s__Sutterella_parvirubra
#> $`2021-03-31.QinN_2014.relative_abundance`
#> dropping rows without rowTree matches:
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Atopobiaceae|g__Olsenella|s__Olsenella_profusa
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Collinsella|s__Collinsella_stercoris
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Enorma|s__[Collinsella]_massiliensis
#>   k__Bacteria|p__Firmicutes|c__Bacilli|o__Bacillales|f__Bacillales_unclassified|g__Gemella|s__Gemella_bergeri
#>   k__Bacteria|p__Firmicutes|c__Bacilli|o__Lactobacillales|f__Carnobacteriaceae|g__Granulicatella|s__Granulicatella_elegans
#>   k__Bacteria|p__Firmicutes|c__Clostridia|o__Clostridiales|f__Ruminococcaceae|g__Ruminococcus|s__Ruminococcus_champanellensis
#>   k__Bacteria|p__Firmicutes|c__Erysipelotrichia|o__Erysipelotrichales|f__Erysipelotrichaceae|g__Bulleidia|s__Bulleidia_extructa
#>   k__Bacteria|p__Proteobacteria|c__Betaproteobacteria|o__Burkholderiales|f__Sutterellaceae|g__Sutterella|s__Sutterella_parvirubra
#>   k__Bacteria|p__Synergistetes|c__Synergistia|o__Synergistales|f__Synergistaceae|g__Cloacibacillus|s__Cloacibacillus_evryensis
#> $`2021-03-31.ThomasAM_2018a.relative_abundance`
#> dropping rows without rowTree matches:
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Atopobiaceae|g__Olsenella|s__Olsenella_profusa
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Collinsella|s__Collinsella_stercoris
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Enorma|s__[Collinsella]_massiliensis
#>   k__Bacteria|p__Firmicutes|c__Bacilli|o__Bacillales|f__Bacillales_unclassified|g__Gemella|s__Gemella_bergeri
#>   k__Bacteria|p__Firmicutes|c__Bacilli|o__Lactobacillales|f__Carnobacteriaceae|g__Granulicatella|s__Granulicatella_elegans
#>   k__Bacteria|p__Firmicutes|c__Clostridia|o__Clostridiales|f__Ruminococcaceae|g__Ruminococcus|s__Ruminococcus_champanellensis
#>   k__Bacteria|p__Firmicutes|c__Erysipelotrichia|o__Erysipelotrichales|f__Erysipelotrichaceae|g__Bulleidia|s__Bulleidia_extructa
#>   k__Bacteria|p__Proteobacteria|c__Betaproteobacteria|o__Burkholderiales|f__Sutterellaceae|g__Sutterella|s__Sutterella_parvirubra
#>   k__Bacteria|p__Synergistetes|c__Synergistia|o__Synergistales|f__Synergistaceae|g__Cloacibacillus|s__Cloacibacillus_evryensis
#> $`2021-03-31.XieH_2016.relative_abundance`
#> dropping rows without rowTree matches:
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Atopobiaceae|g__Olsenella|s__Olsenella_profusa
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Collinsella|s__Collinsella_stercoris
#>   k__Bacteria|p__Actinobacteria|c__Coriobacteriia|o__Coriobacteriales|f__Coriobacteriaceae|g__Enorma|s__[Collinsella]_massiliensis
#>   k__Bacteria|p__Firmicutes|c__Bacilli|o__Lactobacillales|f__Carnobacteriaceae|g__Granulicatella|s__Granulicatella_elegans
#>   k__Bacteria|p__Firmicutes|c__Clostridia|o__Clostridiales|f__Ruminococcaceae|g__Ruminococcus|s__Ruminococcus_champanellensis
#>   k__Bacteria|p__Firmicutes|c__Erysipelotrichia|o__Erysipelotrichales|f__Erysipelotrichaceae|g__Bulleidia|s__Bulleidia_extructa
#>   k__Bacteria|p__Proteobacteria|c__Betaproteobacteria|o__Burkholderiales|f__Sutterellaceae|g__Sutterella|s__Sutterella_parvirubra
#>   k__Bacteria|p__Synergistetes|c__Synergistia|o__Synergistales|f__Synergistaceae|g__Cloacibacillus|s__Cloacibacillus_evryensis
#> Warning: There was 1 warning in `mutate()`.
#> ℹ In argument: `across(.fns = ~replace_na(.x, 0))`.
#> Caused by warning:
#> ! Using `across()` without supplying `.cols` was deprecated in dplyr 1.1.0.
#> ℹ Please supply `.cols` instead.
#> class: TreeSummarizedExperiment 
#> dim: 833 702 
#> metadata(0):
#> assays(1): relative_abundance
#> rownames(833):
#>   k__Bacteria|p__Bacteroidetes|c__Bacteroidia|o__Bacteroidales|f__Prevotellaceae|g__Prevotella|s__Prevotella_copri
#>   k__Bacteria|p__Bacteroidetes|c__Bacteroidia|o__Bacteroidales|f__Prevotellaceae|g__Prevotella|s__Prevotella_sp_CAG_520
#>   ...
#>   k__Bacteria|p__Actinobacteria|c__Actinobacteria|o__Corynebacteriales|f__Corynebacteriaceae|g__Corynebacterium|s__Corynebacterium_aurimucosum
#>   k__Bacteria|p__Actinobacteria|c__Actinobacteria|o__Corynebacteriales|f__Corynebacteriaceae|g__Corynebacterium|s__Corynebacterium_coyleae
#> rowData names(7): superkingdom phylum ... genus species
#> colnames(702): JAS_1 JAS_10 ... YSZC12003_37879 YSZC12003_37880
#> colData names(45): study_name subject_id ... inr zigosity
#> reducedDimNames(0):
#> mainExpName: NULL
#> altExpNames(0):
#> rowLinks: a LinkDataFrame (833 rows)
#> rowTree: 1 phylo tree(s) (10430 leaves)
#> colLinks: NULL
#> colTree: NULL

devtools::session_info()
#> ─ Session info ───────────────────────────────────────────────────────────────
#>  setting  value
#>  version  R version 4.4.2 (2024-10-31)
#>  os       Ubuntu 24.04.1 LTS
#>  system   x86_64, linux-gnu
#>  ui       X11
#>  language (EN)
#>  collate  en_US.UTF-8
#>  ctype    en_US.UTF-8
#>  tz       America/New_York
#>  date     2025-01-27
#>  pandoc   3.2 @ /usr/lib/rstudio-server/bin/quarto/bin/tools/x86_64/ (via rmarkdown)
#> 
#> ─ Packages ───────────────────────────────────────────────────────────────────
#>  package                  * version  date (UTC) lib source
#>  abind                      1.4-8    2024-09-12 [3] CRAN (R 4.4.1)
#>  AnnotationDbi              1.68.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  AnnotationHub              3.14.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  ape                        5.8-1    2024-12-16 [3] CRAN (R 4.4.2)
#>  backports                  1.5.0    2024-05-23 [3] CRAN (R 4.4.0)
#>  base64enc                  0.1-3    2015-07-28 [3] CRAN (R 4.4.0)
#>  beachmat                   2.22.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  beeswarm                   0.4.0    2021-06-01 [3] CRAN (R 4.4.0)
#>  Biobase                  * 2.66.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  BiocFileCache              2.14.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  BiocGenerics             * 0.52.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  BiocManager                1.30.25  2024-08-28 [3] CRAN (R 4.4.1)
#>  BiocNeighbors              2.0.1    2024-11-28 [3] Bioconductor 3.20 (R 4.4.2)
#>  BiocParallel               1.40.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  BiocSingular               1.22.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  BiocVersion                3.20.0   2024-05-01 [3] Bioconductor 3.20 (R 4.4.1)
#>  Biostrings               * 2.74.1   2024-12-16 [3] Bioconductor 3.20 (R 4.4.2)
#>  bit                        4.5.0.1  2024-12-03 [3] CRAN (R 4.4.2)
#>  bit64                      4.6.0-1  2025-01-16 [3] CRAN (R 4.4.2)
#>  blob                       1.2.4    2023-03-17 [3] CRAN (R 4.4.0)
#>  bluster                    1.16.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  boot                       1.3-31   2024-08-28 [4] CRAN (R 4.4.2)
#>  cachem                     1.1.0    2024-05-16 [3] CRAN (R 4.4.0)
#>  cellranger                 1.1.0    2016-07-27 [3] CRAN (R 4.4.0)
#>  checkmate                  2.3.2    2024-07-29 [3] CRAN (R 4.4.1)
#>  cli                        3.6.3    2024-06-21 [3] CRAN (R 4.4.1)
#>  cluster                    2.1.8    2024-12-11 [4] CRAN (R 4.4.2)
#>  coda                       0.19-4.1 2024-01-31 [3] RSPM (R 4.4.0)
#>  codetools                  0.2-20   2024-03-31 [2] CRAN (R 4.4.2)
#>  colorspace                 2.1-1    2024-07-26 [3] CRAN (R 4.4.1)
#>  crayon                     1.5.3    2024-06-20 [3] CRAN (R 4.4.1)
#>  curatedMetagenomicData   * 3.14.0   2024-10-31 [3] Bioconductor 3.20 (R 4.4.1)
#>  curl                       6.2.0    2025-01-23 [3] CRAN (R 4.4.2)
#>  data.table                 1.16.4   2024-12-06 [3] CRAN (R 4.4.2)
#>  DBI                        1.2.3    2024-06-02 [3] CRAN (R 4.4.0)
#>  dbplyr                     2.5.0    2024-03-19 [3] CRAN (R 4.4.0)
#>  DECIPHER                   3.2.0    2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  decontam                   1.26.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  DelayedArray               0.32.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  DelayedMatrixStats         1.28.1   2025-01-09 [3] Bioconductor 3.20 (R 4.4.2)
#>  devtools                   2.4.5    2022-10-11 [3] CRAN (R 4.4.0)
#>  digest                     0.6.37   2024-08-19 [3] CRAN (R 4.4.1)
#>  DirichletMultinomial       1.48.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  dplyr                    * 1.1.4    2023-11-17 [3] CRAN (R 4.4.0)
#>  ellipsis                   0.3.2    2021-04-29 [3] CRAN (R 4.4.0)
#>  emmeans                    1.10.6   2024-12-12 [3] CRAN (R 4.4.2)
#>  estimability               1.5.1    2024-05-12 [3] CRAN (R 4.4.2)
#>  evaluate                   1.0.3    2025-01-10 [3] CRAN (R 4.4.2)
#>  ExperimentHub              2.14.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  farver                     2.1.2    2024-05-13 [3] CRAN (R 4.4.0)
#>  fastmap                    1.2.0    2024-05-15 [3] CRAN (R 4.4.0)
#>  filelock                   1.0.3    2023-12-11 [3] CRAN (R 4.4.0)
#>  fillpattern                1.0.2    2024-06-24 [3] CRAN (R 4.4.2)
#>  foreign                    0.8-87   2024-06-26 [4] CRAN (R 4.4.2)
#>  Formula                    1.2-5    2023-02-24 [3] CRAN (R 4.4.0)
#>  fs                         1.6.5    2024-10-30 [3] CRAN (R 4.4.1)
#>  generics                   0.1.3    2022-07-05 [3] CRAN (R 4.4.0)
#>  GenomeInfoDb             * 1.42.1   2024-11-28 [3] Bioconductor 3.20 (R 4.4.2)
#>  GenomeInfoDbData           1.2.13   2024-10-12 [3] Bioconductor
#>  GenomicRanges            * 1.58.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  ggbeeswarm                 0.7.2    2023-04-29 [3] CRAN (R 4.4.0)
#>  ggnewscale                 0.5.0    2024-07-19 [3] CRAN (R 4.4.1)
#>  ggplot2                    3.5.1    2024-04-23 [3] RSPM (R 4.4.0)
#>  ggrepel                    0.9.6    2024-09-07 [3] CRAN (R 4.4.1)
#>  ggtext                     0.1.2    2022-09-16 [1] CRAN (R 4.4.2)
#>  glue                       1.8.0    2024-09-30 [3] CRAN (R 4.4.1)
#>  gridExtra                  2.3      2017-09-09 [3] CRAN (R 4.4.0)
#>  gridtext                   0.1.5    2022-09-16 [1] CRAN (R 4.4.2)
#>  gtable                     0.3.6    2024-10-25 [3] CRAN (R 4.4.1)
#>  Hmisc                      5.2-2    2025-01-10 [3] CRAN (R 4.4.2)
#>  hms                        1.1.3    2023-03-21 [3] CRAN (R 4.4.0)
#>  htmlTable                  2.4.3    2024-07-21 [3] CRAN (R 4.4.1)
#>  htmltools                  0.5.8.1  2024-04-04 [3] CRAN (R 4.4.0)
#>  htmlwidgets                1.6.4    2023-12-06 [3] CRAN (R 4.4.0)
#>  httpuv                     1.6.15   2024-03-26 [3] RSPM (R 4.4.0)
#>  httr                       1.4.7    2023-08-15 [3] CRAN (R 4.4.0)
#>  igraph                     2.1.4    2025-01-23 [3] CRAN (R 4.4.2)
#>  IRanges                  * 2.40.1   2024-12-05 [3] Bioconductor 3.20 (R 4.4.2)
#>  irlba                      2.3.5.1  2022-10-03 [3] CRAN (R 4.4.0)
#>  jsonlite                   1.8.9    2024-09-20 [3] CRAN (R 4.4.1)
#>  KEGGREST                   1.46.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  knitr                      1.49     2024-11-08 [3] CRAN (R 4.4.2)
#>  later                      1.4.1    2024-11-27 [3] CRAN (R 4.4.2)
#>  lattice                    0.22-6   2024-03-20 [2] CRAN (R 4.4.2)
#>  lazyeval                   0.2.2    2019-03-15 [3] CRAN (R 4.4.0)
#>  lifecycle                  1.0.4    2023-11-07 [3] CRAN (R 4.4.0)
#>  lme4                       1.1-36   2025-01-11 [3] CRAN (R 4.4.2)
#>  lpSolve                    5.6.23   2024-12-14 [3] CRAN (R 4.4.2)
#>  magrittr                   2.0.3    2022-03-30 [3] CRAN (R 4.4.0)
#>  MASS                       7.3-64   2025-01-04 [4] CRAN (R 4.4.2)
#>  Matrix                     1.7-2    2025-01-23 [3] CRAN (R 4.4.2)
#>  MatrixGenerics           * 1.18.1   2025-01-09 [3] Bioconductor 3.20 (R 4.4.2)
#>  matrixStats              * 1.5.0    2025-01-07 [3] CRAN (R 4.4.2)
#>  mediation                  4.5.0    2019-10-08 [3] CRAN (R 4.4.1)
#>  memoise                    2.0.1    2021-11-26 [3] CRAN (R 4.4.0)
#>  mgcv                       1.9-1    2023-12-21 [4] CRAN (R 4.3.2)
#>  mia                        1.14.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  mime                       0.12     2021-09-28 [3] CRAN (R 4.4.0)
#>  miniUI                     0.1.1.1  2018-05-18 [3] CRAN (R 4.4.0)
#>  minqa                      1.2.8    2024-08-17 [3] CRAN (R 4.4.1)
#>  multcomp                   1.4-26   2024-07-18 [3] CRAN (R 4.4.1)
#>  MultiAssayExperiment       1.32.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  munsell                    0.5.1    2024-04-01 [3] RSPM (R 4.4.0)
#>  mvtnorm                    1.3-3    2025-01-10 [3] CRAN (R 4.4.2)
#>  nlme                       3.1-166  2024-08-14 [4] CRAN (R 4.4.2)
#>  nloptr                     2.1.1    2024-06-25 [3] CRAN (R 4.4.1)
#>  nnet                       7.3-20   2025-01-01 [4] CRAN (R 4.4.2)
#>  parallelly                 1.41.0   2024-12-18 [3] CRAN (R 4.4.2)
#>  patchwork                  1.3.0    2024-09-16 [3] CRAN (R 4.4.1)
#>  permute                    0.9-7    2022-01-27 [3] CRAN (R 4.4.0)
#>  pillar                     1.10.1   2025-01-07 [3] CRAN (R 4.4.2)
#>  pkgbuild                   1.4.6    2025-01-16 [3] CRAN (R 4.4.2)
#>  pkgconfig                  2.0.3    2019-09-22 [3] CRAN (R 4.4.0)
#>  pkgload                    1.4.0    2024-06-28 [3] CRAN (R 4.4.1)
#>  plyr                       1.8.9    2023-10-02 [3] CRAN (R 4.4.0)
#>  png                        0.1-8    2022-11-29 [3] CRAN (R 4.4.0)
#>  profvis                    0.4.0    2024-09-20 [3] CRAN (R 4.4.1)
#>  promises                   1.3.2    2024-11-28 [3] CRAN (R 4.4.2)
#>  purrr                      1.0.2    2023-08-10 [3] CRAN (R 4.4.0)
#>  R6                         2.5.1    2021-08-19 [3] CRAN (R 4.4.0)
#>  ragg                       1.3.3    2024-09-11 [3] CRAN (R 4.4.1)
#>  rappdirs                   0.3.3    2021-01-31 [3] CRAN (R 4.4.0)
#>  rbibutils                  2.3      2024-10-04 [3] CRAN (R 4.4.1)
#>  rbiom                      2.0.13   2025-01-24 [3] CRAN (R 4.4.2)
#>  Rcpp                       1.0.14   2025-01-12 [3] CRAN (R 4.4.2)
#>  Rdpack                     2.6.2    2024-11-15 [3] CRAN (R 4.4.2)
#>  readr                      2.1.5    2024-01-10 [3] CRAN (R 4.4.0)
#>  readxl                     1.4.3    2023-07-06 [3] CRAN (R 4.4.0)
#>  reformulas                 0.4.0    2024-11-03 [3] CRAN (R 4.4.1)
#>  remotes                    2.5.0    2024-03-17 [3] CRAN (R 4.4.0)
#>  reprex                     2.1.1    2024-07-06 [3] CRAN (R 4.4.1)
#>  reshape2                   1.4.4    2020-04-09 [3] CRAN (R 4.4.0)
#>  rhdf5                      2.50.2   2025-01-09 [3] Bioconductor 3.20 (R 4.4.2)
#>  rhdf5filters               1.18.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  Rhdf5lib                   1.28.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  rlang                      1.1.5    2025-01-17 [3] CRAN (R 4.4.2)
#>  rmarkdown                  2.29     2024-11-04 [3] CRAN (R 4.4.1)
#>  rpart                      4.1.24   2025-01-07 [2] CRAN (R 4.4.2)
#>  RSQLite                    2.3.9    2024-12-03 [3] CRAN (R 4.4.2)
#>  rstudioapi                 0.17.1   2024-10-22 [3] CRAN (R 4.4.1)
#>  rsvd                       1.0.5    2021-04-16 [3] CRAN (R 4.4.0)
#>  S4Arrays                   1.6.0    2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  S4Vectors                * 0.44.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  sandwich                   3.1-1    2024-09-15 [3] CRAN (R 4.4.1)
#>  ScaledMatrix               1.14.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  scales                     1.3.0    2023-11-28 [3] RSPM (R 4.4.0)
#>  scater                     1.34.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  scuttle                    1.16.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  sessioninfo                1.2.2    2021-12-06 [3] CRAN (R 4.4.0)
#>  shiny                      1.10.0   2024-12-14 [3] CRAN (R 4.4.2)
#>  SingleCellExperiment     * 1.28.1   2024-11-10 [3] Bioconductor 3.20 (R 4.4.2)
#>  slam                       0.1-55   2024-11-13 [3] CRAN (R 4.4.2)
#>  SparseArray                1.6.1    2025-01-19 [3] Bioconductor 3.20 (R 4.4.2)
#>  sparseMatrixStats          1.18.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  stringi                    1.8.4    2024-05-06 [3] CRAN (R 4.4.0)
#>  stringr                    1.5.1    2023-11-14 [3] RSPM (R 4.4.0)
#>  SummarizedExperiment     * 1.36.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  survival                   3.8-3    2024-12-17 [4] CRAN (R 4.4.2)
#>  systemfonts                1.2.1    2025-01-20 [3] CRAN (R 4.4.2)
#>  textshaping                1.0.0    2025-01-20 [3] CRAN (R 4.4.2)
#>  TH.data                    1.1-3    2025-01-17 [3] CRAN (R 4.4.2)
#>  tibble                     3.2.1    2023-03-20 [3] CRAN (R 4.4.0)
#>  tidyr                      1.3.1    2024-01-24 [3] CRAN (R 4.4.0)
#>  tidyselect                 1.2.1    2024-03-11 [3] CRAN (R 4.4.0)
#>  tidytree                   0.4.6    2023-12-12 [3] CRAN (R 4.4.0)
#>  treeio                     1.30.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  TreeSummarizedExperiment * 2.14.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  tzdb                       0.4.0    2023-05-12 [3] CRAN (R 4.4.0)
#>  UCSC.utils                 1.2.0    2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  urlchecker                 1.0.1    2021-11-30 [3] CRAN (R 4.4.0)
#>  usethis                    3.1.0    2024-11-26 [3] CRAN (R 4.4.2)
#>  vctrs                      0.6.5    2023-12-01 [3] CRAN (R 4.4.0)
#>  vegan                      2.6-8    2024-08-28 [3] CRAN (R 4.4.1)
#>  vipor                      0.4.7    2023-12-18 [3] CRAN (R 4.4.0)
#>  viridis                    0.6.5    2024-01-29 [3] CRAN (R 4.4.0)
#>  viridisLite                0.4.2    2023-05-02 [3] RSPM (R 4.4.0)
#>  withr                      3.0.2    2024-10-28 [3] CRAN (R 4.4.1)
#>  xfun                       0.50     2025-01-07 [3] CRAN (R 4.4.2)
#>  xml2                       1.3.6    2023-12-04 [3] CRAN (R 4.4.0)
#>  xtable                     1.8-4    2019-04-21 [3] CRAN (R 4.4.0)
#>  XVector                  * 0.46.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  yaml                       2.3.10   2024-07-26 [3] CRAN (R 4.4.1)
#>  yulab.utils                0.1.9    2025-01-07 [3] CRAN (R 4.4.2)
#>  zlibbioc                   1.52.0   2024-10-29 [3] Bioconductor 3.20 (R 4.4.1)
#>  zoo                        1.8-12   2023-04-13 [3] CRAN (R 4.4.0)
#> 
#>  [1] /mnt/STORE1/bighome/levi/R/x86_64-pc-linux-gnu-library/4.4
#>  [2] /usr/local/lib/R/site-library
#>  [3] /usr/lib/R/site-library
#>  [4] /usr/lib/R/library
#> 
#> ──────────────────────────────────────────────────────────────────────────────

Created on 2025-01-27 with reprex v2.1.1

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees

@schifferl schifferl

Labels
Bug Report
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
3 participants
@lwaldron @schifferl @Guan06