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