我在flexsurvreg
中为生存模型拟合了不同的分布。对数逻辑分布运行良好,但对数正态分布和伽马分布族显示以下警告:
优化可能尚未收敛到最大可能性-黑森州不是正定的
我尝试了不同的优化器来收敛模型:
library(flexsurv)
recsurv<-Surv(data$tree_age,data$dummy)
flexAFT_gengamma=flexsurvreg(recsurv~log(LAR_cm2_g)+log(root_colar_cm),
anc = list(shape = ~ species+health+height_class_cm), method = "L-BFGS-B", lower = -Inf,
upper = Inf, control = list(5), data, dist="gengamma")
或:
flexAFT_genf<-flexsurvreg(recsurv~log(LAR_cm2_g)+log(root_colar_cm),
anc = list(shape = ~ species + height_class_cm+health), method = "Nelder-Mead",
control=list(fnscale = 2500), data, dist="genf")
但没有成功。 我的数据在链接-link to data中。如果有人可以给我一个提示,我将不胜感激。有关flexsurvreg功能的说明不清楚,也没有提供示例。