joey711 / phyloseq

phyloseq is a set of classes, wrappers, and tools (in R) to make it easier to import, store, and analyze phylogenetic sequencing data; and to reproducibly share that data and analysis with others. See the phyloseq front page:
http://joey711.github.io/phyloseq/
567 stars 187 forks source link

rarefy_even_depth() Error en rep_len(x["OTUi"], x["times"]): valor de 'length.out' no válido #1753

Open imonteroo opened 1 month ago

imonteroo commented 1 month ago

I have an Error un rarefy_even_depth(). Never had it before. Thank you

ps.met.aa
phyloseq-class experiment-level object
otu_table()   OTU Table:          [ 28680 taxa and 74 samples ]:
sample_data() Sample Data:        [ 74 samples by 102 sample variables ]:
tax_table()   Taxonomy Table:     [ 28680 taxa by 10 taxonomic ranks ]:
taxa are rows
otMA=data.frame(otu_table(ps.met.aa))
mindeepth = min(sapply(otMA,sum)) 
mindeepth
11156903
ps.met.aar = phyloseq::rarefy_even_depth(ps.met.aa, rngseed=123, sample.size=11156900, replace=FALSE, trimOTUs = TRUE)
`set.seed(123)` was used to initialize repeatable random subsampling.
Please record this for your records so others can reproduce.
Try `set.seed(123); .Random.seed` for the full vector
...
Error en rep_len(x["OTUi"], x["times"]): valor de 'length.out' no válido