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 | ||
pessoais:walmes:cursorcpao [2012/06/27 18:46] walmes |
pessoais:walmes:cursorcpao [2012/06/28 09:03] walmes [section 5] |
||
---|---|---|---|
Linha 44: | Linha 44: | ||
* Apresentação Prezi {{http://prezi.com/nldtrhcqyaql/estatistica-experimental/}}; | * Apresentação Prezi {{http://prezi.com/nldtrhcqyaql/estatistica-experimental/}}; | ||
* Experimental Designs for Research in Education {{http://prezi.com/-nfyorqf2gyw/experimental-designs-for-research-in-education/}}; | * Experimental Designs for Research in Education {{http://prezi.com/-nfyorqf2gyw/experimental-designs-for-research-in-education/}}; | ||
+ | |||
+ | <code R> | ||
+ | #========================================================================================== | ||
+ | # Curso de capacitação em Estatística Experimental com uso do R | ||
+ | # Embrapa Agropecuária Oeste - Dourados/MS - de 25 à 29 de Junho de 2012 | ||
+ | # Aula 4 - Regressão logistica e modelos aditivos generalizados | ||
+ | # Walmes Marques Zeviani - LEG/UFPR | ||
+ | #========================================================================================== | ||
+ | |||
+ | #------------------------------------------------------------------------------------------ | ||
+ | # lendo dados de latência em pêssego | ||
+ | |||
+ | lat <- read.table("http://www.leg.ufpr.br/~walmes/data/chimarrita.txt", | ||
+ | header=TRUE, sep="\t") | ||
+ | str(lat) | ||
+ | |||
+ | #------------------------------------------------------------------------------------------ | ||
+ | # ver | ||
+ | |||
+ | pairs(lat[,c("lesao48","lesao60","cor0","ida0","ms0","ss0","f0","b0")]) | ||
+ | |||
+ | pairs(lat[lat$fer=="sim",c("lesao48","lesao60","cor0","ida0","ms0","ss0","f0","b0")]) | ||
+ | |||
+ | #------------------------------------------------------------------------------------------ | ||
+ | # ajuste do modelo com todas as variáveis | ||
+ | |||
+ | lat0 <- subset(lat, fer=="sim", | ||
+ | select=c("lesao48","lesao60","cor0","ida0","ms0","ss0","f0","b0")) | ||
+ | lat0 <- na.omit(lat0) | ||
+ | |||
+ | m0 <- lm(lesao60~cor0+ida0+ms0+ss0+f0+b0, data=lat0) | ||
+ | par(mfrow=c(2,2)); plot(m0); layout(1) | ||
+ | anova(m0) | ||
+ | |||
+ | m1 <- lm(lesao60~(cor0+ida0+ms0+ss0+f0+b0)^2, data=lat0) | ||
+ | par(mfrow=c(2,2)); plot(m1); layout(1) | ||
+ | anova(m1) | ||
+ | drop1(m1, scope=.~., test="F") | ||
+ | summary(m1) | ||
+ | |||
+ | m2 <- step(m1, k=2) # AIC | ||
+ | anova(m2) | ||
+ | summary(m2) | ||
+ | |||
+ | m3 <- step(m1, k=log(nrow(lat0))) # BIC | ||
+ | anova(m3) | ||
+ | summary(m3) | ||
+ | |||
+ | par(mfrow=c(2,2)); plot(m3); layout(1) | ||
+ | |||
+ | plot(lat0$lesao60~fitted(m3)) | ||
+ | |||
+ | summary(lat0) | ||
+ | pred <- expand.grid(cor0=seq(24,112,l=20), ms0=seq(1,18,l=20), | ||
+ | ss0=seq(3.5,9,l=3), b0=seq(9,15,l=3)) | ||
+ | pred$y <- predict(m3, newdata=pred) | ||
+ | |||
+ | wireframe(y~cor0+ms0|ss0*b0, data=pred) | ||
+ | |||
+ | wire <- wireframe(y~cor0+ms0|ss0*b0, data=pred, drape=TRUE) | ||
+ | |||
+ | require(latticeExtra) | ||
+ | |||
+ | useOuterStrips(wire) | ||
+ | |||
+ | #------------------------------------------------------------------------------------------ | ||
+ | # regressão logística | ||
+ | |||
+ | |||
+ | </code> | ||