我无法绘制适合我的数据的逻辑或高斯分布。
我加载我的数据并创建一个 survreg 对象
gbcs <- read.csv("https://ryanwomack.com/data/gbcs.csv")
recsurv <- Surv(rectime,censrec, gbcs)
如果我尝试拟合指数分布,它不会绘制
fit_exp <- survreg(recsurv~1, dist="exponential")
plot(fit_exp)
Error in xy.coords(x, y, xlabel, ylabel, log) :
'x' is a list, but does not have components 'x' and 'y'
相反,当我使用包 flexsurv 时,它很容易绘图
library(flexsurv)
fit_exp <- flexsurvreg(recsurv~1, dist="exp")
plot(fit_exp)
我的观点是:flexsurv 包没有一些重要的发行版。如何绘制符合 survreg() 的标准分布?