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 02:20] walmes [Gráfico para estudo de medidas de influência em modelos de regressão linear] |
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 299: | Linha 302: | ||
| #------------------------------------------------------------ | #------------------------------------------------------------ | ||
| </code> | </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> | ||
| + | |||
