truemean<- mean(means[,1,1,1,1])
BIAS<- apply(means[,,,,c(2:9,11)],c(2,3,4),mean)-truemean
VARS <- apply(means[,,,,c(2:9,11)],c(2,3,4),function(x) mean((x-mean(x))^2))
MSE <- BIAS^2 + VARS
AVARS <- apply(means[,,,,10],c(2,3),mean)

print("Naive")
print(MSE[1,,1])
print("DR")
print(MSE[,,3])
print("VAR")
print(sqrt(VARS[,,7]))
print(AVARS)