简体   繁体   中英

Phylo correlogram in R bootstrapping error?

I have a 16S rDNA tree in a phylo4d format (called 'tree') and I am trying to create a phylo correlogram to test for phylogenetic signal. I am receiving the following error, and I am unsure how to fix the issue.

    library("phylosignal")
    library("ape")

    pc <- phyloCorrelogram(tree)
    Error in boot::boot(X, function(x, z) moranTest(xr = x[z], Wr = prop.table(Wi[z,  : 
  no data in call to 'boot'

I am a novice in R and internet searches yielded results over my head. Thank you for your help. Below is my attached tree (the dropbox link).

You have no data associated to your tree. Your tree is a phylo4d object which has the "tree" information but no data attached to it. You need something like that

library(phylobase)
g1 <- as(geospiza_raw$tree, "phylo4")
geodata <- geospiza_raw$data
g2 <- phylo4d(g1, geodata)
pc <- phyloCorrelogram(g2)

The technical post webpages of this site follow the CC BY-SA 4.0 protocol. If you need to reprint, please indicate the site URL or the original address.Any question please contact:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM