R学习之统计实验（十）--新乒乓问题（R语言编程）-----数模

问题

盒中有12个乒乓球．其中9个是新的，第一次比赛时从中选出3个，比赛后仍放回盒中，第二次比赛时再从盒中选出3个．(1）求第二次取出的球都是新球的概率；(2）又已知第二次取出的球都是新球，求第一次取出的球都是新球的概率．分别作理论和模拟分析．

直接R模拟（大数定律）

newpingpang<-function(n){
  m<-0;z<-0
  for(i in 1:n){
    u<-runif(3)
    x<-floor(u*12)+1
    used=3
    used=used+sum(x>(12-used))
    u<-runif(3)
    x<-floor(u*12)+1
    if(sum(x>(12-used))==3){m=m+1}
    if(sum(x>(12-used))==3){if(used==6){z=z+1}}
  }
  rt<-c('第二次取出的球都是新球的概率'=m/n,'假设第二次取出球是新球，第一次取出球是新球的概率'=z/m);rt
}
newpingpang(30000)

运行结果

> newpingpang<-function(n){
+   m<-0;z<-0
+   for(i in 1:n){
+     u<-runif(3)
+     x<-floor(u*12)+1
+     used=3
+     used=used+sum(x>(12-used))
+     u<-runif(3)
+     x<-floor(u*12)+1
+     if(sum(x>(12-used))==3){m=m+1}
+     if(sum(x>(12-used))==3){if(used==6){z=z+1}}
+   }
+   rt<-c('第二次取出的球都是新球的概率'=m/n,'假设第二次取出球是新球，第一次取出球是新球的概率'=z/m);rt
+ }
> newpingpang(30000)
                    第二次取出的球都是新球的概率 
                                      0.03636667 
假设第二次取出球是新球，第一次取出球是新球的概率 
                                      0.05591201 

DONE！！！