/****************************************************************/ /* S A S S A M P L E L I B R A R Y */ /* */ /* NAME: SVRGS */ /* TITLE: Getting Started Examples for PROC SURVEYREG */ /* PRODUCT: STAT */ /* SYSTEM: ALL */ /* KEYS: regression, survey sampling */ /* KEYS: stratification, clustering */ /* KEYS: unequal weighting, linear model */ /* PROCS: SURVEYREG */ /* DATA: */ /* */ /* SUPPORT: Tony An */ /* REF: PROC SURVEYREG, Getting Started */ /* */ /* MISC: */ /* */ /****************************************************************/ data IceCream; input Grade Spending Income Kids @@; datalines; 7 7 39 2 7 7 38 1 8 12 47 1 9 10 47 4 7 1 34 4 7 10 43 2 7 3 44 4 8 20 60 3 8 19 57 4 7 2 35 2 7 2 36 1 9 15 51 1 8 16 53 1 7 6 37 4 7 6 41 2 7 6 39 2 9 15 50 4 8 17 57 3 8 14 46 2 9 8 41 2 9 8 41 1 9 7 47 3 7 3 39 3 7 12 50 2 7 4 43 4 9 14 46 3 8 18 58 4 9 9 44 3 7 2 37 1 7 1 37 2 7 4 44 2 7 11 42 2 9 8 41 2 8 10 42 2 8 13 46 1 7 2 40 3 9 6 45 1 9 11 45 4 7 2 36 1 7 9 46 1 ; title1 'Ice Cream Spending Analysis'; title2 'Simple Random Sample Design'; proc surveyreg data=IceCream total=4000; class Kids; model Spending = Income Kids / solution; run; data StudentTotals; input Grade _TOTAL_; datalines; 7 1824 8 1025 9 1151 ; data IceCream; set IceCream; if Grade=7 then Prob=20/1824; if Grade=8 then Prob=9/1025; if Grade=9 then Prob=11/1151; Weight=1/Prob; run; ods graphics on; title1 'Ice Cream Spending Analysis'; title2 'Stratified Sample Design'; proc surveyreg data=IceCream total=StudentTotals; strata Grade /list; model Spending = Income; weight Weight; run; title1 'Ice Cream Spending Analysis'; title2 'Stratified Sample Design'; proc surveyreg data=IceCream total=StudentTotals; strata Grade /list; model Spending = Income; weight Weight; ods output ParameterEstimates = MyParmEst; run; proc print data=MyParmEst; run;