Open chemy-cm opened 8 months ago
Hello Ming,
Thank you for reporting the issue.
I've fixed the problem (closed via 5e9c680). Please pull the latest version of metagMisc
.
It works on a test dataset:
library(phyloseq)
library(metagMisc)
data("GlobalPatterns")
res <- phyloseq_phylo_ses(
physeq = GlobalPatterns,
measures = c("PD", "MPD", "VPD"),
null_model = "richness",
package = "picante",
abundance_weighted = FALSE,
nsim = 5,
swapiter = 5
)
Unfortunately, multi-threading is not implemented yet for the diversity estimation for the randomized datasets, so it may take a while to run.
My Code:
Error information:
I can used this function to calculate "PD" and "MPD", But not "VPD" parameter.
I guess the "VPD" calculate got wrong.
Show "OTU" content:
Show "otu_tree" content:
How can I fix this error? Many Thanks!