我正在尝试为根据年龄分组的树社区创建Rank Abundance Curves。我在8个密封阶段(15(×8))(N = 120)有15个采样点。如果我不得不考虑我的数据,它将看起来像这样:
Age <- factor(c(rep(1,15), rep(2,15), rep(3,15), rep(4,15), rep(5,15), rep(6,15), rep(7,15), rep(8,15)), labels =c("12","15","19","23","27","31","35","38"))
现在,当我使用函数'rankabuncomp'(BiodiversityR包)为每个seral阶段创建RAC曲线时,我收到一条错误消息:
diversitycomp(x,y,factor,index =“richness”)出错: 指定因子1'年龄'不是一个因素 另外:警告信息: 在if((方法%in%METHOD)== F){: 条件的长度> 1,只使用第一个元素
这是我使用的代码:
Trees_2015 <- read.csv(file = "Trees_2015.csv")
Trees<- Trees_2015[,-1]
Reg_age <- read.csv(file = "Age.csv")
RAC_trees <- rankabundance(Trees,y=Reg_age, factor = "Age", level = c("12","15","19","23","27","31","35","38"))
RAC_trees
rankabunplot(RAC_trees,scale='abundance', addit=FALSE, specnames=c(1,2,3))
rankabuncomp(Trees, y=Reg_age, factor='Age',
scale='proportion', legend=TRUE)
为什么R会产生此错误?我怎样才能纠正它?
'Reg_age'数据框(120变量为1变量)看起来像这样:
Age
1 12
2 12
16 15
20 15
120 38
'Trees'数据框有120个obs。 75个变量(即75种不同的物种)
由于 埃琳娜