模型不收敛

时间:2020-02-21 13:59:12

标签: convergence

我在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功能的说明不清楚,也没有提供示例。

0 个答案:

没有答案