Испытание гипотезы по линейной модели n R

Question

Я очень новичок в R, и после прочтения нескольких результатов поиска все еще не удалось выяснить, как проверить различную гипотезу от линейных моделей в R. У меня есть 3 репликации каждого из 3 графиков (лечение). Я хочу проверить гипотезу P2> (P1 + P3)/2 с rep как блокирующий фактор (rcbd) и без rep как фактор блокировки (crd).

Я думаю ЖВВ линейные модели могут соответствовать по rcbd lm(resp ~ factor(plot) + factor(rep), data=pheno) CRD lm(resp ~ factor(plot), data=pheno)

Как проверить гипотезу P2> (P1 + P3)/2 (P означает сюжет)?

pheno 

rep plot resp 
R1 P1 12.6 
R2 P1 11.9 
R3 P1 12.3 
R1 P2 13.0 
R2 P2 12.4 
R3 P2 12.4 
R1 P3 11.3 
R2 P3 11.9 
R3 P3 10.9 

Answer 1

Вы можете использовать glht из multcomp пакета.

У нас есть страница об этом: http://www.ats.ucla.edu/stat/r/faq/testing_contrasts.htm