01>LBL "MLRAN" 02 STO "MX" 03 TRANS 04 X<>Y 05 STO "MY" 06 × 07 RCL "MX" 08 TRANS 09 LASTX 10 × 11 ÷ 12 STO "MC" # Calculate coefficients of regression 13 RCL "MX" 14 DIM? 15 1 16 - 17 STO "K" # Store number of independent variables 18 X<>Y 19 STO "N" # Store number of observations 20 RCL "MY" 21 TRANS 22 RSUM 23 XEQ 00 24 X^2 25 RCL "N" 26 ÷ 27 STO "SY" # Store Sum(y)^2 / N 28 RCL "MC" 29 TRANS 30 RCL "MX" 31 TRANS 32 × 33 RCL "MY" 34 × 35 XEQ 00 36 RCL "SY" 37 - 38 STO "SSR" # Calculate SSR 39 RCL "MY" 40 TRANS 41 RCL "MY" 42 × 43 XEQ 00 44 RCL "SY" 45 - 46 STO "SST" # Calculate SST 47 RCL "SSR" 48 - 49 STO "SSE" # Calculate SSE 50 RCL "SSR" 51 RCL "SST" 52 ÷ 53 STO "R2" # Calculate R2 54 +/- 55 1 56 + 57 RCL "N" 58 1 59 - 60 × 61 RCL "N" 62 1 63 - 64 RCL "K" 65 - 66 ÷ 67 +/- 68 1 69 + 70 STO "R2ADJ" # Calculate R2adj 71 RCL "K" 72 STO "DFR" # Calculate DF residuals 73 +/- 74 1 75 - 76 RCL "N" 77 + 78 STO "DFE" # Calculate DF error 79 RCL "K" 80 + 81 STO "DFT" # Calculate total DF 82 RCL "SSR" 83 RCL "DFR" 84 ÷ 85 STO "MSR" # Calculate MSR 86 RCL "SSE" 87 RCL "DFE" 88 ÷ 89 STO "MSE" # Calculate MSE 90 ÷ 91 STO "FSTAT" # Calculate F stat 92 PRV "MC" # Print results 93 PRV "R2" 94 PRV "R2ADJ" 95 "ANOVA TABLE" 96 PRA 97 PRV "SSR" 98 PRV "SSE" 99 PRV "SST" 100 PRV "DFR" 101 PRV "DFE" 102 PRV "DFT" 103 PRV "MSR" 104 PRV "MSE" 105 PRV "FSTAT" 106 RCL "R2" 107 RTN 108>LBL 00 109 STO "MM" 110 INDEX "MM" 111 RCLEL 112 RTN 113>LBL "CLRANV" 114 CLV "SSR" 115 CLV "MM" 116 CLV "K" 117 CLV "SY" 118 CLV "N" 119 CLV "SSE" 120 CLV "SST" 121 CLV "DFR" 122 CLV "DFE" 123 CLV "DFT" 124 CLV "MSR" 125 CLV "MSE" 126 CLV "FSTAT" 127 CLV "R2" 128 CLV "R2ADJ" 129 .END. 130>LBL "VWRS" 131 RCL "MC" 132 EDIT 133 STOP 134 "R2=" 135 ARCL "R2" 136 PROMPT 137 "R2ADJ=" 138 ARCL "R2ADJ" 139 PROMPT 140 "SSR=" 141 ARCL "SSR" 142 PROMPT 143 "SSE=" 144 ARCL "SSE" 145 PROMPT 146 "SST=" 147 ARCL "SST" 148 PROMPT 149 "DFR=" 150 ARCL "DFR" 151 PROMPT 152 "DFE=" 153 ARCL "DFE" 154 PROMPT 155 "DFT=" 156 ARCL "DFT" 157 PROMPT 158 "MSR=" 159 ARCL "MSR" 160 PROMPT 161 "MSE=" 162 ARCL "MSE" 163 PROMPT 164 "FSTAT=" 165 ARCL "FSTAT" 166 PROMPT 167 BEEP 168 RTN 169 .END. Namir Shammas Program Version 1.0.0 Documentation Version 1.1.0 June 1, 2005 hp41cxcalc@aol.com