Please submit code along with solutions.

### Question 1. (6 points)

Write a `R`

(or python or matlab) function to compute the Breslow partial log likelihood for the Cox proportional hazards model (log of Equation 8.4.1 in textbook) for a one-dimensional parameter \(\beta\). This function should take \(\beta\) as well as the data as arguments. Test your function by finding the MLE on the `kirc_small.RData`

data set using SLC7A11_ex as a predictor. You can use a grid search or an `R`

optimizer such as `optimize`

. Your MLE should match the parameter estimate here:

```
load("kirc_small.RData")
library(survival)
fit <- coxph(Surv(time,death_obs)~SLC7A11_ex,data=kirc_small,method="breslow")
coef(fit)
```

```
## SLC7A11_ex
## 0.002459945
```

You may use the solution from the above code to guide your starting value for the grid search or `optimize`

function.

### Question 2. (8 points=2 points/part)

Complete Question 8.10 in textbook (on page 291).