I'm trying do multilevel regression and poststratification based on this paper: http://aje.oxfordjournals.org/content/179/8/1025.short I stuck the model prediction step now. I found their lecture note online (online note: http://www.cste2.org/webinars/webinarpdfs/GeospatialDataMethods82313.pdf) and try to modify it because I focused on county-level not census track block. My codes as following: proc GLIMMIX data=brfss.b12c_sr noclprint; class FIPS_C state Ragegrp sex Rrace4 Rrurality2 ad_mean; model MH (EVENT='1') = Ragegrp sex Rrace4 Rrurality2 ad_mean /dist=binary solution; random state fips_c(state) / solution; weight rw; ods output parameterEstimates=beta_fixed SolutionR=beta_random; run; %MACRO SIMU(outcome=, year=); %do i=1 %to 1000; data temp1; set beta&year._&outcome; xbeta=int+int_se*rannor(ABS(INT(RAND('UNIFORM')*100000))) +agegrp+age_se*rannor(ABS(INT(RAND('UNIFORM')*100000))) +sexgrp+sex_se*rannor(ABS(INT(RAND('UNIFORM')*100000))) +race4+race4_se*rannor(ABS(INT(RAND('UNIFORM')*100000))) +fips_c+fips_c_se*rannor(ABS(INT(RAND('UNIFORM')*100000))) +state+state_se*rannor(ABS(INT(RAND('UNIFORM')*100000))) +ad_mean*(ad_mean+ad_mean_se*rannor(ABS(INT(RAND('UNIFORM')*100000)))) ; p=exp(xbeta)/(1+exp(xbeta)); keep fips_c pop p; run; *summarize obesity by block (or county, or congressional district); proc sql; create table temp2 as select distinct fips_c, sum(pop) as pop, (sum(p*pop)/sum(pop))*100 as &outcome from temp1 group by fips_C; quit; proc append base=temp3 data=temp2 force;run; %MEND; However, it shows error: ERROR: There were 1 unclosed %DO statements. The macro SIMU will not be compiled. ERROR: A dummy macro will be compiled. Does anyone ever done it before or know how should I modify it? Thank you so much!
... View more