我有两个估算师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,但我不太确定如何使用它。
答案 0 :(得分: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