Пусть
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)