Как я могу построить стандартные дистрибутивы, подходящие с помощью surreg ()? - PullRequest
Новая
       72

Как я могу построить стандартные дистрибутивы, подходящие с помощью surreg ()?

0 голосов
/ 24 января 2019

Я не могу построить логистические или гауссовские распределения, подходящие для моих данных.

Я загружаю свои данные и создаю объект surreg 

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 нет некоторых важных дистрибутивов.Как я могу построить стандартные дистрибутивы, которые подходят при помощи surreg ()?

...