0

我很难找到一个事后测试来破译“会话”或我在会话内接受治疗的时间影响。

这是我的数据：

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

我试图破译何时有治疗影响。有任何想法吗？

4

1 回答 1

0

目前，您正在将会话建模为连续变量，因此对于TR:SESSIONTRT = 2 中的 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.

这是一个警告，但我认为它不适用于这里，因为模型中的变量是相同的，只是结构不同。

于 2017-07-08T18:57:23.723 回答