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 | ||
ridiculas [2011/06/07 19:25] walmes |
ridiculas [2011/07/04 16:52] walmes [section 2] |
||
---|---|---|---|
Linha 7: | Linha 7: | ||
---- | ---- | ||
+ | |||
+ | ==== Gráfico de valores observados e curva de valores preditos para modelo linear generalizado ==== | ||
+ | |||
+ | Temporariamente sem descrição.\\ | ||
+ | palavras-chave: #glm, #poisson, #predict. | ||
+ | |||
+ | <code R> | ||
+ | #------------------------------------------------------------------------------------------ | ||
+ | # por Walmes | ||
+ | #------------------------------------------------------------------------------------------ | ||
+ | |||
+ | da <- expand.grid(trat=gl(2,1), tempo=1:20) | ||
+ | y <- rpois(nrow(da), lambda=da$tempo/5) | ||
+ | |||
+ | g0 <- glm(y~trat/(tempo+I(tempo^2)), data=da, family=poisson) | ||
+ | |||
+ | new <- expand.grid(trat=gl(2,1), tempo=seq(1,20,l=50)) | ||
+ | new$p0 <- predict(g0, newdata=new, type="response") | ||
+ | |||
+ | plot(y~tempo, da, col=da$trat) | ||
+ | with(subset(new, trat=="1"), lines(p0~tempo, col=1)) | ||
+ | with(subset(new, trat=="2"), lines(p0~tempo, col=2)) | ||
+ | |||
+ | #------------------------------------------------------------------------------------------ | ||
+ | </code> | ||
+ | |||
+ | ---- | ||
+ | |||
==== Gráficos de barras com intervalos de confiança para as médias ==== | ==== Gráficos de barras com intervalos de confiança para as médias ==== |