如何从正态分布中提取观察值并使用自定义估计器计算平均值，然后 运行 此过程在 r 中循环
How to draw observations from a normal distribution and compute the mean using a customized estimator, and then running this procedure in a loop in r
我想使用自定义估算器 Yhat=(1/(n-4))*sum(Y_i) 作为 Y_i.
均值的估算器
我如何从分布 N(4,10) 中抽取 20 个观测值并使用 Yhat 计算样本均值的估计值，然后在循环中重复此过程 k 次并保存结果在 r?
的矩阵中
首先，既然你想要均值的估计，你不需要矩阵，向量应该没问题。
Yhat <- function(x) {
if (length(x) <= 4) {
stop("vector x must be at least 5 elements long")
} else {
sum(x) / (length(x) - 1)
}
}
k <- 100
result <- replicate(k, {
simulations <- rnorm(20, 4, 10)
Yhat(simulations)
})
result
我想使用自定义估算器 Yhat=(1/(n-4))*sum(Y_i) 作为 Y_i.
均值的估算器我如何从分布 N(4,10) 中抽取 20 个观测值并使用 Yhat 计算样本均值的估计值，然后在循环中重复此过程 k 次并保存结果在 r?
的矩阵中首先，既然你想要均值的估计，你不需要矩阵，向量应该没问题。
Yhat <- function(x) {
if (length(x) <= 4) {
stop("vector x must be at least 5 elements long")
} else {
sum(x) / (length(x) - 1)
}
}
k <- 100
result <- replicate(k, {
simulations <- rnorm(20, 4, 10)
Yhat(simulations)
})
result