Our OmicsMLRepo project aims to improve the AI/ML-readiness of Omics datasets available through Bioconductor. One of the main activities under this project is metadata harmonization (e.g., remove redundant information) and standardization (i.e., incorporate ontology).
Currently, we released the harmonized version of metadata for two Bioconductor data packages - curatedMetagenomicDatacontaining human microbiome data and cBioPortalData package on cancer genomics data. OmicsMLRepoR is a software package allowing users to easily access the harmonized metadata and to leverage ontology in metadata search.
OmicsMLRepoR package provides the three major functions:
1. Download the harmonized metadata
2. Browse the harmonized metadata using ontology
3. Manipulate the ‘shape’ of the harmonized metadata
You can download the harmonized version of metadata using the
getMetadata
function. Currently, two options are available
- cMD
and cBioPortalData
.
cmd <- getMetadata("cMD")
cmd
#> # A tibble: 22,668 × 87
#> study_name subject_id sample_id curation_id target_condition
#> * <chr> <chr> <chr> <chr> <chr>
#> 1 AsnicarF_2017 MV_FEI1 MV_FEI1_t1Q14 AsnicarF_2017:MV_FEI… human gut micro…
#> 2 AsnicarF_2017 MV_FEI2 MV_FEI2_t1Q14 AsnicarF_2017:MV_FEI… human gut micro…
#> 3 AsnicarF_2017 MV_FEI3 MV_FEI3_t1Q14 AsnicarF_2017:MV_FEI… human gut micro…
#> 4 AsnicarF_2017 MV_FEI4 MV_FEI4_t1Q14 AsnicarF_2017:MV_FEI… human gut micro…
#> 5 AsnicarF_2017 MV_FEI4 MV_FEI4_t2Q15 AsnicarF_2017:MV_FEI… human gut micro…
#> 6 AsnicarF_2017 MV_FEI5 MV_FEI5_t1Q14 AsnicarF_2017:MV_FEI… human gut micro…
#> 7 AsnicarF_2017 MV_FEI5 MV_FEI5_t2Q14 AsnicarF_2017:MV_FEI… human gut micro…
#> 8 AsnicarF_2017 MV_FEI5 MV_FEI5_t3Q15 AsnicarF_2017:MV_FEI… human gut micro…
#> 9 AsnicarF_2017 MV_FEM1 MV_FEM1_t1Q14 AsnicarF_2017:MV_FEM… human gut micro…
#> 10 AsnicarF_2017 MV_FEM2 MV_FEM2_t1Q14 AsnicarF_2017:MV_FEM… human gut micro…
#> # ℹ 22,658 more rows
#> # ℹ 82 more variables: target_condition_ontology_term_id <chr>, control <chr>,
#> # control_ontology_term_id <chr>, country <chr>,
#> # country_ontology_term_id <chr>, body_site <chr>,
#> # body_site_ontology_term_id <chr>, body_site_details <chr>,
#> # body_site_details_ontology_term_id <chr>, age_group <chr>,
#> # age_group_ontology_term_id <chr>, age_max <dbl>, age_min <dbl>, …
cbio <- getMetadata("cBioPortal")
cbio
#> # A tibble: 189,439 × 79
#> studyId patientId sampleId curation_id SAMPLE_COUNT acronym
#> * <chr> <chr> <chr> <chr> <dbl> <chr>
#> 1 acc_tcga TCGA-OR-A5J1 TCGA-OR-A5J1-01 acc_tcga:TCGA-OR-… 1 <NA>
#> 2 acc_tcga TCGA-OR-A5J2 TCGA-OR-A5J2-01 acc_tcga:TCGA-OR-… 1 <NA>
#> 3 acc_tcga TCGA-OR-A5J3 TCGA-OR-A5J3-01 acc_tcga:TCGA-OR-… 1 <NA>
#> 4 acc_tcga TCGA-OR-A5J4 TCGA-OR-A5J4-01 acc_tcga:TCGA-OR-… 1 <NA>
#> 5 acc_tcga TCGA-OR-A5J5 TCGA-OR-A5J5-01 acc_tcga:TCGA-OR-… 1 <NA>
#> 6 acc_tcga TCGA-OR-A5J6 TCGA-OR-A5J6-01 acc_tcga:TCGA-OR-… 1 <NA>
#> 7 acc_tcga TCGA-OR-A5J7 TCGA-OR-A5J7-01 acc_tcga:TCGA-OR-… 1 <NA>
#> 8 acc_tcga TCGA-OR-A5J8 TCGA-OR-A5J8-01 acc_tcga:TCGA-OR-… 1 <NA>
#> 9 acc_tcga TCGA-OR-A5J9 TCGA-OR-A5J9-01 acc_tcga:TCGA-OR-… 1 <NA>
#> 10 acc_tcga TCGA-OR-A5JA TCGA-OR-A5JA-01 acc_tcga:TCGA-OR-… 1 <NA>
#> # ℹ 189,429 more rows
#> # ℹ 73 more variables: acronym_ontology_term_id <chr>, age_at_death <dbl>,
#> # age_at_death_max <dbl>, age_at_death_min <dbl>, age_at_diagnosis <chr>,
#> # age_at_diagnosis_max <chr>, age_at_diagnosis_min <chr>,
#> # age_at_metastasis <dbl>, age_at_metastasis_max <dbl>,
#> # age_at_metastasis_min <dbl>, age_at_procurement <dbl>,
#> # age_at_procurement_max <dbl>, age_at_procurement_min <dbl>, …
Harmonized metadata can be easily searched by dplyr
functions. To fully leverage ontologies incorporated in harmonized
metadata and provide more robust data browsing experience, the package
provides the tree_filter
function. Note, that
tree_filter
can be used on the attributes mapped to the
ontology terms:
#> [1] "target_condition" "control"
#> [3] "country" "body_site"
#> [5] "body_site_details" "age_group"
#> [7] "ancestry_details" "ancestry"
#> [9] "disease_details" "disease"
#> [11] "disease_response_recist" "feces_phenotype"
#> [13] "hla" "neonatal_delivery_procedure"
#> [15] "neonatal_preterm_birth" "obgyn_menopause"
#> [17] "obgyn_pregnancy" "probing_pocket_depth"
#> [19] "sex" "smoker"
#> [21] "treatment"
Compared to the typical searching in the original metadata from the
curatedMetagenomicData
, OmicsMLRepoR enables more robust
data browsing, including case-insensitive, synonyms and descendant
searching capabilities.
Searching the same information in the original, unharmonized metadata
(sampleMetadata
) from the
curatedMetagenomicData
package is much less robust:
## Information spread out in two different columns
nrow(sampleMetadata |> filter(study_condition == "CRC"))
#> [1] 701
nrow(sampleMetadata |> filter(disease == "CRC"))
#> [1] 625
## Case sensitive
nrow(sampleMetadata |> filter(study_condition == "CRC"))
#> [1] 701
nrow(sampleMetadata |> filter(study_condition == "crc"))
#> [1] 0
## Synonyms not covered
nrow(sampleMetadata |> filter(study_condition == "Colorectal Carcinoma"))
#> [1] 0
nrow(sampleMetadata |> filter(study_condition == "Colorectal Cancer"))
#> [1] 0
tree_filter
is not case-sensitive.
tree_filter
includes the synonyms of the queried terms
in its searching.
syn_res1 <- cmd |> tree_filter(disease, "CRC")
syn_res2 <- cmd |> tree_filter(disease, "Colorectal Cancer")
syn_res3 <- cmd |> tree_filter(disease, "Colorectal Carcinoma")
nrow(syn_res1)
#> [1] 701
nrow(syn_res2)
#> [1] 701
nrow(syn_res3)
#> [1] 701
Check that the returned results are identical.
unique(syn_res1$disease)
#> [1] "Colorectal Carcinoma;Hepatic Steatosis;Hypertension;Carcinoma"
#> [2] "Colorectal Carcinoma;Carcinoma"
#> [3] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Hepatic Steatosis;Hypertension;Carcinoma"
#> [4] "Colorectal Carcinoma;Hypertension;Carcinoma"
#> [5] "Colorectal Carcinoma;Hepatic Steatosis;Carcinoma"
#> [6] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Hypertension;Carcinoma"
#> [7] "Colorectal Carcinoma;Adenocarcinoma"
#> [8] "Colorectal Carcinoma"
#> [9] "Colorectal Carcinoma;Hypercholesterolemia;Adenocarcinoma"
#> [10] "Colorectal Carcinoma;Hypertension;Adenocarcinoma"
#> [11] "Colorectal Carcinoma;Hypercholesterolemia;Hypertension;Adenocarcinoma"
#> [12] "Colorectal Carcinoma;Metastatic Malignant Neoplasm;Adenocarcinoma"
#> [13] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Carcinoma"
unique(syn_res2$disease)
#> [1] "Colorectal Carcinoma;Hepatic Steatosis;Hypertension;Carcinoma"
#> [2] "Colorectal Carcinoma;Carcinoma"
#> [3] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Hepatic Steatosis;Hypertension;Carcinoma"
#> [4] "Colorectal Carcinoma;Hypertension;Carcinoma"
#> [5] "Colorectal Carcinoma;Hepatic Steatosis;Carcinoma"
#> [6] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Hypertension;Carcinoma"
#> [7] "Colorectal Carcinoma;Adenocarcinoma"
#> [8] "Colorectal Carcinoma"
#> [9] "Colorectal Carcinoma;Hypercholesterolemia;Adenocarcinoma"
#> [10] "Colorectal Carcinoma;Hypertension;Adenocarcinoma"
#> [11] "Colorectal Carcinoma;Hypercholesterolemia;Hypertension;Adenocarcinoma"
#> [12] "Colorectal Carcinoma;Metastatic Malignant Neoplasm;Adenocarcinoma"
#> [13] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Carcinoma"
unique(syn_res3$disease)
#> [1] "Colorectal Carcinoma;Hepatic Steatosis;Hypertension;Carcinoma"
#> [2] "Colorectal Carcinoma;Carcinoma"
#> [3] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Hepatic Steatosis;Hypertension;Carcinoma"
#> [4] "Colorectal Carcinoma;Hypertension;Carcinoma"
#> [5] "Colorectal Carcinoma;Hepatic Steatosis;Carcinoma"
#> [6] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Hypertension;Carcinoma"
#> [7] "Colorectal Carcinoma;Adenocarcinoma"
#> [8] "Colorectal Carcinoma"
#> [9] "Colorectal Carcinoma;Hypercholesterolemia;Adenocarcinoma"
#> [10] "Colorectal Carcinoma;Hypertension;Adenocarcinoma"
#> [11] "Colorectal Carcinoma;Hypercholesterolemia;Hypertension;Adenocarcinoma"
#> [12] "Colorectal Carcinoma;Metastatic Malignant Neoplasm;Adenocarcinoma"
#> [13] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Carcinoma"
tree_filter
includes all the descendants of the queried
term in its searching.
onto_res <- cmd |> tree_filter(disease, "Intestinal Disorder")
unique(onto_res$disease)
#> [1] "Crohn Disease;Schizophrenia"
#> [2] "Colorectal Carcinoma;Hepatic Steatosis;Hypertension;Carcinoma"
#> [3] "Colorectal Carcinoma;Carcinoma"
#> [4] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Hepatic Steatosis;Hypertension;Carcinoma"
#> [5] "Colorectal Carcinoma;Hypertension;Carcinoma"
#> [6] "Colorectal Carcinoma;Hepatic Steatosis;Carcinoma"
#> [7] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Hypertension;Carcinoma"
#> [8] "Colorectal Carcinoma;Adenocarcinoma"
#> [9] "Inflammatory Bowel Disease;Crohn Disease"
#> [10] "Inflammatory Bowel Disease;Ulcerative Colitis"
#> [11] "Colorectal Carcinoma"
#> [12] "Cytomegaloviral Infection;Celiac Disease;Gestational Diabetes"
#> [13] "Type 1 Diabetes Mellitus;Celiac Disease;Irritable Bowel Syndrome"
#> [14] "Inflammatory Bowel Disease"
#> [15] "Inflammatory Bowel Disease;Fecal Microbiota Transplantation"
#> [16] "Melanoma;Colitis"
#> [17] "Inflammatory Bowel Disease;Anorectal Fistula;Crohn Disease"
#> [18] "Colorectal Carcinoma;Hypercholesterolemia;Adenocarcinoma"
#> [19] "Colorectal Carcinoma;Hypertension;Adenocarcinoma"
#> [20] "Colorectal Carcinoma;Hypercholesterolemia;Hypertension;Adenocarcinoma"
#> [21] "Colorectal Carcinoma;Metastatic Malignant Neoplasm;Adenocarcinoma"
#> [22] "Inflammatory Bowel Disease;Colitis"
#> [23] "Colorectal Carcinoma;Type 2 Diabetes Mellitus;Carcinoma"
#> [24] "Adenoma;Small Intestinal Adenoma"
#> [25] "Adenoma;Colorectal Adenoma"
For example, you can search for any row including a disease related to either “migraine” or “diabetes.”
We can also change the “OR” argument (default) to either “AND” or “NOT” and change the filtering action. “AND” will return any rows including a disease value that is related to both “migraine” and “diabetes,” and “NOT” will return any rows including a disease value that is not related to either “migraine” or “diabetes.”
res_and <- cmd %>% tree_filter(disease, c("migraine", "diabetes"), "AND")
res_not <- cmd %>% tree_filter(disease, c("migraine", "diabetes"), "NOT")
You can combine tree_filter
and dplyr
functions. For example, if you want all rows with a disease value
related to either “migraine” or “diabetes,” as well as with an age_years
value under 30,
Some metadata columns (e.g., biomarker
) contain
multiple, similar attributes separated with a specific delimiter (i.e.,
<;>
). Our harmonization use this structure because
they are related information often looked up together.
cmd_biomarker <- cmd %>%
filter(!is.na(biomarker)) %>%
select(curation_id, biomarker)
wtb <- getWideMetaTb(cmd_biomarker, "biomarker")
head(wtb)
#> # A tibble: 6 × 36
#> curation_id `Adiponectin_in_ug/mL` Alanine_Aminotransfe…¹ `Albumin_in_g/dL`
#> <chr> <chr> <chr> <chr>
#> 1 ChengpingW_20… <NA> 23 57
#> 2 ChengpingW_20… <NA> 13 49.4
#> 3 ChengpingW_20… <NA> 19 50.8
#> 4 ChengpingW_20… <NA> 29 52
#> 5 ChengpingW_20… <NA> 14 47.1
#> 6 ChengpingW_20… <NA> 5 46.1
#> # ℹ abbreviated name: ¹`Alanine_Aminotransferase_in_U/L`
#> # ℹ 32 more variables: `Aspartate_Aminotransferase_in_U/L` <chr>,
#> # `Autoantibody_titer_positive_(finding)` <chr>, `C-peptide_in_ng/ml` <chr>,
#> # `Cholesterol_in_mg/dL` <chr>, `Creatine_in_umol/L` <chr>,
#> # `Creatinine_in_umol/L` <chr>, `Diastolic_Blood_Pressure_in_mm/Hg` <chr>,
#> # `Direct_Bilirubin_in_umol/L` <chr>,
#> # `Erythrocyte_Sedimentation_Rate_in_mm/hr` <chr>, …
cbio_sub <- cbio %>%
getLongMetaTb("treatment_name", "<;>") %>%
filter(treatment_name == "Fluorouracil") %>%
filter(age_at_diagnosis > 50) %>%
filter(sex == "Female") %>%
getShortMetaTb(idCols = "curation_id", targetCols = "treatment_name")
dim(cbio_sub)
#> [1] 56 79
studies <- unique(cbio_sub$studyId)
studies
#> [1] "aml_ohsu_2018" "aml_ohsu_2022" "egc_msk_2017"
A simple for
loop can collect samples from multiple
studies.
cbio_api <- cBioPortal()
resAll <- as.list(vector(length = length(studies)))
for (i in seq_along(studies)) {
study <- studies[i]
samples <- cbio_sub %>%
filter(studyId == study) %>%
pull(sampleId)
res <- cBioPortalData(
api = cbio_api,
by = "hugoGeneSymbol",
studyId = study,
sampleIds = samples,
genePanelId = "IMPACT341"
)
resAll[[i]] <- res
}
sessionInfo()
#> R version 4.4.1 (2024-06-14)
#> Platform: x86_64-pc-linux-gnu
#> Running under: Ubuntu 24.04.1 LTS
#>
#> Matrix products: default
#> BLAS: /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3
#> LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.26.so; LAPACK version 3.12.0
#>
#> locale:
#> [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
#> [3] LC_TIME=en_US.UTF-8 LC_COLLATE=C
#> [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: Etc/UTC
#> tzcode source: system (glibc)
#>
#> attached base packages:
#> [1] stats4 stats graphics grDevices utils datasets methods
#> [8] base
#>
#> other attached packages:
#> [1] cBioPortalData_2.19.0 MultiAssayExperiment_1.31.5
#> [3] AnVIL_1.19.0 AnVILBase_1.1.0
#> [5] curatedMetagenomicData_3.13.0 TreeSummarizedExperiment_2.13.0
#> [7] Biostrings_2.75.0 XVector_0.45.0
#> [9] SingleCellExperiment_1.27.2 SummarizedExperiment_1.35.5
#> [11] Biobase_2.67.0 GenomicRanges_1.57.2
#> [13] GenomeInfoDb_1.41.2 IRanges_2.39.2
#> [15] S4Vectors_0.43.2 BiocGenerics_0.53.0
#> [17] MatrixGenerics_1.17.1 matrixStats_1.4.1
#> [19] dplyr_1.1.4 OmicsMLRepoR_1.1.0
#> [21] BiocStyle_2.35.0
#>
#> loaded via a namespace (and not attached):
#> [1] fs_1.6.4 bitops_1.0-9
#> [3] DirichletMultinomial_1.49.0 lubridate_1.9.3
#> [5] httr_1.4.7 RColorBrewer_1.1-3
#> [7] GenomicDataCommons_1.29.7 tools_4.4.1
#> [9] backports_1.5.0 utf8_1.2.4
#> [11] R6_2.5.1 DT_0.33
#> [13] vegan_2.6-8 lazyeval_0.2.2
#> [15] mgcv_1.9-1 permute_0.9-7
#> [17] withr_3.0.2 TCGAutils_1.25.1
#> [19] gridExtra_2.3 cli_3.6.3
#> [21] formatR_1.14 sandwich_3.1-1
#> [23] slam_0.1-54 sass_0.4.9
#> [25] mvtnorm_1.3-1 readr_2.1.5
#> [27] rapiclient_0.1.8 Rsamtools_2.21.2
#> [29] yulab.utils_0.1.7 foreign_0.8-87
#> [31] scater_1.33.4 decontam_1.27.0
#> [33] rstudioapi_0.17.1 RSQLite_2.3.7
#> [35] visNetwork_2.1.2 generics_0.1.3
#> [37] BiocIO_1.17.0 vroom_1.6.5
#> [39] rbiom_1.0.3 Matrix_1.7-1
#> [41] futile.logger_1.4.3 ggbeeswarm_0.7.2
#> [43] fansi_1.0.6 DECIPHER_3.3.0
#> [45] abind_1.4-8 lifecycle_1.0.4
#> [47] yaml_2.3.10 RaggedExperiment_1.29.5
#> [49] SparseArray_1.5.45 BiocFileCache_2.15.0
#> [51] grid_4.4.1 blob_1.2.4
#> [53] promises_1.3.0 ExperimentHub_2.13.1
#> [55] crayon_1.5.3 miniUI_0.1.1.1
#> [57] lattice_0.22-6 beachmat_2.23.0
#> [59] GenomicFeatures_1.57.1 KEGGREST_1.45.1
#> [61] sys_3.4.3 maketools_1.3.1
#> [63] pillar_1.9.0 knitr_1.48
#> [65] rjson_0.2.23 boot_1.3-31
#> [67] lpSolve_5.6.21 codetools_0.2-20
#> [69] glue_1.8.0 data.table_1.16.2
#> [71] vctrs_0.6.5 png_0.1-8
#> [73] treeio_1.29.2 gtable_0.3.6
#> [75] cachem_1.1.0 xfun_0.48
#> [77] S4Arrays_1.5.11 mime_0.12
#> [79] RTCGAToolbox_2.35.1 DiagrammeR_1.0.11
#> [81] bluster_1.17.0 nlme_3.1-166
#> [83] bit64_4.5.2 filelock_1.0.3
#> [85] data.tree_1.1.0 bslib_0.8.0
#> [87] irlba_2.3.5.1 vipor_0.4.7
#> [89] rpart_4.1.23 colorspace_2.1-1
#> [91] DBI_1.2.3 Hmisc_5.2-0
#> [93] nnet_7.3-19 tidyselect_1.2.1
#> [95] bit_4.5.0 compiler_4.4.1
#> [97] curl_5.2.3 rvest_1.0.4
#> [99] httr2_1.0.5 htmlTable_2.4.3
#> [101] BiocNeighbors_2.1.0 xml2_1.3.6
#> [103] DelayedArray_0.33.1 rtracklayer_1.65.0
#> [105] checkmate_2.3.2 scales_1.3.0
#> [107] rappdirs_0.3.3 stringr_1.5.1
#> [109] digest_0.6.37 minqa_1.2.8
#> [111] rmarkdown_2.28 htmltools_0.5.8.1
#> [113] pkgconfig_2.0.3 base64enc_0.1-3
#> [115] lme4_1.1-35.5 sparseMatrixStats_1.17.2
#> [117] dbplyr_2.5.0 fastmap_1.2.0
#> [119] rlang_1.1.4 htmlwidgets_1.6.4
#> [121] UCSC.utils_1.1.0 shiny_1.9.1
#> [123] DelayedMatrixStats_1.29.0 jquerylib_0.1.4
#> [125] zoo_1.8-12 jsonlite_1.8.9
#> [127] BiocParallel_1.41.0 BiocSingular_1.23.0
#> [129] RCurl_1.98-1.16 magrittr_2.0.3
#> [131] Formula_1.2-5 scuttle_1.15.5
#> [133] GenomeInfoDbData_1.2.13 munsell_0.5.1
#> [135] Rcpp_1.0.13 ape_5.8
#> [137] viridis_0.6.5 stringi_1.8.4
#> [139] RJSONIO_1.3-1.9 zlibbioc_1.51.2
#> [141] MASS_7.3-61 AnnotationHub_3.15.0
#> [143] plyr_1.8.9 mediation_4.5.0
#> [145] parallel_4.4.1 ggrepel_0.9.6
#> [147] splines_4.4.1 hms_1.1.3
#> [149] igraph_2.1.1 buildtools_1.0.0
#> [151] reshape2_1.4.4 ScaledMatrix_1.13.0
#> [153] futile.options_1.0.1 BiocVersion_3.21.1
#> [155] XML_3.99-0.17 evaluate_1.0.1
#> [157] RcppParallel_5.1.9 lambda.r_1.2.4
#> [159] BiocManager_1.30.25 nloptr_2.1.1
#> [161] tzdb_0.4.0 httpuv_1.6.15
#> [163] rols_3.1.0 tidyr_1.3.1
#> [165] purrr_1.0.2 ggplot2_3.5.1
#> [167] BiocBaseUtils_1.9.0 rsvd_1.0.5
#> [169] xtable_1.8-4 restfulr_0.0.15
#> [171] tidytree_0.4.6 later_1.3.2
#> [173] viridisLite_0.4.2 tibble_3.2.1
#> [175] GenomicAlignments_1.41.0 memoise_2.0.1
#> [177] beeswarm_0.4.0 AnnotationDbi_1.69.0
#> [179] cluster_2.1.6 timechange_0.3.0
#> [181] mia_1.13.47