Я пытаюсь оценить совокупные вероятности распространения твида ie с пакетом твида ie. Мне кажется, что функции ptweed ie и qtweed ie работают для параметра значения дисперсии = 1. После использования значения> 1 я получаю результаты, которые мне не нужны.
Здесь Вот небольшой пример, иллюстрирующий проблему:
library(tweedie)
power <- 1
mu <- 1000
phi <- 1
Fn <- ecdf(rtweedie(1000000, phi=phi, xi=var.power, mu=mu))
Fn(mu)
ptweedie(mu,phi=phi, xi=var.power, mu=mu)
qtweedie(0.5,phi=phi, xi=var.power, mu=mu)
Запустив этот код, я получаю следующий результат, который имеет смысл:
> Fn(mu)
[1] 0.50869
> ptweedie(mu,phi=phi, xi=var.power, mu=mu)
[1] 0.5084094
> qtweedie(0.5,phi=phi, xi=var.power, mu=mu)
[1] 1000
Но однажды я использую значение> 1, результаты кумулятивная вероятность, оцененная по эмпирическому распределению (полученному с помощью rtweed ie), отличается от тех, которые получены непосредственно от ptweed ie и qtweed ie.
Например, выполнение:
library(tweedie)
power <- 1
mu <- 1000
phi <- 2
Fn <- ecdf(rtweedie(1000000, phi=phi, xi=var.power, mu=mu))
Fn(mu)
ptweedie(mu,phi=phi, xi=var.power, mu=mu)
qtweedie(0.5,phi=phi, xi=var.power, mu=mu)
дает следующий результат:
> Fn(mu)
[1] 0.512228
> ptweedie(mu,phi=phi, xi=var.power, mu=mu)
[1] 1
> qtweedie(0.5,phi=phi, xi=var.power, mu=mu)
[1] 500
Среднее значение, по-видимому, делится на фи, что не относится к эмпирическому распределение.
Может ли кто-нибудь объяснить, почему теоретическое распределение, сгенерированное ptweed ie и / или qtweed ie, не соответствует эмпирическому распределению, оцененному благодаря rtweed ie (думаю, правильные ) пока параметры точно такие же?
Большое спасибо.