Você está aqui: start » pessoais » walmes » ridiculas

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 anteriorRevisão anterior
Próxima revisão		Revisão anterior
pessoais:walmes:ridiculas [2010/12/21 14:25] – walmespessoais:walmes:ridiculas [2011/02/11 17:01] (atual) – walmes
Linha 1: Linha 1:
 ====== Ridiculas - dicas curtas em R ====== ====== Ridiculas - dicas curtas em R ======
  
-===== Ridiculas - dicas curtas em R ===== +===== Regressão na análise de variância =====
- +
-== Regressão na análise de variância ==+
  
 <code R> <code R>
Linha 111: 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)