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"
r
generalized-linear-model
survival
gamlss
NiroshaR
fonte
fonte
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.A
glm()
função não suporta a distribuição de Weibull em R infelizmente. Você pode tentar?family
ver quais distribuições estão disponíveis. Eu tentaria usar asurvreg()
partir dosurvival
pacote.fonte
Eu usei o
brms
pacote, que é bayesiano. Ele suporta as famílias Weibull, exponencial, lognormal, Frechet e outras famílias e a censura (esquerda / direita / intervalo) para implementar os modelos AFT. Ele também inclui efeitos aleatórios que são conhecidos nos modelos de sobrevivência como "fragilidade" e uma série de outras opções de regressão, como smoothers estilo gam.Como as abordagens bayesianas usam a amostragem MCMC, é mais lenta que
glm
,gamlss
ousurvreg
, mas também é uma solução abrangente de regressão, e ser bayesiano tem outras vantagens. (Eu amo issostanplot
, que fornece uma série de gráficos de diagnóstico esclarecedores.)fonte