“参数'env'的种类(list)不对和参数的长度不是一错误”，跪求大神解答

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位老师，请教个问题：
使用dbstats包中的dblm函数做基于距离的线性模型，spe为物种数据，env为环境因子数据，代码如下：
library(vegan)
library(dbstats)
> spe=read.csv("spe.csv",row.names=1)
> env=read.csv("env.csv",row.names=1)
> spe=log1p(spe)
> env=decostand(env,"standardize")
> spe.dist=vegdist(spe,"bray")
> spe.dblm=dblm(spe.dist~env,data=env)

报错：Error in model.frame.default(formula = formula, data = data, drop.unused.levels = TRUE) : 参数'env'的种类(list)不对




使用BiodiversityR包中的CAPdiscrim函数做典范主坐标排序，spe为物种数据，env为环境因子数据，代码如下：
library(vegan)
library(BiodiversityR)
> spe=read.csv("spe.csv",row.names=1)
> env=read.csv("env.csv",row.names=1)
> spe=log1p(spe)
> env=decostand(env,"standardize")

环境因子以env表示：
> spe.cap=CAPdiscrim(spe~env,data=env,dist="bray",axes=2,m=0,permutations=1000)
报错：Error in `[.data.frame`(y, , group) : undefined columns selected


环境因子以各自名字表示：
>spe.cap=CAPdiscrim(spe~X1+X2+X3+X4+X5+X6+X7+X8+X9+X10+X11+X12+X13+X14,data=env,dist="bray",axes=2,m=0,permutations=1000)
报错：Error in eval(predvars, data, env) : 数值'envir'参数的长度不是一

跪求大神帮忙解答，万分感谢！附原始CSV数据。
原始CSV数据.rar (4.52 KB) 本附件包括：

• env.csv
• spe.csv

2016-12-7 12:14:47 上传

原始CSV数据

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：list IST Biodiversity Permutation standardize library Error 种类 模型