#60.
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")