sampleSize <- 30;
lambda <- 1;

numIter <- 1000;
aveMSE <- rep(0,100);

for (alpha in 1:100) {

  mse <- rep(0,1000);
  for (i in 1:numIter) {
    x <- rpois(sampleSize,lambda);
    mse[i] <- (lambda - 
	((alpha/100)*mean(x)) + ((1-(alpha/100))*var(x)))^2
  }
  aveMSE[alpha] <- mean(mse);
}

plot(1:100,aveMSE)