Question

我真的很难看到我从同事那里得到的一些代码。

我有以下模型，其中因子1具有ABC选项而因子2具有ab选项

model<-aov(response~factor1*factor2,data=mydata)
coef(model)

(Intercept)                factor1-A                     factor1-B 
  1.82404                        -0.29790                  -0.80398 
 factor2-b                factor1-A:factor2-b     factor1-B:factor2-b 
 -0.37816                        0.04560                   0.07672

然后，我想确定此模型中的主效应和成对比较。所以我的同事给了我这段代码：

lambda <- rbind(c(0, 1, 0, 0,1/2, 0),   #Main effect: factor1, A-C
                c(0, 0, 1, 0, 0,1/2),   #Main effect: factor1, B-C
                c(0, 0, 0, 1,1/3, 1/3), #Main effect factor2: a-b
esticon(model,lambda)

我无法弄清楚如何操纵lambda代码以便我可以获得“factor1，A-B”的主要效果有谁知道怎么做？我在网上看过，我发现对于如何使用这些矩阵没有很好的解释。如果有人能指出我对这些lambda矩阵如何工作的一个很好的参考，那也会很棒。另外，我已经安装了libary（doBy）。这就是esticon的来源。

谢谢！

Answer 1

无法访问您的数据，但听起来您正在寻求进行因子方差分析。

首先，您需要运行TukeyHSD(model, ...)，然后使用$factor1和$factor2获取有关因子1和2的更多信息。如果你想看看两者之间的相互作用，你可以使用

$`factor1:factor2`

Here is a link了解您可能会发现有用的更多信息。

在r中使用lambda（）表示矩阵的系数

1 个答案: