R1 Lecture 09 Class Notes

R1 Lecture 09 Class Notes

苹果盘子问题解答

apple_dish <- function(N){
    x <- sample(3,4*N,replace=TRUE)
    apple <- matrix(x,ncol=4)
    result <- apply(apple,1,f1)
    table(result)/N
}

f1 <- function(x){
    # x中有几个盘子中有两个苹果
    length(which(table(x)==2))
}
# 5000次实验结果
apple_dish(5000)

## result
##      0      1      2 
## 0.3280 0.4442 0.2278

中心极限定理的验证

$x_{1},x_{2},...,x_{n},...$ 独立同分布

则根据中心极限定理,当 $(n\rightarrow\infty)$ 时

$$\frac{\sum(x_{i}-\mu)}{\sqrt{n}\sigma}\sim N(0,1)$$

验证步骤

1.找一个分布 distribution,计算它的 $\mu$,$\sigma$

选择分布为 $Binominal(10,0.2)$ ,$\mu=2,\sigma=1.6$
二项分布:伯努利试验独立的重复n次,其试验成功的次数满足参数为(n,p)的二项分布

2.生成n个满足上述分布的随机数,计算其对应的统计量

r_binom10 <- function(n,p){
    mat <- matrix(rbinom(10*n,1,p),ncol=10)
    xi <- apply(mat,1,sum)
    mu <- 10*p
    sigma <- 10*p*(1-p)
    sum(xi-mu)/(sqrt(n)*sigma)
}

3.重复2多次,得到这个统计量的分布

r_binom10_clt <- function(N=100,n,p){
    x <- numeric(length=N)
    for(i in 1:N){
        x[i] <- r_binom10(n,p)
    }
}

4.在图形上表现出这个统计量的分布,及其与正态分布的差距

作业

1. 完善课堂上关于用 $binominal(n,p)$ 验证CLT的程序，用一个函数（包含多个参数）直接生成图形结果

2. 用一个其他的分布来验证CLT