R - 需要帮助模拟估算器

时间:2015-05-24 12:32:50

标签: r statistics normal-distribution

我有两个估算师T1 = ((X1+X2)/2)*(Y1+Y2)/2)和。{ T2 = (X1*Y1+X2*Y2)/2

这用于查找估算区域,例如T1 = ((503+505)/2)*((334+330)/2)T2 = ((503*334 + 505*330)/2)

X1, X2通常以平均µ1和方差σ^2分发 Y1,Y2通常以平均µ2和差异σ^2分发。现在我想用10 ^ 4个样本模拟这个数据。我想知道我需要使用哪些相关命令。我在想rnorm,但我不太确定如何使用它。

1 个答案:

答案 0 :(得分:1)

由于T1T2都依赖于X1X2Y1Y2,因此您应首先模拟这四个随机变量:

X1 <- rnorm(1e4, mu1, sigma)
X2 <- rnorm(1e4, mu1, sigma)
Y1 <- rnorm(1e4, mu2, sigma)
Y2 <- rnorm(1e4, mu2, sigma)

然后,您可以运行代码以获取T1T2的所有模拟值:

T1 <- ((X1+X2)/2)*(Y1+Y2)/2)
T2 <- (X1*Y1+X2*Y2)/2