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 Próxima revisão Ambos lados da revisão seguinte | ||
|
pessoais:walmes:cursorcpao [2012/06/26 21:10] walmes [Materiais do curso] |
pessoais:walmes:cursorcpao [2012/06/28 09:03] walmes [section 5] |
||
|---|---|---|---|
| Linha 1: | Linha 1: | ||
| ====== Curso em análise de experimentos com afastamento das pressuposições ====== | ====== Curso em análise de experimentos com afastamento das pressuposições ====== | ||
| + | |||
| + | {{:pessoais:walmes:dsc05114_2_.jpg?nolink&800|}} | ||
| + | |||
| + | ---- | ||
| ==== Descrição ==== | ==== Descrição ==== | ||
| Linha 40: | 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> | ||
