[QUOTE=Arges1984;2131423]OPPURE QUESTE?????
library(MASS)
mu <- 120
sigma <- 12
n1 <- 4
n2 <- 7
n3 <- 8
n.rip <- 10000
F <- seq(1:n.rip)
for (i in 1:n.rip) {
c1 <- rnorm(n1, mu, sigma)
c2 <- rnorm(n2, mu, sigma)
c3 <- rnorm(n3, mu, sigma)
y <- c(c1, c2, c3)
m1 <- mean(c1)
m2 <- mean(c2)
m3 <- mean(c3)
10.75
Simulazione
# stima della varianza entro
m.groups <- c(rep(m1, n1), rep(m2, n2), rep(m3, n3))
dev.entro <- sum( (y - m.groups)^2 )
gdl.entro <- (n1+n2+n3 - 3)
var.entro <- dev.entro / gdl.entro
# stima della varianza tra
m.tot <- mean(y)
dev.tra<-n1*(m1-m.tot)^2+n2*(m2-m.tot)^2+n3*(m3-m.tot)^2
gdl.tra <- 3-1
var.tra <- dev.tra / gdl.tra
F[i] <- var.tra / var.entro[/Q
queste sono per l'appendice r
per i grafici - anova una via sono :
lineplot.CI
e poi quelle per le assunzioni di normalità, omoschedasticità, indipendenza...
nn me le ricordo ...l'ho ricopiate...cmq ci sono sulla tesina del prof.







Rispondi citando

crepi il lupo,è il mio ultimo esame...


