R - 需要帮助模拟估算器答案

【问题标题】：R - Need help simulating estimatorsR - 需要帮助模拟估算器
【发布时间】：2015-05-24 12:32:50
【问题描述】：

我有两个估算器，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，但我不太确定如何使用它。

【问题讨论】：

是的，你应该使用rnorm。
@Roland 我应该如何在 rnorm 中使用这些估计器？对于T1、(((rnorm(10^4,0,1)+rnorm(10^4,0,1))/2)*((rnorm(10^4,0,1)+rnorm(10^4,0,1))/2) 这样的事情是否正确？我想比较这些估计器的标准差
你为什么怀疑你的代码是正确的？
@Roland 我是 R 新手，不确定我是否这样做正确 - 例如，我只是将平均值设为 0，标准差设为 1。我不完全确定我做的是否正确。
阅读?rnorm 的示例。创建 x1、x2、y1 和 y2，假设 µ 和 σ^2 有一些值。他们计算 T1 和 T2 并使用 summary()、sd() 或包中的其他函数进行分析，例如 ?psych::describe()

【解决方案1】：

由于T1 和T2 都依赖于X1、X2、Y1 和Y2，您应该首先模拟这四个随机变量：

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

然后你可以运行你的代码来获取T1和T2的所有模拟值：

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

【讨论】：