I performed the demo with the CRC data with no hassle but now while loading my data I have a problem with the dimensions of the data
MMUPHinTax <- read_tsv("./METADATA/MMUPHin_TAXA.csv") MMUPHinmeta <- read_tsv("./METADATA/MMUPHin_METAa.csv") MMUPHinTaxa <- MMUPHinTax %>% column_to_rownames("ID") %>% as.matrix() MMUPHinMETA <- MMUPHinmeta %>% column_to_rownames("ID") fit_adjust_batch <- adjust_batch(feature_abd = Taxa, batch = "study", covariates = NULL, data = META, control = list(verbose = FALSE))
The dimensions seem to match at least with this checking. The labels (sample or ID) match so I don’t know what can be causing trouble here.
> class(MMUPHinMETA)  "data.frame" > class(MMUPHinTaxa)  "matrix" "array" > dim(MMUPHinTaxa)  807 344 > dim(MMUPHinMETA)  344 4
 vegan_2.5-7 lattice_0.20-41 permute_0.9-5
 magrittr_2.0.1 MMUPHin_1.4.2 lpsymphony_1.17.0
 compositions_2.0-2 forcats_0.5.1 stringr_1.4.0
 dplyr_1.0.7 purrr_0.3.4 readr_2.0.1
 tidyr_1.1.3 tibble_3.1.3 ggplot2_3.3.5
 tidyverse_1.3.1 nlme_3.1-152
I already checked this previous issue but maybe I am missing something obvious. Thanks in advance for any help on this issue.