Как мне изменить значения по оси x в кривой выживания Каплана-Мейера?

#r

#r

Вопрос:

Пожалуйста, посмотрите на мой R-код и помогите мне, как редактировать метки оси x. Мои попытки изменить метки оси x полностью меняют график. Заранее спасибо.

 attach(Pigments)
summary(Pigments)


library(survival)
library(survminer)

km.model <- survfit(Surv(Time.hours, Status)~1, data = Pigments)
plot(km.model,conf.int = T, xlab = "Time in minutes", ylab = "Cum Survival", main= "Cumulative Survival vs Time")
fit=survfit(Surv(Time.hours,Status) ~ Treatment )
survdiff(Surv(Time.hours,Status) ~ Treatment)
res <- pairwise_survdiff(Surv(Time.hours, Status) ~ Treatment,
                         data = Pigments)
res
plot(fit, conf.int = F)
model1<- survfit(Surv(Time.hours, Status) ~ rep(30,20), data =Pigments[1:20,], type= "kaplan-meier")
model2<- survfit(Surv(Time.hours, Status) ~ rep(40,21), data = Pigments[21:41,], type= "kaplan-meier")
model3<- survfit(Surv(Time.hours, Status) ~ rep(50,30), data = Pigments[42:71,], type= "kaplan-meier")
plot(model1, xlab = "Time in hours", ylab = "Cum Survival",conf.int = F)
lines(model1, xlab = "Time in hours", ylab = "cum survival", conf.int = F, col = 5)
lines(model2, xlab = "Time in hours", ylab = "Cum Survival",conf.int = F, col = 2)
lines(model3, xlab = "Time in hours", ylab = "Cum Survival",conf.int = F, col = 3,xlim = c(0,4)) 
legend("topright", cex = 0.65,bty ="n",lty =1, col =1:4)
       legend = c("Trained","Naive","Fed")
plot(1:300,
     xaxt = "n",
     xlab = "My X-Axis")
axis(1,
     at = c(10,50,100,150,200,250),
     labels = c(10,50,100,150,200,250))
 

Комментарии:

1. Не используйте вторую plot команду. Если вы хотите использовать axis() , установите axes = FALSE в исходном plot() вызове и используйте axis() впоследствии, чтобы добавить их к графику.

2. Привет, мы можем получить образец данных?

3. файл:///C:/Users/admin/Desktop/Micro cage Refined data 02-09-2020.csv

4. Грегор Томас и Франческо Гроссетти. Пожалуйста, найдите гиперссылку на образец моего файла в предыдущем сообщении. Заранее спасибо.