R version 2.8.1 (2008-12-22) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > y <- c(255 + ,280.2 + ,299.9 + ,339.2 + ,374.2 + ,393.5 + ,389.2 + ,381.7 + ,375.2 + ,369 + ,357.4 + ,352.1 + ,346.5 + ,342.9 + ,340.3 + ,328.3 + ,322.9 + ,314.3 + ,308.9 + ,294 + ,285.6 + ,281.2 + ,280.3 + ,278.8 + ,274.5 + ,270.4 + ,263.4 + ,259.9 + ,258 + ,262.7 + ,284.7 + ,311.3 + ,322.1 + ,327 + ,331.3 + ,333.3 + ,321.4 + ,327 + ,320 + ,314.7 + ,316.7 + ,314.4 + ,321.3 + ,318.2 + ,307.2 + ,301.3 + ,287.5 + ,277.7 + ,274.4 + ,258.8 + ,253.3 + ,251 + ,248.4 + ,249.5 + ,246.1 + ,244.5 + ,243.6 + ,244 + ,240.8 + ,249.8 + ,248 + ,259.4 + ,260.5 + ,260.8 + ,261.3 + ,259.5 + ,256.6 + ,257.9 + ,256.5 + ,254.2 + ,253.3 + ,253.8 + ,255.5 + ,257.1 + ,257.3 + ,253.2 + ,252.8 + ,252 + ,250.7 + ,252.2 + ,250 + ,251 + ,253.4 + ,251.2 + ,255.6 + ,261.1 + ,258.9 + ,259.9 + ,261.2 + ,264.7 + ,267.1 + ,266.4 + ,267.7 + ,268.6 + ,267.5 + ,268.5 + ,268.5 + ,270.5 + ,270.9 + ,270.1 + ,269.3 + ,269.8 + ,270.1 + ,264.9 + ,263.7 + ,264.8 + ,263.7 + ,255.9 + ,276.2 + ,360.1 + ,380.5 + ,373.7 + ,369.8 + ,366.6 + ,359.3 + ,345.8 + ,326.2 + ,324.5 + ,328.1 + ,327.5 + ,324.4 + ,316.5 + ,310.9 + ,301.5 + ,291.7 + ,290.4 + ,287.4 + ,277.7 + ,281.6 + ,288 + ,276 + ,272.9 + ,283 + ,283.3 + ,276.8 + ,284.5 + ,282.7 + ,281.2 + ,287.4 + ,283.1 + ,284 + ,285.5 + ,289.2 + ,292.5 + ,296.4 + ,305.2 + ,303.9 + ,311.5 + ,316.3 + ,316.7 + ,322.5 + ,317.1 + ,309.8 + ,303.8 + ,290.3 + ,293.7 + ,291.7 + ,296.5 + ,289.1 + ,288.5 + ,293.8 + ,297.7 + ,305.4 + ,302.7 + ,302.5 + ,303 + ,294.5 + ,294.1 + ,294.5 + ,297.1 + ,289.4 + ,292.4 + ,287.9 + ,286.6 + ,280.5 + ,272.4 + ,269.2 + ,270.6 + ,267.3 + ,262.5 + ,266.8 + ,268.8 + ,263.1 + ,261.2 + ,266 + ,262.5 + ,265.2 + ,261.3 + ,253.7 + ,249.2 + ,239.1 + ,236.4 + ,235.2 + ,245.2 + ,246.2 + ,247.7 + ,251.4 + ,253.3 + ,254.8 + ,250 + ,249.3 + ,241.5 + ,243.3 + ,248 + ,253 + ,252.9 + ,251.5 + ,251.6 + ,253.5 + ,259.8 + ,334.1 + ,448 + ,445.8 + ,445 + ,448.2 + ,438.2 + ,439.8 + ,423.4 + ,410.8 + ,408.4 + ,406.7 + ,405.9 + ,402.7 + ,405.1 + ,399.6 + ,386.5 + ,381.4 + ,375.2 + ,357.7 + ,359 + ,355 + ,352.7 + ,344.4 + ,343.8 + ,338 + ,339 + ,333.3 + ,334.4 + ,328.3 + ,330.7 + ,330 + ,331.6 + ,351.2 + ,389.4 + ,410.9 + ,442.8 + ,462.8 + ,466.9 + ,461.7 + ,439.2 + ,430.3 + ,416.1 + ,402.5 + ,397.3 + ,403.3 + ,395.9 + ,387.8 + ,378.6 + ,377.1 + ,370.4 + ,362 + ,350.3 + ,348.2 + ,344.6 + ,343.5 + ,342.8 + ,347.6 + ,346.6 + ,349.5 + ,342.1 + ,342 + ,342.8 + ,339.3 + ,348.2 + ,333.7 + ,334.7 + ,354 + ,367.7 + ,363.3 + ,358.4 + ,353.1 + ,343.1 + ,344.6 + ,344.4 + ,333.9 + ,331.7 + ,324.3 + ,321.2 + ,322.4 + ,321.7 + ,320.5 + ,312.8 + ,309.7 + ,315.6 + ,309.7 + ,304.6 + ,302.5 + ,301.5 + ,298.8 + ,291.3 + ,293.6 + ,294.6 + ,285.9 + ,297.6 + ,301.1 + ,293.8 + ,297.7 + ,292.9 + ,292.1 + ,287.2 + ,288.2 + ,283.8 + ,299.9 + ,292.4 + ,293.3 + ,300.8 + ,293.7 + ,293.1 + ,294.4 + ,292.1 + ,291.9 + ,282.5 + ,277.9 + ,287.5 + ,289.2 + ,285.6 + ,293.2 + ,290.8 + ,283.1 + ,275 + ,287.8 + ,287.8 + ,287.4 + ,284 + ,277.8 + ,277.6 + ,304.9 + ,294 + ,300.9 + ,324 + ,332.9 + ,341.6 + ,333.4 + ,348.2 + ,344.7 + ,344.7 + ,329.3 + ,323.5 + ,323.2 + ,317.4 + ,330.1 + ,329.2 + ,334.9 + ,315.8 + ,315.4 + ,319.6 + ,317.3 + ,313.8 + ,315.8 + ,311.3) > x <- c(87.28 + ,87.28 + ,87.09 + ,86.92 + ,87.59 + ,90.72 + ,90.69 + ,90.3 + ,89.55 + ,88.94 + ,88.41 + ,87.82 + ,87.07 + ,86.82 + ,86.4 + ,86.02 + ,85.66 + ,85.32 + ,85 + ,84.67 + ,83.94 + ,82.83 + ,81.95 + ,81.19 + ,80.48 + ,78.86 + ,69.47 + ,68.77 + ,70.06 + ,73.95 + ,75.8 + ,77.79 + ,81.57 + ,83.07 + ,84.34 + ,85.1 + ,85.25 + ,84.26 + ,83.63 + ,86.44 + ,85.3 + ,84.1 + ,83.36 + ,82.48 + ,81.58 + ,80.47 + ,79.34 + ,82.13 + ,81.69 + ,80.7 + ,79.88 + ,79.16 + ,78.38 + ,77.42 + ,76.47 + ,75.46 + ,74.48 + ,78.27 + ,80.7 + ,79.91 + ,78.75 + ,77.78 + ,81.14 + ,81.08 + ,80.03 + ,78.91 + ,78.01 + ,76.9 + ,75.97 + ,81.93 + ,80.27 + ,78.67 + ,77.42 + ,76.16 + ,74.7 + ,76.39 + ,76.04 + ,74.65 + ,73.29 + ,71.79 + ,74.39 + ,74.91 + ,74.54 + ,73.08 + ,72.75 + ,71.32 + ,70.38 + ,70.35 + ,70.01 + ,69.36 + ,67.77 + ,69.26 + ,69.8 + ,68.38 + ,67.62 + ,68.39 + ,66.95 + ,65.21 + ,66.64 + ,63.45 + ,60.66 + ,62.34 + ,60.32 + ,58.64 + ,60.46 + ,58.59 + ,61.87 + ,61.85 + ,67.44 + ,77.06 + ,91.74 + ,93.15 + ,94.15 + ,93.11 + ,91.51 + ,89.96 + ,88.16 + ,86.98 + ,88.03 + ,86.24 + ,84.65 + ,83.23 + ,81.7 + ,80.25 + ,78.8 + ,77.51 + ,76.2 + ,75.04 + ,74 + ,75.49 + ,77.14 + ,76.15 + ,76.27 + ,78.19 + ,76.49 + ,77.31 + ,76.65 + ,74.99 + ,73.51 + ,72.07 + ,70.59 + ,71.96 + ,76.29 + ,74.86 + ,74.93 + ,71.9 + ,71.01 + ,77.47 + ,75.78 + ,76.6 + ,76.07 + ,74.57 + ,73.02 + ,72.65 + ,73.16 + ,71.53 + ,69.78 + ,67.98 + ,69.96 + ,72.16 + ,70.47 + ,68.86 + ,67.37 + ,65.87 + ,72.16 + ,71.34 + ,69.93 + ,68.44 + ,67.16 + ,66.01 + ,67.25 + ,70.91 + ,69.75 + ,68.59 + ,67.48 + ,66.31 + ,64.81 + ,66.58 + ,65.97 + ,64.7 + ,64.7 + ,60.94 + ,59.08 + ,58.42 + ,57.77 + ,57.11 + ,53.31 + ,49.96 + ,49.4 + ,48.84 + ,48.3 + ,47.74 + ,47.24 + ,46.76 + ,46.29 + ,48.9 + ,49.23 + ,48.53 + ,48.03 + ,54.34 + ,53.79 + ,53.24 + ,52.96 + ,52.17 + ,51.7 + ,58.55 + ,78.2 + ,77.03 + ,76.19 + ,77.15 + ,75.87 + ,95.47 + ,109.67 + ,112.28 + ,112.01 + ,107.93 + ,105.96 + ,105.06 + ,102.98 + ,102.2 + ,105.23 + ,101.85 + ,99.89 + ,96.23 + ,94.76 + ,91.51 + ,91.63 + ,91.54 + ,85.23 + ,87.83 + ,87.38 + ,84.44 + ,85.19 + ,84.03 + ,86.73 + ,102.52 + ,104.45 + ,106.98 + ,107.02 + ,99.26 + ,94.45 + ,113.44 + ,157.33 + ,147.38 + ,171.89 + ,171.95 + ,132.71 + ,126.02 + ,121.18 + ,115.45 + ,110.48 + ,117.85 + ,117.63 + ,124.65 + ,109.59 + ,111.27 + ,99.78 + ,98.21 + ,99.2 + ,97.97 + ,89.55 + ,87.91 + ,93.34 + ,94.42 + ,93.2 + ,90.29 + ,91.46 + ,89.98 + ,88.35 + ,88.41 + ,82.44 + ,79.89 + ,75.69 + ,75.66 + ,84.5 + ,96.73 + ,87.48 + ,82.39 + ,83.48 + ,79.31 + ,78.16 + ,72.77 + ,72.45 + ,68.46 + ,67.62 + ,68.76 + ,70.07 + ,68.55 + ,65.3 + ,58.96 + ,59.17 + ,62.37 + ,66.28 + ,55.62 + ,55.23 + ,55.85 + ,56.75 + ,50.89 + ,53.88 + ,52.95 + ,55.08 + ,53.61 + ,58.78 + ,61.85 + ,55.91 + ,53.32 + ,46.41 + ,44.57 + ,50 + ,50 + ,53.36 + ,46.23 + ,50.45 + ,49.07 + ,45.85 + ,48.45 + ,49.96 + ,46.53 + ,50.51 + ,47.58 + ,48.05 + ,46.84 + ,47.67 + ,49.16 + ,55.54 + ,55.82 + ,58.22 + ,56.19 + ,57.77 + ,63.19 + ,54.76 + ,55.74 + ,62.54 + ,61.39 + ,69.6 + ,79.23 + ,80 + ,93.68 + ,107.63 + ,100.18 + ,97.3 + ,90.45 + ,80.64 + ,80.58 + ,75.82 + ,85.59 + ,89.35 + ,89.42 + ,104.73 + ,95.32 + ,89.27 + ,90.44 + ,86.97 + ,79.98 + ,81.22 + ,87.35 + ,83.64 + ,82.22 + ,94.4 + ,102.18) > par2 = '36' > par1 = '0' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: > par1 <- as.numeric(par1) > par2 <- as.numeric(par2) > x <- as.ts(x) > y <- as.ts(y) > mylm <- lm(y~x) > cbind(mylm$resid) [,1] 1 -7.129220e+01 2 -4.609220e+01 3 -2.604228e+01 4 1.357080e+01 5 4.733687e+01 6 6.087241e+01 7 5.662766e+01 8 4.984591e+01 9 4.472718e+01 10 3.965060e+01 11 2.902669e+01 12 2.481329e+01 13 2.059455e+01 14 1.745497e+01 15 1.562848e+01 16 4.328317e+00 17 -4.086770e-01 18 -8.382505e+00 19 -1.319317e+01 20 -2.748541e+01 21 -3.454098e+01 22 -3.689671e+01 23 -3.617603e+01 24 -3.627635e+01 25 -3.926876e+01 26 -4.038523e+01 27 -3.009183e+01 28 -3.230265e+01 29 -3.657842e+01 30 -3.904257e+01 31 -2.044968e+01 32 2.485369e+00 33 6.323807e+00 34 8.461282e+00 35 1.042234e+01 36 1.102267e+01 37 -1.153587e+00 38 6.269679e+00 39 4.299398e-01 40 -1.004519e+01 41 -5.945671e+00 42 -6.035651e+00 43 2.227194e+00 44 7.478754e-01 45 -8.594610e+00 46 -1.245034e+01 47 -2.416924e+01 48 -3.910754e+01 49 -4.159719e+01 50 -5.537393e+01 51 -5.936375e+01 52 -6.033774e+01 53 -6.150122e+01 54 -5.863321e+01 55 -6.028361e+01 56 -6.002351e+01 57 -5.911866e+01 58 -6.569864e+01 59 -7.337393e+01 60 -6.291900e+01 61 -6.258265e+01 62 -4.939621e+01 63 -5.448427e+01 64 -5.407377e+01 65 -5.164000e+01 66 -5.137732e+01 67 -5.261980e+01 68 -4.927553e+01 69 -4.896277e+01 70 -6.223920e+01 71 -6.008200e+01 72 -5.663531e+01 73 -5.263321e+01 74 -4.871269e+01 75 -4.582383e+01 76 -5.303627e+01 77 -5.279169e+01 78 -5.103175e+01 79 -4.982706e+01 80 -4.556453e+01 81 -5.255291e+01 82 -5.251058e+01 83 -4.942916e+01 84 -4.894030e+01 85 -4.393255e+01 86 -3.579894e+01 87 -3.626776e+01 88 -3.521251e+01 89 -3.328634e+01 90 -2.858924e+01 91 -2.326097e+01 92 -2.670507e+01 93 -2.639958e+01 94 -2.288439e+01 95 -2.258471e+01 96 -2.300281e+01 97 -2.035079e+01 98 -1.514626e+01 99 -1.737986e+01 100 -1.230489e+01 101 -7.966598e+00 102 -1.056063e+01 103 -6.540426e+00 104 -8.646398e+00 105 -1.319826e+01 106 -8.654314e+00 107 -1.579503e+01 108 -2.355820e+01 109 -1.355321e+01 110 5.262980e+01 111 4.599389e+01 112 3.659712e+01 113 3.085543e+01 114 2.957078e+01 115 2.521748e+01 116 1.457209e+01 117 -1.712885e+00 118 -1.239699e+00 119 4.265339e-01 120 3.123147e+00 121 2.951423e+00 122 -2.333387e+00 123 -5.115612e+00 124 -1.184517e+01 125 -1.897473e+01 126 -1.789896e+01 127 -1.848635e+01 128 -2.605000e+01 129 -2.023465e+01 130 -1.657876e+01 131 -3.161754e+01 132 -3.289427e+01 133 -2.301527e+01 134 -2.625130e+01 135 -2.962044e+01 136 -2.343062e+01 137 -2.401511e+01 138 -2.245792e+01 139 -1.353223e+01 140 -1.518020e+01 141 -1.155451e+01 142 -1.257762e+01 143 -1.685211e+01 144 -1.091850e+01 145 -7.147417e+00 146 7.232883e+00 147 7.571981e+00 148 3.274708e+00 149 1.118715e+01 150 1.007697e+01 151 1.685306e+01 152 1.421559e+01 153 9.770198e+00 154 4.451621e+00 155 -9.987638e+00 156 -3.585694e+00 157 -2.362748e+00 158 5.752281e+00 159 -5.294251e+00 160 -9.945954e+00 161 -1.533510e+00 162 5.331600e+00 163 1.577571e+01 164 1.583823e+01 165 4.054046e+00 166 6.064226e+00 167 1.609991e-01 168 2.505107e+00 169 5.262461e+00 170 9.980397e+00 171 -3.290047e-03 172 -3.743850e+00 173 -6.107498e+00 174 -5.271145e+00 175 -9.326877e+00 176 -1.527211e+01 177 -1.570958e+01 178 -1.756936e+01 179 -1.974594e+01 180 -2.220700e+01 181 -1.790700e+01 182 -8.982269e+00 183 -1.125674e+01 184 -1.194123e+01 185 -5.944134e+00 186 -8.228623e+00 187 1.469773e+00 188 3.739412e+00 189 -2.829246e+00 190 -6.297903e+00 191 -1.540339e+01 192 -1.707205e+01 193 -1.735121e+01 194 -6.467202e+00 195 -4.601611e+00 196 -7.908404e+00 197 -4.816160e+00 198 -1.626981e+00 199 7.938603e-01 200 -1.562716e+01 201 -1.531423e+01 202 -2.210131e+01 203 -1.978564e+01 204 -1.363071e+01 205 -7.765117e+00 206 -2.048065e+01 207 -5.806972e+01 208 -5.581495e+01 209 -5.236794e+01 210 -4.783595e+01 211 2.882140e+01 212 1.066244e+02 213 7.827251e+01 214 7.266572e+01 215 7.636297e+01 216 7.387704e+01 217 7.910515e+01 218 6.436267e+01 219 5.559337e+01 220 5.462988e+01 221 4.734958e+01 222 5.277447e+01 223 5.318417e+01 224 6.232473e+01 225 5.953201e+01 226 5.241748e+01 227 4.709647e+01 228 4.106223e+01 229 3.518325e+01 230 3.169487e+01 231 2.852363e+01 232 3.163818e+01 233 2.195691e+01 234 2.349327e+01 235 1.272072e+01 236 -1.535946e+01 237 -2.461390e+01 238 -2.817336e+01 239 -3.434703e+01 240 -1.765557e+01 241 -9.497072e+00 242 -4.287064e+01 243 -1.041021e+02 244 -4.757736e+01 245 -7.121702e+01 246 -3.942752e+01 247 5.284013e+01 248 6.926099e+01 249 7.297474e+01 250 6.102758e+01 251 6.128075e+01 252 3.350754e+01 253 2.031271e+01 254 2.184097e+00 255 3.591985e+01 256 2.542582e+01 257 3.848676e+01 258 3.217820e+01 259 2.885493e+01 260 2.442020e+01 261 3.152718e+01 262 2.284754e+01 263 1.074720e+01 264 5.158179e+00 265 6.305032e+00 266 1.096433e+01 267 1.360956e+01 268 1.533525e+01 269 2.123720e+01 270 1.372669e+01 271 2.462154e+01 272 3.011783e+01 273 3.435290e+01 274 4.330815e+01 275 1.252768e+01 276 -8.996110e+00 277 2.733946e+01 278 5.041363e+01 279 4.400619e+01 280 4.678601e+01 281 4.360395e+01 282 4.353062e+01 283 4.561996e+01 284 5.276827e+01 285 4.381529e+01 286 3.951577e+01 287 2.970316e+01 288 2.940252e+01 289 3.658799e+01 290 4.756426e+01 291 4.597751e+01 292 3.238412e+01 293 2.208314e+01 294 4.761549e+01 295 4.243374e+01 296 3.619190e+01 297 3.243438e+01 298 4.222665e+01 299 3.402001e+01 300 2.823278e+01 301 2.660999e+01 302 3.031727e+01 303 1.209577e+01 304 1.814180e+01 305 3.258140e+01 306 3.005136e+01 307 4.667739e+01 308 4.526608e+01 309 3.446574e+01 310 2.956574e+01 311 2.437769e+01 312 3.310889e+01 313 4.143699e+01 314 3.647851e+01 315 4.330873e+01 316 4.602035e+01 317 3.613941e+01 318 4.185639e+01 319 3.582649e+01 320 3.892262e+01 321 3.785703e+01 322 3.068546e+01 323 2.455687e+01 324 3.141276e+01 325 2.136282e+01 326 1.724715e+01 327 2.042711e+01 328 2.176573e+01 329 1.115587e+01 330 -6.926056e+00 331 2.139933e+01 332 1.959448e+01 333 6.671038e+00 334 5.388973e+00 335 -1.593125e+01 336 -3.386665e+01 337 -7.984750e+00 338 -4.407898e+01 339 -6.287046e+01 340 -2.604992e+01 341 -1.184587e+01 342 9.469661e+00 343 1.933657e+01 344 3.424707e+01 345 3.951349e+01 346 2.152024e+01 347 -8.044879e-01 348 -6.733406e+00 349 -3.522957e+01 350 -2.369934e+01 351 1.428468e-01 352 -2.911923e+00 353 9.178718e+00 354 2.952083e+00 355 2.683962e-01 356 -6.821122e+00 357 -2.288477e+00 358 -3.173287e+00 359 -2.360499e+01 360 -4.243328e+01 > library(lattice) > postscript(file="/var/www/rcomp/tmp/1bvm11257515031.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(y,type='l',main='Run Sequence Plot of Y[t]',xlab='time or index',ylab='value') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/2st1r1257515031.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x,type='l',main='Run Sequence Plot of X[t]',xlab='time or index',ylab='value') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/3alzd1257515031.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x,y,main='Scatter Plot',xlab='X[t]',ylab='Y[t]') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/4huk01257515031.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mylm$resid,type='l',main='Run Sequence Plot of e[t]',xlab='time or index',ylab='value') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5q1601257515031.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mylm$resid,main='Histogram of e[t]') > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6n39l1257515031.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > if (par1 > 0) + { + densityplot(~mylm$resid,col='black',main=paste('Density Plot of e[t] bw = ',par1),bw=par1) + } else { + densityplot(~mylm$resid,col='black',main='Density Plot of e[t]') + } > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7ljvz1257515031.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(mylm$resid,main='QQ plot of e[t]') > qqline(mylm$resid) > grid() > dev.off() null device 1 > if (par2 > 0) + { + postscript(file="/var/www/rcomp/tmp/85oj11257515031.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + acf(mylm$resid,lag.max=par2,main='Residual Autocorrelation Function') + grid() + dev.off() + } null device 1 > summary(x) Min. 1st Qu. Median Mean 3rd Qu. Max. 44.57 66.87 76.78 77.54 87.00 172.00 > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Model: Y[t] = c + b X[t] + e[t]',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'c',1,TRUE) > a<-table.element(a,mylm$coeff[[1]]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'b',1,TRUE) > a<-table.element(a,mylm$coeff[[2]]) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/9va291257515032.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Descriptive Statistics about e[t]',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'# observations',header=TRUE) > a<-table.element(a,length(mylm$resid)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'minimum',header=TRUE) > a<-table.element(a,min(mylm$resid)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Q1',header=TRUE) > a<-table.element(a,quantile(mylm$resid,0.25)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'median',header=TRUE) > a<-table.element(a,median(mylm$resid)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'mean',header=TRUE) > a<-table.element(a,mean(mylm$resid)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Q3',header=TRUE) > a<-table.element(a,quantile(mylm$resid,0.75)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'maximum',header=TRUE) > a<-table.element(a,max(mylm$resid)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/108rjk1257515032.tab") > > system("convert tmp/1bvm11257515031.ps tmp/1bvm11257515031.png") > system("convert tmp/2st1r1257515031.ps tmp/2st1r1257515031.png") > system("convert tmp/3alzd1257515031.ps tmp/3alzd1257515031.png") > system("convert tmp/4huk01257515031.ps tmp/4huk01257515031.png") > system("convert tmp/5q1601257515031.ps tmp/5q1601257515031.png") > system("convert tmp/6n39l1257515031.ps tmp/6n39l1257515031.png") > system("convert tmp/7ljvz1257515031.ps tmp/7ljvz1257515031.png") > system("convert tmp/85oj11257515031.ps tmp/85oj11257515031.png") > > > proc.time() user system elapsed 2.27 1.85 3.11