#58.
Y<-read.table("C:/Users/Marija/Desktop/OPMS/prihodi.txt")
Ys<-sort(Y$V1)
s<-0
for(i in 1:30){
  s<-s+Ys[i]*((2*i-1)/30-1)
}
Tk<-s/sum(Ys)
Tk #ocena Djini indexa
#ocena standardne greske
Tm<-c()
for (i in 1:30){
  y<-Ys[-i]
  s<-0
  for(j in 1:29){
    s<-s+y[j]*((2*j-1)/29-1)
  }
  Tm[i]<-s/sum(y)
}
Tm
t<-mean(Tm)
t
Tjack<-30*Tk-29*t #djeknajf ocena 
S<-0
for(i in 1:30){
  S<-S+(Tm[i]-t)^2
}
S2<-29*S/30
Se<-sqrt(S2)
Se