如何计算期望值和方差，然后从R中的此分布模拟500个样本

Question

假设我们有一个离散的随机变量，其概率函数为

x = c(0, 3, 5, 10)

和

px = c(0.2, 0.1, 0.1, 0.6)

如何计算期望值和方差，然后从R中的这种分布模拟500个样本？

Answer 1

让我们

x <- c(0, 3, 5, 10)
px <- c(0.2, 0.1, 0.1, 0.6)

您可以如下计算x的期望值和方差：

Ex <- sum(x * px)
Vx <- sum(((x - Ex) ^ 2) * px)

然后使用sample模拟数据：

sample(x, size = 500, prob = px, replace = TRUE)