В настоящее время я использую пакет для выживания, чтобы посмотреть на группу пациентов, поступивших в отделение интенсивной терапии.
Я использовал информацию на этой странице, чтобы помочь: https://cran.r-project.org/web/packages/survivalAnalysis/vignettes/multivariate.html
Мой код:
covariate_names <- c(age="Age at Dx",
sex="Sex",
comorb_no = "Number of Comorbdities")
map(vars(dataset_onlyfirstadmit_hospdischarge$age, dataset_onlyfirstadmit_hospdischarge$sex, dataset_onlyfirstadmit_hospdischarge$comorb_no), function(by)
{
analyse_multivariate(dataset_onlyfirstadmit_hospdischarge,
vars(dataset_onlyfirstadmit_hospdischarge$survivalyrs_2, dataset_onlyfirstadmit_hospdischarge$status),
covariates = list(by), # covariates expects a list
covariate_name_dict = covariate_names)
}) %>%
forest_plot(factor_labeller = covariate_names,
endpoint_labeller = c(time="OS"),
orderer = ~order(HR),
labels_displayed = c("endpoint", "factor", "n"),
ggtheme = ggplot2::theme_bw(base_size = 10))
Однако, когда я ввожу свои переменные и изменяю набор данных, я получаю следующую ошибку: Ошибка: индексы больше, чем количество элементов: ( затем следуют многие числа в этом формате: 57,5, 75,699997, 40,700001, 76,800003, 61.200001, 91,900002, 69,199997, 78,5, 78,050003, 66,199997, 70,300003, 84,699997, 78,5, 83, 64,5, 70,699997, 61,299999, 65,999 , 69.800003, 85.199997, 77.099998, 62.200001, 89.25 и т. Д.)
Ниже этой ошибки есть еще одна ошибка: в create_indexes_from_evaluated_names (curr_names, variable_names) (за этим также следуют числа, как указано выше - но разные числа сами по себе ).
Я дважды проверил, все ли переменные введены правильно, но мне не повезло.