Estou tentando ajustar o modelo linear generalizado para a família weibull, mas quando o tento em R, isso gera um erro. Sei que o weibull não se encaixa na família exponencial, mas li alguns artigos de pesquisa sobre como ajustar o GLM para a família weibull. Se alguém puder me ajudar com isso, eu realmente aprecio. Dá o seguinte erro.
> 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
suporta a distribuição Weibull viaWEI
,WEI2
eWEI3
, todos os 2 parâmetros, embora obviamente diferentes. Porém, não tenho certeza se ele suporta censura, o que seria um elemento-chave de um modelo de sobrevivência da AFT.