Contrastes de médias

# o objetivo aqui e demonstrar testes de contrastes de grupos de medias

# vamos utilizar exemplo de Kaps and Lamberson(2009, pg 254), contido no pacote easyanova
data(data1)

data1

# realizando a analise de variância

r<-ea1(data1, design=1)

# onservando o resultado

r

# desdobrando em contrastes ortogonais utilizando a funçao ec() do pacote easyanova

# primeiro contraste

mg1=312;mg2=c(278,280); sdg1=7.7028;sdg2=c(7.7028,7.7028); df=12
ec(mg1,mg2,sdg1,sdg2,df)

# segundo contraste

mg1=280;mg2=278; sdg1=7.7028;sdg2=7.7028; df=12
ec(mg1,mg2,sdg1,sdg2,df)

# ou entao fazendo via funções R base

# matriz de contrastes ortogonais

c=matrix(c(-1,-1,2,1,-1,0),ncol=2)

c=t(c)

c

# inversa generalizada da matriz de contrastes ortogonais

c1=ginv(c)

# modelo com os contrastes e o resultado

contrasts(data1$Diet)=c1

m=lm(Gain~Diet, data=data1)

summary(m)