Hi all, I am performing survival analysis with time-dependent covariates. Essentially, the values for methadone_1-methadone_12 are either OTP, OBOT, both, or none. Same for otp_methadone_bup_1-otp_methadone_bup_12. Non_med_1-non_med_12 are either 0 or 1. There were no missing values when I ran frequency tables for all of these! I continue getting the note that 606 event observations have missing values in the time-dependent covariates. Here is my proc phreg proc phreg data = outlib.EXAMPLE;
class age (ref='25-34') sex (ref='F') race (ref='White') ethnicity (ref='Not Hispanic')
rural (ref='Urban')
cancer (ref='0') chronic_kidney_disease (ref='0') chronic_lung_disease (ref='0')
chronic_liver_disease (ref='0') diabetes_mellitus (ref='0') heart_disease (ref='0')
obesity (ref='0') smoking (ref='0') meth_use_disorder (ref='0')
cocaine_use_disorder (ref='0');
model dx_survival_months*covid_dx(0)=
/*demographics*/
age race ethnicity rural
/*comorbidities/predisposing conditions*/
cancer chronic_kidney_disease chronic_lung_disease chronic_liver_disease diabetes_mellitus
heart_disease obesity smoking meth_use_disorder cocaine_use_disorder
/*time dependent vars*/
m otp non
/ ties=efron risklimits;
ARRAY methadone_(*) methadone_1-methadone_12;
m = methadone_[dx_survival_months];
ARRAY otp_methadone_bup_(*) otp_methadone_bup_1-otp_methadone_bup_12;
otp = otp_methadone_bup_[dx_survival_months];
ARRAY non_med_(*) non_med_1-non_med_12;
non = non_med_[dx_survival_months];
format age age.;
run; Any help is appreciated, and I can try to share more if this it is not clear. I'm not sure if it's an issue with how the time-dependent covariates are formatted in the data set, or what. Thanks!
... View more