#66. X<-rnorm(100) T<-median(X) Tb<-c() for(i in 1:1000){ Xs<-sample(X,100,replace=T) Tb[i]<-median(Xs) } t<-mean(Tb) s2<-var(Tb) se<-sqrt(s2) se #drugi nacin f<-function(x){ n<-length(x) Xs<-sample(x,1000*n,replace = TRUE) Xs<-array(Xs,c(n,1000)) #pravi matricu od n vrsta i 1000 kolona Tb<-apply(Xs,2,median) #primenjuje funkciju median na kolone od Xs se<-sqrt(var(Tb)) return(se) }