Hello everyone,   I am working on survival data with incidence of cancer and some food consumption data in quartiles. Exposition : quartile of sugared food Confusion factor : age, BMI, sex. Outcome : cancer   I would like to draw a curve of cumulative hazard rate function (for example with Nelson Aalen estimator, a non parametric estimator of the cumulative hazard rate function) like this graph below : where we would see the cumulative risk of cancer in function of the 4 quartiles of food consumption. Caption of the graph : "Nelson-Aalen curves showing the outcome of new-onset overweight depending on ultraprocessed food consumption at baseline. Adjusted for sex, age, marital status, educational status"   If someone had an idea of a beginning of script, i would be really grateful. i am doing Cox model with this type of script : proc phreg data=mydata; class sexe; model survival_time*censure(1)=sexe age IMC; strata quartile_ultraprocessed ; run; Thanks in advance, have a nice day, ... View more

Hello everyone,   I am doing a Cox's model with delayed entry.     proc phreg data=food; class sexe_cl IPAQ1(ref="1") tabac_2classes; model(age_inclus,age_fin)*censure(1)=ultraprocessed sexe_cl IMC taille nb_enq Alc kcal IPAQ1 diplome_3cl tabac_2classes atcd_K0 / risklimits ties=efron; run ; age_inclus is the age at the inclusion of the participant. age_fin is the age at the minimum between the date of censure or the date of death or the date of last follow-up. ultraprocessed is my variable of interest. The other variables are explicative.   I would like to draw the graph of Martingale or schoenfeld residuals to see if the variables follow the proportional hazards hypothesis and i don't know how to dot it. Do you have an idea please ? Thanks in advance, Best regards,       ... View more