Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anterior Revisão anterior Próxima revisão | Revisão anterior | ||
pessoais:walmes:ridiculas [2010/12/21 14:25] walmes |
pessoais: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 113: | 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> |