如何在R中进行蒙特卡罗模拟

Question

我有以下信息：

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

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

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

Answer 1

对于粗略的蒙特卡罗，您可以运行n实验：

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