Do you know how to plot a case-cohort study where we perform a Cox regression to estimate HR of different covariates, in which I am specially interested in potential molecular biomarkers and its prognostic efficiency. The thing is I tried the available functions in the package ‘survminer’. It requires an object class “survfit” which doesn’t seem to support ‘cch’objects.

As I mentioned I wuld like to assess the prognostic value of several variables. The object comes from a databse with 189 controls and 55 cases. I assess multiple molecular markers, but I will paste the object of one particular and the generic formula applied

Thank you so much as always!

#case-cohort: cch function from the survival package

library(survival)

result <- cch(Surv(timetoevent, event) ~ age+ sex+ smoke+ box + weight + cholesterol+ cholesterol_HDL +  molecular_marker, data = dat, subcoh =  ~ casecohort2, id = ~parti, cohort.size = 5404, method = "LinYing", robust = TRUE)

The object “cch” in example