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 | ||
|
dicas:rstudio [2011/04/25 01:42] walmes |
dicas:rstudio [2011/04/25 10:44] (atual) walmes |
||
|---|---|---|---|
| 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. | + | * ✔ 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); | ||
| + | * ✔ Taxa de erro tipo I e tipo II (Walmes). | ||
| Dicas sobre o editor: | Dicas sobre o editor: | ||
| Linha 241: | Linha 244: | ||
| dx=slider(-rx+0.001, rx, initial=0), | dx=slider(-rx+0.001, rx, initial=0), | ||
| dy=slider(-ry, ry, 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> | ||
| + | |||
| + | ==== Taxa de erro tipo I e tipo II ==== | ||
| + | |||
| + | <code R> | ||
| + | # por Walmes ------------------------------------------------ | ||
| + | |||
| + | require(manipulate) | ||
| + | |||
| + | n <- 70 | ||
| + | xx <- seq(qnorm(0.9), 4, l=n) | ||
| + | yy <- dnorm(xx, 0, 1) | ||
| + | |||
| + | area <- function(x,y){ | ||
| + | da <- rbind(cbind(x,y), c(x[1],y[1])) | ||
| + | DA <- sapply(1:length(x), | ||
| + | function(o){ | ||
| + | dir <- da[o,1]*da[o+1,2] | ||
| + | esp <- da[o,2]*da[o+1,1] | ||
| + | c(dir, -esp) | ||
| + | } | ||
| + | ) | ||
| + | abs(sum(apply(DA, 1, sum)/2)) | ||
| + | } | ||
| + | |||
| + | manipulate({ | ||
| + | curve(dnorm(x, 0, 1), -4, 9, col=1, | ||
| + | ylim=c(-0.15,0.5), yaxt="n", ylab="f(x)") | ||
| + | axis(2, at=seq(0,0.4,0.1)) | ||
| + | curve(dnorm(x, i, 1), col=2, add=TRUE) | ||
| + | xx2 <- seq(-4+i, qnorm(0.9, i, 1)+2, l=n) | ||
| + | yy2 <- dnorm(xx2, i, 1) | ||
| + | xx2 <- c(min(xx2), xx2, max(xx2)) | ||
| + | yy2 <- c(0,yy2,0) | ||
| + | yy2 <- pmin(yy2, dnorm(xx2, 0, 1)) | ||
| + | polygon(xx2, yy2, col="gray90") | ||
| + | polygon(c(xx, rev(xx)), | ||
| + | c(yy, rep(0, length(yy))), den=10) | ||
| + | segments(i, 0, i, 0.4, col=2) | ||
| + | segments(0, 0, 0, 0.4, col=1) | ||
| + | tipo2 <- round(area(xx2,yy2), 3) | ||
| + | text(2, -0.075, pos=1, label=expression(alpha==0.10)) | ||
| + | arrows(2, -0.075, 2, -0.01, length=0.1) | ||
| + | text(i/2, 0.45, pos=3, label=substitute(beta==b, list(b=tipo2))) | ||
| + | arrows(i/2, 0.45, i/2, max(yy2)+0.01, length=0.1) | ||
| + | }, | ||
| + | i=slider(0, 6, step=0.01, initial=0) | ||
| + | ) | ||
| #------------------------------------------------------------ | #------------------------------------------------------------ | ||
| </code> | </code> | ||
