我正在尝试使用 vegan 包执行简单的 RDA,以使用以下数据框测试深度、盆地和扇区对遗传种群结构的影响。
数据文件.
“ALL”变量是遗传种群分配(结构)。
如果指向我的数据的链接不能正常工作,我将在此处粘贴我的数据框片段。
我以这种方式读取数据:
RDAmorph_Oct6 <- read.csv("RDAmorph_Oct6.csv")
我的问题有两个方面:1)我似乎无法正确读取我的遗传变量。我尝试了三件事来解决这个问题。
gen=rda(ALL ~ Depth + Basin + Sector, data=RDAmorph_Oct6, na.action="na.exclude")
Error in eval(specdata, environment(formula), enclos = globalenv()) :
object 'ALL' not found
In addition: There were 12 warnings (use warnings() to see them)
所以,我尝试了这样的事情:
> gen=rda("ALL ~ Depth + Basin + Sector", data=RDAmorph_Oct6, na.action="na.exclude")
Error in colMeans(x, na.rm = TRUE) : 'x' must be numeric
所以我指定了数字
> RDAmorph_Oct6$ALL = as.numeric(RDAmorph_Oct6$ALL)
> gen=rda("ALL ~ Depth + Basin + Sector", data=RDAmorph_Oct6, na.action="na.exclude")
Error in colMeans(x, na.rm = TRUE) : 'x' must be numeric
我真的很困惑。我也试过用 指定每个变量dataset$variable
,但这也不起作用。
奇怪的是,如果我查看环境变量对不同的复合变量的影响,我可以让 rda 工作
MC = RDAmorph_Oct6[,5:6]
H_morph_var=rda(MC ~ Depth + Basin + Sector, data=RDAmorph_Oct6, na.action="na.exclude")
请注意,我确实尝试仅提取上面遗传 rda 的 ALL 列。这也不起作用。无论如何,这导致了我的第二个问题。
当我尝试绘制 rda 时,我得到了一个非常奇怪的情节。注意三个地方的五个点。我不知道这些是从哪里来的。
我将不得不绘制遗传 rda,我想我会提出同样的问题,所以我想我现在就问。
我已经阅读了几个教程,并尝试了每个问题的多次迭代。我在这里提供的是我认为最好的总结。如果有人能给我一些线索,我将不胜感激。