r - R:用蒙特卡罗计算 pi 近似的近似误差
问题描述
我是计算机模拟和 R 编程的新手,这是我必须在大学学习的一个新模块 我在理解练习中提出的问题时遇到了一些问题
我正在尝试的代码:
runs <- 100000
#runif samples from a uniform distribution
xs <- runif(runs,min=-0.5,max=0.5)
ys <- runif(runs,min=-0.5,max=0.5)
in.circle <- xs^2 + ys^2 <= 0.5^2
mc.pi <- (sum(in.circle)/runs)*4
# absdif <- abs(mc.pi - 3.14159265)
# print(absdif)
plot(xs,ys,pch=".",col=ifelse(in.circle,"blue","red"),xlab='',ylab='',asp=1, main=paste("MC Approximation of Pi =",mc.pi))
变量“运行”是否被视为练习中的样本量?在此先感谢您的帮助
解决方案
是的,runs
变量配置了抽取的样本数量,所以这是你在这段代码中配置样本大小的方式。如果您查看?runif
帮助页面,第一个参数是要绘制的样本数,因此您的代码从均匀分布中runif(runs, ...)
绘制样本点。runs
推荐阅读
- r - 使用 Case_When 跳过 NA
- javascript - 在 Javascript 中等待服务器响应
- java - 发布时出现 404 错误,但我的服务器可以正常获取数据但无法正确返回
- json - 仅从 scala 中的 HDFS 目录获取所需的文件详细信息
- amazon-web-services - aws 负载均衡器未向实例注册
- bluetooth - 如何解决缺少 printk 输出的问题?
- python - Django迁移大文件
- swift - AVPlayerItem 的 selectedMediaOption(in:) 与 seek(to:) 调用后的可用选择选项不一致
- php - 'key' meta_query 不起作用自定义高级搜索(WordPress)
- php - 将特定价格添加到产品 Prestashop