-
Notifications
You must be signed in to change notification settings - Fork 4
/
Copy pathhomework6_SAScode.sas
219 lines (209 loc) · 6.17 KB
/
homework6_SAScode.sas
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
* ==================================
* Homework 06 - logistic regression
*
* Melinda Higgins, PhD
* dated 11/25/2018
* ==================================
* ==================================
* we're be working with the
* helpmkh dataset
* ==================================;
* remember to CHANGE to the correct directory on your computer;
libname library 'C:\MyGithub\N736Fall2017_HELPdataset\' ;
* apply values;
proc format library = library ;
value TREAT
0 = 'usual care'
1 = 'HELP clinic' ;
value FEMALE
0 = 'Male'
1 = 'Female' ;
value HOMELESS
0 = 'no'
1 = 'yes' ;
value G1B
0 = 'no'
1 = 'yes' ;
value F1A
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1B
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1C
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1D
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1E
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1F
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1G
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1H
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1I
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1J
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1K
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1L
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1M
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1N
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1O
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1P
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1Q
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1R
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1S
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value F1T
0 = 'Not at all or less than 1 day'
1 = '1-2 days'
2 = '3-4 days'
3 = '5-7 days or nearly every day for 2 weeks' ;
value SATREAT
0 = 'no'
1 = 'yes' ;
value DRINKSTATUS
0 = 'no'
1 = 'yes' ;
value ANYSUBSTATUS
0 = 'no'
1 = 'yes' ;
value LINKSTATUS
0 = 'no'
1 = 'yes' ;
proc datasets library = library;
modify helpmkh / correctencoding="WLATIN1";
format treat TREAT.;
format female FEMALE.;
format homeless HOMELESS.;
format g1b G1B.;
format f1a F1A.;
format f1b F1B.;
format f1c F1C.;
format f1d F1D.;
format f1e F1E.;
format f1f F1F.;
format f1g F1G.;
format f1h F1H.;
format f1i F1I.;
format f1j F1J.;
format f1k F1K.;
format f1l F1L.;
format f1m F1M.;
format f1n F1N.;
format f1o F1O.;
format f1p F1P.;
format f1q F1Q.;
format f1r F1R.;
format f1s F1S.;
format f1t F1T.;
format satreat SATREAT.;
format drinkstatus DRINKSTATUS.;
format anysubstatus ANYSUBSTATUS.;
format linkstatus LINKSTATUS.;
quit;
* make a copy to WORK;
data helpmkh;
set library.helpmkh;
run;
* ============================================
* For this homework using the helpmkh dataset,
* You will be working with the e2b variable
* Number of times in past 6 months entered a
* detox program (collected at Baseline)
*
* For this logistic regression homework 6,
* I've provided the code below to capture the
* individuals who did NOT say they had entered a detox
* program in the 6 months preceeding baseline
*
* nodetox = 0 for individuals who did enter
a detox program and nodetox = 1 for individuals
who did NOT say they entered a detox program
* ============================================;
data help2;
set helpmkh;
nodetox = missing(e2b);
run;
* check results - the amount of missing for e2b
should equal the number of 1's for nodetox;
proc freq data=help2;
tables e2b nodetox;
run;
* alternatively you can use the / missing option
to get a table showing missing values for the 2 variables;
proc freq data=help2;
tables e2b * nodetox / missing;
run;
* ============================================
* For homework 6 focus on nodetox as the main outcome variable
* which is dichotomous coded 0 and 1. We'll use
* logistic regression to look at predicting whether someone
* was in detox or not using these variables
* age, female, pss_fr, pcs, mcs, and cesd
* ============================================;