# File "Rjack5". numreps = 100 n = 100 estlist = jacklist = NULL for (rep in 1:numreps) { # Generate some data. x = runif(n) estimate = max(x) # Compute the jackknife estimate of bias. thetam = rep(0,n) for (i in 1:n) { if (i==1) xm = x[2:n] else if (i==n) xm = x[1:(n-1)] else xm = c(x[1:(i-1)], x[(i+1):n]) thetam[i] = max(xm) } thetadot = mean(thetam) jackbiasest = (n-1) * (thetadot - estimate) jacklist = c(jacklist, jackbiasest) estlist = c(estlist, estimate) } truebiasest = mean(estlist - 1) print(jacklist) print(mean(jacklist)) print(truebiasest)