1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
70
71 /* MathReg2.sas */
72 %include '/home/brunner0/441s20/readmath2b.sas'; /* Creates data table mathex */
NOTE: Format YNFMT is already on the library WORK.FORMATS.
NOTE: Format YNFMT has been output.
NOTE: Format CRSFMT is already on the library WORK.FORMATS.
NOTE: Format CRSFMT has been output.
NOTE: Format NFMT is already on the library WORK.FORMATS.
NOTE: Format NFMT has been output.
NOTE: PROCEDURE FORMAT used (Total process time):
real time 0.00 seconds
user cpu time 0.00 seconds
system cpu time 0.00 seconds
memory 244.06k
OS Memory 31652.00k
Timestamp 02/01/2020 09:52:11 PM
Step Count 43 Switch Count 0
Page Faults 0
Page Reclaims 25
Page Swaps 0
Voluntary Context Switches 0
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 32
176 title2 'Residual analysis';
177
NOTE: The infile '/home/brunner0/441s20/exploremath.data.txt' is:
Filename=/home/brunner0/441s20/exploremath.data.txt,
Owner Name=brunner0,Group Name=oda,
Access Permission=-rw-r--r--,
Last Modified=26Jan2020:18:49:34,
File Size (bytes)=44583
NOTE: 579 records were read from the infile '/home/brunner0/441s20/exploremath.data.txt'.
The minimum record length was 75.
The maximum record length was 75.
NOTE: Missing values were generated as a result of performing an operation on missing values.
Each place is given by: (Number of times) at (Line):(Column).
99 at 95:24 99 at 134:13
NOTE: The data set WORK.MATHEX has 579 observations and 34 variables.
NOTE: DATA statement used (Total process time):
real time 0.00 seconds
user cpu time 0.00 seconds
system cpu time 0.00 seconds
memory 880.62k
OS Memory 32168.00k
Timestamp 02/01/2020 09:52:11 PM
Step Count 44 Switch Count 3
Page Faults 0
Page Reclaims 133
Page Swaps 0
Voluntary Context Switches 25
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 520
178 proc reg plots(only) = ResidualPlot data=mathex;
179 title3 'Model H: hsgpa hscalc hsengl totscore';
180 model grade = hsgpa hscalc hsengl totscore;
181 output out = Explor predicted = yhat
182 residual = resid
183 rstudent = delstud;
184 /* Deleted Studentized Residual */
185
186 /* Could have included LCL and UCL for upper and lower limits of a
187 95% prediction interval for each case in the file */
188
189
190 /* What is a big (Studentized deleted) residual? If the model is correct,
191 each one has a t distribution with n-p-1 = 283 df (practically standard
192 normal), so the Studentized deleted residual can be treated directly as
193 a t-test statistic. Values that are too big in absolute value will cause
194 rejection of the null hypothesis that x_i*beta is the same for this case
195 and the other n-1 cases. Tests are NOT independent, but use a Bonferroni
196 correction for n = 289 tests. Get the critical value from proc iml. */
197
NOTE: The data set WORK.EXPLOR has 579 observations and 37 variables.
NOTE: PROCEDURE REG used (Total process time):
real time 0.37 seconds
user cpu time 0.16 seconds
system cpu time 0.02 seconds
memory 19480.59k
OS Memory 47556.00k
Timestamp 02/01/2020 09:52:11 PM
Step Count 45 Switch Count 4
Page Faults 0
Page Reclaims 4195
Page Swaps 0
Voluntary Context Switches 273
Involuntary Context Switches 1
Block Input Operations 0
Block Output Operations 1600
198 proc iml;
NOTE: IML Ready
199 title3 'Critical value for Joint t-test on Studentized Residuals';
200 Alpha = 0.05/289;
200 ! print Alpha;
201 Critval = tinv(1-Alpha/2,283);
201 ! print Critval;
202
NOTE: Exiting IML.
NOTE: PROCEDURE IML used (Total process time):
real time 0.03 seconds
user cpu time 0.04 seconds
system cpu time 0.00 seconds
memory 1010.21k
OS Memory 45952.00k
Timestamp 02/01/2020 09:52:12 PM
Step Count 46 Switch Count 1
Page Faults 0
Page Reclaims 267
Page Swaps 0
Voluntary Context Switches 10
Involuntary Context Switches 1
Block Input Operations 0
Block Output Operations 176
203 proc univariate normal plot;
204 var delstud;
205
206 /* Tests for normality indicate residuals are not normal. No st resids
207 greater than critical value. Next, a few more plots. */
208
NOTE: PROCEDURE UNIVARIATE used (Total process time):
real time 0.27 seconds
user cpu time 0.15 seconds
system cpu time 0.00 seconds
memory 3807.84k
OS Memory 47452.00k
Timestamp 02/01/2020 09:52:12 PM
Step Count 47 Switch Count 0
Page Faults 0
Page Reclaims 567
Page Swaps 0
Voluntary Context Switches 292
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 528
209 proc sgplot;
210 title3 'Plot of Y-hat by Y';
211 scatter y=grade x=yhat;
212
213 proc sgplot;
NOTE: PROCEDURE SGPLOT used (Total process time):
real time 0.15 seconds
user cpu time 0.06 seconds
system cpu time 0.01 seconds
memory 1531.93k
OS Memory 47408.00k
Timestamp 02/01/2020 09:52:12 PM
Step Count 48 Switch Count 2
Page Faults 0
Page Reclaims 330
Page Swaps 0
Voluntary Context Switches 197
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 536
NOTE: There were 579 observations read from the data set WORK.EXPLOR.
214 title3 'Calculus sub-test by deleted studentized residual';
215 scatter x=calc y=delstud;
216
217 proc sgplot;
NOTE: PROCEDURE SGPLOT used (Total process time):
real time 0.14 seconds
user cpu time 0.05 seconds
system cpu time 0.00 seconds
memory 2551.56k
OS Memory 47920.00k
Timestamp 02/01/2020 09:52:12 PM
Step Count 49 Switch Count 2
Page Faults 0
Page Reclaims 392
Page Swaps 0
Voluntary Context Switches 196
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 560
NOTE: There were 579 observations read from the data set WORK.EXPLOR.
218 title3 'Pre-calculus sub-test by deleted studentized residual';
219 scatter x=precalc y=delstud;
220
221 proc sgplot;
NOTE: PROCEDURE SGPLOT used (Total process time):
real time 0.12 seconds
user cpu time 0.06 seconds
system cpu time 0.01 seconds
memory 2281.15k
OS Memory 47920.00k
Timestamp 02/01/2020 09:52:12 PM
Step Count 50 Switch Count 2
Page Faults 0
Page Reclaims 305
Page Swaps 0
Voluntary Context Switches 198
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 552
NOTE: There were 579 observations read from the data set WORK.EXPLOR.
222 title3 'Mother tongue by deleted studentized residual';
223 scatter x=mtongue y=delstud;
224
225 quit;
NOTE: PROCEDURE SGPLOT used (Total process time):
real time 0.14 seconds
user cpu time 0.05 seconds
system cpu time 0.00 seconds
memory 2277.37k
OS Memory 47920.00k
Timestamp 02/01/2020 09:52:12 PM
Step Count 51 Switch Count 2
Page Faults 0
Page Reclaims 289
Page Swaps 0
Voluntary Context Switches 197
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 536
NOTE: There were 579 observations read from the data set WORK.EXPLOR.
226
227 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
238