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 20:43] walmes [Conteúdo] |
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 37: | Linha 41: | ||
==== Materiais do curso ==== | ==== Materiais do curso ==== | ||
- | * Diretório com todos os {{http://www.leg.ufpr.br/~walmes/cursoR/cnpaf/|arquivos do Curso}}; | + | * Diretório com todos os {{http://www.leg.ufpr.br/~walmes/cursoR/cpao/|arquivos do Curso}}; |
+ | * 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/}}; | ||
+ | <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> | ||
+ | |||
+ | |||
+ | ---- | ||
==== Links úteis ==== | ==== Links úteis ==== | ||