Saya mencoba menyesuaikan model linear umum untuk keluarga weibull, tetapi ketika saya mencobanya di R, itu memberikan kesalahan. Saya tahu bahwa weibull tidak cocok dengan keluarga eksponensial, tetapi saya telah membaca beberapa artikel penelitian tentang pemasangan GLM untuk keluarga weibull. Jika ada yang bisa membantu saya dengan ini, saya sangat menghargai. Ini memberikan kesalahan berikut.
> 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
sumber
sumber
gamlss
mendukung distribusi Weibull viaWEI
,WEI2
danWEI3
, semua 2-parameter meskipun jelas berbeda. Tidak yakin apakah itu mendukung sensor, yang akan menjadi elemen kunci dari model survival AFT.The
glm()
Fungsi tidak mendukung distribusi Weibull di R sayangnya. Anda dapat mencoba?family
melihat distribusi mana yang tersedia. Saya akan mencoba menggunakansurvreg()
darisurvival
paket sebagai gantinya.sumber
Saya telah menggunakan
brms
paket itu, yaitu Bayesian. Ini mendukung Weibull, eksponensial, lognormal, Frechet, dan keluarga lainnya dan sensor (kiri / kanan / interval) sehingga menerapkan model AFT. Ini juga termasuk efek acak yang dikenal dalam model survival sebagai "kelemahan", dan sejumlah opsi regresi lainnya seperti smoothers bergaya gam.Karena pendekatan Bayesian menggunakan pengambilan sampel MCMC, itu lebih lambat daripada
glm
,,gamlss
atausurvreg
, tetapi juga merupakan solusi regresi yang komprehensif, dan menjadi Bayesian memiliki keunggulan lain. (Saya suka itustanplot
, yang menyediakan sejumlah plot diagnostik yang menerangi.)sumber