0

我正在尝试使用 field.ID 使用部分 db-RDA 来纠正样本的重复测量特征。但是，包括 Condition(field.ID) 会导致主要关注因素的质心从图中消失（下图左）。

设计：连续两年重复对12个田地进行物种数据抽样。此外，每年都会从参考字段中抽取 3 个样本。由于前一个字段不可用，这三个字段在第二年发生了变化。此外，还对一些环境变量进行了采样（氮、土壤水分、温度）。每个字段都有一个标识符（field.ID）。使用 field.ID 作为 Condition 似乎错误地删除了 F1 因素。但是，使用抽样活动 (SC) 作为 Condition 不会。后者是纠正部分 db-RDA 中重复测量的正确方法吗？

设置种子（1234）

df.exp <- data.frame(field.ID = factor(c(1:12,13,14,15,1:12,16,17,18)),
                     SC = factor(rep(c(1,2), each=15)),
                     F1 = factor(rep(rep(c("A","B","C","D","E"),each=3),2)),
                     Nitrogen  = rnorm(30,mean=0.16, sd=0.07),
                     Temp  = rnorm(30,mean=13.5, sd=3.9),
                     Moist  = rnorm(30,mean=19.4, sd=5.8))

df.rsp <- data.frame(Spec1 = rpois(30, 5),
                     Spec2 = rpois(30,1),
                     Spec3 = rpois(30,4.5),
                     Spec4 = rpois(30,3),
                     Spec5 = rpois(30,7),
                     Spec6 = rpois(30,7),
                     Spec7 = rpois(30,5))

data=cbind(df.exp, df.rsp)                  



dbRDA <- capscale(df.rsp ~ F1 + Nitrogen + Temp + Moist + Condition(SC), df.exp); ordiplot(dbRDA)
dbRDA <- capscale(df.rsp ~ F1 + Nitrogen + Temp + Moist + Condition(field.ID), df.exp); ordiplot(dbRDA)

在此处输入图像描述

4

1 回答 1

1

您因变量而被部分排除ID，然后您尝试解释别名为 this 的变量ID，但它已经被部分排除。打印输出中的关键行是这样的：

Some constraints were aliased because they were collinear (redundant)

事实上，当你询问细节时，你会得到

> alias(dbRDA, names=TRUE)
[1] "F1B" "F1C" "F1D" "F1E"

F1?变量是恒定的， 其中ID已经被部分排除了，没有什么可以解释的了。

于 2016-03-30T18:44:35.053 回答