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 | ||
|
dicas:rstudio [2011/04/25 01:42] walmes |
dicas:rstudio [2011/04/25 02:20] walmes [Gráfico para estudo de medidas de influência em modelos de regressão linear] |
||
|---|---|---|---|
| Linha 245: | Linha 245: | ||
| </code> | </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> | ||
