用R

Question

在R中使用mvabund时，是否有人有过提取计算估计相关矩阵的经验？我不清楚它在任何分析创建的对象中的位置。下面的代码是运行该函数的示例。我想知道如何将这个矩阵从名为glm.test的对象中拉出来。

我使用的是R版本3.2.1（2015-06-18）

平台：x86_64-apple-darwin13.4.0（64位）

在OS X 10.10.5（Yosemite）下运行

library("mvabund")

#generating some fake data
species_data<-data.frame()

for(i in 1:100){
    species_data<-rbind(species_data, rpois(100,round(runif(1,1,100)))) 
}
#generating a vector of treatments
treatments<-c(rep("A",50),rep("B",50))

#turning species data into an mvabund object
species_data<-mvabund(species_data)

#running the manyglm() function
glm.test<-manyglm(species_data~treatments,family="negative.binomial")