【问题标题】:How to do Monte Carlo simulations in R [closed]如何在 R 中进行蒙特卡洛模拟 [关闭]
【发布时间】:2015-04-29 03:52:02
【问题描述】:

我有以下信息:

A = B/C
B has mean of 50 with SD of 10
C has mean of 3 with SD of 0.5

假定所有变量的分布都是正态的。

如何针对这种情况运行蒙特卡罗模拟?感谢您的帮助。

【问题讨论】:

  • ?rnorm?replicate
  • 我可以通过这些获得正态分布和重复。但是如何进行模拟呢?你会推荐一些包或一些功能吗?
  • 您希望 B 和 C 以某种方式相关或相关,还是相互独立?
  • A、B 和 C 由给定的公式关联,尽管我们可以认为 B 和 C 是独立的。
  • 两个正态随机变量的比率服从柯西分布,可以由rcauchy()生成。如果提到具体目标,回答会容易得多。

标签: r montecarlo


【解决方案1】:

对于粗略的蒙特卡洛,您可以运行n 实验:

set.seed(1)
n = 1000
A = rnorm(n, 50, 10)/rnorm(n, 3, 0.5)
mean(A)
#[1] 17.09732

【讨论】:

    猜你喜欢
    • 2018-09-24
    • 1970-01-01
    • 2016-08-19
    • 2021-01-06
    • 1970-01-01
    • 1970-01-01
    • 2023-01-27
    • 2018-04-03
    • 2023-03-05
    相关资源
    最近更新 更多