Я имитирую некоторые гамма-случайные числа
plot(density(rgamma(10000,8.1,rate=0.00510)),lwd=2,las=1,cex.axis=0.75,
main=expression(paste("Gamma Distribution with",' scale ',alpha," and rate
",beta)))
plot(density(rgamma(10000,2.1,rate=0.00110)),lwd=2,las=1,cex.axis=0.75,
main=expression(paste("Gamma Distribution with",' scale ',alpha," and rate
",beta)))
plot(density(rgamma(10000,2.1,rate=110)),lwd=2,las=1,cex.axis=0.75,
main=expression(paste("Gamma Distribution with",' scale ',alpha," and rate
",beta)))
Мне нужно смоделировать этот тип гаммы с некоторым хвостом, в среднем около 1200. Я выбирал случайные числа, чтобы получить эти значения, учитывая определение ожидания и дисперсии для гамма-распределения, но в первом случае я получаю отрицательные числа, Я не хочу этого. Во втором случае то же самое, но также на обоих графиках вероятность по оси y настолько мала, что я хотел бы увеличить эту вероятность, но я не знаю, как выбрать адекватные параметры для этого.
С другой стороны, параметры, указанные на третьем графике, дают странную плотность, потому что вероятность по оси y больше единицы. Я могу получить значения больше 1. Я этого не понимаю.
dgamma
вместоrgamma
. - person MichaelChirico   schedule 10.03.2018x <- seq(0,5000, length = 10000); y <- dgamma(x, 8.1, .0051); plot(x,y, type = "l")
- person hpesoj626   schedule 10.03.2018