Estoy tratando de ajustar el modelo lineal generalizado para la familia weibull, pero cuando lo intento en R, aparece un error. Sé que weibull no encaja en la familia exponencial, pero he leído algunos artículos de investigación sobre cómo ajustar GLM para la familia weibull. Si alguien me puede ayudar con esto, realmente lo aprecio. Da el siguiente error.
> data(lung)
> glm(time ~ age+sex+ph.ecog+ wt.loss, family = weibull(link='log'), data = lung)
Error in glm(time ~ age + sex + ph.ecog + wt.loss, family = weibull(link = "log"), :
could not find function "weibull"
gamlss
admite la distribución de Weibull a través deWEI
,WEI2
yWEI3
, todos los 2 parámetros, aunque obviamente diferentes. Sin embargo, no estoy seguro de si es compatible con la censura, que sería un elemento clave de un modelo de supervivencia AFT.