#PRIHODI
Y<-c(3841,7084,7254,15228,18042,19089,22588,23972,25694,27592,27927,31576,32528,32921,33724,36887,37776,37992,39464,40506,44516,46538,51088,51955,54339,57935,75137,82612,83381,84741)
Ys<-sort(Y)
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 #džeknajf ocena 
S<-0
for(i in 1:30){
  S<-S+(Tm[i]-t)^2
}
S2<-29*S/30
Se<-sqrt(S2)
Se

#BAZA NERV
X<-read.table("C:/Users/Marija/Desktop/OPS/nerv.txt")
x<-c(X$V1,X$V2,X$V3,X$V4,X$V5,X$V6)
x1<-x[!is.na(x)]
n<-length(x1)
T<-sqrt(n)*sum((x1-mean(x1))^3)/(sum((x1-mean(x1))^2))^(3/2)
T #ocena koeficijenta asimetrije

#ocena standardne greske ocene
Tt<-c()
for (i in 1:n){
  x2<-x1[-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