我尝试使用varpart
功能在纯素包中执行RDA分析。
我在不同的站点(行)中有一个社区结构数据矩阵,我想用3个环境数据矩阵来解释。 3个矩阵是:
在每个矩阵中,行数相同(站点),并且列号不同(测量变量)。
我使用以下脚本:
WQ_RDA <- read.csv("F:/Sher_sophi_new3/RDA/WQ_RDA_1.csv", row.names = 1)
LU2_loc_RDA <- read.csv("F:/Sher_sophi_new3/RDA/LU2_loc_RDA_1.csv", row.names = 1)
LU2_bas_RDA <- read.csv("F:/Sher_sophi_new3/RDA/LU2_basin_RDA_1.csv", row.names = 1)
comm <- read.csv("F:/Sher_sophi_new3/RDA/Final_true_OTUs_97_1.csv", row.names = 1)
x1 = as.matrix(WQ_RDA)
x2 = as.matrix(LU2_loc_RDA)
x3 = as.matrix(LU2_bas_RDA)
comm = as.matrix(t(comm))
RDA_Ger = varpart(comm, x1, x2, x3, transfo="hellinger", scale = FALSE,
na.action = na.omit)
我收到此错误消息:
&#39;&#39; as.vector(x,mode)出错:无法强制类型&#39;关闭&#39;至 矢量类型&#39;任何&#39;
任何人都可以帮我弄清楚我做错了什么吗?