Я использую модель со смешанными эффектами в R, и я хочу знать, есть ли код, который я могу запустить, чтобы получить контрасты для взаимодействия для последующего c тестирования.
A , Ie. Тест Тьюки, который я сейчас выполняю, таков:
Есть 2 области (ML и HG) и 3 обработки (Low_control, Frog и High_control).
L.S <- pairs(lsmeans(juvweightmodelfeb6th, ~ Region*Treatment_three))
test(L.S, adjust = "tukey")
B. Что дает мне эти сравнения в выводе:
HG,Low_control - ML,Low_control
HG,Low_control - HG,Frog
HG,Low_control - ML,Frog
HG,Low_control - HG,High_control
HG,Low_control - ML,High_control
ML,Low_control - HG,Frog
ML,Low_control - ML,Frog
ML,Low_control - HG,High_control
ML,Low_control - ML,High_control
HG,Frog - ML,Frog
HG,Frog - HG,High_control
HG,Frog - ML,High_control
ML,Frog - HG,High_control
ML,Frog - ML,High_control
HG,High_control - ML,High_control
C. Но что я действительно хочу, так это три сравнения, которые по сути являются контрастами:
1. HG, Low_control - HG, Frog - ML, Low_control - ML, Frog
2. HG, High_control - HG, Frog - ML, High_control - ML, Frog
3. HG, High_control - HG, Low_control - ML, High_control - ML, Low_control
D. Как мне написать код для трех приведенных выше контрастов в тесте post-ho c для модели со смешанными эффектами в R?
Спасибо !!