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 Ambos lados da revisão seguinte | ||
cursos:ruel:sessao0 [2007/12/19 17:29] uel |
cursos:ruel:sessao0 [2007/12/20 21:56] paulojus |
||
---|---|---|---|
Linha 352: | Linha 352: | ||
hist(resid(lmo.h)) | hist(resid(lmo.h)) | ||
</code> | </code> | ||
+ | |||
+ | Regressão múltipla, fórmulas, transformação e comparação e seleção de modelos. | ||
+ | <code R> | ||
+ | ## | ||
+ | head(hills) | ||
+ | |||
+ | lm2.h <- lm(time ~ dist + climb, data=hills) | ||
+ | anova(lm2.h) | ||
+ | summary(lm2.h) | ||
+ | par(mfrow=c(2,2)) | ||
+ | plot(lm2.h) | ||
+ | |||
+ | shapiro.test(resid(lm2.h)) | ||
+ | |||
+ | par(mfrow=c(1,1)) | ||
+ | boxcox(time ~dist+climb, data=hills) | ||
+ | |||
+ | |||
+ | lm2r.h <- lm(sqrt(time) ~ dist + climb, data=hills) | ||
+ | coef(lm2r.h) | ||
+ | ## ou... aproveitando o modelo previamente definido... | ||
+ | lm2r.h <- update(lm2.h, sqrt(time) ~ ., data=hills) | ||
+ | coef(lm2r.h) | ||
+ | |||
+ | summary(lm2r.h) | ||
+ | |||
+ | ## modelo retirando variavel climb | ||
+ | lm3r.h <- update(lm2r.h, . ~ . - climb) | ||
+ | coef(lm3r.h) | ||
+ | |||
+ | anova(lm3r.h, lm2r.h) | ||
+ | |||
+ | stepAIC(lm(time ~ dist*climb, data=hills)) | ||
+ | <code> | ||
+ | |||
+ | * Link de material sobre fórmulas e declaração de modelos | ||
+ | * Análise de esperimentos inteiramente casualisados, contrastes etc | ||
+ | * experimentos fatoriais | ||
+ |