На данных (запись) я хочу проверить гипотезу о том, что времена заражения стафилококками одинаковы для двух дезинфицирующих групп с использованием модели, которая корректируется с учетом каждого из факторов (пол, раса, общая площадь сожженной поверхности, и тип ожога (химический, ожоговый, электрический, пламенный)).
Для этого мне нужно создать набор постоянных ковариат для каждого фактора (Z2, Z3, Z4 и Z11). Как я могу это сделать?
Я пытался:
library(KMsurv)
data(burn)
burn
library(survival)
fit.burn<-coxph(Surv(T3,D3)~Z1+Z2+Z3+Z4+Z11, method="breslow", data=burn)
fit.burn
summary(fit.burn)
Это не дает результатов, которые я ищу.