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 Última revisão Ambos lados da revisão seguinte | ||
|
dicas:rstudio [2011/04/19 22:52] paulojus |
dicas:rstudio [2011/04/25 02:21] walmes [RStudio ferramenta didática] |
||
|---|---|---|---|
| Linha 17: | Linha 17: | ||
| * ✔ Outras Aproximações pela normal; | * ✔ Outras Aproximações pela normal; | ||
| * Convergência da média de realizações binomial, Poisson, beta, etc, para uma distribuição normal controlando tamanho da amostra; | * Convergência da média de realizações binomial, Poisson, beta, etc, para uma distribuição normal controlando tamanho da amostra; | ||
| + | * ✔ Gráfico para estudo de medidas de influência em modelos de regressão linear (Walmes); | ||
| + | * ✔ Teste de normalidade aplicado aos dados e aos resíduos (Walmes); | ||
| + | * ✔ Teste de correlação para dados e resíduos de experimentos (Walmes). | ||
| Dicas sobre o editor: | Dicas sobre o editor: | ||
| Linha 202: | Linha 205: | ||
| ) | ) | ||
| | | ||
| + | #------------------------------------------------------------ | ||
| + | </code> | ||
| + | |||
| + | ==== Gráfico para estudo de medidas de influência em modelos de regressão linear ==== | ||
| + | |||
| + | <code R> | ||
| + | # por Walmes ------------------------------------------------ | ||
| + | |||
| + | require(manipulate) | ||
| + | data(anscombe) | ||
| + | ans0 <- anscombe[,c("x1","y1")] | ||
| + | ans0 <- ans0[order(ans0$x1),] | ||
| + | rownames(ans0) <- NULL | ||
| + | ans1 <- ans0 | ||
| + | rx <- 2*diff(range(ans1$x1)) | ||
| + | ry <- 2*diff(range(ans1$y1)) | ||
| + | cols <- rep(1,nrow(ans0)) | ||
| + | |||
| + | layout(matrix(c(1,2,1,3,4,5),2,3)) | ||
| + | |||
| + | manipulate({ | ||
| + | ans1[po,] <- ans0[po,]+c(dx,dy) | ||
| + | plot(ans1) | ||
| + | points(ans1[po,], col="red", pch=19) | ||
| + | abline(a=3, b=0.5, col="gray50", lty=2) | ||
| + | m1 <- lm(y1~x1, data=ans1) | ||
| + | abline(m1, col=2) | ||
| + | h <- hatvalues(m1)[po] | ||
| + | r <- residuals(m1)[po] | ||
| + | legend("topleft", bty="n", | ||
| + | legend=c(substitute(h==ha, list(ha=round(h,3))), | ||
| + | substitute(r==re, list(re=round(r,3))))) | ||
| + | legend("bottomright", legend=po, bty="n") | ||
| + | plot(m1) | ||
| + | }, | ||
| + | po=slider(1,nrow(ans0), step=1, initial=5), | ||
| + | dx=slider(-rx+0.001, rx, initial=0), | ||
| + | dy=slider(-ry, ry, initial=0)) | ||
| + | |||
| + | #------------------------------------------------------------ | ||
| + | </code> | ||
| + | |||
| + | ==== Teste de normalidade aplicado aos dados e aos resíduos ==== | ||
| + | |||
| + | <code R> | ||
| + | # por Walmes ------------------------------------------------ | ||
| + | |||
| + | par(mfrow=c(2,1)) | ||
| + | |||
| + | manipulate({ | ||
| + | m <- rep(seq(0,by=h1,length.out=nlev), nrep) | ||
| + | x <- rnorm(m, m, sd) | ||
| + | xp <- qqnorm(x); qqline(x) | ||
| + | rug(xp$x); rug(xp$y, side=2) | ||
| + | legend("topleft", legend=shapiro.test(x)$p, bty="n") | ||
| + | m0 <- lm(x~factor(m)) | ||
| + | xp <- qqnorm(residuals(m0)); qqline(residuals(m0)) | ||
| + | rug(xp$x);rug(xp$y, side=2) | ||
| + | legend("topleft", bty="n", | ||
| + | legend=shapiro.test(residuals(m0))$p) | ||
| + | }, | ||
| + | h1=slider(0.001, 10, initial=1), | ||
| + | nlev=slider(2, 15, initial=5), | ||
| + | nrep=slider(2, 25, initial=5), | ||
| + | sd=slider(0.01, 10, initial=1)) | ||
| + | | ||
| + | #------------------------------------------------------------ | ||
| + | </code> | ||
| + | |||
| + | ==== Teste de correlação para dados e resíduos de experimentos ==== | ||
| + | |||
| + | <code R> | ||
| + | # por Walmes ------------------------------------------------ | ||
| + | |||
| + | require(MASS) | ||
| + | par(mfrow=c(2,1)) | ||
| + | |||
| + | manipulate({ | ||
| + | m <- rep(seq(0,by=h1,length.out=nlev), nrep) | ||
| + | x <- mvrnorm(length(m), mu=c(0,0), | ||
| + | Sigma=matrix(c(1,cor,cor,1),2,2)) | ||
| + | x[,1] <- x[,1]+m; x[,2] <- x[,2]+m | ||
| + | plot(x) | ||
| + | legend("topleft", legend=cor.test(x[,1], x[,2])$est, bty="n") | ||
| + | m0 <- aov(x~factor(m)) | ||
| + | r <- residuals(m0) | ||
| + | plot(r) | ||
| + | legend("topleft", legend=cor.test(r[,1], r[,2])$est, bty="n") | ||
| + | }, | ||
| + | h1=slider(0,19.99,initial=0.01), | ||
| + | nrep=slider(10,300,initial=20), | ||
| + | nlev=slider(2,15,initial=5), | ||
| + | cor=slider(-0.99,0.99,initial=0)) | ||
| + | |||
| #------------------------------------------------------------ | #------------------------------------------------------------ | ||
| </code> | </code> | ||
