Hello everyone,
I wonder anyone can give me some guidance about db-RDA analysis. Thank you! The following is my code.
data(B) #B is the genetic distance matrix for bulk soil#
data(env_B) # env_B is the soil characteristics data containing SOM, TN, TP, TK#
data(SD) # SD is the geographic distance matrix#
SD=pcnm(SD)
#condition geographic distance#
b=capscale(B~SOM+TN+TK+TP+Condition(scores(SD)),data=env_B)
But there is a error message "error in rownames(X)：can not find X".
I recommend start using capscale by going through the examples, and the documentation provided by Jari Oksanen (e.g. https://cran.r-project.org/web/packages/vegan/vignettes/intro-vegan.pdf)
You cant feed a distance matrix into capscale! for including spatial distance, a converions of your distance matrix into pcnm eigenvectors is an option, there is plenty of literature on that , see again documentation for pcnm in vegan...
Hi, Robert and Malte,
Thank you for your consideration! I have already solved this problem. Thanks
- Badges
- Science topic
- Similar topics
- Computer Science and Engineering
- Bioinformatics