R попарно пост-хо c чтобы получить контрасты для взаимодействия? - PullRequest
0 голосов
/ 10 февраля 2020

Я использую модель со смешанными эффектами в 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?

Спасибо !!

...