Мой вопрос должен быть очень простым.
У меня есть объект flexsurvreg, и я хочу нанести его на график, чтобы увидеть соответствие моих данных о выживании некоторым распределениям.
График в основном работает, но по оси Y данные идут только от 0,6 до 1. Поэтому мне нужно ограничить эту ось. Мне нужно избегать его начала, чтобы я мог правильно видеть соответствие.
Я использую функцию и параметр следующим образом, но это не работает. Ось Y продолжает отображаться от 0 до 1 и даже не возвращает никакого предупреждения.
plot(fit, ylim = c(0.6, 1))
Параметр работает для оси X, но не для оси Y, например:
plot(fit, xlim = c(300, 1300), ylim = c(0.6, 1))
Можно воспроизвести ошибку (с другим набором данных) следующим образом:
library(survival)
library(flexsurv)
gbcs <- read.csv("https://ryanwomack.com/data/gbcs.csv")
recsurv <- Surv(gbcs$rectime, gbcs$censrec)
fit <- flexsurvreg(recsurv ~ 1, dist="exp")
plot(fit, ylim = c(0.6, 1))
Этот же сбой распространяется на все остальные дистрибутивы: dist = "weibull", "gamma", "gengamma "," genf "," lnorm "," gompertz "