我无法绘制符合我的数据的逻辑或高斯分布。
我加载数据并创建一个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()的标准分布?