Как вы управляете выводом функции R cuminc?

Question

Например, в консоли,

> x=cuminc(ftime=c(1,2,3,4,5,6,5,4),c(1,0,0,1,1,0,0,0),0)
> x
Estimates and Variances:
$`est`
        1     2     3   4         5
0 1 0.125 0.125 0.125 0.3 0.5333333

$var
       1        2        3        4      5
0 1 0.015625 0.015625 0.015625 0.040625 0.0725

Как бы я получил 'est' данные?

Answer 1

Если мы проверяем str(x), это вложенный список

library(cmprsk)
str(x)
#List of 1
# $ 0 1:List of 3
#  ..$ time: num [1:8] 0 1 1 4 4 5 5 6
#  ..$ est : num [1:8] 0 0 0.125 0.125 0.3 ...
#  ..$ var : num [1:8] 0 0 0.0156 0.0156 0.0406 ...
# - attr(*, "class")= chr "cuminc"

Извлеките list с помощью [[, а затем используйте $ или [[ для извлечения элемента 'est'

x[[1]]$est
#[1] 0.0000000 0.0000000 0.1250000 0.1250000 0.3000000 0.3000000 0.5333333 0.5333333

Если мы хотим получить тот же результат, что и в распечатанной смете

d1 <- data.frame(x[[1]][1:2])
library(dplyr)
library(tidyr)
d1 %>% 
  distinct(est, .keep_all = TRUE) %>% 
  filter(time > 0) %>%
  complete(time = full_seq(c(min(time), max(time)), 1)) %>% 
  fill(est) %>% 
  pull(est)
#[1] 0.1250000 0.1250000 0.1250000 0.3000000 0.5333333