Computing the survival curve from a Cox model

In this exercise, we will reproduce the example from the video following the steps:

Compute Cox model
Decide on "imaginary patients"
Compute survival curves
Create data.frame with survival curve information
Plot

We will focus now on the first three steps in this exercise and do the next two steps in the upcoming exercise.

The survival and survminer packages and the GBSG2 data are loaded for you in this exercise.

This exercise is part of the course

Survival Analysis in R

View Course

Hands-on interactive exercise

Have a go at this exercise by completing this sample code.

# Cox model
cxmod <- ___(Surv(time, cens) ~ horTh + tsize, data = GBSG2)

# Imaginary patients
newdat <- expand.grid(
  horTh = levels(GBSG2$___),
  tsize = quantile(GBSG2$___, probs = c(0.25, 0.5, 0.75)))
___(newdat) <- letters[1:6]

# Inspect newdat
newdat

Edit and Run Code