R - 需要帮助模拟估算器
R - Need help simulating estimators
我有两个估算器，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，但我不太确定我将如何使用它。
由于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
我有两个估算器，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，但我不太确定我将如何使用它。
由于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