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

最新推荐文章于 2022-11-06 14:19:03 发布

默默小艺

最新推荐文章于 2022-11-06 14:19:03 发布

阅读量7.5k

点赞数

分类专栏： R 数模文章标签： R 数模

本文链接：https://blog.csdn.net/qq_36306781/article/details/81275774

版权

数模同时被 2 个专栏收录

18 篇文章 5 订阅

订阅专栏

17 篇文章 11 订阅

订阅专栏

问题

问题一颗股子投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