#part1
Data <- read.table("mediandata.dat")
X <- Data[,1]
results <- matrix(0,1000,1)

for (i in 1:1000){
 boot.samp <- sample(X,20,replace=TRUE)
 results[i] <- median(boot.samp)
}
results <- results - mean(results)
results <- sort(results)

median(X) - results[25]
median(X) - results[975]


#part 2

Y <- rnorm(20,mean=16,sd=4)

results2 <- matrix(0,10000,1)
for (i in 1:10000){
  boot.samp <- sample(Y,20,replace=TRUE)
  results2[i] <- var(boot.samp)
}


results2 <- results2 - mean(results2)
results2 <- sort(results2)

var(Y) - results2[50]
var(Y) - results2[9950]
19*var(Y)/38.580
19*var(Y)/6.843