-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathGMASQAA.htm
155 lines (97 loc) · 2.73 KB
/
GMASQAA.htm
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
LIBNAME IN 'C:\Users\';
DATA TEMPX; SET IN.data1;
KEEP ID Gender G ZSHS_TOT_I MASQ_AA_I;
PROC SORT; BY ID;
RUN;
PROC STANDARD DATA=TEMPX OUT=STDTEMPX MEAN=0 STD=1;
VAR MASQ_AA_I G;
RUN;
PROC SORT; BY ID;
RUN;
DATA FU0; SET IN.data1;
KEEP ID TIME MASQ_AA_I ZSHS_TOT_I;
TIME = 0;
RENAME MASQ_AA_I= CESDRESSION;
RENAME ZSHS_TOT_I = STRESS;
PROC SORT; BY ID;
RUN;
DATA FU1; SET IN.data1;
KEEP ID TIME MASQ_AA_1 ZSHS_TOT_1;
TIME = 1;
RENAME MASQ_AA_1 = CESDRESSION;
RENAME ZSHS_TOT_1 = STRESS;
PROC SORT; BY ID;
RUN;
DATA FU2; SET IN.data1;
KEEP ID TIME MASQ_AA_2 ZSHS_TOT_2;
TIME = 2;
RENAME MASQ_AA_2 = CESDRESSION;
RENAME ZSHS_TOT_2 = STRESS;
PROC SORT; BY ID;
RUN;
DATA FU3; SET IN.data1;
KEEP ID TIME MASQ_AA_3 ZSHS_TOT_3;
TIME = 3;
RENAME MASQ_AA_3 = CESDRESSION;
RENAME ZSHS_TOT_3 = STRESS;
PROC SORT; BY ID;
RUN;
DATA FU4; SET IN.data1;
KEEP ID TIME MASQ_AA_4 ZSHS_TOT_4;
TIME = 4;
RENAME MASQ_AA_4 = CESDRESSION;
RENAME ZSHS_TOT_4 = STRESS;
PROC SORT; BY ID;
RUN;
DATA FU5; SET IN.data1;
KEEP ID TIME MASQ_AA_5 ZSHS_TOT_5;
TIME = 5;
RENAME MASQ_AA_5 = CESDRESSION;
RENAME ZSHS_TOT_5 = STRESS;
PROC SORT; BY ID;
RUN;
DATA FU6; SET IN.data1;
KEEP ID TIME MASQ_AA_6 ZSHS_TOT_6;
TIME = 6;
RENAME MASQ_AA_6 = CESDRESSION;
RENAME ZSHS_TOT_6 = STRESS;
PROC SORT; BY ID;
RUN;
DATA TEMPY; MERGE FU1 FU2 FU3 FU4 FU5 FU6; BY ID TIME;
PROC SORT; BY ID TIME;
RUN;
DATA TEMPZ; MERGE STDTEMPX TEMPY; BY ID;
IF CESDRESSION = . THEN DELETE;
IF STRESS = . THEN DELETE;
PROC SORT; BY ID TIME;
RUN;
PROC MEANS MEAN DATA=TEMPZ NOPRINT; BY ID;
VAR CESDRESSION STRESS;
OUTPUT OUT=TOTAL MEAN= MN_CESDRESSION MN_STRESS;
RUN;
PROC STANDARD DATA=TOTAL MEAN=0 STD=1 OUT=STDTOTAL;
VAR MN_CESDRESSION MN_STRESS;
RUN;
DATA FINAL; MERGE TEMPZ STDTOTAL; BY ID;
PROC STANDARD MEAN=0 DATA=FINAL OUT=FINALY; BY ID;
VAR STRESS;
RUN;
PROC MIXED NOCLPRINT METHOD=ML COVTEST NOITPRINT;
CLASS ID TIME;
MODEL CESDRESSION = Gender MASQ_AA_I STRESS G STRESS*G
/SOLUTION DDFM=BW;
RANDOM INTERCEPT STRESS/SUB=ID TYPE=UN;
repeated time /type=ARH(1) subject=id rcorr;
estimate 'mean slope' STRESS 1 ;
estimate 'slope for high G' STRESS 1 G*STRESS 1.5;
estimate 'slope for low G' STRESS 1 G*STRESS -1.5;
estimate 'slope h vs slope l' G*STRESS 3;
estimate ' h G, high STRESS' int 1 MASQ_AA_I .48 STRESS 0.48 G 1.5
G*STRESS 0.72;
estimate ' h G, low STRESS' int 1 MASQ_AA_I .48 STRESS -0.48 G 1.5
G*STRESS -0.72;
estimate ' l G, high STRESS' int 1 MASQ_AA_I -.47 STRESS 0.48 G -1.5
G*STRESS -0.72;
estimate ' l G, low STRESS' int 1 MASQ_AA_I -.47 STRESS -0.48 G -1.5
G*STRESS 0.72;
RUN;