1 The SAS System 11:14 Sunday, November 21, 2004
NOTE: Copyright (c) 1999-2001 by SAS Institute Inc., Cary, NC, USA.
NOTE: SAS (r) Proprietary Software Release 8.2 (TS2M0)
Licensed to AGENCY FOR HEALTHCARE RESEARCH & QUALITY, Site 0040776001.
NOTE: This session is executing on the WIN_PRO platform.
NOTE: This installation is running Base SAS hot fix bundle 82BX08.
NOTE: AUTOEXEC processing beginning; file is C:\PROGRA~1\SASINS~1\SAS\V8\autoexec.sas.
1 ** C:\PROGRAM FILES\SAS INSTITUTE\SAS\V8\AUTOEXEC.SAS ;
2
3 OPTIONS LS= 132 PS= 70 MERGENOBY=ERROR MSGLEVEL=i;
NOTE: AUTOEXEC processing completed.
1 /****************************************************************\
2
3 PROGRAM: C:\MEPS\PROG\EXAMPLE_EM1.SAS
4
5 DESCRIPTION: THIS EXAMPLE SHOWS HOW TO BUILD AN ANALYTIC FILE
6 AND CREATE NEW VARIABLES TO EXAMINE THE
7 RELATIONSHIP BETWEEN PERCEIVED HEALTH STATUS
8 AND A PERSON'S WEEKLY EARNINGS, WHERE WEEKLY
9 EARNINGS ARE THOSE OF THE CURRENT MAIN JOB.
10
11 PERSON-LEVEL RECORDS ARE DIVIDED INTO QUARTILES
12 BASED ON WEEKLY EARNINGS AND PERSON-LEVEL WEIGHTS.
13 THE RESULT IS 4 EQUALLY WEIGHTED QUARTILES FOR
14 WEEKLY EARNINGS.
15
16 INPUT FILE: C:\MEPS\DATA\H62.SAS7BDAT (2002 FULL-YEAR POPULATION)
17 -- PERSON-LEVEL FILE
18
19 \****************************************************************/
20
21 LIBNAME CMEPS V8 'C:\MEPS\DATA' ;
NOTE: Libref CMEPS was successfully assigned as follows:
Engine: V8
Physical Name: C:\MEPS\DATA
22
23 FOOTNOTE 'PROGRAM: C:\MEPS\PROG\EXAMPLE_EM1.SAS';
24
25 TITLE1 'AHRQ MEPS DATA USERS WORKSHOP (EMPLOYMENT) -- NOV/DEC 2004';
26 TITLE2 'PERCEIVED HEALTH STATUS AND WEEKLY EARNINGS - 2002 DATA';
27 TITLE3 ' ';
28
29 PROC FORMAT;
30 VALUE GTZF
31 -10 = '-10'
32 0<-HIGH = '>0';
NOTE: Format GTZF has been output.
33 VALUE HLTHF
34 1 = '1 EXCELLENT'
35 2 = '2 VERY GOOD'
36 3 = '3 GOOD'
37 4 = '4 FAIR'
38 5 = '5 POOR';
NOTE: Format HLTHF has been output.
39 VALUE QUARTF
40 1 = ' 1 LOWEST '
41 2 = ' 2 '
42 3 = ' 3 '
43 4 = ' 4 HIGHEST';
NOTE: Format QUARTF has been output.
44 VALUE WKLYF
2 The SAS System 11:14 Sunday, November 21, 2004
45 0.13-<49.99 = ' 0.13 - 49.99'
46 50-499.99 = ' 50.00 - 499.99'
47 500-999.99 = ' 500.00 - 999.99'
48 1000-HIGH = '1,000.00+ ';
NOTE: Format WKLYF has been output.
49 RUN;
50
51 /***** THIS DATA STEP READS IN THE REQUIRED VARIABLES FROM THE *****/
52 /***** FULL-YEAR POPULATION FILE (HC-062). *****/
53
54 /***** THE 'WHERE' STATEMENT SUBSETS TO THE DESIRED SUBSET OF *****/
55 /***** PERSONS: THOSE WITH A CURRENT MAIN JOB (CMJ) IN ROUNDS *****/
56 /***** 4/2 THAT REPORTED EARNINGS, WEEKLY HOURS AND PERCEIVED *****/
57 /***** HEALTH STATUS. *****/
58
59 /***** A VALUE OF "-2" FOR THE '42' VARIABLES INDICATES THAT *****/
60 /***** THE ACTUAL VALUE CAN BE OBTAINED FROM THE PRIOR ROUND, *****/
61 /***** THE '31' VARIABLES. *****/
62
63 /***** A WEEKLY EARNINGS VARIABLE "WKLYEARN' IS CREATED AS *****/
64 /***** THE PRODUCE OF HOURLY WAGE ("HRLYWAGE") TIMES NUMBER *****/
65 /***** OF WEEKLY WORK HOURS ("HOURS"). *****/
66
67 DATA H62A;
68 SET CMEPS.H62 (KEEP= DUPERSID PERWT02P VARSTR VARPSU EMPST42
69 SELFCM42 SELFCM31 HRWG42X HRWG31X HOUR42 HOUR31
70 RTHLTH42 TEMPJB42 AGE42X);
71 WHERE (EMPST42 IN (1,2)) AND
72 (SELFCM42 = 2 OR (SELFCM42 = -2 AND SELFCM31 = 2)) AND
73 ((HRWG42X = -10) OR (HRWG42X > 0) OR
74 (HRWG42X = -2 AND (HRWG31X = -10 OR HRWG31X > 0))) AND
75 (HOUR42 > 0 OR (HOUR42 = -2 AND HOUR31 > 0)) AND
76 (RTHLTH42 IN (1,2,3,4,5)) AND
77 (TEMPJB42 = 2) AND
78 (AGE42X > 24);
79 IF HRWG42X = -2
80 THEN HRLYWAGE = HRWG31X;
81 ELSE HRLYWAGE = HRWG42X;
82 IF HRLYWAGE = -10
83 THEN HRLYWAGE = 61.98;
84 IF HOUR42 > 0
85 THEN HOURS = HOUR42;
86 ELSE IF HOUR42 = -2
87 THEN HOURS = HOUR31;
88 WKLYEARN = HRLYWAGE*HOURS;
89 RUN;
NOTE: There were 11882 observations read from the data set CMEPS.H62.
WHERE EMPST42 in (1, 2) and ((SELFCM42=2) or ((SELFCM42=-2) and (SELFCM31=2))) and ((HRWG42X=-10) or (HRWG42X>0) or
((HRWG42X=-2) and ((HRWG31X=-10) or (HRWG31X>0)))) and ((HOUR42>0) or ((HOUR42=-2) and (HOUR31>0))) and RTHLTH42 in (1, 2, 3,
4, 5) and (TEMPJB42=2) and (AGE42X>24);
NOTE: The data set WORK.H62A has 11882 observations and 17 variables.
90
91 PROC SORT DATA= H62A;
92 BY WKLYEARN ;
93 RUN;
NOTE: SAS sort was used.
NOTE: There were 11882 observations read from the data set WORK.H62A.
NOTE: The data set WORK.H62A has 11882 observations and 17 variables.
94
95 /***** THIS DATA STEP READS IN H62A TWICE. THE FIRST TIME IS *****/
96 /***** TO COMPUTE THE TOTAL SUM OF PERWT02P (THE PERSON-LEVEL *****/
97 /***** WEIGHT VARIABLE). THE SECOND TIME IS TO CREATE 4 GROUPS *****/
3 The SAS System 11:14 Sunday, November 21, 2004
98 /***** OF EQUAL WEIGHT: QUARTILE = 1-4. SEE THE TABLE ON LST *****/
99 /***** P. 5. *****/
100
101
102 DATA H62B (KEEP= RTHLTH42 WKLYEARN QUARTILE PERWT02P VARSTR VARPSU);
103 SET H62A (IN= A) H62A ;
104 IF A
105 THEN TOT_WT+PERWT02P;
106 ELSE DO;
107 SUM_WT+PERWT02P;
108 IF SUM_WT <= TOT_WT*0.25
109 THEN QUARTILE = 1;
110 ELSE IF SUM_WT <= TOT_WT*0.5
111 THEN QUARTILE = 2;
112 ELSE IF SUM_WT <= TOT_WT*0.75
113 THEN QUARTILE = 3;
114 ELSE QUARTILE = 4;
115 OUTPUT;
116 END;
117 LABEL QUARTILE = 'WEEKLY EARNINGS QUARTILE'
118 WKLYEARN = 'WEEKLY EARNINGS';
119 RUN;
NOTE: There were 11882 observations read from the data set WORK.H62A.
NOTE: There were 11882 observations read from the data set WORK.H62A.
NOTE: The data set WORK.H62B has 11882 observations and 6 variables.
120
121 TITLE4 'H62B: WEIGHTED FREQUENCY OF QUARTILES AND WKLYEARN';
122
123 PROC FREQ DATA= H62B;
124 TABLES QUARTILE WKLYEARN;
125 WEIGHT PERWT02P;
126 FORMAT WKLYEARN WKLYF. ;
127 RUN;
NOTE: There were 11882 observations read from the data set WORK.H62B.
NOTE: The PROCEDURE FREQ printed page 1.
128
129 TITLE4 'H62B: UNWEIGHTED FREQUENCY OF WKLYEARN';
130
131 PROC FREQ DATA= H62B;
132 TABLES WKLYEARN;
133 FORMAT WKLYEARN WKLYF. ;
134 RUN;
NOTE: There were 11882 observations read from the data set WORK.H62B.
NOTE: The PROCEDURE FREQ printed page 2.
135
136 TITLE4 ' ';
137
138 PROC FREQ DATA= H62B;
139 TABLES QUARTILE*RTHLTH42
140 / LIST MISSING ;
141 FORMAT RTHLTH42 HLTHF. ;
142 RUN;
NOTE: There were 11882 observations read from the data set WORK.H62B.
NOTE: The PROCEDURE FREQ printed page 3.
143
144 PROC TABULATE DATA= H62B;
145 CLASS QUARTILE RTHLTH42;
146 TABLE QUARTILE,RTHLTH42*(PCTN<RTHLTH42>);
147 FORMAT QUARTILE QUARTF. RTHLTH42 HLTHF. ;
148 FREQ PERWT02P;
4 The SAS System 11:14 Sunday, November 21, 2004
149 WEIGHT PERWT02P;
150 RUN;
NOTE: There were 11882 observations read from the data set WORK.H62B.
NOTE: The PROCEDURE TABULATE printed page 4.
151
152 TITLE4 'MEAN REPORTED HEALTH STATUS BY WEEKLY EARNINGS QUARTILE';
153 TITLE5 '1=EXCELLENT, 2=VERY GOOD, 3=GOOD, 4=FAIR, 5=POOR';
154
155 PROC SURVEYMEANS DATA= H62B NOBS SUMWGT MEAN STDERR;
156 VAR RTHLTH42;
157 STRATA VARSTR;
158 CLUSTER VARPSU;
159 WEIGHT PERWT02P;
160 DOMAIN QUARTILE;
161 FORMAT QUARTILE QUARTF. ;
162 RUN;
NOTE: Due to nonpositive weights, 429 observation(s) were deleted.
NOTE: Only one cluster in a stratum in domain QUARTILE for variable(s) RTHLTH42. The variance of RTHLTH42 in that stratum is
estimated by zero.
NOTE: There were 11882 observations read from the data set WORK.H62B.
NOTE: The PROCEDURE SURVEYMEANS printed page 5.
NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414
AHRQ MEPS DATA USERS WORKSHOP (EMPLOYMENT) -- NOV/DEC 2004 1
PERCEIVED HEALTH STATUS AND WEEKLY EARNINGS - 2002 DATA 11:14 Sunday, November 21, 2004
H62B: WEIGHTED FREQUENCY OF QUARTILES AND WKLYEARN
The FREQ Procedure
WEEKLY EARNINGS QUARTILE
Cumulative Cumulative
QUARTILE Frequency Percent Frequency Percent
----------------------------------------------------------------
1 24,691,520 25.0 24,691,520 25.0
2 24,691,126 25.0 49,382,646 50.0
3 24,687,804 25.0 74,070,450 75.0
4 24,696,042 25.0 98,766,492 100.0
WEEKLY EARNINGS
Cumulative Cumulative
WKLYEARN Frequency Percent Frequency Percent
-------------------------------------------------------------------------
0.13 - 49.99 1,557,396 1.6 1,557,396 1.6
50.00 - 499.99 38,403,334 38.9 39,960,730 40.5
500.00 - 999.99 36,501,086 37.0 76,461,816 77.4
1,000.00+ 22,304,676 22.6 98,766,492 100.0
PROGRAM: C:\MEPS\PROG\EXAMPLE_EM1.SAS
AHRQ MEPS DATA USERS WORKSHOP (EMPLOYMENT)
-- NOV/DEC 2004 2
PERCEIVED HEALTH STATUS AND WEEKLY EARNINGS - 2002 DATA 11:14 Sunday, November 21, 2004
H62B: UNWEIGHTED FREQUENCY OF WKLYEARN
The FREQ Procedure
WEEKLY EARNINGS
Cumulative Cumulative
WKLYEARN Frequency Percent Frequency Percent
-------------------------------------------------------------------------
0.13 - 49.99 199 1.7 199 1.7
50.00 - 499.99 5,308 44.7 5,507 46.3
500.00 - 999.99 4,124 34.7 9,631 81.1
1,000.00+ 2,251 18.9 11,882 100.0
PROGRAM: C:\MEPS\PROG\EXAMPLE_EM1.SAS
AHRQ MEPS DATA USERS WORKSHOP (EMPLOYMENT) -- NOV/DEC 2004 3
PERCEIVED HEALTH STATUS AND WEEKLY EARNINGS - 2002 DATA 11:14 Sunday, November 21, 2004
The FREQ Procedure
Cumulative Cumulative
QUARTILE RTHLTH42 Frequency Percent Frequency Percent
----------------------------------------------------------------------------
1 1 EXCELLENT 733 6.17 733 6.17
1 2 VERY GOOD 1191 10.02 1924 16.19
1 3 GOOD 1139 9.59 3063 25.78
1 4 FAIR 398 3.35 3461 29.13
1 5 POOR 50 0.42 3511 29.55
2 1 EXCELLENT 706 5.94 4217 35.49
2 2 VERY GOOD 1159 9.75 5376 45.24
2 3 GOOD 956 8.05 6332 53.29
2 4 FAIR 269 2.26 6601 55.55
2 5 POOR 36 0.30 6637 55.86
3 1 EXCELLENT 742 6.24 7379 62.10
3 2 VERY GOOD 1092 9.19 8471 71.29
3 3 GOOD 735 6.19 9206 77.48
3 4 FAIR 170 1.43 9376 78.91
3 5 POOR 20 0.17 9396 79.08
4 1 EXCELLENT 782 6.58 10178 85.66
4 2 VERY GOOD 990 8.33 11168 93.99
4 3 GOOD 582 4.90 11750 98.89
4 4 FAIR 116 0.98 11866 99.87
4 5 POOR 16 0.13 11882 100.00
PROGRAM: C:\MEPS\PROG\EXAMPLE_EM1.SAS
AHRQ MEPS DATA USERS WORKSHOP (EMPLOYMENT) -- NOV/DEC 2004 4
PERCEIVED HEALTH STATUS AND WEEKLY EARNINGS - 2002 DATA 11:14 Sunday, November 21, 2004
--------------------------------------------------------------------------------------------------
| | PERCEIVED HEALTH STATUS - RD 4/2 |
| |----------------------------------------------------------------|
| |1 EXCELLENT |2 VERY GOOD | 3 GOOD | 4 FAIR | 5 POOR |
| |------------+------------+------------+------------+------------|
| | PctN | PctN | PctN | PctN | PctN |
|-------------------------------+------------+------------+------------+------------+------------|
|WEEKLY EARNINGS QUARTILE | | | | | |
|-------------------------------| | | | | |
|1 LOWEST | 21.80| 36.18| 31.06| 9.64| 1.33|
|-------------------------------+------------+------------+------------+------------+------------|
|2 | 22.96| 38.52| 29.34| 8.18| 0.99|
|-------------------------------+------------+------------+------------+------------+------------|
|3 | 27.51| 40.09| 26.20| 5.67| 0.54|
|-------------------------------+------------+------------+------------+------------+------------|
|4 HIGHEST | 32.31| 39.99| 22.61| 4.43| 0.66|
--------------------------------------------------------------------------------------------------
PROGRAM: C:\MEPS\PROG\EXAMPLE_EM1.SAS
AHRQ MEPS DATA USERS WORKSHOP (EMPLOYMENT) -- NOV/DEC 2004 5
PERCEIVED HEALTH STATUS AND WEEKLY EARNINGS - 2002 DATA 11:14 Sunday, November 21, 2004
MEAN REPORTED HEALTH STATUS BY WEEKLY EARNINGS QUARTILE
1=EXCELLENT, 2=VERY GOOD, 3=GOOD, 4=FAIR, 5=POOR
The SURVEYMEANS Procedure
Data Summary
Number of Strata 203
Number of Clusters 445
Number of Observations 11882
Number of Observations Used 11453
Number of Obs with Nonpositive Weights 429
Sum of Weights 98766491.9
Statistics
Sum of Std Error
Variable Label N Weights Mean of Mean
------------------------------------------------------------------------------------------------------------
RTHLTH42 PERCEIVED HEALTH STATUS - RD 4/2 11453 98766492 2.177565 0.012317
------------------------------------------------------------------------------------------------------------
Domain Analysis: QUARTILE
Sum of Std Error
QUARTILE Variable Label N Weights Mean of Mean
----------------------------------------------------------------------------------------------------------------------------
1 LOWEST RTHLTH42 PERCEIVED HEALTH STATUS - RD 4/2 3394 24691520 2.325160 0.021478
2 RTHLTH42 PERCEIVED HEALTH STATUS - RD 4/2 2990 24691126 2.257211 0.022077
3 RTHLTH42 PERCEIVED HEALTH STATUS - RD 4/2 2660 24687804 2.116474 0.019022
4 HIGHEST RTHLTH42 PERCEIVED HEALTH STATUS - RD 4/2 2409 24696042 2.011436 0.019907
----------------------------------------------------------------------------------------------------------------------------
PROGRAM: C:\MEPS\PROG\EXAMPLE_EM1.SAS