[英]Post-Hoc Statistical Analysis for Repeated Measures ANOVA Treatment within Time Effect
我很难找到事后测试来破译我在会话影响范围内的“会话”或治疗时间。
这是我的数据:
TR SESSION MAC FISHD ID
1 1 1 3.285714286 1
2 1 2 0.571428571 2
2 1 3 3.571428571 3
1 1 4 4 4
1 2 1 4 5
2 2 2 6.571428571 6
2 2 3 3.142857143 7
1 2 4 8.857142857 8
1 3 1 0.714285714 9
2 3 2 1.714285714 10
2 3 3 4.428571429 11
1 3 4 0.714285714 12
这就是我反复测量的方法:
model.b = lme(FISHD ~ TR + SESSION + TR*SESSION,
random = ~1|MAC,
data=TTDall2)
> ACF(model.b)
lag ACF
1 0 1.0000000
2 1 -0.7547232
3 2 0.4852727
> model2 = lme(FISHD ~ TR + SESSION + TR*SESSION,
+ random = ~1|MAC,
+ correlation = corAR1(form = ~ SESSION | MAC,
+ value = -0.7547232),
+ data=TTDall2,
+ method="REML")
> Anova(model2)
Analysis of Deviance Table (Type II tests)
Response: FISHD
Chisq Df Pr(>Chisq)
TR 0.2014 1 0.6536
SESSION 25.0418 1 5.61e-07 ***
TR:SESSION 103.9113 1 < 2.2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
我试图破译何时有治疗效果。 有任何想法吗?
目前,您正在将会话建模为连续变量,因此,对于TRT = 2中的SESSION值,报告的TR:SESSION
系数有点像“趋势测试”,范围为1:3。您希望发现什么? 让我们以正确的方式进行处理,尤其要创建ordered
因素:
model2 = lme(FISHD ~ factor(TR)*ordered(SESSION), # same as TRT+SESSION+TRT*ORDERED
random = ~1|MAC,
correlation = corAR1(form = ~ SESSION | MAC,
value = -0.7547232),
data=dat,
method="REML")
anova(model2)
#------
numDF denDF F-value p-value
(Intercept) 1 4 37.48966 0.0036
factor(TR) 1 2 0.20124 0.6976
ordered(SESSION) 2 4 13.94379 0.0157
factor(TR):ordered(SESSION) 2 4 52.02864 0.0014
因此,这就是说,将TRT与SESSION序列之间的相互作用考虑在内的分析发现,TRT的效果存在显着差异。 如果不对交互进行建模(使用公式: FISHD ~ factor(TR)+ordered(SESSION)
,则这两个分量都不重要:
model1 = lme(FISHD ~ factor(TR)+ordered(SESSION),
random = ~1|MAC,
correlation = corAR1(form = ~ SESSION | MAC, value = -0.7547232),
data=dat,
method="REML")
anova(model1)
#-------------
numDF denDF F-value p-value
(Intercept) 1 6 44.53272 0.0005
factor(TR) 1 2 0.14143 0.7430
ordered(SESSION) 2 6 2.55154 0.1578
和:
> anova(model1, model2)
Model df AIC BIC logLik Test L.Ratio p-value
model1 1 7 54.74813 55.30422 -20.37406
model2 2 9 44.04435 42.17018 -13.02217 1 vs 2 14.70378 6e-04
Warning message:
In anova.lme(model1, model2) :
fitted objects with different fixed effects. REML comparisons are not meaningful.
这是一个警告,但我认为它并不适用于此,因为模型中的变量是相同的,只是结构不同。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.