R中的部分正则对应分析

R中的部分正则对应分析,r,vegan,R,Vegan,在R中进行部分规范对应分析(pCCA)时,我遇到了一些问题。与代码相关的数据非常广泛,因此我无法将其包括在这里 下面的代码在下面生成错误。在pCCA模型中,我试图在解释物种矩阵时考虑环境和空间变量。空间变量是纬度和经度值。Env2变量是一系列连续的和少量的二进制(0,1)环境变量 mod2 <-cca(species ~ env2 + spatial) Error in model.frame.default(~env2 + spatial, na.action = na.pass, xl

在R中进行部分规范对应分析(pCCA)时,我遇到了一些问题。与代码相关的数据非常广泛,因此我无法将其包括在这里

下面的代码在下面生成错误。在pCCA模型中,我试图在解释物种矩阵时考虑环境和空间变量。空间变量是纬度和经度值。Env2变量是一系列连续的和少量的二进制(0,1)环境变量

mod2 <-cca(species ~ env2 + spatial)
Error in model.frame.default(~env2 + spatial, na.action = na.pass, xlev = NULL) : invalid type (list) for variable 'env2'

mod2公式右侧必须有变量,但似乎有多个变量的数据框。这将不起作用,但会给出与您的帖子中类似的错误消息(这是有文档记录的)。此外,您的公式不会定义部分CCA,因为该公式不包含定义部分调用的术语的函数
Condition()

如果将数值矩阵作为术语,公式接口可能会工作,但它不适用于
unlist()
变量

mod2 <-cca(species ~ env2 + spatial)
Error in model.frame.default(~env2 + spatial, na.action = na.pass, xlev = NULL) : invalid type (list) for variable 'env2'
如果您使用的是素食主义者2.5-1或更新版本,您可以将不带配方界面的部分CCA定义为

cca(物种、环境2、空间)

数据帧
env2
spatial
会自动展开为模型矩阵,在分析
env2
术语的影响之前,会将
spatial
术语分割出来。

您不需要提供实际数据,但如果您想知道代码不起作用的原因,提供一个代表您实际数据的玩具示例数据集是更快获得更好答案的好方法。对于找到你的帖子并想从中学习的人来说,这也会让你的问题变得更清楚。更多信息,请参阅。