Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças

Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

Link para esta página de comparações

Ambos lados da revisão anterior Revisão anterior
Próxima revisão
Revisão anterior
pessoais:walmes:ridiculas [2010/12/21 16:03]
walmes
pessoais:walmes:ridiculas [2011/02/11 17:01] (atual)
walmes
Linha 109: Linha 109:
 </​code>​ </​code>​
  
 +===== Experimento com dois fatores de efeito aditivo e perda de muitas parcelas =====
  
 +<code R>
 +#​------------------------------------------------------------------------------------------
 +# dados
  
 +da <- expand.grid(rept=1:​5,​ ep=factor(1:​5),​ tr=factor(1:​4))
 +da$y <- c(58.4, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
 +          68.4, NA, NA, NA, NA, 258.8, 265.6, NA, NA, NA, NA, NA, 250, NA, 278.8,
 +          268.8, NA, NA, NA, 309.6, NA, NA, NA, NA, NA, NA, NA, NA, NA, 254, 598.8,
 +          NA, NA, NA, NA, 250, 399.6, 260, NA, NA, NA, 288.4, NA, NA, NA, 397.2, NA,
 +          NA, 337.6, NA, 415.2, NA, 450.8, NA, NA, NA, NA, 393.2, NA, NA, NA, NA,
 +          NA, NA, NA, 380.4, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 634, 417.2,
 +          NA, NA, NA, NA, NA)
 +
 +#​------------------------------------------------------------------------------------------
 +# ajuste do modelo aditivo com teste F marginal
 +
 +m0 <- lm(y~ep+tr, data=da)
 +drop1(m0, test="​F"​)
 +
 +#​------------------------------------------------------------------------------------------
 +# análise de resíduos
 +
 +par(mfrow=c(2,​2))
 +plot(m0)
 +layout(1)
 +
 +#​------------------------------------------------------------------------------------------
 +# estimativas dos efeitos sob a restrição do R
 +
 +summary(m0)
 +
 +#​------------------------------------------------------------------------------------------
 +# obtenção das médias ajustadas para os níveis de tratamento
 +
 +require(contrast)
 +lapply(levels(da$tr),​
 +       ​function(i){
 +         ​contrast(m0,​ type="​average",​ list(tr=i, ep=levels(da$ep)))
 +       }
 +       )
 +
 +#​------------------------------------------------------------------------------------------
 +# comparação multipla de médias
 +
 +require(multcomp)
 +summary(glht(m0,​ linfct=mcp(tr="​Tukey"​)))
 +
 +#​------------------------------------------------------------------------------------------
 +</​code>​

QR Code
QR Code pessoais:walmes:ridiculas (generated for current page)