Hi all. I'm running a logistic regression for odds of receiving a skeletal survey in children less than 1 year of age admitted to the hospital for an accidental fall. My database is the National Trauma Data Bank, and there are over 5000 facilities included in my cohort. My model is the following: proc logistic data = logreg.alldata;
model sksurvey (event = '1')= mnatam_pi masian mblack moth_mix
mhisp maid mother mloc mISS/rsquare;
run; The first 5 variables are for race, maid = primary payment by medicaid, mother = other primary payment method, mloc = injury occurred outside the home, mISS = injury severity greater than 10. My Rsquare is very low (0.041), and I think this might be due to variation of the outcome (skeletal survey) by facility, for which I have a variable. I've never created a fixed effects model; would someone who is more familiar with this coding help point me in the right direction?
... View more