Code: Select all
n<-length(mtcars$mpg)
r<-n/5;r
kvant<-qnorm(1:(r-1)/r,mean(mtcars$mpg),sd(mtcars$mpg)); kvant
sk<-c()
for (i in 1:(r)) sk[i]<-length(mtcars$mpg[mtcars$mpg<kvant[i]]); sk
Ni<-c()
Ni[1]<-sk[1]; Ni[2:(n/5)]<-diff(sk); Ni
T<-sum(((Ni-5)^2)/5); T
alpha<-0.05
kritiska_vertiba1<-qchisq(1-alpha,r-3); kritiska_vertiba1