unbiase application plot

set.seed(1)x1 <- rpois (100,lambda = 5)x2 <- rpois (100,lambda = 5)x3 <- rpois (100,lambda = 5)x4 <- rpois (100,lambda = 5)x5 <- rpois (100,lambda = 5)x6 <- rpois (100,lambda = 5)x7 <- rpois (100,lambda = 5)x8 <- rpois (100,lambda = 5)x9 <- rpois (100,lambda = 5)x10 <- rpois (100,lambda = 5)
x11 <- rpois (100,lambda = 5)x12 <- rpois (100,lambda = 5)x13 <- rpois (100,lambda = 5)x14 <- rpois (100,lambda = 5)x15 <- rpois (100,lambda = 5)x16 <- rpois (100,lambda = 5)x17 <- rpois (100,lambda = 5)x18 <- rpois (100,lambda = 5)x19 <- rpois (100,lambda = 5)x20 <- rpois (100,lambda = 5)
x21 <- rpois (100,lambda = 5)x22 <- rpois (100,lambda = 5)x23 <- rpois (100,lambda = 5)x24 <- rpois (100,lambda = 5)x25 <- rpois (100,lambda = 5)x26 <- rpois (100,lambda = 5)x27 <- rpois (100,lambda = 5)x28 <- rpois (100,lambda = 5)x29 <- rpois (100,lambda = 5)x30 <- rpois (100,lambda = 5)
x31 <- rpois (100,lambda = 5)x32 <- rpois (100,lambda = 5)x33 <- rpois (100,lambda = 5)x34 <- rpois (100,lambda = 5)x35 <- rpois (100,lambda = 5)x36 <- rpois (100,lambda = 5)x37 <- rpois (100,lambda = 5)x38 <- rpois (100,lambda = 5)x39 <- rpois (100,lambda = 5)x40 <- rpois (100,lambda = 5)
x41 <- rpois (100,lambda = 5)x42 <- rpois (100,lambda = 5)x43 <- rpois (100,lambda = 5)x44 <- rpois (100,lambda = 5)x45 <- rpois (100,lambda = 5)x46 <- rpois (100,lambda = 5)x47 <- rpois (100,lambda = 5)x48 <- rpois (100,lambda = 5)x49 <- rpois (100,lambda = 5)x50 <- rpois (100,lambda = 5)
a <- data.frame(x1,x2,x3,x4,x5,x6,x7,x8,x9,x10, x11,x12,x13,x14,x15,x16,x17,x18,x19,x20, x21,x22,x23,x24,x25,x26,x27,x28,x29,x30, x31,x32,x33,x34,x35,x36,x37,x38,x39,x40, x41,x42,x43,x44,x45,x46,x47,x48,x49,x50)
N <- sapply(a[,1:50], FUN = mean) ＃生成一个向量
m1 <- mean(N[1])m5 <- mean(N[1:5])m10 <- mean(N[1:10])m15 <- mean(N[1:15])m20 <- mean(N[1:20])m25 <- mean(N[1:25])m30 <- mean(N[1:30])m35 <- mean(N[1:35])m40 <- mean(N[1:40])m45 <- mean(N[1:45])m50 <- mean(N[1:50])
b <-c(m1,m5,m10,m15,m20,m25,m30,m35,m40,m45,m50)c <-c(1,5,10,15,20,25,30,35,40,45,50)
par(family='STKaiti')            ＃写中文。。不懂
plot(c,b,xlim = c(0,50),ylim = c(4.85,5.2),xlab = "观测次数",ylab = "均值",main = "无偏性在泊松分布lamda=5的应用")abline(h=c(5),lty=2,lwd=2,col="red")