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 | |||
|
pessoais:walmes:cursorcpao [2012/06/29 23:57] walmes [section 7] |
pessoais:walmes:cursorcpao [2012/06/29 23:57] (atual) walmes [Materiais do curso] |
||
|---|---|---|---|
| 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> | ||
| - | #------------------------------------------------------------------------------------------ | ||
| - | # análise de dados binomial | ||
| - | |||
| - | trat <- gl(3, 8) | ||
| - | X <- model.matrix(~trat) | ||
| - | pars <- c(2,1,1) | ||
| - | nu <- X%*%pars | ||
| - | p <- exp(nu)/(1+exp(nu)) | ||
| - | y <- rbinom(length(p), size=10, prob=p) | ||
| - | |||
| - | dados <- data.frame(trat=trat, y=y) | ||
| - | |||
| - | g0 <- glm(cbind(y, 10-y)~trat, data=trat, family=binomial) | ||
| - | summary(g0) # se Residual Deviance/graus de liberade próximo de 1, use test="Chisq" | ||
| - | anova(g0, test="Chisq") | ||
| - | |||
| - | #------------------------------------------------------------------------------------------ | ||
| - | # como fazer testes de "médias" | ||
| - | |||
| - | require(doBy) | ||
| - | popMeans(g0, effect="trat") # estimativas dos parâmetros do preditor | ||
| - | |||
| - | require(glht) | ||
| - | summary(glht(g0, linfct=mcp(trat="Tukey"))) # diferença nas estimativas | ||
| - | |||
| - | #------------------------------------------------------------------------------------------ | ||
| - | # experimento em bloco | ||
| - | |||
| - | dados <- expand.grid(bloco=gl(3,1), trat=gl(5,1)) | ||
| - | X <- model.matrix(~bloco+trat, dados) | ||
| - | pars <- c(1, 1.1,1.2, -1,-0.5,0.5,0.3) | ||
| - | nu <- X%*%pars | ||
| - | |||
| - | p <- exp(nu)/(1+exp(nu)) | ||
| - | dados$y <- rbinom(nrow(dados), size=50, prob=p) | ||
| - | |||
| - | g0 <- glm(cbind(y, 50-y)~bloco+trat, data=dados, family=binomial) | ||
| - | summary(g0) | ||
| - | |||
| - | popMeans(g0, effect="trat") | ||
| - | |||
| - | summary(glht(g0, linfct=mcp(trat="Tukey"))) | ||
| - | |||
| - | #------------------------------------------------------------------------------------------ | ||
| - | </code> | ||
| - | |||
| ---- | ---- | ||
