#59. X<-read.table("C:/Users/Marija/Desktop/OPMS/nerv.txt") x<-X$V1 n<-length(x) T<-sqrt(n)*sum((x-mean(x))^3)/(sum((x-mean(x))^2))^(3/2) T #ocena koeficijenta asimetrije #ocena standardne greske ocene Tt<-c() for (i in 1:n){ x2<-x[-i] Tt[i]<-sqrt(n-1)*sum((x2-mean(x2))^3)/(sum((x2-mean(x2))^2))^(3/2) } s2<-(n-1)/n*sum((Tt-mean(Tt))^2) se<-sqrt(s2) se