#odredjivanje trake poverenja za histogram
X<-read.table("C:/Users/Marija/Desktop/OPMS/SDSS1.txt")
x<-X[,2]
n<-length(x)
#formira se histogram sa 308 podeoka 
m<-308 
h<-diff(range(x)/m)
breaks<-min(x)+h*0:m
#crta se histogram, ali sa ocenjenim gustinama na y-osi (zato je freq=FALSE)
p<-hist(x, breaks = breaks, freq = FALSE, ylim = c(0,80), main = "Traka poverenja")

alfa<-0.05
c<-qnorm(alfa/(2*m))*sqrt(m/n)/2
c 
f<-p$density #ocenjene vrednosti gustine
#c je negativna vrednost pa ce u biti donja, a l gornja granica
l<-(sqrt(f)-c)^2
u<-(pmax(sqrt(f)+c,0))^2
#dodajemo po jednu vrednost na kraj identicnu poslednjoj vrednost da bi se nacrtao i poslednji deo trake
l[length(l)+1]<-l[length(l)]
u[length(u)+1]<-u[length(u)]
lines(breaks, l, type = "s", col="blue")
lines(breaks, u, type = "s", col="red")