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")