我正在尝试使用afex包和mixed（）函数（在R中）为混合效果模型建模，但仅接收每个因子的一般p值。如何计算主要效果的成对比较？特别是，我对每只动物进行了3次测试，并对这些测试之间的差异感兴趣。
这是我的代码：
m1<-mixed(stand~trt*test+(1|calf),
method="KR",
check_contrasts = TRUE,
data
print(m1) #summary(m1) only reports correlations of the fixed effects
输出： 混合模型Anova表（类型3测试，KR方法） 型号：站立〜trt *测试+（1 |小腿） 数据：数据 效果df F p.value 1 trt 1，22 0.78 .39 2测试2，44 7.93 ** .001 3 trt：测试2，44 1.19 .31
不幸的是，我不知道如何准备数据示例。因此，我复制了自己的部分数据。我希望这会有所帮助。
小牛测试trt新郎站
母亲3075的1 30918
2 30919 of control 58 313
母亲的3 30920 31 397
母亲4 30921 0 239
5 30922 OF控件16674
6 30923 of control 1 434