#74. X<-read.table("C:/Users/Marija/Desktop/OPMS/suicide.txt") x<-X$V1 par(mfrow=c(2,2)) n<-length(x) hist(x, xlim=c(-200,1000)) #Sturdzisovo prvilo: h=R/(1+log2(n)) broj<-ceiling(1+log2(n)) #ceo deo + 1 sirina<-diff(range(x)/broj) breaks<-min(x)+sirina*0:broj hist(x, breaks = breaks, xlim=c(-200,1000), main="Sturdzesovo pravilo") #Skotovo pravilo: h=3.5*sd*n^(-1/3) h1<-3.5*sd(x)*n^(-1/3) m<-min(x) M<-max(x) broj1<-ceiling((M-m)/h1) breaks<-m+h1*0:broj1 hist(x, breaks = breaks, xlim=c(-200,1000), main="Skotovo pravilo") #Fridman-Dijakonis: h=2*IQR*n^(-1/3) h2<-2*IQR(x)*n^(-1/3) broj2<-ceiling((M-m)/h2) breaks<-m+h2*0:broj2 hist(x, breaks = breaks,xlim=c(-200,1000), main="Fridman-Diakonisovo pravilo") #malo pomeren histogram u odnosu na poslednji breaks<-m+h2*0:broj2-30 hist(x, breaks = breaks,xlim=c(-200,1000), main="Fridman-Diakonisovo pravilo")