My population is divided into 10 groups. I want to run the cox-proportional hazard model separately for each of these 10 groups so that in the end i can compare which groups has highest hazard ratio depending uopn the covariates. How can i do this in R? Earlier i was using this formula.
Model_1 <- coxph(serve_object ~ GROUP + age+ Smoking+ Diabetes, data = DF summary(Model_1)
But the output is very confusing. it gives hazard ratio of 9 out of 10 groups. I think it dealt group as covariate.