Я хочу использовать смесь гамма-распределения в качестве параметрической модели для анализа выживания на цензурированных данных с использованием R. В пакете "flexsurv" есть разные распределения, но я не смог найти распределение гамма-смеси. В этом пакете говорится, что:
«Может быть установлено любое параметрическое распределение, определенное пользователем, с учетом, по крайней мере, функции R, определяющей плотность вероятности или опасность».
https://cran.r -project.org / веб / пакеты / flexsurv / flexsurv.pdf
Есть ли способ напрямую определить распределение гамма-смеси (с заранее заданным количеством компонентов) параметрическим способом, чтобы напрямую использовать этот пакет для оценки максимального правдоподобия?
data <- Surv(ages, censored)
fit_gammamixture <- flexsurvreg(data~1, dist=???)
Я нашел эту статью об анализе выживания со смесью гамма-распределений, но трудно понять и реализовать алгоритм, представленный здесь.
Моделирование цензурированных данных времени жизни с использованием смеси базовых гамма-характеристик
https://projecteuclid.org/download/pdf_1/euclid.ba/1340371053