I am running a median regression and a few of my independent variables such as severity of stroke is categorized into 5 categories. When I run proc quantreg, I do not get medians, IQR, number of observations for each separate subgroup(severity of stroke as none/minor/moder....). I would like to report this as there are missing values in my data set due to its retrospective nature, so I need to get information specific to the model. The question I am trying to answer is looking at a measure of stroke severity at discharge and how that depends on admission severity of stroke and other variables like age and so on. I am wondering how to get this data. My current code is below. Btw this is my first time using SAS 9.4.
proc quantreg data=filtereddata ci=resampling; class corticalinvolvement sex alcohol age_cat admission_nihss_score; model last_nihss_score = age_cat sex alcohol prophylaxis admission_nihss_score infarctsize corticalinvolvement / quantile = .5; estimate 'Diff in Medians' sex 1 -1 / CL; run;
It's not clear to me whether you want the univariate statistics (median, IQR,...) or whether you are asking for the dependent variable of the quantile regression.
If you want the univariate statistics, use PROC UNIVARIATE. You can use the BY statement if you want information for subgroups.
For the quantile regression, you can get a predicted (conditional) median for each value of the explanatory variables. Similarly, you can look at the 25th and 75th predicted quantiles to get a conditional IQR for each value of the explanatory variables. The value of the median and IQR depend on the value for the explanatory variables.
Secure your spot at the must-attend AI and analytics event of 2024: SAS Innovate 2024! Get ready for a jam-packed agenda featuring workshops, super demos, breakout sessions, roundtables, inspiring keynotes and incredible networking events.
Register by March 1 to snag the Early Bird rate of just $695! Don't miss out on this exclusive offer.