Я работаю над дизайном повторных измерений GEE и хотел бы получить результат post-Ho c теста , но я не знаю, как это сделать .
Я просто использовал эти методы, чтобы получить результат теста post-Ho c , но он кажется неправильным.
Моя таблица данных выглядит следующим образом:
id times y dose doTimes
1 1 200 250 3
1 2 300 250 3
1 3 280 210 3
1 4 280 125 3
2 1 248 254 3
2 2 345 148 3
2 3 2654.4 73 3
3 1 912 223 3
3 2 523.5 98 3
......
Затем я получил подходящие значения:
gee <- geeglm(y ~ doTimes + dose, data=myData, id=id, family = gaussian,corstr = "exchangeable")
myData$fitted.Values <- gee$fitted.values
myData$fitted.Values
id times y dose doTimes fitted.Values
6 2 313.5 100 3 1036.6992
6 3 1101.8 100 3 1036.6992
6 4 2501.8 100 3 1036.6992
9 1 246 350 3 954.9299
9 2 676.5 350 3 954.9299
9 3 1125.6 350 3 954.9299
9 4 2123.8 350 3 954.9299
10 1 794 300 3 971.2838
10 2 598.5 300 3 971.2838
10 3 1625.4 300 3 971.2838
10 4 4379.2 300 3 971.2838
11 1 200 200 3 1003.9915
После этого я запускаю TukeyHSD:
f <- aov(myData$fitted.Values ~ myData$doTimes)
HSD <- TukeyHSD(f, conf.level = 0.95)
HSD$`myData$doTimes`
diff lwr upr p adj
1-0 546.4474 535.9451 556.9497 0.0000000003426653
2-0 414.4079 405.7935 423.0224 0.0000000003426653
3-0 161.3648 150.0971 172.6325 0.0000000003426653
2-1 -132.0395 -140.8439 -123.2350 0.0000000003426653
3-1 -385.0826 -396.4962 -373.6690 0.0000000003426653
3-2 -253.0432 -262.7478 -243.3385 0.0000000003426653
Весь код выглядит так:
gee <- geeglm(y ~ doTimes + dose, data=myData, id=id, family = gaussian,corstr = "exchangeable")
myData$fitted.Values <- gee$fitted.values
f <- aov(myData$fitted.Values ~ myData$doTimes)
HSD <- TukeyHSD(f, conf.level = 0.95)
HSD$`myData$doTimes`
Значение p.adj кажется неправильным, я знаю, может быть, я допустил очень низкую ошибку ... но если кто-нибудь сможет мне помочь, это будет оценено! Спасибо.