R学习之统计实验（二）--掷色子问题（R语言编程）-----数模

问题

问题一颗股子投4次至少得到一个六点与两颗般子投24次至少得到 一个双六点，这两个事件哪一个概率更大？

直接R模拟（大数定律）

注意括号各种英文输入，R语言编写时候看不出明显区别

sice<-function(n){
  p1=0;p2=0;
  for (i in 1:n) {
    u1<-runif(4);x1<-floor(u1*6)+1;
    if(max(x1)==6){p1=p1+1}
    x2<-floor(6*runif(24))+1;
    x3<-floor(6*runif(24))+1;
    if(max(x2+x3)==12){p2=p2+1}
  }
  rt<-c('P(A)'= p1/n,'P(B)'= p2/n);rt;
}
sice(10000);

runif等函数之前介绍过，这里不介绍
max
这里表示找到向量得最大值

> max(c(1,2,1.11,1,2,4,5))
[1] 5

用法基本都围绕最大值详见

运行结果

> source('C:/Users/默默/Desktop/R/sice.R', echo=TRUE)

> sice<-function(n){
+   p1=0;p2=0;
+   for (i in 1:n) {
+     u1<-runif(4);x1<-floor(u1*6)+1;
+     if(max(x1)==6){p1=p1+1}
+     x2<-floor(6*runif(2 .... [TRUNCATED] 

> sice(10000);
  P(A)   P(B) 
0.5173 0.4974 

结论第一种几率大