CentriIntervalli <- seq(57,127,7)
freq <- c(4,39,122,197,192,126,93,33,16,2,1)
interv <- seq(57,134,7)-3.5; interv
source("http://macosa.dima.unige.it/R/daticlas.txt")
# se hai già scaricato daticlas.txt modifica il comando precedente
# con l'opportuno indirizzo (usa eventualmente Change Dir)
z <- function(x) dnorm(x,mean=mean(XxXx),sd=sd(XxXx))
curve(z,add=TRUE)
abline(v=mean(XxXx),col="red",lty=3)
# Il valore di chi^2:
f <- function(i) integrate(z,53.5+i*7,53.5+(i+1)*7)$value
Prob <- c(f(0),f(1),f(2),f(3),f(4),f(5),f(6),f(7),f(8),f(9),f(10))
FrOs <- freq
n <- sum(FrOs)
FrAt <- Prob*n
sum( (FrOs-FrAt)^2/FrAt )