Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

--- dicas:rstudio [2011/04/25 02:21] – [RStudio ferramenta didática] walmes
+++ dicas:rstudio [2011/04/25 10:44] (atual) – walmes
@@ Linha 17: / Linha 17: @@
   * ✔ 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;
-  * ✔ 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;
+  * ✔ 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:
@@ Linha 300: / Linha 302: @@
 #------------------------------------------------------------
 </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>

Diferenças

Navegação

Busca

Ferramentas

QR Code