R : ggplot2 / plotly를 사용하여 "계단"플로팅

R 프로그래밍 언어를 사용하고 있습니다. 이 튜토리얼을 따르려고합니다 : https://rviews.rstudio.com/2017/09/25/survival-analysis-with-r/ (페이지 하단).

이 자습서의 코드를 약간 수정하고 데이터의 3 가지 관측치에 해당하는 "계단"(예 : "생존 함수", 아래 그림에서 "빨강", "파란색", "녹색"))을 플로팅했습니다.

 library(survival)
    library(dplyr)
    library(ranger)
    library(data.table)
library(ggplot2)
library(plotly)
    
a = na.omit(lung)
a$ID <- seq_along(a[,1])

r_fit <- ranger(Surv(time,status) ~ age + sex + ph.ecog + ph.karno + pat.karno + meal.cal + wt.loss, data = a, mtry = 4, 
importance = "permutation", splitrule = "extratrees", verbose = TRUE)

death_times <- r_fit$unique.death.times
surv_prob  <-data.frame(r_fit$survival)
avg_prob <- sapply(surv_prob, mean)

plot(r_fit$unique.death.times, r_fit$survival[1,], type = "l", ylim = c(0,1), col = "red", xlab = "Days", ylab = "survival", main = "Survival Curves")

new = a[1:3,]

pred <- predict(r_fit, new, type = 'response')$survival
pred <- data.table(pred)
colnames(pred) <- as.character(r_fit$unique.death.times)

plot(r_fit$unique.death.times, pred[1,], type = "l", col = "red")

lines(r_fit$unique.death.times, r_fit$survival[2,], type = "l", col = "green")

lines(r_fit$unique.death.times, r_fit$survival[3,], type = "l", col = "blue")

여기에서 위의 플롯을 "대화 형"으로 만들고 싶습니다. 커브 중 하나 위로 마우스를 움직일 때 다음과 같이 만들고 싶습니다.

1. 해당 곡선에 속하는 "속성"(객체 "a"에서) 호버 (예 : ID, 나이, 성별, ph.ecog 등)
2. 1의 동일한 "호버 상자"에서 마우스가 (주어진 곡선에 대해) 마우스를 올려 놓은 각 위치에 대한 x 좌표 (r_fit $ unique) 및 y 좌표 ( "pred"에서)도 표시합니다.

내 계획은 먼저 "grob"객체를 "ggplot"객체로 변환 한 다음 "ggplot"객체를 "plotly"객체로 변환하는 것이 었습니다.

 grob= plot(r_fit$unique.death.times, pred[1,], type = "l", col = "red")
basic_plot = ggpubr::as_ggplot(grob)

그러나 "basic_plot"을 검사하려고하면 "NULL"로 표시됩니다.

 ggplot(f)
Error: `data` must be a data frame, or other object coercible by `fortify()`, not an S3 object with class gg/ggplot

이것이 효과가 있었다면 결국 ggplot 객체를 플롯으로 변환했을 것입니다.

plotly_plot = ggplotly(final_plot)

누군가이 대화 형 플롯을 만드는 방법을 보여줄 수 있습니까?

나는 이것과 비슷한 것을 얻기 위해 노력하고있다 : https://plotly.com/python/v3/ipython-notebooks/survival-analysis-r-vs-python/ (페이지 하단에 "수명"이라는 제목으로 플롯 서로 다른 종양 DNA 프로필 ")

(참고 : 저는 USB 포트 나 인터넷 연결이없는 컴퓨터로 작업하고 있습니다. 몇 개의 사전 설치된 라이브러리가있는 R 만 있습니다. "ggplotify"또는 "survminer"가 없습니다.)

감사