R version 2.9.0 (2009-04-17) Copyright (C) 2009 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. > x <- array(list(1.3866 + ,126.64 + ,40.7819 + ,1.3582 + ,126.81 + ,39.5915 + ,1.3332 + ,125.84 + ,38.8859 + ,1.3595 + ,126.77 + ,39.9068 + ,1.3617 + ,124.34 + ,41.47 + ,1.3684 + ,124.4 + ,41.5613 + ,1.3394 + ,120.48 + ,41.6005 + ,1.3262 + ,118.54 + ,41.4113 + ,1.3173 + ,117.66 + ,41.84 + ,1.3085 + ,116.97 + ,42.2892 + ,1.327 + ,120.11 + ,43.1521 + ,1.3182 + ,119.16 + ,43.5998 + ,1.293 + ,116.9 + ,43.116 + ,1.291 + ,116.11 + ,42.4185 + ,1.2984 + ,114.98 + ,42.3687 + ,1.2795 + ,113.65 + ,42.2975 + ,1.299 + ,115.82 + ,42.8528 + ,1.3174 + ,117.59 + ,43.535 + ,1.326 + ,118.57 + ,44.7265 + ,1.3111 + ,118.07 + ,45.7293 + ,1.2816 + ,114.98 + ,45.7585 + ,1.276 + ,114.04 + ,46.1685 + ,1.2849 + ,115.02 + ,46.5075 + ,1.2818 + ,114.28 + ,46.527 + ,1.2829 + ,115.04 + ,46.601 + ,1.2796 + ,116.7 + ,46.4607 + ,1.3008 + ,119.21 + ,46.7135 + ,1.2967 + ,118.39 + ,46.4113 + ,1.2938 + ,116.5 + ,45.55 + ,1.2833 + ,115.46 + ,44.6081 + ,1.2823 + ,117.59 + ,44.4395 + ,1.2765 + ,117.33 + ,44.9847 + ,1.2634 + ,116.2 + ,45.7558 + ,1.2596 + ,116.83 + ,45.3942 + ,1.2705 + ,118.99 + ,45.697 + ,1.2591 + ,118.62 + ,45.5664 + ,1.2798 + ,121.09 + ,46.0205 + ,1.2763 + ,122.4 + ,45.9195 + ,1.2795 + ,123.76 + ,45.8005 + ,1.2782 + ,125.33 + ,45.535 + ,1.2644 + ,123.23 + ,45.4977 + ,1.2596 + ,122.52 + ,45.5782 + ,1.2615 + ,123.64 + ,45.7697 + ,1.2555 + ,124.67 + ,45.2445 + ,1.2555 + ,124.71 + ,45.0615 + ,1.2658 + ,122.53 + ,45.2865 + ,1.2565 + ,124.4 + ,44.791 + ,1.2783 + ,125.45 + ,44.7625 + ,1.2786 + ,125.35 + ,44.7644 + ,1.2782 + ,124.3 + ,44.9973 + ,1.2905 + ,127.03 + ,44.7265 + ,1.3042 + ,128.51 + ,45.1465 + ,1.2942 + ,128.1 + ,44.7465 + ,1.313 + ,128.94 + ,45.1795 + ,1.3671 + ,129.67 + ,45.6515 + ,1.3549 + ,129.87 + ,45.492 + ,1.3558 + ,131.12 + ,45.2775 + ,1.3507 + ,132.68 + ,45.2115 + ,1.3494 + ,132.24 + ,45.411 + ,1.3607 + ,133.63 + ,45.4005 + ,1.3295 + ,129.91 + ,44.7692 + ,1.3193 + ,127.93 + ,44.8913 + ,1.3308 + ,131.17 + ,45.032 + ,1.3246 + ,130.86 + ,44.879 + ,1.3392 + ,133.48 + ,44.833 + ,1.3425 + ,134.08 + ,44.8257 + ,1.3496 + ,136.02 + ,44.7815 + ,1.3255 + ,132.8 + ,44.479 + ,1.3231 + ,132.37 + ,44.6317 + ,1.3273 + ,133.05 + ,44.5043 + ,1.3276 + ,132.57 + ,44.3217 + ,1.3173 + ,130.7 + ,44.1005 + ,1.3196 + ,130.5 + ,44.047 + ,1.3058 + ,129.67 + ,43.6835 + ,1.2966 + ,127.8 + ,43.7864 + ,1.2932 + ,126.82 + ,44.1807 + ,1.2947 + ,126.85 + ,43.9595 + ,1.305 + ,128.28 + ,43.937 + ,1.3232 + ,128.3 + ,43.991 + ,1.3125 + ,126.82 + ,43.865 + ,1.2992 + ,125.08 + ,43.671 + ,1.3266 + ,128.53 + ,43.93 + ,1.3275 + ,130.34 + ,43.863 + ,1.3223 + ,131.52 + ,43.7095 + ,1.3403 + ,132.59 + ,43.9435 + ,1.3322 + ,131.17 + ,43.736 + ,1.3363 + ,132.72 + ,43.6295 + ,1.3425 + ,133.36 + ,43.598 + ,1.3574 + ,132.82 + ,43.8726 + ,1.3683 + ,132.9 + ,43.8935 + ,1.3623 + ,130.9 + ,43.5957 + ,1.3563 + ,129.41 + ,43.7155 + ,1.3518 + ,128.67 + ,43.528 + ,1.3494 + ,129.28 + ,43.3415 + ,1.3612 + ,130.91 + ,43.3374 + ,1.369 + ,131.06 + ,43.332 + ,1.3771 + ,130.84 + ,43.3869 + ,1.3972 + ,131.41 + ,43.5016 + ,1.401 + ,133.22 + ,43.4875 + ,1.3908 + ,132.06 + ,43.6023 + ,1.3901 + ,132.48 + ,43.3886 + ,1.3856 + ,134.38 + ,43.3105 + ,1.4098 + ,135.22 + ,43.4455 + ,1.422 + ,134.89 + ,43.5185 + ,1.4238 + ,136.09 + ,43.5755 + ,1.4207 + ,136.33 + ,43.6217 + ,1.4095 + ,136.32 + ,43.644 + ,1.4177 + ,137.48 + ,43.5789 + ,1.3866 + ,136.53 + ,43.5215 + ,1.3959 + ,136.8 + ,43.5033 + ,1.4102 + ,138.03 + ,43.632 + ,1.3969 + ,137.39 + ,43.263 + ,1.4004 + ,137.55 + ,43.3717 + ,1.385 + ,136.08 + ,43.2745 + ,1.389 + ,134.78 + ,43.2647 + ,1.384 + ,133.28 + ,43.324 + ,1.392 + ,133.57 + ,43.4455 + ,1.3932 + ,134.84 + ,43.4098 + ,1.3858 + ,133.02 + ,43.41 + ,1.3978 + ,133.49 + ,43.93 + ,1.4029 + ,133.77 + ,43.8104 + ,1.394 + ,134.34 + ,43.54 + ,1.4096 + ,134.5 + ,43.858 + ,1.4058 + ,134.03 + ,43.8375 + ,1.4134 + ,135.51 + ,43.881 + ,1.4096 + ,136.53 + ,43.887 + ,1.4049 + ,135.95 + ,43.8009 + ,1.4009 + ,134.32 + ,43.7877 + ,1.3897 + ,132.44 + ,43.811 + ,1.4019 + ,133.61 + ,44.0625 + ,1.3901 + ,131.02 + ,44.125 + ,1.399 + ,130.05 + ,44.52 + ,1.3901 + ,128.21 + ,45.4005 + ,1.3975 + ,129.03 + ,45.89 + ,1.3991 + ,130.34 + ,45.189 + ,1.4089 + ,131.57 + ,44.9035 + ,1.413 + ,132.63 + ,44.9351 + ,1.409 + ,132.06 + ,44.801 + ,1.4217 + ,134.44 + ,43.98 + ,1.4223 + ,134.1 + ,44.11 + ,1.4191 + ,132.49 + ,44.2661 + ,1.4229 + ,134.23 + ,44.361 + ,1.4227 + ,134.92 + ,44.099 + ,1.4269 + ,135.61 + ,43.8435 + ,1.4229 + ,134.53 + ,43.8914 + ,1.4104 + ,133.86 + ,44.217 + ,1.4053 + ,133.89 + ,44.506 + ,1.4138 + ,135.33 + ,44.54 + ,1.4303 + ,135.86 + ,44.4465 + ,1.4384 + ,136.22 + ,44.842 + ,1.441 + ,137.38 + ,44.8946 + ,1.437 + ,137.31 + ,44.951 + ,1.4357 + ,136.89 + ,45.445 + ,1.4202 + ,138.01 + ,45.0035 + ,1.4166 + ,136.72 + ,45.769 + ,1.417 + ,135.77 + ,46.09 + ,1.4293 + ,137.52 + ,45.412 + ,1.4294 + ,135.61 + ,45.12 + ,1.4072 + ,132.94 + ,45.48 + ,1.4101 + ,134.12 + ,45.105 + ,1.4112 + ,132.55 + ,45.056 + ,1.4243 + ,134.11 + ,45.22 + ,1.433 + ,134.19 + ,45.39 + ,1.4323 + ,135.57 + ,45.041 + ,1.4324 + ,135.05 + ,44.9399 + ,1.427 + ,134.32 + ,44.9315 + ,1.4268 + ,133.61 + ,45.1935 + ,1.4364 + ,134.75 + ,45.3466 + ,1.4272 + ,133.1 + ,45.4645 + ,1.4314 + ,133.26 + ,45.5685 + ,1.422 + ,131.63 + ,45.3921 + ,1.4335 + ,132.47 + ,45.34 + ,1.4262 + ,132.45 + ,45.1308 + ,1.433 + ,133.33 + ,45.1005 + ,1.4473 + ,133.57 + ,45.37 + ,1.4522 + ,134.13 + ,45.2 + ,1.4545 + ,133.92 + ,44.9614 + ,1.4594 + ,132.62 + ,44.8015 + ,1.4561 + ,132.3 + ,44.9152 + ,1.4611 + ,133.26 + ,45.095 + ,1.4671 + ,132.6 + ,44.9271 + ,1.4712 + ,134.38 + ,44.6026 + ,1.4705 + ,134.17 + ,44.5 + ,1.4658 + ,135.46 + ,44.54 + ,1.478 + ,135.09 + ,44.5532 + ,1.4783 + ,134.96 + ,44.407 + ,1.4768 + ,133.85 + ,44.259 + ,1.467 + ,132.59 + ,44.1365 + ,1.465 + ,131.15 + ,44.112 + ,1.4549 + ,130.91 + ,43.8814 + ,1.4643 + ,131.07 + ,43.98 + ,1.4539 + ,130.78 + ,43.7294 + ,1.4537 + ,129.95 + ,43.9119 + ,1.4616 + ,131.41 + ,43.955 + ,1.4722 + ,131.21 + ,43.9 + ,1.4694 + ,130.68 + ,43.7065 + ,1.4763 + ,130.46 + ,43.6939 + ,1.475 + ,131.12 + ,43.6587 + ,1.4765 + ,132.99 + ,43.5885 + ,1.4864 + ,133.02 + ,43.8885 + ,1.4881 + ,133.39 + ,43.8216 + ,1.4864 + ,134.07 + ,43.751 + ,1.4869 + ,135.6 + ,43.699 + ,1.4918 + ,135.66 + ,43.7425 + ,1.4971 + ,135.53 + ,43.639 + ,1.4921 + ,135.82 + ,43.589 + ,1.5 + ,136.9 + ,43.606 + ,1.502 + ,137.97 + ,43.5325 + ,1.5019 + ,138.09 + ,43.385 + ,1.4874 + ,136.91 + ,43.3745 + ,1.4785 + ,134.76 + ,43.236 + ,1.4788 + ,135.13 + ,43.1957 + ,1.48 + ,134.66 + ,43.01 + ,1.4772 + ,132.95 + ,43.1401 + ,1.4658 + ,132.25 + ,43.0487 + ,1.4761 + ,134.3 + ,43.1972 + ,1.4867 + ,134.3 + ,43.2461 + ,1.4862 + ,134.76 + ,43.0866 + ,1.4984 + ,134.81 + ,43.0865 + ,1.4966 + ,134.51 + ,43.0194 + ,1.5037 + ,135.11 + ,43.08 + ,1.4922 + ,134.32 + ,43.007 + ,1.4868 + ,133.51 + ,42.9278 + ,1.4965 + ,134.02 + ,42.9545 + ,1.4875 + ,132.76 + ,42.7995 + ,1.4957 + ,133.39 + ,42.9048 + ,1.4863 + ,132.05 + ,42.9468 + ,1.4815 + ,131.87 + ,43.08 + ,1.4968 + ,133.03 + ,43.1274 + ,1.4969 + ,132.57 + ,43.1625 + ,1.5083 + ,132.1 + ,43.45 + ,1.5071 + ,130.7 + ,43.831 + ,1.4918 + ,129.2 + ,43.7769 + ,1.5023 + ,129.77 + ,43.98 + ,1.5074 + ,131.02 + ,43.92 + ,1.509 + ,131.55 + ,44.11 + ,1.512 + ,133.17 + ,44.03 + ,1.5068 + ,133.08 + ,44.1582 + ,1.4787 + ,133.24 + ,44.14 + ,1.4774 + ,130.74 + ,45.07 + ,1.4768 + ,129.91 + ,44.8737 + ,1.473 + ,130.03 + ,44.8505 + ,1.4757 + ,131.13 + ,44.373 + ,1.4647 + ,129.55 + ,44.075 + ,1.4541 + ,130.22 + ,43.9725 + ,1.456 + ,130.61 + ,44.094 + ,1.4343 + ,129.27 + ,44.191 + ,1.4337 + ,129.68 + ,43.9685 + ,1.4368 + ,130.1 + ,43.79 + ,1.4279 + ,130.83 + ,43.6041 + ,1.4276 + ,130.95 + ,43.1707 + ,1.4398 + ,131.73 + ,42.71 + ,1.4405 + ,131.86 + ,42.755 + ,1.4433 + ,132.44 + ,43.3316 + ,1.4338 + ,132.35 + ,43.5 + ,1.4406 + ,133.16 + ,43.154 + ,1.4389 + ,133.62 + ,43.16 + ,1.4442 + ,132.54 + ,43.1 + ,1.435 + ,132.69 + ,42.85 + ,1.4304 + ,133.5 + ,42.6175 + ,1.4273 + ,133.36 + ,42.5 + ,1.4528 + ,134.23 + ,42.6285 + ,1.4481 + ,132.41 + ,42.6974 + ,1.4563 + ,133.02 + ,43.04 + ,1.4486 + ,132.88 + ,42.673 + ,1.4374 + ,130.76 + ,42.5015 + ,1.4369 + ,130.33 + ,42.538 + ,1.4279 + ,129.79 + ,42.3735 + ,1.4132 + ,128.65 + ,42.014 + ,1.4064 + ,129.14 + ,41.8618 + ,1.4135 + ,127.35 + ,42.1824 + ,1.4151 + ,127.74 + ,42.605 + ,1.4085 + ,126.31 + ,42.7345 + ,1.4072 + ,125.95 + ,42.615 + ,1.3999 + ,126.36 + ,42.465 + ,1.3966 + ,126.15 + ,42.34 + ,1.3913 + ,125.6 + ,42.251 + ,1.3937 + ,126.2 + ,42.0475 + ,1.3984 + ,126.73 + ,41.86 + ,1.3847 + ,125.68 + ,41.685 + ,1.3691 + ,122.49 + ,41.735 + ,1.3675 + ,122.07 + ,41.706 + ,1.376 + ,123.4 + ,41.764 + ,1.374 + ,123.01 + ,41.58 + ,1.3718 + ,123.03 + ,41.373 + ,1.3572 + ,122.33 + ,41.088 + ,1.3607 + ,122.42 + ,41.137 + ,1.3649 + ,122.68 + ,41.1587 + ,1.3726 + ,124.69 + ,41.185 + ,1.3567 + ,123.3 + ,40.819 + ,1.3519 + ,124.17 + ,40.633 + ,1.3626 + ,124.38 + ,40.858 + ,1.3577 + ,123.19 + ,40.794 + ,1.3547 + ,122.16 + ,40.69 + ,1.3489 + ,120.66 + ,40.595 + ,1.357 + ,120.92 + ,40.7305 + ,1.3525 + ,120.67 + ,40.5471 + ,1.3548 + ,120.68 + ,40.5145 + ,1.3641 + ,121.1 + ,40.7 + ,1.3668 + ,120.86 + ,40.7 + ,1.3582 + ,121.48 + ,40.522 + ,1.3662 + ,123.48 + ,40.6165 + ,1.3557 + ,121.72 + ,40.3985 + ,1.361 + ,123.16 + ,40.2815 + ,1.3657 + ,123.84 + ,40.245 + ,1.3765 + ,124.57 + ,40.3055 + ,1.3705 + ,124.3 + ,40.2696 + ,1.3723 + ,124.22 + ,40.251 + ,1.3756 + ,124.43 + ,40.127 + ,1.366 + ,123.33 + ,39.95 + ,1.3548 + ,122.86 + ,39.675 + ,1.3471 + ,121.25 + ,39.954 + ,1.3519 + ,122.16 + ,39.8828 + ,1.3338 + ,122.62 + ,39.62 + ,1.3356 + ,123.44 + ,39.5415 + ,1.3353 + ,124 + ,39.525 + ,1.3471 + ,124.75 + ,39.8145 + ,1.3482 + ,124.8 + ,39.6675 + ,1.3479 + ,125.93 + ,39.695 + ,1.3468 + ,126.28 + ,39.5985 + ,1.3396 + ,126.04 + ,39.2735 + ,1.334 + ,125.04 + ,39.1435 + ,1.3296 + ,123.76 + ,39.1742 + ,1.3384 + ,125.34 + ,39.2025 + ,1.3585 + ,126.99 + ,39.3946 + ,1.3583 + ,126.34 + ,39.5025 + ,1.3615 + ,127.42 + ,39.4845 + ,1.3544 + ,126.18 + ,39.33 + ,1.3535 + ,125.3 + ,39.295 + ,1.3432 + ,123.5 + ,39.2675 + ,1.3486 + ,125.32 + ,39.2535 + ,1.3373 + ,124.65 + ,38.9845 + ,1.3339 + ,124.03 + ,38.9285 + ,1.3311 + ,125.11 + ,38.8592 + ,1.3321 + ,125.46 + ,38.77 + ,1.329 + ,124.7 + ,38.79 + ,1.3245 + ,124.48 + ,38.8205 + ,1.3256 + ,124.76 + ,38.7577 + ,1.3315 + ,125.81 + ,38.839 + ,1.3238 + ,124.95 + ,38.78 + ,1.3089 + ,123.66 + ,38.54 + ,1.2924 + ,122.66 + ,38.511 + ,1.2727 + ,119.34 + ,38.615 + ,1.2746 + ,117.84 + ,38.898 + ,1.2969 + ,120.97 + ,38.8691 + ,1.2698 + ,117.38 + ,38.384 + ,1.2686 + ,118.06 + ,38.0277 + ,1.2587 + ,116.99 + ,37.72 + ,1.2492 + ,115.55 + ,37.7325 + ,1.2349 + ,114.17 + ,37.626 + ,1.2428 + ,115.32 + ,37.603 + ,1.227 + ,112.49 + ,37.78 + ,1.2334 + ,111.93 + ,38.559 + ,1.2497 + ,112.08 + ,39.0459 + ,1.236 + ,111.63 + ,38.45 + ,1.2223 + ,109.53 + ,38.505 + ,1.2309 + ,111.35 + ,38.2885 + ,1.2255 + ,110.79 + ,37.795 + ,1.2384 + ,113.06 + ,37.92 + ,1.2307 + ,112.62 + ,38.034 + ,1.2155 + ,110.65 + ,38.029 + ,1.2218 + ,112.36 + ,38.063 + ,1.2268 + ,113.74 + ,37.9828 + ,1.206 + ,111.73 + ,37.745 + ,1.1959 + ,109.86 + ,37.969 + ,1.1942 + ,109.32 + ,38.007 + ,1.201 + ,109.99 + ,38.0615 + ,1.2045 + ,109.84 + ,38.0912 + ,1.2127 + ,111.13 + ,38.091 + ,1.2249 + ,112.43 + ,38.431 + ,1.2258 + ,111.77 + ,38.48 + ,1.2277 + ,112.15 + ,38.35 + ,1.2363 + ,112.89 + ,38.214 + ,1.2372 + ,112.12 + ,38.384 + ,1.2391 + ,113.1 + ,38.1375 + ,1.2258 + ,111.09 + ,38.0075 + ,1.2271 + ,110.76 + ,38.0524 + ,1.2262 + ,109.59 + ,38.235 + ,1.2294 + ,109.99 + ,38.31 + ,1.2339 + ,110.25 + ,38.2615 + ,1.2198 + ,108.31 + ,38.13 + ,1.2271 + ,108.79 + ,38.282 + ,1.2328 + ,108.14 + ,38.581 + ,1.2548 + ,109.88 + ,39.0801 + ,1.2531 + ,109.93 + ,39.0387 + ,1.2579 + ,110.46 + ,39.1015 + ,1.2567 + ,109.56 + ,39.1503 + ,1.266 + ,111.49 + ,39.14 + ,1.2637 + ,111.85 + ,39.0275 + ,1.2572 + ,111.35 + ,38.7665 + ,1.2569 + ,110.95 + ,38.691 + ,1.2703 + ,112.49 + ,38.849 + ,1.2828 + ,113.11 + ,39.1644 + ,1.3 + ,112.54 + ,39.4907 + ,1.2957 + ,112.84 + ,39.5095 + ,1.2844 + ,111.5 + ,39.2795 + ,1.2817 + ,111.52 + ,39.0437 + ,1.285 + ,111.57 + ,39.1355 + ,1.2897 + ,112.48 + ,39.143 + ,1.2931 + ,112.31 + ,39.185 + ,1.3033 + ,113.79 + ,39.355 + ,1.2992 + ,114.01 + ,39.297 + ,1.3069 + ,113.64 + ,39.4514 + ,1.3028 + ,112.62 + ,39.4173 + ,1.3073 + ,113.27 + ,39.4305 + ,1.3221 + ,113.51 + ,39.384 + ,1.3206 + ,112.92 + ,39.3261 + ,1.3184 + ,113.66 + ,39.301 + ,1.3176 + ,113.14 + ,39.35 + ,1.3253 + ,113.48 + ,39.64 + ,1.3133 + ,113.23 + ,39.4723 + ,1.3016 + ,110.56 + ,39.3685 + ,1.279 + ,109.5 + ,39.1906 + ,1.2799 + ,109.78 + ,39.1183 + ,1.282 + ,109.49 + ,39.1325 + ,1.286 + ,109.66 + ,39.1144 + ,1.288 + ,109.93 + ,39.1614 + ,1.2836 + ,109.82 + ,39.0908 + ,1.2711 + ,108.54 + ,38.9199 + ,1.2704 + ,108.23 + ,38.913 + ,1.2611 + ,106.19 + ,38.9655 + ,1.2613 + ,106.49 + ,39.029 + ,1.2693 + ,107.15 + ,39.089 + ,1.2713 + ,107.74 + ,39.07 + ,1.27 + ,107.54 + ,39.0046 + ,1.268 + ,107.07 + ,39.1038 + ,1.28 + ,107.54 + ,39.3572 + ,1.2818 + ,107.81 + ,39.388 + ,1.2834 + ,108.38 + ,39.382 + ,1.2874 + ,108.42 + ,39.4398 + ,1.2744 + ,106.86 + ,39.2537 + ,1.2697 + ,106.41 + ,39.2301 + ,1.2715 + ,106.46 + ,39.2763 + ,1.2725 + ,106.84 + ,39.282 + ,1.2801 + ,107.69 + ,39.3325 + ,1.285 + ,107.04 + ,39.557 + ,1.2989 + ,111.04 + ,40.1 + ,1.3078 + ,111.93 + ,40.5875 + ,1.306 + ,111.98 + ,40.485 + ,1.3074 + ,112.07 + ,40.55 + ,1.312 + ,112.05 + ,40.7955 + ,1.3364 + ,113.14 + ,41.456 + ,1.3323 + ,112.49 + ,41.3557 + ,1.3412 + ,113.2 + ,41.374 + ,1.3477 + ,113.52 + ,41.2235 + ,1.346 + ,113.22 + ,41.15 + ,1.3611 + ,113.85 + ,41.3725 + ,1.3648 + ,113.68 + ,41.6923 + ,1.3726 + ,114.26 + ,41.8 + ,1.3705 + ,114.1 + ,41.8045 + ,1.378 + ,114.8 + ,41.64 + ,1.3856 + ,114.98 + ,41.36 + ,1.397 + ,115.1 + ,41.5745 + ,1.3874 + ,114.21 + ,41.593 + ,1.3936 + ,114.24 + ,41.575 + ,1.3833 + ,113.35 + ,41.68 + ,1.3958 + ,114.23 + ,42.0055 + ,1.4101 + ,114.43 + ,42.3188 + ,1.4089 + ,114.28 + ,42.565 + ,1.3896 + ,113 + ,42.3575 + ,1.3859 + ,113.16 + ,42.29 + ,1.3861 + ,112.59 + ,42.695 + ,1.4016 + ,113.65 + ,43.0028 + ,1.3934 + ,113.18 + ,42.4507 + ,1.4031 + ,113.21 + ,42.4705 + ,1.3912 + ,113.11 + ,42.2875 + ,1.3803 + ,112.78 + ,42.3172 + ,1.3857 + ,112.57 + ,42.55 + ,1.3857 + ,111.87 + ,42.7523 + ,1.3926 + ,111.94 + ,42.8993 + ,1.4018 + ,113.18 + ,43.1555 + ,1.4014 + ,113.67 + ,43.1885 + ,1.4244 + ,115.15 + ,43.43 + ,1.4084 + ,114.41 + ,43.31 + ,1.3917 + ,112.88 + ,42.815 + ,1.3945 + ,112.44 + ,42.7017 + ,1.377 + ,113.48 + ,42.28 + ,1.37 + ,112.78 + ,41.922 + ,1.3711 + ,112.59 + ,42.17 + ,1.3626 + ,113.31 + ,42.1962 + ,1.3612 + ,113.21 + ,42.3215 + ,1.3481 + ,112.5 + ,42.3173 + ,1.3647 + ,113.72 + ,42.391 + ,1.3674 + ,114.09 + ,42.463 + ,1.3647 + ,113.97 + ,42.4125 + ,1.3496 + ,112.5 + ,42.304 + ,1.3339 + ,111.28 + ,41.813 + ,1.3321 + ,111.35 + ,41.651 + ,1.3225 + ,110.92 + ,41.539 + ,1.3146 + ,110.73 + ,41.1575 + ,1.2998 + ,109 + ,40.9545) + ,dim=c(3 + ,491) + ,dimnames=list(c('Dollar' + ,'Yen' + ,'Roebel') + ,1:491)) > y <- array(NA,dim=c(3,491),dimnames=list(c('Dollar','Yen','Roebel'),1:491)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '3' > #'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: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Roebel Dollar Yen 1 40.7819 1.3866 126.64 2 39.5915 1.3582 126.81 3 38.8859 1.3332 125.84 4 39.9068 1.3595 126.77 5 41.4700 1.3617 124.34 6 41.5613 1.3684 124.40 7 41.6005 1.3394 120.48 8 41.4113 1.3262 118.54 9 41.8400 1.3173 117.66 10 42.2892 1.3085 116.97 11 43.1521 1.3270 120.11 12 43.5998 1.3182 119.16 13 43.1160 1.2930 116.90 14 42.4185 1.2910 116.11 15 42.3687 1.2984 114.98 16 42.2975 1.2795 113.65 17 42.8528 1.2990 115.82 18 43.5350 1.3174 117.59 19 44.7265 1.3260 118.57 20 45.7293 1.3111 118.07 21 45.7585 1.2816 114.98 22 46.1685 1.2760 114.04 23 46.5075 1.2849 115.02 24 46.5270 1.2818 114.28 25 46.6010 1.2829 115.04 26 46.4607 1.2796 116.70 27 46.7135 1.3008 119.21 28 46.4113 1.2967 118.39 29 45.5500 1.2938 116.50 30 44.6081 1.2833 115.46 31 44.4395 1.2823 117.59 32 44.9847 1.2765 117.33 33 45.7558 1.2634 116.20 34 45.3942 1.2596 116.83 35 45.6970 1.2705 118.99 36 45.5664 1.2591 118.62 37 46.0205 1.2798 121.09 38 45.9195 1.2763 122.40 39 45.8005 1.2795 123.76 40 45.5350 1.2782 125.33 41 45.4977 1.2644 123.23 42 45.5782 1.2596 122.52 43 45.7697 1.2615 123.64 44 45.2445 1.2555 124.67 45 45.0615 1.2555 124.71 46 45.2865 1.2658 122.53 47 44.7910 1.2565 124.40 48 44.7625 1.2783 125.45 49 44.7644 1.2786 125.35 50 44.9973 1.2782 124.30 51 44.7265 1.2905 127.03 52 45.1465 1.3042 128.51 53 44.7465 1.2942 128.10 54 45.1795 1.3130 128.94 55 45.6515 1.3671 129.67 56 45.4920 1.3549 129.87 57 45.2775 1.3558 131.12 58 45.2115 1.3507 132.68 59 45.4110 1.3494 132.24 60 45.4005 1.3607 133.63 61 44.7692 1.3295 129.91 62 44.8913 1.3193 127.93 63 45.0320 1.3308 131.17 64 44.8790 1.3246 130.86 65 44.8330 1.3392 133.48 66 44.8257 1.3425 134.08 67 44.7815 1.3496 136.02 68 44.4790 1.3255 132.80 69 44.6317 1.3231 132.37 70 44.5043 1.3273 133.05 71 44.3217 1.3276 132.57 72 44.1005 1.3173 130.70 73 44.0470 1.3196 130.50 74 43.6835 1.3058 129.67 75 43.7864 1.2966 127.80 76 44.1807 1.2932 126.82 77 43.9595 1.2947 126.85 78 43.9370 1.3050 128.28 79 43.9910 1.3232 128.30 80 43.8650 1.3125 126.82 81 43.6710 1.2992 125.08 82 43.9300 1.3266 128.53 83 43.8630 1.3275 130.34 84 43.7095 1.3223 131.52 85 43.9435 1.3403 132.59 86 43.7360 1.3322 131.17 87 43.6295 1.3363 132.72 88 43.5980 1.3425 133.36 89 43.8726 1.3574 132.82 90 43.8935 1.3683 132.90 91 43.5957 1.3623 130.90 92 43.7155 1.3563 129.41 93 43.5280 1.3518 128.67 94 43.3415 1.3494 129.28 95 43.3374 1.3612 130.91 96 43.3320 1.3690 131.06 97 43.3869 1.3771 130.84 98 43.5016 1.3972 131.41 99 43.4875 1.4010 133.22 100 43.6023 1.3908 132.06 101 43.3886 1.3901 132.48 102 43.3105 1.3856 134.38 103 43.4455 1.4098 135.22 104 43.5185 1.4220 134.89 105 43.5755 1.4238 136.09 106 43.6217 1.4207 136.33 107 43.6440 1.4095 136.32 108 43.5789 1.4177 137.48 109 43.5215 1.3866 136.53 110 43.5033 1.3959 136.80 111 43.6320 1.4102 138.03 112 43.2630 1.3969 137.39 113 43.3717 1.4004 137.55 114 43.2745 1.3850 136.08 115 43.2647 1.3890 134.78 116 43.3240 1.3840 133.28 117 43.4455 1.3920 133.57 118 43.4098 1.3932 134.84 119 43.4100 1.3858 133.02 120 43.9300 1.3978 133.49 121 43.8104 1.4029 133.77 122 43.5400 1.3940 134.34 123 43.8580 1.4096 134.50 124 43.8375 1.4058 134.03 125 43.8810 1.4134 135.51 126 43.8870 1.4096 136.53 127 43.8009 1.4049 135.95 128 43.7877 1.4009 134.32 129 43.8110 1.3897 132.44 130 44.0625 1.4019 133.61 131 44.1250 1.3901 131.02 132 44.5200 1.3990 130.05 133 45.4005 1.3901 128.21 134 45.8900 1.3975 129.03 135 45.1890 1.3991 130.34 136 44.9035 1.4089 131.57 137 44.9351 1.4130 132.63 138 44.8010 1.4090 132.06 139 43.9800 1.4217 134.44 140 44.1100 1.4223 134.10 141 44.2661 1.4191 132.49 142 44.3610 1.4229 134.23 143 44.0990 1.4227 134.92 144 43.8435 1.4269 135.61 145 43.8914 1.4229 134.53 146 44.2170 1.4104 133.86 147 44.5060 1.4053 133.89 148 44.5400 1.4138 135.33 149 44.4465 1.4303 135.86 150 44.8420 1.4384 136.22 151 44.8946 1.4410 137.38 152 44.9510 1.4370 137.31 153 45.4450 1.4357 136.89 154 45.0035 1.4202 138.01 155 45.7690 1.4166 136.72 156 46.0900 1.4170 135.77 157 45.4120 1.4293 137.52 158 45.1200 1.4294 135.61 159 45.4800 1.4072 132.94 160 45.1050 1.4101 134.12 161 45.0560 1.4112 132.55 162 45.2200 1.4243 134.11 163 45.3900 1.4330 134.19 164 45.0410 1.4323 135.57 165 44.9399 1.4324 135.05 166 44.9315 1.4270 134.32 167 45.1935 1.4268 133.61 168 45.3466 1.4364 134.75 169 45.4645 1.4272 133.10 170 45.5685 1.4314 133.26 171 45.3921 1.4220 131.63 172 45.3400 1.4335 132.47 173 45.1308 1.4262 132.45 174 45.1005 1.4330 133.33 175 45.3700 1.4473 133.57 176 45.2000 1.4522 134.13 177 44.9614 1.4545 133.92 178 44.8015 1.4594 132.62 179 44.9152 1.4561 132.30 180 45.0950 1.4611 133.26 181 44.9271 1.4671 132.60 182 44.6026 1.4712 134.38 183 44.5000 1.4705 134.17 184 44.5400 1.4658 135.46 185 44.5532 1.4780 135.09 186 44.4070 1.4783 134.96 187 44.2590 1.4768 133.85 188 44.1365 1.4670 132.59 189 44.1120 1.4650 131.15 190 43.8814 1.4549 130.91 191 43.9800 1.4643 131.07 192 43.7294 1.4539 130.78 193 43.9119 1.4537 129.95 194 43.9550 1.4616 131.41 195 43.9000 1.4722 131.21 196 43.7065 1.4694 130.68 197 43.6939 1.4763 130.46 198 43.6587 1.4750 131.12 199 43.5885 1.4765 132.99 200 43.8885 1.4864 133.02 201 43.8216 1.4881 133.39 202 43.7510 1.4864 134.07 203 43.6990 1.4869 135.60 204 43.7425 1.4918 135.66 205 43.6390 1.4971 135.53 206 43.5890 1.4921 135.82 207 43.6060 1.5000 136.90 208 43.5325 1.5020 137.97 209 43.3850 1.5019 138.09 210 43.3745 1.4874 136.91 211 43.2360 1.4785 134.76 212 43.1957 1.4788 135.13 213 43.0100 1.4800 134.66 214 43.1401 1.4772 132.95 215 43.0487 1.4658 132.25 216 43.1972 1.4761 134.30 217 43.2461 1.4867 134.30 218 43.0866 1.4862 134.76 219 43.0865 1.4984 134.81 220 43.0194 1.4966 134.51 221 43.0800 1.5037 135.11 222 43.0070 1.4922 134.32 223 42.9278 1.4868 133.51 224 42.9545 1.4965 134.02 225 42.7995 1.4875 132.76 226 42.9048 1.4957 133.39 227 42.9468 1.4863 132.05 228 43.0800 1.4815 131.87 229 43.1274 1.4968 133.03 230 43.1625 1.4969 132.57 231 43.4500 1.5083 132.10 232 43.8310 1.5071 130.70 233 43.7769 1.4918 129.20 234 43.9800 1.5023 129.77 235 43.9200 1.5074 131.02 236 44.1100 1.5090 131.55 237 44.0300 1.5120 133.17 238 44.1582 1.5068 133.08 239 44.1400 1.4787 133.24 240 45.0700 1.4774 130.74 241 44.8737 1.4768 129.91 242 44.8505 1.4730 130.03 243 44.3730 1.4757 131.13 244 44.0750 1.4647 129.55 245 43.9725 1.4541 130.22 246 44.0940 1.4560 130.61 247 44.1910 1.4343 129.27 248 43.9685 1.4337 129.68 249 43.7900 1.4368 130.10 250 43.6041 1.4279 130.83 251 43.1707 1.4276 130.95 252 42.7100 1.4398 131.73 253 42.7550 1.4405 131.86 254 43.3316 1.4433 132.44 255 43.5000 1.4338 132.35 256 43.1540 1.4406 133.16 257 43.1600 1.4389 133.62 258 43.1000 1.4442 132.54 259 42.8500 1.4350 132.69 260 42.6175 1.4304 133.50 261 42.5000 1.4273 133.36 262 42.6285 1.4528 134.23 263 42.6974 1.4481 132.41 264 43.0400 1.4563 133.02 265 42.6730 1.4486 132.88 266 42.5015 1.4374 130.76 267 42.5380 1.4369 130.33 268 42.3735 1.4279 129.79 269 42.0140 1.4132 128.65 270 41.8618 1.4064 129.14 271 42.1824 1.4135 127.35 272 42.6050 1.4151 127.74 273 42.7345 1.4085 126.31 274 42.6150 1.4072 125.95 275 42.4650 1.3999 126.36 276 42.3400 1.3966 126.15 277 42.2510 1.3913 125.60 278 42.0475 1.3937 126.20 279 41.8600 1.3984 126.73 280 41.6850 1.3847 125.68 281 41.7350 1.3691 122.49 282 41.7060 1.3675 122.07 283 41.7640 1.3760 123.40 284 41.5800 1.3740 123.01 285 41.3730 1.3718 123.03 286 41.0880 1.3572 122.33 287 41.1370 1.3607 122.42 288 41.1587 1.3649 122.68 289 41.1850 1.3726 124.69 290 40.8190 1.3567 123.30 291 40.6330 1.3519 124.17 292 40.8580 1.3626 124.38 293 40.7940 1.3577 123.19 294 40.6900 1.3547 122.16 295 40.5950 1.3489 120.66 296 40.7305 1.3570 120.92 297 40.5471 1.3525 120.67 298 40.5145 1.3548 120.68 299 40.7000 1.3641 121.10 300 40.7000 1.3668 120.86 301 40.5220 1.3582 121.48 302 40.6165 1.3662 123.48 303 40.3985 1.3557 121.72 304 40.2815 1.3610 123.16 305 40.2450 1.3657 123.84 306 40.3055 1.3765 124.57 307 40.2696 1.3705 124.30 308 40.2510 1.3723 124.22 309 40.1270 1.3756 124.43 310 39.9500 1.3660 123.33 311 39.6750 1.3548 122.86 312 39.9540 1.3471 121.25 313 39.8828 1.3519 122.16 314 39.6200 1.3338 122.62 315 39.5415 1.3356 123.44 316 39.5250 1.3353 124.00 317 39.8145 1.3471 124.75 318 39.6675 1.3482 124.80 319 39.6950 1.3479 125.93 320 39.5985 1.3468 126.28 321 39.2735 1.3396 126.04 322 39.1435 1.3340 125.04 323 39.1742 1.3296 123.76 324 39.2025 1.3384 125.34 325 39.3946 1.3585 126.99 326 39.5025 1.3583 126.34 327 39.4845 1.3615 127.42 328 39.3300 1.3544 126.18 329 39.2950 1.3535 125.30 330 39.2675 1.3432 123.50 331 39.2535 1.3486 125.32 332 38.9845 1.3373 124.65 333 38.9285 1.3339 124.03 334 38.8592 1.3311 125.11 335 38.7700 1.3321 125.46 336 38.7900 1.3290 124.70 337 38.8205 1.3245 124.48 338 38.7577 1.3256 124.76 339 38.8390 1.3315 125.81 340 38.7800 1.3238 124.95 341 38.5400 1.3089 123.66 342 38.5110 1.2924 122.66 343 38.6150 1.2727 119.34 344 38.8980 1.2746 117.84 345 38.8691 1.2969 120.97 346 38.3840 1.2698 117.38 347 38.0277 1.2686 118.06 348 37.7200 1.2587 116.99 349 37.7325 1.2492 115.55 350 37.6260 1.2349 114.17 351 37.6030 1.2428 115.32 352 37.7800 1.2270 112.49 353 38.5590 1.2334 111.93 354 39.0459 1.2497 112.08 355 38.4500 1.2360 111.63 356 38.5050 1.2223 109.53 357 38.2885 1.2309 111.35 358 37.7950 1.2255 110.79 359 37.9200 1.2384 113.06 360 38.0340 1.2307 112.62 361 38.0290 1.2155 110.65 362 38.0630 1.2218 112.36 363 37.9828 1.2268 113.74 364 37.7450 1.2060 111.73 365 37.9690 1.1959 109.86 366 38.0070 1.1942 109.32 367 38.0615 1.2010 109.99 368 38.0912 1.2045 109.84 369 38.0910 1.2127 111.13 370 38.4310 1.2249 112.43 371 38.4800 1.2258 111.77 372 38.3500 1.2277 112.15 373 38.2140 1.2363 112.89 374 38.3840 1.2372 112.12 375 38.1375 1.2391 113.10 376 38.0075 1.2258 111.09 377 38.0524 1.2271 110.76 378 38.2350 1.2262 109.59 379 38.3100 1.2294 109.99 380 38.2615 1.2339 110.25 381 38.1300 1.2198 108.31 382 38.2820 1.2271 108.79 383 38.5810 1.2328 108.14 384 39.0801 1.2548 109.88 385 39.0387 1.2531 109.93 386 39.1015 1.2579 110.46 387 39.1503 1.2567 109.56 388 39.1400 1.2660 111.49 389 39.0275 1.2637 111.85 390 38.7665 1.2572 111.35 391 38.6910 1.2569 110.95 392 38.8490 1.2703 112.49 393 39.1644 1.2828 113.11 394 39.4907 1.3000 112.54 395 39.5095 1.2957 112.84 396 39.2795 1.2844 111.50 397 39.0437 1.2817 111.52 398 39.1355 1.2850 111.57 399 39.1430 1.2897 112.48 400 39.1850 1.2931 112.31 401 39.3550 1.3033 113.79 402 39.2970 1.2992 114.01 403 39.4514 1.3069 113.64 404 39.4173 1.3028 112.62 405 39.4305 1.3073 113.27 406 39.3840 1.3221 113.51 407 39.3261 1.3206 112.92 408 39.3010 1.3184 113.66 409 39.3500 1.3176 113.14 410 39.6400 1.3253 113.48 411 39.4723 1.3133 113.23 412 39.3685 1.3016 110.56 413 39.1906 1.2790 109.50 414 39.1183 1.2799 109.78 415 39.1325 1.2820 109.49 416 39.1144 1.2860 109.66 417 39.1614 1.2880 109.93 418 39.0908 1.2836 109.82 419 38.9199 1.2711 108.54 420 38.9130 1.2704 108.23 421 38.9655 1.2611 106.19 422 39.0290 1.2613 106.49 423 39.0890 1.2693 107.15 424 39.0700 1.2713 107.74 425 39.0046 1.2700 107.54 426 39.1038 1.2680 107.07 427 39.3572 1.2800 107.54 428 39.3880 1.2818 107.81 429 39.3820 1.2834 108.38 430 39.4398 1.2874 108.42 431 39.2537 1.2744 106.86 432 39.2301 1.2697 106.41 433 39.2763 1.2715 106.46 434 39.2820 1.2725 106.84 435 39.3325 1.2801 107.69 436 39.5570 1.2850 107.04 437 40.1000 1.2989 111.04 438 40.5875 1.3078 111.93 439 40.4850 1.3060 111.98 440 40.5500 1.3074 112.07 441 40.7955 1.3120 112.05 442 41.4560 1.3364 113.14 443 41.3557 1.3323 112.49 444 41.3740 1.3412 113.20 445 41.2235 1.3477 113.52 446 41.1500 1.3460 113.22 447 41.3725 1.3611 113.85 448 41.6923 1.3648 113.68 449 41.8000 1.3726 114.26 450 41.8045 1.3705 114.10 451 41.6400 1.3780 114.80 452 41.3600 1.3856 114.98 453 41.5745 1.3970 115.10 454 41.5930 1.3874 114.21 455 41.5750 1.3936 114.24 456 41.6800 1.3833 113.35 457 42.0055 1.3958 114.23 458 42.3188 1.4101 114.43 459 42.5650 1.4089 114.28 460 42.3575 1.3896 113.00 461 42.2900 1.3859 113.16 462 42.6950 1.3861 112.59 463 43.0028 1.4016 113.65 464 42.4507 1.3934 113.18 465 42.4705 1.4031 113.21 466 42.2875 1.3912 113.11 467 42.3172 1.3803 112.78 468 42.5500 1.3857 112.57 469 42.7523 1.3857 111.87 470 42.8993 1.3926 111.94 471 43.1555 1.4018 113.18 472 43.1885 1.4014 113.67 473 43.4300 1.4244 115.15 474 43.3100 1.4084 114.41 475 42.8150 1.3917 112.88 476 42.7017 1.3945 112.44 477 42.2800 1.3770 113.48 478 41.9220 1.3700 112.78 479 42.1700 1.3711 112.59 480 42.1962 1.3626 113.31 481 42.3215 1.3612 113.21 482 42.3173 1.3481 112.50 483 42.3910 1.3647 113.72 484 42.4630 1.3674 114.09 485 42.4125 1.3647 113.97 486 42.3040 1.3496 112.50 487 41.8130 1.3339 111.28 488 41.6510 1.3321 111.35 489 41.5390 1.3225 110.92 490 41.1575 1.3146 110.73 491 40.9545 1.2998 109.00 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Dollar Yen 19.0193 4.0538 0.1431 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -3.5966 -1.1512 -0.1688 1.1625 5.9636 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 19.01933 1.45278 13.092 <2e-16 *** Dollar 4.05376 1.59126 2.548 0.0112 * Yen 0.14305 0.01298 11.020 <2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 1.803 on 488 degrees of freedom Multiple R-squared: 0.437, Adjusted R-squared: 0.4346 F-statistic: 189.4 on 2 and 488 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 0.0026666094 5.333219e-03 9.973334e-01 [2,] 0.0022358579 4.471716e-03 9.977641e-01 [3,] 0.0008667966 1.733593e-03 9.991332e-01 [4,] 0.0001689667 3.379333e-04 9.998310e-01 [5,] 0.0001024587 2.049174e-04 9.998975e-01 [6,] 0.0036008232 7.201646e-03 9.963992e-01 [7,] 0.0127634132 2.552683e-02 9.872366e-01 [8,] 0.0082973085 1.659462e-02 9.917027e-01 [9,] 0.0039323630 7.864726e-03 9.960676e-01 [10,] 0.0027291378 5.458276e-03 9.972709e-01 [11,] 0.0017294849 3.458970e-03 9.982705e-01 [12,] 0.0007830174 1.566035e-03 9.992170e-01 [13,] 0.0006139249 1.227850e-03 9.993861e-01 [14,] 0.0029020568 5.804114e-03 9.970979e-01 [15,] 0.0374410160 7.488203e-02 9.625590e-01 [16,] 0.0816696461 1.633393e-01 9.183304e-01 [17,] 0.1267665762 2.535332e-01 8.732334e-01 [18,] 0.2082555905 4.165112e-01 7.917444e-01 [19,] 0.2613762998 5.227526e-01 7.386237e-01 [20,] 0.3335484144 6.670968e-01 6.664516e-01 [21,] 0.4172763406 8.345527e-01 5.827237e-01 [22,] 0.6014830368 7.970339e-01 3.985170e-01 [23,] 0.6763156277 6.473687e-01 3.236844e-01 [24,] 0.6681362901 6.637274e-01 3.318637e-01 [25,] 0.6402187650 7.195625e-01 3.597812e-01 [26,] 0.6077662251 7.844675e-01 3.922338e-01 [27,] 0.5784053068 8.431894e-01 4.215947e-01 [28,] 0.5629884286 8.740231e-01 4.370116e-01 [29,] 0.5475220654 9.049559e-01 4.524779e-01 [30,] 0.5360876025 9.278248e-01 4.639124e-01 [31,] 0.5208458141 9.583084e-01 4.791542e-01 [32,] 0.5437860690 9.124279e-01 4.562139e-01 [33,] 0.5460628578 9.078743e-01 4.539371e-01 [34,] 0.5414919183 9.170162e-01 4.585081e-01 [35,] 0.5216067525 9.567865e-01 4.783932e-01 [36,] 0.5103876237 9.792248e-01 4.896124e-01 [37,] 0.5070187740 9.859625e-01 4.929812e-01 [38,] 0.5013802974 9.972394e-01 4.986197e-01 [39,] 0.5036172726 9.927655e-01 4.963827e-01 [40,] 0.5064472221 9.871056e-01 4.935528e-01 [41,] 0.5025820047 9.948360e-01 4.974180e-01 [42,] 0.5100093258 9.799813e-01 4.899907e-01 [43,] 0.4925378159 9.850756e-01 5.074622e-01 [44,] 0.4768712092 9.537424e-01 5.231288e-01 [45,] 0.4703816771 9.407634e-01 5.296183e-01 [46,] 0.4621884348 9.243769e-01 5.378116e-01 [47,] 0.5105653658 9.788693e-01 4.894346e-01 [48,] 0.5040778840 9.918442e-01 4.959221e-01 [49,] 0.5704218791 8.591562e-01 4.295781e-01 [50,] 0.8753432654 2.493135e-01 1.246567e-01 [51,] 0.9465531264 1.068937e-01 5.344687e-02 [52,] 0.9706143018 5.877140e-02 2.938570e-02 [53,] 0.9787291389 4.254172e-02 2.127086e-02 [54,] 0.9849254792 3.014904e-02 1.507452e-02 [55,] 0.9894471568 2.110569e-02 1.055284e-02 [56,] 0.9889029483 2.219410e-02 1.109705e-02 [57,] 0.9889986022 2.200280e-02 1.100140e-02 [58,] 0.9889827999 2.203440e-02 1.101720e-02 [59,] 0.9885491003 2.290180e-02 1.145090e-02 [60,] 0.9876520309 2.469594e-02 1.234797e-02 [61,] 0.9866284423 2.674312e-02 1.337156e-02 [62,] 0.9850702846 2.985943e-02 1.492972e-02 [63,] 0.9838171234 3.236575e-02 1.618288e-02 [64,] 0.9829462482 3.410750e-02 1.705375e-02 [65,] 0.9816725978 3.665480e-02 1.832740e-02 [66,] 0.9802542950 3.949141e-02 1.974571e-02 [67,] 0.9799077463 4.018451e-02 2.009225e-02 [68,] 0.9793262784 4.134744e-02 2.067372e-02 [69,] 0.9816088546 3.678229e-02 1.839115e-02 [70,] 0.9846270958 3.074581e-02 1.537290e-02 [71,] 0.9872541070 2.549179e-02 1.274589e-02 [72,] 0.9895298238 2.094035e-02 1.047018e-02 [73,] 0.9906547615 1.869048e-02 9.345238e-03 [74,] 0.9906307771 1.873845e-02 9.369223e-03 [75,] 0.9912343596 1.753128e-02 8.765640e-03 [76,] 0.9928077560 1.438449e-02 7.192244e-03 [77,] 0.9928328466 1.433431e-02 7.167153e-03 [78,] 0.9927679899 1.446402e-02 7.232010e-03 [79,] 0.9930000529 1.399989e-02 6.999947e-03 [80,] 0.9926231930 1.475361e-02 7.376807e-03 [81,] 0.9924406735 1.511865e-02 7.559327e-03 [82,] 0.9920540989 1.589180e-02 7.945901e-03 [83,] 0.9913858316 1.722834e-02 8.614168e-03 [84,] 0.9908759627 1.824807e-02 9.124037e-03 [85,] 0.9906166201 1.876676e-02 9.383380e-03 [86,] 0.9897625917 2.047482e-02 1.023741e-02 [87,] 0.9893586812 2.128264e-02 1.064132e-02 [88,] 0.9886380859 2.272383e-02 1.136191e-02 [89,] 0.9875990343 2.480193e-02 1.240097e-02 [90,] 0.9860262235 2.794755e-02 1.397378e-02 [91,] 0.9842342478 3.153150e-02 1.576575e-02 [92,] 0.9826906818 3.461864e-02 1.730932e-02 [93,] 0.9829032995 3.419340e-02 1.709670e-02 [94,] 0.9823531437 3.529371e-02 1.764686e-02 [95,] 0.9815922530 3.681549e-02 1.840775e-02 [96,] 0.9794626985 4.107460e-02 2.053730e-02 [97,] 0.9762050476 4.758990e-02 2.379495e-02 [98,] 0.9744372814 5.112544e-02 2.556272e-02 [99,] 0.9744792545 5.104149e-02 2.552075e-02 [100,] 0.9739078426 5.218431e-02 2.609216e-02 [101,] 0.9725476228 5.490475e-02 2.745238e-02 [102,] 0.9697462351 6.050753e-02 3.025376e-02 [103,] 0.9664176964 6.716461e-02 3.358230e-02 [104,] 0.9617004966 7.659901e-02 3.829950e-02 [105,] 0.9561391535 8.772169e-02 4.386085e-02 [106,] 0.9507611493 9.847770e-02 4.923885e-02 [107,] 0.9431214213 1.137572e-01 5.687858e-02 [108,] 0.9347699472 1.304601e-01 6.523005e-02 [109,] 0.9265946410 1.468107e-01 7.340536e-02 [110,] 0.9172330131 1.655340e-01 8.276699e-02 [111,] 0.9081108594 1.837783e-01 9.188914e-02 [112,] 0.8987150025 2.025700e-01 1.012850e-01 [113,] 0.8877439434 2.245121e-01 1.122561e-01 [114,] 0.8773289696 2.453421e-01 1.226710e-01 [115,] 0.8752589666 2.494821e-01 1.247410e-01 [116,] 0.8706063081 2.587874e-01 1.293937e-01 [117,] 0.8596740301 2.806519e-01 1.403260e-01 [118,] 0.8557572538 2.884855e-01 1.442427e-01 [119,] 0.8507921855 2.984156e-01 1.492078e-01 [120,] 0.8458143397 3.083713e-01 1.541857e-01 [121,] 0.8384512138 3.230976e-01 1.615488e-01 [122,] 0.8290778032 3.418444e-01 1.709222e-01 [123,] 0.8207505689 3.584989e-01 1.792494e-01 [124,] 0.8152075426 3.695849e-01 1.847925e-01 [125,] 0.8155197762 3.689604e-01 1.844802e-01 [126,] 0.8215112176 3.569776e-01 1.784888e-01 [127,] 0.8457223152 3.085554e-01 1.542777e-01 [128,] 0.9087881920 1.824236e-01 9.121181e-02 [129,] 0.9628856250 7.422875e-02 3.711437e-02 [130,] 0.9777417126 4.451657e-02 2.225829e-02 [131,] 0.9843579497 3.128410e-02 1.564205e-02 [132,] 0.9890639764 2.187205e-02 1.093602e-02 [133,] 0.9919388261 1.612235e-02 8.061174e-03 [134,] 0.9914152825 1.716943e-02 8.584717e-03 [135,] 0.9912003039 1.759939e-02 8.799696e-03 [136,] 0.9915636310 1.687274e-02 8.436369e-03 [137,] 0.9920015804 1.599684e-02 7.998420e-03 [138,] 0.9916709449 1.665811e-02 8.329055e-03 [139,] 0.9906853796 1.862924e-02 9.314620e-03 [140,] 0.9898067902 2.038642e-02 1.019321e-02 [141,] 0.9901158886 1.976822e-02 9.884111e-03 [142,] 0.9915772402 1.684552e-02 8.422760e-03 [143,] 0.9927091474 1.458171e-02 7.290853e-03 [144,] 0.9931518207 1.369636e-02 6.848179e-03 [145,] 0.9945038316 1.099234e-02 5.496168e-03 [146,] 0.9956254598 8.749080e-03 4.374540e-03 [147,] 0.9966631808 6.673638e-03 3.336819e-03 [148,] 0.9981979624 3.604075e-03 1.802038e-03 [149,] 0.9988163249 2.367350e-03 1.183675e-03 [150,] 0.9996307238 7.385524e-04 3.692762e-04 [151,] 0.9999349536 1.300927e-04 6.504636e-05 [152,] 0.9999748470 5.030596e-05 2.515298e-05 [153,] 0.9999877712 2.445760e-05 1.222880e-05 [154,] 0.9999974264 5.147289e-06 2.573644e-06 [155,] 0.9999992267 1.546536e-06 7.732680e-07 [156,] 0.9999997706 4.588499e-07 2.294250e-07 [157,] 0.9999999379 1.241880e-07 6.209398e-08 [158,] 0.9999999859 2.815717e-08 1.407859e-08 [159,] 0.9999999955 9.095481e-09 4.547740e-09 [160,] 0.9999999984 3.155722e-09 1.577861e-09 [161,] 0.9999999996 8.974587e-10 4.487294e-10 [162,] 0.9999999999 1.482803e-10 7.414017e-11 [163,] 1.0000000000 1.998267e-11 9.991337e-12 [164,] 1.0000000000 1.264649e-12 6.323245e-13 [165,] 1.0000000000 5.309827e-14 2.654913e-14 [166,] 1.0000000000 1.740517e-15 8.702586e-16 [167,] 1.0000000000 8.204818e-17 4.102409e-17 [168,] 1.0000000000 3.163380e-18 1.581690e-18 [169,] 1.0000000000 1.298885e-19 6.494423e-20 [170,] 1.0000000000 4.385982e-21 2.192991e-21 [171,] 1.0000000000 2.515075e-22 1.257538e-22 [172,] 1.0000000000 2.975891e-23 1.487946e-23 [173,] 1.0000000000 8.186113e-24 4.093057e-24 [174,] 1.0000000000 1.298464e-24 6.492318e-25 [175,] 1.0000000000 1.248876e-25 6.244381e-26 [176,] 1.0000000000 3.058699e-26 1.529350e-26 [177,] 1.0000000000 1.498817e-26 7.494087e-27 [178,] 1.0000000000 8.662831e-27 4.331416e-27 [179,] 1.0000000000 2.588883e-27 1.294441e-27 [180,] 1.0000000000 1.586309e-27 7.931544e-28 [181,] 1.0000000000 1.287313e-27 6.436567e-28 [182,] 1.0000000000 1.392717e-27 6.963585e-28 [183,] 1.0000000000 1.353819e-27 6.769094e-28 [184,] 1.0000000000 1.465372e-27 7.326858e-28 [185,] 1.0000000000 1.416660e-27 7.083300e-28 [186,] 1.0000000000 1.725975e-27 8.629875e-28 [187,] 1.0000000000 1.898431e-27 9.492156e-28 [188,] 1.0000000000 1.811785e-27 9.058924e-28 [189,] 1.0000000000 1.886157e-27 9.430783e-28 [190,] 1.0000000000 3.016189e-27 1.508094e-27 [191,] 1.0000000000 5.205113e-27 2.602556e-27 [192,] 1.0000000000 9.938801e-27 4.969400e-27 [193,] 1.0000000000 1.843520e-26 9.217599e-27 [194,] 1.0000000000 3.282721e-26 1.641360e-26 [195,] 1.0000000000 6.148164e-26 3.074082e-26 [196,] 1.0000000000 1.175475e-25 5.877377e-26 [197,] 1.0000000000 2.185099e-25 1.092550e-25 [198,] 1.0000000000 3.889469e-25 1.944734e-25 [199,] 1.0000000000 7.239788e-25 3.619894e-25 [200,] 1.0000000000 1.370518e-24 6.852588e-25 [201,] 1.0000000000 2.546686e-24 1.273343e-24 [202,] 1.0000000000 4.762270e-24 2.381135e-24 [203,] 1.0000000000 8.692595e-24 4.346297e-24 [204,] 1.0000000000 1.474532e-23 7.372662e-24 [205,] 1.0000000000 2.550253e-23 1.275127e-23 [206,] 1.0000000000 4.225570e-23 2.112785e-23 [207,] 1.0000000000 6.901136e-23 3.450568e-23 [208,] 1.0000000000 1.001631e-22 5.008153e-23 [209,] 1.0000000000 1.591400e-22 7.957001e-23 [210,] 1.0000000000 2.472186e-22 1.236093e-22 [211,] 1.0000000000 4.103358e-22 2.051679e-22 [212,] 1.0000000000 6.732447e-22 3.366224e-22 [213,] 1.0000000000 9.955987e-22 4.977993e-22 [214,] 1.0000000000 1.141790e-21 5.708952e-22 [215,] 1.0000000000 1.212189e-21 6.060947e-22 [216,] 1.0000000000 1.159949e-21 5.799746e-22 [217,] 1.0000000000 1.321668e-21 6.608341e-22 [218,] 1.0000000000 1.462406e-21 7.312032e-22 [219,] 1.0000000000 1.263339e-21 6.316695e-22 [220,] 1.0000000000 1.006988e-21 5.034938e-22 [221,] 1.0000000000 7.257141e-22 3.628571e-22 [222,] 1.0000000000 6.936748e-22 3.468374e-22 [223,] 1.0000000000 9.044975e-22 4.522487e-22 [224,] 1.0000000000 8.237857e-22 4.118928e-22 [225,] 1.0000000000 7.275740e-22 3.637870e-22 [226,] 1.0000000000 5.509873e-22 2.754937e-22 [227,] 1.0000000000 6.040611e-22 3.020306e-22 [228,] 1.0000000000 8.940658e-22 4.470329e-22 [229,] 1.0000000000 1.173401e-21 5.867005e-22 [230,] 1.0000000000 1.358342e-21 6.791711e-22 [231,] 1.0000000000 1.794668e-21 8.973341e-22 [232,] 1.0000000000 2.309903e-21 1.154952e-21 [233,] 1.0000000000 3.596346e-21 1.798173e-21 [234,] 1.0000000000 5.483559e-21 2.741779e-21 [235,] 1.0000000000 2.007639e-21 1.003820e-21 [236,] 1.0000000000 1.207150e-21 6.035752e-22 [237,] 1.0000000000 5.927225e-22 2.963612e-22 [238,] 1.0000000000 6.883045e-22 3.441522e-22 [239,] 1.0000000000 9.274847e-22 4.637424e-22 [240,] 1.0000000000 8.699361e-22 4.349680e-22 [241,] 1.0000000000 6.654666e-22 3.327333e-22 [242,] 1.0000000000 1.044569e-22 5.222846e-23 [243,] 1.0000000000 2.263476e-23 1.131738e-23 [244,] 1.0000000000 7.999671e-24 3.999835e-24 [245,] 1.0000000000 1.610170e-24 8.050852e-25 [246,] 1.0000000000 7.444748e-25 3.722374e-25 [247,] 1.0000000000 8.707281e-25 4.353640e-25 [248,] 1.0000000000 1.013344e-24 5.066718e-25 [249,] 1.0000000000 6.647955e-25 3.323977e-25 [250,] 1.0000000000 1.297107e-25 6.485537e-26 [251,] 1.0000000000 7.063908e-26 3.531954e-26 [252,] 1.0000000000 2.741389e-26 1.370694e-26 [253,] 1.0000000000 2.018854e-26 1.009427e-26 [254,] 1.0000000000 1.129540e-26 5.647699e-27 [255,] 1.0000000000 5.210447e-27 2.605223e-27 [256,] 1.0000000000 2.235476e-27 1.117738e-27 [257,] 1.0000000000 2.598034e-27 1.299017e-27 [258,] 1.0000000000 3.051996e-27 1.525998e-27 [259,] 1.0000000000 3.364119e-27 1.682060e-27 [260,] 1.0000000000 3.707273e-27 1.853637e-27 [261,] 1.0000000000 3.875951e-27 1.937976e-27 [262,] 1.0000000000 4.095949e-27 2.047975e-27 [263,] 1.0000000000 3.567012e-27 1.783506e-27 [264,] 1.0000000000 2.391915e-27 1.195957e-27 [265,] 1.0000000000 1.169636e-27 5.848179e-28 [266,] 1.0000000000 8.957943e-28 4.478972e-28 [267,] 1.0000000000 3.918948e-28 1.959474e-28 [268,] 1.0000000000 1.058671e-28 5.293354e-29 [269,] 1.0000000000 3.464938e-29 1.732469e-29 [270,] 1.0000000000 6.496167e-30 3.248084e-30 [271,] 1.0000000000 1.161961e-30 5.809803e-31 [272,] 1.0000000000 1.681770e-31 8.408849e-32 [273,] 1.0000000000 3.652320e-32 1.826160e-32 [274,] 1.0000000000 1.333003e-32 6.665016e-33 [275,] 1.0000000000 2.966988e-33 1.483494e-33 [276,] 1.0000000000 4.606454e-34 2.303227e-34 [277,] 1.0000000000 7.453966e-35 3.726983e-35 [278,] 1.0000000000 1.312741e-35 6.563705e-36 [279,] 1.0000000000 3.303718e-36 1.651859e-36 [280,] 1.0000000000 9.936349e-37 4.968175e-37 [281,] 1.0000000000 1.865340e-37 9.326700e-38 [282,] 1.0000000000 4.287650e-38 2.143825e-38 [283,] 1.0000000000 1.241672e-38 6.208358e-39 [284,] 1.0000000000 2.974514e-39 1.487257e-39 [285,] 1.0000000000 6.035814e-40 3.017907e-40 [286,] 1.0000000000 7.708733e-41 3.854367e-41 [287,] 1.0000000000 1.604135e-41 8.020677e-42 [288,] 1.0000000000 4.108030e-42 2.054015e-42 [289,] 1.0000000000 1.485494e-42 7.427469e-43 [290,] 1.0000000000 7.009080e-43 3.504540e-43 [291,] 1.0000000000 4.561400e-43 2.280700e-43 [292,] 1.0000000000 2.933381e-43 1.466691e-43 [293,] 1.0000000000 2.166003e-43 1.083001e-43 [294,] 1.0000000000 2.080251e-43 1.040126e-43 [295,] 1.0000000000 2.271119e-43 1.135560e-43 [296,] 1.0000000000 1.871002e-43 9.355010e-44 [297,] 1.0000000000 1.357544e-43 6.787718e-44 [298,] 1.0000000000 1.080823e-43 5.404117e-44 [299,] 1.0000000000 7.836016e-44 3.918008e-44 [300,] 1.0000000000 5.780875e-44 2.890437e-44 [301,] 1.0000000000 4.684979e-44 2.342490e-44 [302,] 1.0000000000 3.871814e-44 1.935907e-44 [303,] 1.0000000000 3.309872e-44 1.654936e-44 [304,] 1.0000000000 2.328580e-44 1.164290e-44 [305,] 1.0000000000 1.440965e-44 7.204824e-45 [306,] 1.0000000000 6.934087e-45 3.467043e-45 [307,] 1.0000000000 6.362572e-45 3.181286e-45 [308,] 1.0000000000 5.256576e-45 2.628288e-45 [309,] 1.0000000000 3.149915e-45 1.574958e-45 [310,] 1.0000000000 1.763467e-45 8.817337e-46 [311,] 1.0000000000 9.402959e-46 4.701480e-46 [312,] 1.0000000000 6.527663e-46 3.263831e-46 [313,] 1.0000000000 4.708904e-46 2.354452e-46 [314,] 1.0000000000 2.873639e-46 1.436820e-46 [315,] 1.0000000000 1.657856e-46 8.289279e-47 [316,] 1.0000000000 8.408259e-47 4.204129e-47 [317,] 1.0000000000 4.549811e-47 2.274906e-47 [318,] 1.0000000000 3.033652e-47 1.516826e-47 [319,] 1.0000000000 1.974115e-47 9.870573e-48 [320,] 1.0000000000 1.411539e-47 7.057696e-48 [321,] 1.0000000000 1.230383e-47 6.151914e-48 [322,] 1.0000000000 1.031939e-47 5.159694e-48 [323,] 1.0000000000 8.008696e-48 4.004348e-48 [324,] 1.0000000000 5.719865e-48 2.859932e-48 [325,] 1.0000000000 4.762025e-48 2.381013e-48 [326,] 1.0000000000 3.803296e-48 1.901648e-48 [327,] 1.0000000000 2.492313e-48 1.246157e-48 [328,] 1.0000000000 1.557104e-48 7.785522e-49 [329,] 1.0000000000 1.041584e-48 5.207920e-49 [330,] 1.0000000000 5.690901e-49 2.845451e-49 [331,] 1.0000000000 3.377480e-49 1.688740e-49 [332,] 1.0000000000 2.511240e-49 1.255620e-49 [333,] 1.0000000000 1.536516e-49 7.682581e-50 [334,] 1.0000000000 8.907507e-50 4.453753e-50 [335,] 1.0000000000 5.385109e-50 2.692555e-50 [336,] 1.0000000000 2.997672e-50 1.498836e-50 [337,] 1.0000000000 2.859463e-50 1.429731e-50 [338,] 1.0000000000 4.324619e-50 2.162310e-50 [339,] 1.0000000000 8.084544e-50 4.042272e-50 [340,] 1.0000000000 1.273156e-49 6.365779e-50 [341,] 1.0000000000 1.918702e-49 9.593510e-50 [342,] 1.0000000000 1.511763e-49 7.558817e-50 [343,] 1.0000000000 7.553471e-50 3.776736e-50 [344,] 1.0000000000 6.078556e-50 3.039278e-50 [345,] 1.0000000000 7.452663e-50 3.726332e-50 [346,] 1.0000000000 5.785702e-50 2.892851e-50 [347,] 1.0000000000 1.205607e-49 6.028033e-50 [348,] 1.0000000000 2.757959e-49 1.378979e-49 [349,] 1.0000000000 6.436311e-49 3.218155e-49 [350,] 1.0000000000 1.810813e-48 9.054066e-49 [351,] 1.0000000000 3.529473e-48 1.764736e-48 [352,] 1.0000000000 1.005324e-47 5.026619e-48 [353,] 1.0000000000 2.604745e-47 1.302373e-47 [354,] 1.0000000000 5.731836e-47 2.865918e-47 [355,] 1.0000000000 1.611712e-46 8.058558e-47 [356,] 1.0000000000 4.214047e-46 2.107024e-46 [357,] 1.0000000000 1.101382e-45 5.506912e-46 [358,] 1.0000000000 2.948109e-45 1.474055e-45 [359,] 1.0000000000 6.717737e-45 3.358869e-45 [360,] 1.0000000000 7.369079e-45 3.684540e-45 [361,] 1.0000000000 6.323156e-45 3.161578e-45 [362,] 1.0000000000 6.458530e-45 3.229265e-45 [363,] 1.0000000000 7.578894e-45 3.789447e-45 [364,] 1.0000000000 1.185967e-44 5.929834e-45 [365,] 1.0000000000 1.406122e-44 7.030608e-45 [366,] 1.0000000000 1.618911e-44 8.094554e-45 [367,] 1.0000000000 2.497651e-44 1.248826e-44 [368,] 1.0000000000 6.682431e-44 3.341215e-44 [369,] 1.0000000000 1.610826e-43 8.054128e-44 [370,] 1.0000000000 4.760651e-43 2.380326e-43 [371,] 1.0000000000 1.295860e-42 6.479298e-43 [372,] 1.0000000000 3.634071e-42 1.817036e-42 [373,] 1.0000000000 8.975619e-42 4.487810e-42 [374,] 1.0000000000 2.062887e-41 1.031444e-41 [375,] 1.0000000000 5.967443e-41 2.983722e-41 [376,] 1.0000000000 1.450724e-40 7.253621e-41 [377,] 1.0000000000 3.571322e-40 1.785661e-40 [378,] 1.0000000000 7.854555e-40 3.927278e-40 [379,] 1.0000000000 1.705054e-39 8.525272e-40 [380,] 1.0000000000 3.401216e-39 1.700608e-39 [381,] 1.0000000000 6.967066e-39 3.483533e-39 [382,] 1.0000000000 1.420324e-38 7.101620e-39 [383,] 1.0000000000 3.339879e-38 1.669939e-38 [384,] 1.0000000000 7.298977e-38 3.649488e-38 [385,] 1.0000000000 1.850032e-37 9.250161e-38 [386,] 1.0000000000 5.476579e-37 2.738289e-37 [387,] 1.0000000000 1.766455e-36 8.832274e-37 [388,] 1.0000000000 5.800661e-36 2.900330e-36 [389,] 1.0000000000 2.224177e-35 1.112088e-35 [390,] 1.0000000000 8.136546e-35 4.068273e-35 [391,] 1.0000000000 2.999178e-34 1.499589e-34 [392,] 1.0000000000 1.131037e-33 5.655185e-34 [393,] 1.0000000000 4.241987e-33 2.120994e-33 [394,] 1.0000000000 1.569523e-32 7.847617e-33 [395,] 1.0000000000 5.573554e-32 2.786777e-32 [396,] 1.0000000000 1.979327e-31 9.896635e-32 [397,] 1.0000000000 7.171423e-31 3.585712e-31 [398,] 1.0000000000 2.433935e-30 1.216968e-30 [399,] 1.0000000000 7.955272e-30 3.977636e-30 [400,] 1.0000000000 2.409614e-29 1.204807e-29 [401,] 1.0000000000 2.390336e-29 1.195168e-29 [402,] 1.0000000000 1.546067e-29 7.730335e-30 [403,] 1.0000000000 1.128152e-29 5.640760e-30 [404,] 1.0000000000 6.891378e-30 3.445689e-30 [405,] 1.0000000000 4.097730e-30 2.048865e-30 [406,] 1.0000000000 3.239951e-30 1.619975e-30 [407,] 1.0000000000 4.251843e-30 2.125921e-30 [408,] 1.0000000000 1.378128e-29 6.890642e-30 [409,] 1.0000000000 3.750121e-29 1.875061e-29 [410,] 1.0000000000 9.173311e-29 4.586655e-29 [411,] 1.0000000000 1.603323e-28 8.016615e-29 [412,] 1.0000000000 2.533209e-28 1.266604e-28 [413,] 1.0000000000 4.126375e-28 2.063188e-28 [414,] 1.0000000000 9.777191e-28 4.888595e-28 [415,] 1.0000000000 2.306934e-27 1.153467e-27 [416,] 1.0000000000 8.627358e-27 4.313679e-27 [417,] 1.0000000000 3.272919e-26 1.636460e-26 [418,] 1.0000000000 1.074630e-25 5.373148e-26 [419,] 1.0000000000 3.029765e-25 1.514883e-25 [420,] 1.0000000000 7.748488e-25 3.874244e-25 [421,] 1.0000000000 2.425362e-24 1.212681e-24 [422,] 1.0000000000 6.826373e-24 3.413186e-24 [423,] 1.0000000000 1.771510e-23 8.857550e-24 [424,] 1.0000000000 3.894320e-23 1.947160e-23 [425,] 1.0000000000 7.186012e-23 3.593006e-23 [426,] 1.0000000000 1.539434e-22 7.697172e-23 [427,] 1.0000000000 3.282689e-22 1.641344e-22 [428,] 1.0000000000 5.286182e-22 2.643091e-22 [429,] 1.0000000000 4.959243e-22 2.479622e-22 [430,] 1.0000000000 1.216395e-22 6.081975e-23 [431,] 1.0000000000 1.116152e-24 5.580761e-25 [432,] 1.0000000000 1.181772e-24 5.908861e-25 [433,] 1.0000000000 4.595804e-24 2.297902e-24 [434,] 1.0000000000 1.569703e-23 7.848515e-24 [435,] 1.0000000000 5.529235e-23 2.764617e-23 [436,] 1.0000000000 2.248515e-22 1.124258e-22 [437,] 1.0000000000 9.014920e-22 4.507460e-22 [438,] 1.0000000000 3.930941e-21 1.965471e-21 [439,] 1.0000000000 1.709589e-20 8.547945e-21 [440,] 1.0000000000 6.408618e-20 3.204309e-20 [441,] 1.0000000000 1.921927e-19 9.609634e-20 [442,] 1.0000000000 5.791555e-19 2.895777e-19 [443,] 1.0000000000 2.256463e-18 1.128232e-18 [444,] 1.0000000000 9.112242e-18 4.556121e-18 [445,] 1.0000000000 3.659606e-17 1.829803e-17 [446,] 1.0000000000 1.181028e-16 5.905138e-17 [447,] 1.0000000000 9.376474e-17 4.688237e-17 [448,] 1.0000000000 4.054078e-17 2.027039e-17 [449,] 1.0000000000 1.641145e-17 8.205725e-18 [450,] 1.0000000000 1.114286e-18 5.571430e-19 [451,] 1.0000000000 2.373841e-19 1.186920e-19 [452,] 1.0000000000 5.882430e-20 2.941215e-20 [453,] 1.0000000000 1.225813e-20 6.129063e-21 [454,] 1.0000000000 1.413604e-20 7.068020e-21 [455,] 1.0000000000 4.595829e-20 2.297915e-20 [456,] 1.0000000000 1.195438e-19 5.977189e-20 [457,] 1.0000000000 8.048818e-19 4.024409e-19 [458,] 1.0000000000 5.464866e-18 2.732433e-18 [459,] 1.0000000000 1.824283e-17 9.121417e-18 [460,] 1.0000000000 1.681314e-17 8.406569e-18 [461,] 1.0000000000 9.281967e-18 4.640983e-18 [462,] 1.0000000000 3.277106e-17 1.638553e-17 [463,] 1.0000000000 2.383790e-16 1.191895e-16 [464,] 1.0000000000 1.884364e-15 9.421818e-16 [465,] 1.0000000000 1.225616e-14 6.128079e-15 [466,] 1.0000000000 6.090087e-14 3.045043e-14 [467,] 1.0000000000 2.438464e-13 1.219232e-13 [468,] 1.0000000000 2.236338e-12 1.118169e-12 [469,] 1.0000000000 9.124981e-12 4.562490e-12 [470,] 1.0000000000 4.311212e-11 2.155606e-11 [471,] 1.0000000000 7.661610e-11 3.830805e-11 [472,] 0.9999999996 7.097622e-10 3.548811e-10 [473,] 0.9999999993 1.361633e-09 6.808164e-10 [474,] 0.9999999984 3.107004e-09 1.553502e-09 [475,] 0.9999999966 6.777357e-09 3.388679e-09 [476,] 0.9999999760 4.803148e-08 2.401574e-08 [477,] 0.9999999441 1.117465e-07 5.587323e-08 [478,] 0.9999992772 1.445661e-06 7.228304e-07 [479,] 0.9999875743 2.485138e-05 1.242569e-05 [480,] 0.9998203254 3.593491e-04 1.796746e-04 > postscript(file="/var/www/html/rcomp/tmp/1vz2k1292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2682n1292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3682n1292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4682n1292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5682n1292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(mysum$resid, main='Residual Normal Q-Q Plot') > qqline(mysum$resid) > grid() > dev.off() null device 1 > (myerror <- as.ts(mysum$resid)) Time Series: Start = 1 End = 491 Frequency = 1 1 2 3 4 5 6 -1.974472429 -3.074064414 -3.539560020 -2.758312228 -0.856414381 -0.800857700 7 8 9 10 11 12 -0.083335145 0.058495203 0.649159359 1.232738277 1.571460707 2.190733119 13 14 15 16 17 18 2.132385220 1.556003747 1.637854559 1.833529690 1.999358703 2.353767616 19 20 21 22 23 24 3.370214399 4.504941403 5.095757759 5.662927614 5.825658268 5.963583338 25 26 27 28 29 30 5.924404754 5.560016010 5.367815981 5.199538965 4.620362971 3.869801450 31 32 33 34 35 36 3.400554659 4.006459973 4.992312904 4.555994501 4.505616386 4.474158476 37 38 39 40 41 42 4.491007404 4.216797579 3.890274953 3.405453356 3.724504266 3.926029174 43 44 45 46 47 48 3.949608896 3.301388007 3.112665931 3.607765331 2.882458259 2.615381751 49 50 51 52 53 54 2.630370813 3.015096816 2.303903850 2.456650495 2.155839395 2.392465072 55 56 57 58 59 60 2.540728663 2.402074178 2.005110914 1.736624132 2.004336859 1.749187205 61 62 63 64 65 66 1.776517655 2.223208794 1.853802367 1.770281781 1.290300871 1.183792315 67 68 69 70 71 72 0.833289914 1.089112704 1.313054050 1.071352956 0.956201741 1.044262547 73 74 75 76 77 78 1.010049275 0.821224269 1.228925937 1.777199592 1.545627392 1.276809423 79 80 81 82 83 84 1.254169918 1.383261987 1.492087331 1.146485190 0.816912861 0.515691177 85 86 87 88 89 90 0.523657927 0.552127100 0.207276227 0.059089685 0.350536660 0.315806503 91 92 93 94 95 96 0.328432880 0.685702787 0.622303123 0.358270491 0.073161499 0.014684371 97 98 99 100 101 102 0.068220318 0.019900119 -0.268528121 0.053560458 -0.217383708 -0.549040394 103 104 105 106 107 108 -0.632305031 -0.561553798 -0.683512853 -0.659078648 -0.589945995 -0.854227050 109 110 111 112 113 114 -0.649655747 -0.744179747 -0.849402382 -1.072934131 -1.001310602 -0.825796372 115 116 117 118 119 120 -0.665843946 -0.371697282 -0.324112429 -0.546352860 -0.255800559 0.148319904 121 122 123 124 125 126 -0.032008815 -0.347869918 -0.115996909 -0.053858219 -0.252883626 -0.377392271 127 128 129 130 131 132 -0.361469486 -0.125279837 0.212359873 0.247033252 0.727872070 1.225553935 133 134 135 136 137 138 2.405347917 2.747547519 1.852663507 1.351482800 1.214827359 1.178481991 139 140 141 142 143 144 -0.034464314 0.141741076 0.541126678 0.371712072 0.011817011 -0.359414602 145 146 147 148 149 150 -0.140803499 0.331313300 0.636695928 0.430244212 0.194039632 0.505205475 151 152 153 154 155 156 0.381325487 0.463954168 1.023305858 0.484421036 1.449051534 1.904329337 157 158 159 160 161 162 0.926127235 0.906950993 1.738893088 1.183335934 1.354468283 1.242203033 163 164 165 166 167 168 1.365491152 0.821917160 0.794798773 0.912716977 1.277094580 1.228199296 169 170 171 172 173 174 1.619429546 1.683515441 1.778395405 1.559513544 1.382767044 1.199015789 175 176 177 178 179 180 1.376214537 1.106242038 0.888359285 0.894563325 1.067417348 1.089618712 181 182 183 184 185 186 0.991810396 0.396057586 0.326336118 0.200851845 0.217525154 0.088705772 187 188 189 190 191 192 0.105574027 0.203046291 0.392648555 0.237324007 0.274930341 0.107974517 193 194 195 196 197 198 0.410018348 0.212237851 0.142878355 0.036546397 0.027446858 -0.096897507 199 200 201 202 203 204 -0.440685208 -0.185109008 -0.311829607 -0.472813506 -0.745709798 -0.730656346 205 206 207 208 209 210 -0.837044536 -0.908260779 -1.077781553 -1.312454612 -1.476715465 -1.259634672 211 212 213 214 215 216 -1.054494600 -1.148939932 -1.272270052 -0.886200765 -0.831251547 -1.017761696 217 218 219 220 221 222 -1.011831572 -1.235108566 -1.291817057 -1.308704715 -1.362717566 -1.276088301 223 224 225 226 227 228 -1.217525945 -1.303103906 -1.241374651 -1.259438198 -0.987643286 -0.809235886 229 230 231 232 233 234 -0.989798651 -0.889300152 -0.580778643 0.005358535 0.227858946 0.306854861 235 236 237 238 239 240 0.047365797 0.155062270 -0.168843098 -0.006688865 0.066133540 1.359033188 241 242 243 244 245 246 1.283898524 1.258936591 0.613134341 0.585747728 0.430372830 0.488380440 247 248 249 250 251 252 0.865036623 0.586317600 0.335169139 0.080919731 -0.368430368 -0.990166748 253 254 255 256 257 258 -0.966601129 -0.484321766 -0.264536356 -0.753973978 -0.806886458 -0.733875342 259 260 261 262 263 264 -0.968038517 -1.297763254 -1.382669325 -1.481995410 -1.133688266 -0.911590774 265 266 267 268 269 270 -1.227349541 -1.050177374 -0.950138175 -1.000906290 -1.137736820 -1.332466672 271 272 273 274 275 276 -0.784585475 -0.424261737 -0.063442687 -0.126174111 -0.305232930 -0.386814616 277 278 279 280 281 282 -0.375651131 -0.674711301 -0.957081491 -0.926340455 -0.356766199 -0.319198381 283 284 285 286 287 288 -0.485914388 -0.606016622 -0.806959383 -0.932638128 -0.910700965 -0.943220260 289 290 291 292 293 294 -1.235668552 -1.338371592 -1.629368690 -1.477784843 -1.351689645 -1.296184899 295 296 297 298 299 300 -1.153095226 -1.087624192 -1.217019288 -1.260373459 -1.172655244 -1.149267945 301 302 303 304 305 306 -1.381097772 -1.605131673 -1.528795824 -1.873275502 -2.026103477 -2.113811995 307 308 309 310 311 312 -2.086765409 -2.101218029 -2.268636343 -2.249363135 -2.411726607 -1.871199077 313 314 315 316 317 318 -2.092034366 -2.347265150 -2.550364482 -2.645757419 -2.511380737 -2.669992470 319 320 321 322 323 324 -2.802924992 -2.945034019 -3.206514477 -3.170761507 -2.939118519 -3.172513631 325 326 327 328 329 330 -3.297929885 -3.096235396 -3.281703489 -3.230037420 -3.135503360 -2.863756187 331 332 333 334 335 336 -3.160000964 -3.287348680 -3.240873709 -3.453319231 -3.596641159 -3.455355051 337 338 339 340 341 342 -3.375141704 -3.482455375 -3.575277068 -3.480038465 -3.475100457 -3.294161483 343 344 345 346 347 348 -2.635370056 -2.145494350 -2.712545696 -2.574232417 -3.022943197 -3.137445418 349 350 351 352 353 354 -2.880439940 -2.731559518 -2.951093926 -2.305207603 -1.472042613 -1.072676718 355 356 357 358 359 360 -1.548666823 -1.137721289 -1.649438105 -2.040938725 -2.292960073 -2.084803269 361 362 363 364 365 366 -1.746373839 -1.982531293 -2.280411729 -2.146359156 -1.613909103 -1.491769680 367 368 369 370 371 372 -1.560680036 -1.523710417 -1.741688219 -1.637111588 -1.497345718 -1.689407589 373 374 375 376 377 378 -1.966128349 -1.689626770 -2.084019782 -1.872570424 -1.785733187 -1.432114075 379 380 381 382 383 384 -1.427306874 -1.531242298 -1.328063563 -1.274320939 -0.905443645 -0.744436718 385 386 387 388 389 390 -0.786097918 -0.818573483 -0.636162256 -0.960252415 -1.114927447 -1.278052043 391 392 393 394 395 396 -1.295115154 -1.411735494 -1.235699697 -0.897584818 -0.904269212 -0.896772153 397 398 399 400 401 402 -1.124488034 -1.053218043 -1.194947956 -1.142411923 -1.225477111 -1.298328105 403 404 405 406 407 408 -1.122212868 -0.993779503 -1.091805169 -1.232633302 -1.200052036 -1.322092168 409 410 411 412 413 414 -1.195462169 -0.985313783 -1.068605664 -0.743028069 -0.677678033 -0.793680952 415 416 417 418 419 420 -0.746508800 -0.805142671 -0.804874209 -0.841901947 -0.779023487 -0.738739764 421 422 423 424 425 426 -0.356713897 -0.336940220 -0.403784571 -0.515292718 -0.546812446 -0.372270528 427 428 429 430 431 432 -0.234750066 -0.249870851 -0.343896455 -0.308033578 -0.218273705 -0.158447668 433 434 435 436 437 438 -0.126697034 -0.179410519 -0.281313227 0.016307076 -0.069247826 0.254857500 439 440 441 442 443 444 0.152501676 0.198951738 0.428665471 0.834327107 0.843631267 0.724285935 445 446 447 448 449 450 0.501659874 0.477966840 0.549132336 0.878252241 0.871362794 0.907263999 451 452 453 454 455 456 0.612224452 0.275666519 0.426787405 0.611519713 0.564094832 0.838164773 457 458 459 460 461 462 0.987107075 1.213827899 1.486350198 1.540194239 1.464804854 1.950533686 463 464 465 466 467 468 2.043865359 1.592240601 1.568427553 1.447972511 1.569065643 1.810016229 469 470 471 472 473 474 2.112452560 2.221467970 2.262989001 2.227515073 2.164061733 2.214780332 475 476 477 478 479 480 2.006347567 1.944639871 1.445106726 1.215619391 1.486340115 1.443999721 481 482 483 484 485 486 1.589280178 1.739751310 1.571635541 1.579761180 1.557372565 1.720370667 487 488 489 490 491 1.467538051 1.302821189 1.291249621 0.968954202 1.073429670 > postscript(file="/var/www/html/rcomp/tmp/6hz181292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > dum <- cbind(lag(myerror,k=1),myerror) > dum Time Series: Start = 0 End = 491 Frequency = 1 lag(myerror, k = 1) myerror 0 -1.974472429 NA 1 -3.074064414 -1.974472429 2 -3.539560020 -3.074064414 3 -2.758312228 -3.539560020 4 -0.856414381 -2.758312228 5 -0.800857700 -0.856414381 6 -0.083335145 -0.800857700 7 0.058495203 -0.083335145 8 0.649159359 0.058495203 9 1.232738277 0.649159359 10 1.571460707 1.232738277 11 2.190733119 1.571460707 12 2.132385220 2.190733119 13 1.556003747 2.132385220 14 1.637854559 1.556003747 15 1.833529690 1.637854559 16 1.999358703 1.833529690 17 2.353767616 1.999358703 18 3.370214399 2.353767616 19 4.504941403 3.370214399 20 5.095757759 4.504941403 21 5.662927614 5.095757759 22 5.825658268 5.662927614 23 5.963583338 5.825658268 24 5.924404754 5.963583338 25 5.560016010 5.924404754 26 5.367815981 5.560016010 27 5.199538965 5.367815981 28 4.620362971 5.199538965 29 3.869801450 4.620362971 30 3.400554659 3.869801450 31 4.006459973 3.400554659 32 4.992312904 4.006459973 33 4.555994501 4.992312904 34 4.505616386 4.555994501 35 4.474158476 4.505616386 36 4.491007404 4.474158476 37 4.216797579 4.491007404 38 3.890274953 4.216797579 39 3.405453356 3.890274953 40 3.724504266 3.405453356 41 3.926029174 3.724504266 42 3.949608896 3.926029174 43 3.301388007 3.949608896 44 3.112665931 3.301388007 45 3.607765331 3.112665931 46 2.882458259 3.607765331 47 2.615381751 2.882458259 48 2.630370813 2.615381751 49 3.015096816 2.630370813 50 2.303903850 3.015096816 51 2.456650495 2.303903850 52 2.155839395 2.456650495 53 2.392465072 2.155839395 54 2.540728663 2.392465072 55 2.402074178 2.540728663 56 2.005110914 2.402074178 57 1.736624132 2.005110914 58 2.004336859 1.736624132 59 1.749187205 2.004336859 60 1.776517655 1.749187205 61 2.223208794 1.776517655 62 1.853802367 2.223208794 63 1.770281781 1.853802367 64 1.290300871 1.770281781 65 1.183792315 1.290300871 66 0.833289914 1.183792315 67 1.089112704 0.833289914 68 1.313054050 1.089112704 69 1.071352956 1.313054050 70 0.956201741 1.071352956 71 1.044262547 0.956201741 72 1.010049275 1.044262547 73 0.821224269 1.010049275 74 1.228925937 0.821224269 75 1.777199592 1.228925937 76 1.545627392 1.777199592 77 1.276809423 1.545627392 78 1.254169918 1.276809423 79 1.383261987 1.254169918 80 1.492087331 1.383261987 81 1.146485190 1.492087331 82 0.816912861 1.146485190 83 0.515691177 0.816912861 84 0.523657927 0.515691177 85 0.552127100 0.523657927 86 0.207276227 0.552127100 87 0.059089685 0.207276227 88 0.350536660 0.059089685 89 0.315806503 0.350536660 90 0.328432880 0.315806503 91 0.685702787 0.328432880 92 0.622303123 0.685702787 93 0.358270491 0.622303123 94 0.073161499 0.358270491 95 0.014684371 0.073161499 96 0.068220318 0.014684371 97 0.019900119 0.068220318 98 -0.268528121 0.019900119 99 0.053560458 -0.268528121 100 -0.217383708 0.053560458 101 -0.549040394 -0.217383708 102 -0.632305031 -0.549040394 103 -0.561553798 -0.632305031 104 -0.683512853 -0.561553798 105 -0.659078648 -0.683512853 106 -0.589945995 -0.659078648 107 -0.854227050 -0.589945995 108 -0.649655747 -0.854227050 109 -0.744179747 -0.649655747 110 -0.849402382 -0.744179747 111 -1.072934131 -0.849402382 112 -1.001310602 -1.072934131 113 -0.825796372 -1.001310602 114 -0.665843946 -0.825796372 115 -0.371697282 -0.665843946 116 -0.324112429 -0.371697282 117 -0.546352860 -0.324112429 118 -0.255800559 -0.546352860 119 0.148319904 -0.255800559 120 -0.032008815 0.148319904 121 -0.347869918 -0.032008815 122 -0.115996909 -0.347869918 123 -0.053858219 -0.115996909 124 -0.252883626 -0.053858219 125 -0.377392271 -0.252883626 126 -0.361469486 -0.377392271 127 -0.125279837 -0.361469486 128 0.212359873 -0.125279837 129 0.247033252 0.212359873 130 0.727872070 0.247033252 131 1.225553935 0.727872070 132 2.405347917 1.225553935 133 2.747547519 2.405347917 134 1.852663507 2.747547519 135 1.351482800 1.852663507 136 1.214827359 1.351482800 137 1.178481991 1.214827359 138 -0.034464314 1.178481991 139 0.141741076 -0.034464314 140 0.541126678 0.141741076 141 0.371712072 0.541126678 142 0.011817011 0.371712072 143 -0.359414602 0.011817011 144 -0.140803499 -0.359414602 145 0.331313300 -0.140803499 146 0.636695928 0.331313300 147 0.430244212 0.636695928 148 0.194039632 0.430244212 149 0.505205475 0.194039632 150 0.381325487 0.505205475 151 0.463954168 0.381325487 152 1.023305858 0.463954168 153 0.484421036 1.023305858 154 1.449051534 0.484421036 155 1.904329337 1.449051534 156 0.926127235 1.904329337 157 0.906950993 0.926127235 158 1.738893088 0.906950993 159 1.183335934 1.738893088 160 1.354468283 1.183335934 161 1.242203033 1.354468283 162 1.365491152 1.242203033 163 0.821917160 1.365491152 164 0.794798773 0.821917160 165 0.912716977 0.794798773 166 1.277094580 0.912716977 167 1.228199296 1.277094580 168 1.619429546 1.228199296 169 1.683515441 1.619429546 170 1.778395405 1.683515441 171 1.559513544 1.778395405 172 1.382767044 1.559513544 173 1.199015789 1.382767044 174 1.376214537 1.199015789 175 1.106242038 1.376214537 176 0.888359285 1.106242038 177 0.894563325 0.888359285 178 1.067417348 0.894563325 179 1.089618712 1.067417348 180 0.991810396 1.089618712 181 0.396057586 0.991810396 182 0.326336118 0.396057586 183 0.200851845 0.326336118 184 0.217525154 0.200851845 185 0.088705772 0.217525154 186 0.105574027 0.088705772 187 0.203046291 0.105574027 188 0.392648555 0.203046291 189 0.237324007 0.392648555 190 0.274930341 0.237324007 191 0.107974517 0.274930341 192 0.410018348 0.107974517 193 0.212237851 0.410018348 194 0.142878355 0.212237851 195 0.036546397 0.142878355 196 0.027446858 0.036546397 197 -0.096897507 0.027446858 198 -0.440685208 -0.096897507 199 -0.185109008 -0.440685208 200 -0.311829607 -0.185109008 201 -0.472813506 -0.311829607 202 -0.745709798 -0.472813506 203 -0.730656346 -0.745709798 204 -0.837044536 -0.730656346 205 -0.908260779 -0.837044536 206 -1.077781553 -0.908260779 207 -1.312454612 -1.077781553 208 -1.476715465 -1.312454612 209 -1.259634672 -1.476715465 210 -1.054494600 -1.259634672 211 -1.148939932 -1.054494600 212 -1.272270052 -1.148939932 213 -0.886200765 -1.272270052 214 -0.831251547 -0.886200765 215 -1.017761696 -0.831251547 216 -1.011831572 -1.017761696 217 -1.235108566 -1.011831572 218 -1.291817057 -1.235108566 219 -1.308704715 -1.291817057 220 -1.362717566 -1.308704715 221 -1.276088301 -1.362717566 222 -1.217525945 -1.276088301 223 -1.303103906 -1.217525945 224 -1.241374651 -1.303103906 225 -1.259438198 -1.241374651 226 -0.987643286 -1.259438198 227 -0.809235886 -0.987643286 228 -0.989798651 -0.809235886 229 -0.889300152 -0.989798651 230 -0.580778643 -0.889300152 231 0.005358535 -0.580778643 232 0.227858946 0.005358535 233 0.306854861 0.227858946 234 0.047365797 0.306854861 235 0.155062270 0.047365797 236 -0.168843098 0.155062270 237 -0.006688865 -0.168843098 238 0.066133540 -0.006688865 239 1.359033188 0.066133540 240 1.283898524 1.359033188 241 1.258936591 1.283898524 242 0.613134341 1.258936591 243 0.585747728 0.613134341 244 0.430372830 0.585747728 245 0.488380440 0.430372830 246 0.865036623 0.488380440 247 0.586317600 0.865036623 248 0.335169139 0.586317600 249 0.080919731 0.335169139 250 -0.368430368 0.080919731 251 -0.990166748 -0.368430368 252 -0.966601129 -0.990166748 253 -0.484321766 -0.966601129 254 -0.264536356 -0.484321766 255 -0.753973978 -0.264536356 256 -0.806886458 -0.753973978 257 -0.733875342 -0.806886458 258 -0.968038517 -0.733875342 259 -1.297763254 -0.968038517 260 -1.382669325 -1.297763254 261 -1.481995410 -1.382669325 262 -1.133688266 -1.481995410 263 -0.911590774 -1.133688266 264 -1.227349541 -0.911590774 265 -1.050177374 -1.227349541 266 -0.950138175 -1.050177374 267 -1.000906290 -0.950138175 268 -1.137736820 -1.000906290 269 -1.332466672 -1.137736820 270 -0.784585475 -1.332466672 271 -0.424261737 -0.784585475 272 -0.063442687 -0.424261737 273 -0.126174111 -0.063442687 274 -0.305232930 -0.126174111 275 -0.386814616 -0.305232930 276 -0.375651131 -0.386814616 277 -0.674711301 -0.375651131 278 -0.957081491 -0.674711301 279 -0.926340455 -0.957081491 280 -0.356766199 -0.926340455 281 -0.319198381 -0.356766199 282 -0.485914388 -0.319198381 283 -0.606016622 -0.485914388 284 -0.806959383 -0.606016622 285 -0.932638128 -0.806959383 286 -0.910700965 -0.932638128 287 -0.943220260 -0.910700965 288 -1.235668552 -0.943220260 289 -1.338371592 -1.235668552 290 -1.629368690 -1.338371592 291 -1.477784843 -1.629368690 292 -1.351689645 -1.477784843 293 -1.296184899 -1.351689645 294 -1.153095226 -1.296184899 295 -1.087624192 -1.153095226 296 -1.217019288 -1.087624192 297 -1.260373459 -1.217019288 298 -1.172655244 -1.260373459 299 -1.149267945 -1.172655244 300 -1.381097772 -1.149267945 301 -1.605131673 -1.381097772 302 -1.528795824 -1.605131673 303 -1.873275502 -1.528795824 304 -2.026103477 -1.873275502 305 -2.113811995 -2.026103477 306 -2.086765409 -2.113811995 307 -2.101218029 -2.086765409 308 -2.268636343 -2.101218029 309 -2.249363135 -2.268636343 310 -2.411726607 -2.249363135 311 -1.871199077 -2.411726607 312 -2.092034366 -1.871199077 313 -2.347265150 -2.092034366 314 -2.550364482 -2.347265150 315 -2.645757419 -2.550364482 316 -2.511380737 -2.645757419 317 -2.669992470 -2.511380737 318 -2.802924992 -2.669992470 319 -2.945034019 -2.802924992 320 -3.206514477 -2.945034019 321 -3.170761507 -3.206514477 322 -2.939118519 -3.170761507 323 -3.172513631 -2.939118519 324 -3.297929885 -3.172513631 325 -3.096235396 -3.297929885 326 -3.281703489 -3.096235396 327 -3.230037420 -3.281703489 328 -3.135503360 -3.230037420 329 -2.863756187 -3.135503360 330 -3.160000964 -2.863756187 331 -3.287348680 -3.160000964 332 -3.240873709 -3.287348680 333 -3.453319231 -3.240873709 334 -3.596641159 -3.453319231 335 -3.455355051 -3.596641159 336 -3.375141704 -3.455355051 337 -3.482455375 -3.375141704 338 -3.575277068 -3.482455375 339 -3.480038465 -3.575277068 340 -3.475100457 -3.480038465 341 -3.294161483 -3.475100457 342 -2.635370056 -3.294161483 343 -2.145494350 -2.635370056 344 -2.712545696 -2.145494350 345 -2.574232417 -2.712545696 346 -3.022943197 -2.574232417 347 -3.137445418 -3.022943197 348 -2.880439940 -3.137445418 349 -2.731559518 -2.880439940 350 -2.951093926 -2.731559518 351 -2.305207603 -2.951093926 352 -1.472042613 -2.305207603 353 -1.072676718 -1.472042613 354 -1.548666823 -1.072676718 355 -1.137721289 -1.548666823 356 -1.649438105 -1.137721289 357 -2.040938725 -1.649438105 358 -2.292960073 -2.040938725 359 -2.084803269 -2.292960073 360 -1.746373839 -2.084803269 361 -1.982531293 -1.746373839 362 -2.280411729 -1.982531293 363 -2.146359156 -2.280411729 364 -1.613909103 -2.146359156 365 -1.491769680 -1.613909103 366 -1.560680036 -1.491769680 367 -1.523710417 -1.560680036 368 -1.741688219 -1.523710417 369 -1.637111588 -1.741688219 370 -1.497345718 -1.637111588 371 -1.689407589 -1.497345718 372 -1.966128349 -1.689407589 373 -1.689626770 -1.966128349 374 -2.084019782 -1.689626770 375 -1.872570424 -2.084019782 376 -1.785733187 -1.872570424 377 -1.432114075 -1.785733187 378 -1.427306874 -1.432114075 379 -1.531242298 -1.427306874 380 -1.328063563 -1.531242298 381 -1.274320939 -1.328063563 382 -0.905443645 -1.274320939 383 -0.744436718 -0.905443645 384 -0.786097918 -0.744436718 385 -0.818573483 -0.786097918 386 -0.636162256 -0.818573483 387 -0.960252415 -0.636162256 388 -1.114927447 -0.960252415 389 -1.278052043 -1.114927447 390 -1.295115154 -1.278052043 391 -1.411735494 -1.295115154 392 -1.235699697 -1.411735494 393 -0.897584818 -1.235699697 394 -0.904269212 -0.897584818 395 -0.896772153 -0.904269212 396 -1.124488034 -0.896772153 397 -1.053218043 -1.124488034 398 -1.194947956 -1.053218043 399 -1.142411923 -1.194947956 400 -1.225477111 -1.142411923 401 -1.298328105 -1.225477111 402 -1.122212868 -1.298328105 403 -0.993779503 -1.122212868 404 -1.091805169 -0.993779503 405 -1.232633302 -1.091805169 406 -1.200052036 -1.232633302 407 -1.322092168 -1.200052036 408 -1.195462169 -1.322092168 409 -0.985313783 -1.195462169 410 -1.068605664 -0.985313783 411 -0.743028069 -1.068605664 412 -0.677678033 -0.743028069 413 -0.793680952 -0.677678033 414 -0.746508800 -0.793680952 415 -0.805142671 -0.746508800 416 -0.804874209 -0.805142671 417 -0.841901947 -0.804874209 418 -0.779023487 -0.841901947 419 -0.738739764 -0.779023487 420 -0.356713897 -0.738739764 421 -0.336940220 -0.356713897 422 -0.403784571 -0.336940220 423 -0.515292718 -0.403784571 424 -0.546812446 -0.515292718 425 -0.372270528 -0.546812446 426 -0.234750066 -0.372270528 427 -0.249870851 -0.234750066 428 -0.343896455 -0.249870851 429 -0.308033578 -0.343896455 430 -0.218273705 -0.308033578 431 -0.158447668 -0.218273705 432 -0.126697034 -0.158447668 433 -0.179410519 -0.126697034 434 -0.281313227 -0.179410519 435 0.016307076 -0.281313227 436 -0.069247826 0.016307076 437 0.254857500 -0.069247826 438 0.152501676 0.254857500 439 0.198951738 0.152501676 440 0.428665471 0.198951738 441 0.834327107 0.428665471 442 0.843631267 0.834327107 443 0.724285935 0.843631267 444 0.501659874 0.724285935 445 0.477966840 0.501659874 446 0.549132336 0.477966840 447 0.878252241 0.549132336 448 0.871362794 0.878252241 449 0.907263999 0.871362794 450 0.612224452 0.907263999 451 0.275666519 0.612224452 452 0.426787405 0.275666519 453 0.611519713 0.426787405 454 0.564094832 0.611519713 455 0.838164773 0.564094832 456 0.987107075 0.838164773 457 1.213827899 0.987107075 458 1.486350198 1.213827899 459 1.540194239 1.486350198 460 1.464804854 1.540194239 461 1.950533686 1.464804854 462 2.043865359 1.950533686 463 1.592240601 2.043865359 464 1.568427553 1.592240601 465 1.447972511 1.568427553 466 1.569065643 1.447972511 467 1.810016229 1.569065643 468 2.112452560 1.810016229 469 2.221467970 2.112452560 470 2.262989001 2.221467970 471 2.227515073 2.262989001 472 2.164061733 2.227515073 473 2.214780332 2.164061733 474 2.006347567 2.214780332 475 1.944639871 2.006347567 476 1.445106726 1.944639871 477 1.215619391 1.445106726 478 1.486340115 1.215619391 479 1.443999721 1.486340115 480 1.589280178 1.443999721 481 1.739751310 1.589280178 482 1.571635541 1.739751310 483 1.579761180 1.571635541 484 1.557372565 1.579761180 485 1.720370667 1.557372565 486 1.467538051 1.720370667 487 1.302821189 1.467538051 488 1.291249621 1.302821189 489 0.968954202 1.291249621 490 1.073429670 0.968954202 491 NA 1.073429670 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -3.074064414 -1.974472429 [2,] -3.539560020 -3.074064414 [3,] -2.758312228 -3.539560020 [4,] -0.856414381 -2.758312228 [5,] -0.800857700 -0.856414381 [6,] -0.083335145 -0.800857700 [7,] 0.058495203 -0.083335145 [8,] 0.649159359 0.058495203 [9,] 1.232738277 0.649159359 [10,] 1.571460707 1.232738277 [11,] 2.190733119 1.571460707 [12,] 2.132385220 2.190733119 [13,] 1.556003747 2.132385220 [14,] 1.637854559 1.556003747 [15,] 1.833529690 1.637854559 [16,] 1.999358703 1.833529690 [17,] 2.353767616 1.999358703 [18,] 3.370214399 2.353767616 [19,] 4.504941403 3.370214399 [20,] 5.095757759 4.504941403 [21,] 5.662927614 5.095757759 [22,] 5.825658268 5.662927614 [23,] 5.963583338 5.825658268 [24,] 5.924404754 5.963583338 [25,] 5.560016010 5.924404754 [26,] 5.367815981 5.560016010 [27,] 5.199538965 5.367815981 [28,] 4.620362971 5.199538965 [29,] 3.869801450 4.620362971 [30,] 3.400554659 3.869801450 [31,] 4.006459973 3.400554659 [32,] 4.992312904 4.006459973 [33,] 4.555994501 4.992312904 [34,] 4.505616386 4.555994501 [35,] 4.474158476 4.505616386 [36,] 4.491007404 4.474158476 [37,] 4.216797579 4.491007404 [38,] 3.890274953 4.216797579 [39,] 3.405453356 3.890274953 [40,] 3.724504266 3.405453356 [41,] 3.926029174 3.724504266 [42,] 3.949608896 3.926029174 [43,] 3.301388007 3.949608896 [44,] 3.112665931 3.301388007 [45,] 3.607765331 3.112665931 [46,] 2.882458259 3.607765331 [47,] 2.615381751 2.882458259 [48,] 2.630370813 2.615381751 [49,] 3.015096816 2.630370813 [50,] 2.303903850 3.015096816 [51,] 2.456650495 2.303903850 [52,] 2.155839395 2.456650495 [53,] 2.392465072 2.155839395 [54,] 2.540728663 2.392465072 [55,] 2.402074178 2.540728663 [56,] 2.005110914 2.402074178 [57,] 1.736624132 2.005110914 [58,] 2.004336859 1.736624132 [59,] 1.749187205 2.004336859 [60,] 1.776517655 1.749187205 [61,] 2.223208794 1.776517655 [62,] 1.853802367 2.223208794 [63,] 1.770281781 1.853802367 [64,] 1.290300871 1.770281781 [65,] 1.183792315 1.290300871 [66,] 0.833289914 1.183792315 [67,] 1.089112704 0.833289914 [68,] 1.313054050 1.089112704 [69,] 1.071352956 1.313054050 [70,] 0.956201741 1.071352956 [71,] 1.044262547 0.956201741 [72,] 1.010049275 1.044262547 [73,] 0.821224269 1.010049275 [74,] 1.228925937 0.821224269 [75,] 1.777199592 1.228925937 [76,] 1.545627392 1.777199592 [77,] 1.276809423 1.545627392 [78,] 1.254169918 1.276809423 [79,] 1.383261987 1.254169918 [80,] 1.492087331 1.383261987 [81,] 1.146485190 1.492087331 [82,] 0.816912861 1.146485190 [83,] 0.515691177 0.816912861 [84,] 0.523657927 0.515691177 [85,] 0.552127100 0.523657927 [86,] 0.207276227 0.552127100 [87,] 0.059089685 0.207276227 [88,] 0.350536660 0.059089685 [89,] 0.315806503 0.350536660 [90,] 0.328432880 0.315806503 [91,] 0.685702787 0.328432880 [92,] 0.622303123 0.685702787 [93,] 0.358270491 0.622303123 [94,] 0.073161499 0.358270491 [95,] 0.014684371 0.073161499 [96,] 0.068220318 0.014684371 [97,] 0.019900119 0.068220318 [98,] -0.268528121 0.019900119 [99,] 0.053560458 -0.268528121 [100,] -0.217383708 0.053560458 [101,] -0.549040394 -0.217383708 [102,] -0.632305031 -0.549040394 [103,] -0.561553798 -0.632305031 [104,] -0.683512853 -0.561553798 [105,] -0.659078648 -0.683512853 [106,] -0.589945995 -0.659078648 [107,] -0.854227050 -0.589945995 [108,] -0.649655747 -0.854227050 [109,] -0.744179747 -0.649655747 [110,] -0.849402382 -0.744179747 [111,] -1.072934131 -0.849402382 [112,] -1.001310602 -1.072934131 [113,] -0.825796372 -1.001310602 [114,] -0.665843946 -0.825796372 [115,] -0.371697282 -0.665843946 [116,] -0.324112429 -0.371697282 [117,] -0.546352860 -0.324112429 [118,] -0.255800559 -0.546352860 [119,] 0.148319904 -0.255800559 [120,] -0.032008815 0.148319904 [121,] -0.347869918 -0.032008815 [122,] -0.115996909 -0.347869918 [123,] -0.053858219 -0.115996909 [124,] -0.252883626 -0.053858219 [125,] -0.377392271 -0.252883626 [126,] -0.361469486 -0.377392271 [127,] -0.125279837 -0.361469486 [128,] 0.212359873 -0.125279837 [129,] 0.247033252 0.212359873 [130,] 0.727872070 0.247033252 [131,] 1.225553935 0.727872070 [132,] 2.405347917 1.225553935 [133,] 2.747547519 2.405347917 [134,] 1.852663507 2.747547519 [135,] 1.351482800 1.852663507 [136,] 1.214827359 1.351482800 [137,] 1.178481991 1.214827359 [138,] -0.034464314 1.178481991 [139,] 0.141741076 -0.034464314 [140,] 0.541126678 0.141741076 [141,] 0.371712072 0.541126678 [142,] 0.011817011 0.371712072 [143,] -0.359414602 0.011817011 [144,] -0.140803499 -0.359414602 [145,] 0.331313300 -0.140803499 [146,] 0.636695928 0.331313300 [147,] 0.430244212 0.636695928 [148,] 0.194039632 0.430244212 [149,] 0.505205475 0.194039632 [150,] 0.381325487 0.505205475 [151,] 0.463954168 0.381325487 [152,] 1.023305858 0.463954168 [153,] 0.484421036 1.023305858 [154,] 1.449051534 0.484421036 [155,] 1.904329337 1.449051534 [156,] 0.926127235 1.904329337 [157,] 0.906950993 0.926127235 [158,] 1.738893088 0.906950993 [159,] 1.183335934 1.738893088 [160,] 1.354468283 1.183335934 [161,] 1.242203033 1.354468283 [162,] 1.365491152 1.242203033 [163,] 0.821917160 1.365491152 [164,] 0.794798773 0.821917160 [165,] 0.912716977 0.794798773 [166,] 1.277094580 0.912716977 [167,] 1.228199296 1.277094580 [168,] 1.619429546 1.228199296 [169,] 1.683515441 1.619429546 [170,] 1.778395405 1.683515441 [171,] 1.559513544 1.778395405 [172,] 1.382767044 1.559513544 [173,] 1.199015789 1.382767044 [174,] 1.376214537 1.199015789 [175,] 1.106242038 1.376214537 [176,] 0.888359285 1.106242038 [177,] 0.894563325 0.888359285 [178,] 1.067417348 0.894563325 [179,] 1.089618712 1.067417348 [180,] 0.991810396 1.089618712 [181,] 0.396057586 0.991810396 [182,] 0.326336118 0.396057586 [183,] 0.200851845 0.326336118 [184,] 0.217525154 0.200851845 [185,] 0.088705772 0.217525154 [186,] 0.105574027 0.088705772 [187,] 0.203046291 0.105574027 [188,] 0.392648555 0.203046291 [189,] 0.237324007 0.392648555 [190,] 0.274930341 0.237324007 [191,] 0.107974517 0.274930341 [192,] 0.410018348 0.107974517 [193,] 0.212237851 0.410018348 [194,] 0.142878355 0.212237851 [195,] 0.036546397 0.142878355 [196,] 0.027446858 0.036546397 [197,] -0.096897507 0.027446858 [198,] -0.440685208 -0.096897507 [199,] -0.185109008 -0.440685208 [200,] -0.311829607 -0.185109008 [201,] -0.472813506 -0.311829607 [202,] -0.745709798 -0.472813506 [203,] -0.730656346 -0.745709798 [204,] -0.837044536 -0.730656346 [205,] -0.908260779 -0.837044536 [206,] -1.077781553 -0.908260779 [207,] -1.312454612 -1.077781553 [208,] -1.476715465 -1.312454612 [209,] -1.259634672 -1.476715465 [210,] -1.054494600 -1.259634672 [211,] -1.148939932 -1.054494600 [212,] -1.272270052 -1.148939932 [213,] -0.886200765 -1.272270052 [214,] -0.831251547 -0.886200765 [215,] -1.017761696 -0.831251547 [216,] -1.011831572 -1.017761696 [217,] -1.235108566 -1.011831572 [218,] -1.291817057 -1.235108566 [219,] -1.308704715 -1.291817057 [220,] -1.362717566 -1.308704715 [221,] -1.276088301 -1.362717566 [222,] -1.217525945 -1.276088301 [223,] -1.303103906 -1.217525945 [224,] -1.241374651 -1.303103906 [225,] -1.259438198 -1.241374651 [226,] -0.987643286 -1.259438198 [227,] -0.809235886 -0.987643286 [228,] -0.989798651 -0.809235886 [229,] -0.889300152 -0.989798651 [230,] -0.580778643 -0.889300152 [231,] 0.005358535 -0.580778643 [232,] 0.227858946 0.005358535 [233,] 0.306854861 0.227858946 [234,] 0.047365797 0.306854861 [235,] 0.155062270 0.047365797 [236,] -0.168843098 0.155062270 [237,] -0.006688865 -0.168843098 [238,] 0.066133540 -0.006688865 [239,] 1.359033188 0.066133540 [240,] 1.283898524 1.359033188 [241,] 1.258936591 1.283898524 [242,] 0.613134341 1.258936591 [243,] 0.585747728 0.613134341 [244,] 0.430372830 0.585747728 [245,] 0.488380440 0.430372830 [246,] 0.865036623 0.488380440 [247,] 0.586317600 0.865036623 [248,] 0.335169139 0.586317600 [249,] 0.080919731 0.335169139 [250,] -0.368430368 0.080919731 [251,] -0.990166748 -0.368430368 [252,] -0.966601129 -0.990166748 [253,] -0.484321766 -0.966601129 [254,] -0.264536356 -0.484321766 [255,] -0.753973978 -0.264536356 [256,] -0.806886458 -0.753973978 [257,] -0.733875342 -0.806886458 [258,] -0.968038517 -0.733875342 [259,] -1.297763254 -0.968038517 [260,] -1.382669325 -1.297763254 [261,] -1.481995410 -1.382669325 [262,] -1.133688266 -1.481995410 [263,] -0.911590774 -1.133688266 [264,] -1.227349541 -0.911590774 [265,] -1.050177374 -1.227349541 [266,] -0.950138175 -1.050177374 [267,] -1.000906290 -0.950138175 [268,] -1.137736820 -1.000906290 [269,] -1.332466672 -1.137736820 [270,] -0.784585475 -1.332466672 [271,] -0.424261737 -0.784585475 [272,] -0.063442687 -0.424261737 [273,] -0.126174111 -0.063442687 [274,] -0.305232930 -0.126174111 [275,] -0.386814616 -0.305232930 [276,] -0.375651131 -0.386814616 [277,] -0.674711301 -0.375651131 [278,] -0.957081491 -0.674711301 [279,] -0.926340455 -0.957081491 [280,] -0.356766199 -0.926340455 [281,] -0.319198381 -0.356766199 [282,] -0.485914388 -0.319198381 [283,] -0.606016622 -0.485914388 [284,] -0.806959383 -0.606016622 [285,] -0.932638128 -0.806959383 [286,] -0.910700965 -0.932638128 [287,] -0.943220260 -0.910700965 [288,] -1.235668552 -0.943220260 [289,] -1.338371592 -1.235668552 [290,] -1.629368690 -1.338371592 [291,] -1.477784843 -1.629368690 [292,] -1.351689645 -1.477784843 [293,] -1.296184899 -1.351689645 [294,] -1.153095226 -1.296184899 [295,] -1.087624192 -1.153095226 [296,] -1.217019288 -1.087624192 [297,] -1.260373459 -1.217019288 [298,] -1.172655244 -1.260373459 [299,] -1.149267945 -1.172655244 [300,] -1.381097772 -1.149267945 [301,] -1.605131673 -1.381097772 [302,] -1.528795824 -1.605131673 [303,] -1.873275502 -1.528795824 [304,] -2.026103477 -1.873275502 [305,] -2.113811995 -2.026103477 [306,] -2.086765409 -2.113811995 [307,] -2.101218029 -2.086765409 [308,] -2.268636343 -2.101218029 [309,] -2.249363135 -2.268636343 [310,] -2.411726607 -2.249363135 [311,] -1.871199077 -2.411726607 [312,] -2.092034366 -1.871199077 [313,] -2.347265150 -2.092034366 [314,] -2.550364482 -2.347265150 [315,] -2.645757419 -2.550364482 [316,] -2.511380737 -2.645757419 [317,] -2.669992470 -2.511380737 [318,] -2.802924992 -2.669992470 [319,] -2.945034019 -2.802924992 [320,] -3.206514477 -2.945034019 [321,] -3.170761507 -3.206514477 [322,] -2.939118519 -3.170761507 [323,] -3.172513631 -2.939118519 [324,] -3.297929885 -3.172513631 [325,] -3.096235396 -3.297929885 [326,] -3.281703489 -3.096235396 [327,] -3.230037420 -3.281703489 [328,] -3.135503360 -3.230037420 [329,] -2.863756187 -3.135503360 [330,] -3.160000964 -2.863756187 [331,] -3.287348680 -3.160000964 [332,] -3.240873709 -3.287348680 [333,] -3.453319231 -3.240873709 [334,] -3.596641159 -3.453319231 [335,] -3.455355051 -3.596641159 [336,] -3.375141704 -3.455355051 [337,] -3.482455375 -3.375141704 [338,] -3.575277068 -3.482455375 [339,] -3.480038465 -3.575277068 [340,] -3.475100457 -3.480038465 [341,] -3.294161483 -3.475100457 [342,] -2.635370056 -3.294161483 [343,] -2.145494350 -2.635370056 [344,] -2.712545696 -2.145494350 [345,] -2.574232417 -2.712545696 [346,] -3.022943197 -2.574232417 [347,] -3.137445418 -3.022943197 [348,] -2.880439940 -3.137445418 [349,] -2.731559518 -2.880439940 [350,] -2.951093926 -2.731559518 [351,] -2.305207603 -2.951093926 [352,] -1.472042613 -2.305207603 [353,] -1.072676718 -1.472042613 [354,] -1.548666823 -1.072676718 [355,] -1.137721289 -1.548666823 [356,] -1.649438105 -1.137721289 [357,] -2.040938725 -1.649438105 [358,] -2.292960073 -2.040938725 [359,] -2.084803269 -2.292960073 [360,] -1.746373839 -2.084803269 [361,] -1.982531293 -1.746373839 [362,] -2.280411729 -1.982531293 [363,] -2.146359156 -2.280411729 [364,] -1.613909103 -2.146359156 [365,] -1.491769680 -1.613909103 [366,] -1.560680036 -1.491769680 [367,] -1.523710417 -1.560680036 [368,] -1.741688219 -1.523710417 [369,] -1.637111588 -1.741688219 [370,] -1.497345718 -1.637111588 [371,] -1.689407589 -1.497345718 [372,] -1.966128349 -1.689407589 [373,] -1.689626770 -1.966128349 [374,] -2.084019782 -1.689626770 [375,] -1.872570424 -2.084019782 [376,] -1.785733187 -1.872570424 [377,] -1.432114075 -1.785733187 [378,] -1.427306874 -1.432114075 [379,] -1.531242298 -1.427306874 [380,] -1.328063563 -1.531242298 [381,] -1.274320939 -1.328063563 [382,] -0.905443645 -1.274320939 [383,] -0.744436718 -0.905443645 [384,] -0.786097918 -0.744436718 [385,] -0.818573483 -0.786097918 [386,] -0.636162256 -0.818573483 [387,] -0.960252415 -0.636162256 [388,] -1.114927447 -0.960252415 [389,] -1.278052043 -1.114927447 [390,] -1.295115154 -1.278052043 [391,] -1.411735494 -1.295115154 [392,] -1.235699697 -1.411735494 [393,] -0.897584818 -1.235699697 [394,] -0.904269212 -0.897584818 [395,] -0.896772153 -0.904269212 [396,] -1.124488034 -0.896772153 [397,] -1.053218043 -1.124488034 [398,] -1.194947956 -1.053218043 [399,] -1.142411923 -1.194947956 [400,] -1.225477111 -1.142411923 [401,] -1.298328105 -1.225477111 [402,] -1.122212868 -1.298328105 [403,] -0.993779503 -1.122212868 [404,] -1.091805169 -0.993779503 [405,] -1.232633302 -1.091805169 [406,] -1.200052036 -1.232633302 [407,] -1.322092168 -1.200052036 [408,] -1.195462169 -1.322092168 [409,] -0.985313783 -1.195462169 [410,] -1.068605664 -0.985313783 [411,] -0.743028069 -1.068605664 [412,] -0.677678033 -0.743028069 [413,] -0.793680952 -0.677678033 [414,] -0.746508800 -0.793680952 [415,] -0.805142671 -0.746508800 [416,] -0.804874209 -0.805142671 [417,] -0.841901947 -0.804874209 [418,] -0.779023487 -0.841901947 [419,] -0.738739764 -0.779023487 [420,] -0.356713897 -0.738739764 [421,] -0.336940220 -0.356713897 [422,] -0.403784571 -0.336940220 [423,] -0.515292718 -0.403784571 [424,] -0.546812446 -0.515292718 [425,] -0.372270528 -0.546812446 [426,] -0.234750066 -0.372270528 [427,] -0.249870851 -0.234750066 [428,] -0.343896455 -0.249870851 [429,] -0.308033578 -0.343896455 [430,] -0.218273705 -0.308033578 [431,] -0.158447668 -0.218273705 [432,] -0.126697034 -0.158447668 [433,] -0.179410519 -0.126697034 [434,] -0.281313227 -0.179410519 [435,] 0.016307076 -0.281313227 [436,] -0.069247826 0.016307076 [437,] 0.254857500 -0.069247826 [438,] 0.152501676 0.254857500 [439,] 0.198951738 0.152501676 [440,] 0.428665471 0.198951738 [441,] 0.834327107 0.428665471 [442,] 0.843631267 0.834327107 [443,] 0.724285935 0.843631267 [444,] 0.501659874 0.724285935 [445,] 0.477966840 0.501659874 [446,] 0.549132336 0.477966840 [447,] 0.878252241 0.549132336 [448,] 0.871362794 0.878252241 [449,] 0.907263999 0.871362794 [450,] 0.612224452 0.907263999 [451,] 0.275666519 0.612224452 [452,] 0.426787405 0.275666519 [453,] 0.611519713 0.426787405 [454,] 0.564094832 0.611519713 [455,] 0.838164773 0.564094832 [456,] 0.987107075 0.838164773 [457,] 1.213827899 0.987107075 [458,] 1.486350198 1.213827899 [459,] 1.540194239 1.486350198 [460,] 1.464804854 1.540194239 [461,] 1.950533686 1.464804854 [462,] 2.043865359 1.950533686 [463,] 1.592240601 2.043865359 [464,] 1.568427553 1.592240601 [465,] 1.447972511 1.568427553 [466,] 1.569065643 1.447972511 [467,] 1.810016229 1.569065643 [468,] 2.112452560 1.810016229 [469,] 2.221467970 2.112452560 [470,] 2.262989001 2.221467970 [471,] 2.227515073 2.262989001 [472,] 2.164061733 2.227515073 [473,] 2.214780332 2.164061733 [474,] 2.006347567 2.214780332 [475,] 1.944639871 2.006347567 [476,] 1.445106726 1.944639871 [477,] 1.215619391 1.445106726 [478,] 1.486340115 1.215619391 [479,] 1.443999721 1.486340115 [480,] 1.589280178 1.443999721 [481,] 1.739751310 1.589280178 [482,] 1.571635541 1.739751310 [483,] 1.579761180 1.571635541 [484,] 1.557372565 1.579761180 [485,] 1.720370667 1.557372565 [486,] 1.467538051 1.720370667 [487,] 1.302821189 1.467538051 [488,] 1.291249621 1.302821189 [489,] 0.968954202 1.291249621 [490,] 1.073429670 0.968954202 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -3.074064414 -1.974472429 2 -3.539560020 -3.074064414 3 -2.758312228 -3.539560020 4 -0.856414381 -2.758312228 5 -0.800857700 -0.856414381 6 -0.083335145 -0.800857700 7 0.058495203 -0.083335145 8 0.649159359 0.058495203 9 1.232738277 0.649159359 10 1.571460707 1.232738277 11 2.190733119 1.571460707 12 2.132385220 2.190733119 13 1.556003747 2.132385220 14 1.637854559 1.556003747 15 1.833529690 1.637854559 16 1.999358703 1.833529690 17 2.353767616 1.999358703 18 3.370214399 2.353767616 19 4.504941403 3.370214399 20 5.095757759 4.504941403 21 5.662927614 5.095757759 22 5.825658268 5.662927614 23 5.963583338 5.825658268 24 5.924404754 5.963583338 25 5.560016010 5.924404754 26 5.367815981 5.560016010 27 5.199538965 5.367815981 28 4.620362971 5.199538965 29 3.869801450 4.620362971 30 3.400554659 3.869801450 31 4.006459973 3.400554659 32 4.992312904 4.006459973 33 4.555994501 4.992312904 34 4.505616386 4.555994501 35 4.474158476 4.505616386 36 4.491007404 4.474158476 37 4.216797579 4.491007404 38 3.890274953 4.216797579 39 3.405453356 3.890274953 40 3.724504266 3.405453356 41 3.926029174 3.724504266 42 3.949608896 3.926029174 43 3.301388007 3.949608896 44 3.112665931 3.301388007 45 3.607765331 3.112665931 46 2.882458259 3.607765331 47 2.615381751 2.882458259 48 2.630370813 2.615381751 49 3.015096816 2.630370813 50 2.303903850 3.015096816 51 2.456650495 2.303903850 52 2.155839395 2.456650495 53 2.392465072 2.155839395 54 2.540728663 2.392465072 55 2.402074178 2.540728663 56 2.005110914 2.402074178 57 1.736624132 2.005110914 58 2.004336859 1.736624132 59 1.749187205 2.004336859 60 1.776517655 1.749187205 61 2.223208794 1.776517655 62 1.853802367 2.223208794 63 1.770281781 1.853802367 64 1.290300871 1.770281781 65 1.183792315 1.290300871 66 0.833289914 1.183792315 67 1.089112704 0.833289914 68 1.313054050 1.089112704 69 1.071352956 1.313054050 70 0.956201741 1.071352956 71 1.044262547 0.956201741 72 1.010049275 1.044262547 73 0.821224269 1.010049275 74 1.228925937 0.821224269 75 1.777199592 1.228925937 76 1.545627392 1.777199592 77 1.276809423 1.545627392 78 1.254169918 1.276809423 79 1.383261987 1.254169918 80 1.492087331 1.383261987 81 1.146485190 1.492087331 82 0.816912861 1.146485190 83 0.515691177 0.816912861 84 0.523657927 0.515691177 85 0.552127100 0.523657927 86 0.207276227 0.552127100 87 0.059089685 0.207276227 88 0.350536660 0.059089685 89 0.315806503 0.350536660 90 0.328432880 0.315806503 91 0.685702787 0.328432880 92 0.622303123 0.685702787 93 0.358270491 0.622303123 94 0.073161499 0.358270491 95 0.014684371 0.073161499 96 0.068220318 0.014684371 97 0.019900119 0.068220318 98 -0.268528121 0.019900119 99 0.053560458 -0.268528121 100 -0.217383708 0.053560458 101 -0.549040394 -0.217383708 102 -0.632305031 -0.549040394 103 -0.561553798 -0.632305031 104 -0.683512853 -0.561553798 105 -0.659078648 -0.683512853 106 -0.589945995 -0.659078648 107 -0.854227050 -0.589945995 108 -0.649655747 -0.854227050 109 -0.744179747 -0.649655747 110 -0.849402382 -0.744179747 111 -1.072934131 -0.849402382 112 -1.001310602 -1.072934131 113 -0.825796372 -1.001310602 114 -0.665843946 -0.825796372 115 -0.371697282 -0.665843946 116 -0.324112429 -0.371697282 117 -0.546352860 -0.324112429 118 -0.255800559 -0.546352860 119 0.148319904 -0.255800559 120 -0.032008815 0.148319904 121 -0.347869918 -0.032008815 122 -0.115996909 -0.347869918 123 -0.053858219 -0.115996909 124 -0.252883626 -0.053858219 125 -0.377392271 -0.252883626 126 -0.361469486 -0.377392271 127 -0.125279837 -0.361469486 128 0.212359873 -0.125279837 129 0.247033252 0.212359873 130 0.727872070 0.247033252 131 1.225553935 0.727872070 132 2.405347917 1.225553935 133 2.747547519 2.405347917 134 1.852663507 2.747547519 135 1.351482800 1.852663507 136 1.214827359 1.351482800 137 1.178481991 1.214827359 138 -0.034464314 1.178481991 139 0.141741076 -0.034464314 140 0.541126678 0.141741076 141 0.371712072 0.541126678 142 0.011817011 0.371712072 143 -0.359414602 0.011817011 144 -0.140803499 -0.359414602 145 0.331313300 -0.140803499 146 0.636695928 0.331313300 147 0.430244212 0.636695928 148 0.194039632 0.430244212 149 0.505205475 0.194039632 150 0.381325487 0.505205475 151 0.463954168 0.381325487 152 1.023305858 0.463954168 153 0.484421036 1.023305858 154 1.449051534 0.484421036 155 1.904329337 1.449051534 156 0.926127235 1.904329337 157 0.906950993 0.926127235 158 1.738893088 0.906950993 159 1.183335934 1.738893088 160 1.354468283 1.183335934 161 1.242203033 1.354468283 162 1.365491152 1.242203033 163 0.821917160 1.365491152 164 0.794798773 0.821917160 165 0.912716977 0.794798773 166 1.277094580 0.912716977 167 1.228199296 1.277094580 168 1.619429546 1.228199296 169 1.683515441 1.619429546 170 1.778395405 1.683515441 171 1.559513544 1.778395405 172 1.382767044 1.559513544 173 1.199015789 1.382767044 174 1.376214537 1.199015789 175 1.106242038 1.376214537 176 0.888359285 1.106242038 177 0.894563325 0.888359285 178 1.067417348 0.894563325 179 1.089618712 1.067417348 180 0.991810396 1.089618712 181 0.396057586 0.991810396 182 0.326336118 0.396057586 183 0.200851845 0.326336118 184 0.217525154 0.200851845 185 0.088705772 0.217525154 186 0.105574027 0.088705772 187 0.203046291 0.105574027 188 0.392648555 0.203046291 189 0.237324007 0.392648555 190 0.274930341 0.237324007 191 0.107974517 0.274930341 192 0.410018348 0.107974517 193 0.212237851 0.410018348 194 0.142878355 0.212237851 195 0.036546397 0.142878355 196 0.027446858 0.036546397 197 -0.096897507 0.027446858 198 -0.440685208 -0.096897507 199 -0.185109008 -0.440685208 200 -0.311829607 -0.185109008 201 -0.472813506 -0.311829607 202 -0.745709798 -0.472813506 203 -0.730656346 -0.745709798 204 -0.837044536 -0.730656346 205 -0.908260779 -0.837044536 206 -1.077781553 -0.908260779 207 -1.312454612 -1.077781553 208 -1.476715465 -1.312454612 209 -1.259634672 -1.476715465 210 -1.054494600 -1.259634672 211 -1.148939932 -1.054494600 212 -1.272270052 -1.148939932 213 -0.886200765 -1.272270052 214 -0.831251547 -0.886200765 215 -1.017761696 -0.831251547 216 -1.011831572 -1.017761696 217 -1.235108566 -1.011831572 218 -1.291817057 -1.235108566 219 -1.308704715 -1.291817057 220 -1.362717566 -1.308704715 221 -1.276088301 -1.362717566 222 -1.217525945 -1.276088301 223 -1.303103906 -1.217525945 224 -1.241374651 -1.303103906 225 -1.259438198 -1.241374651 226 -0.987643286 -1.259438198 227 -0.809235886 -0.987643286 228 -0.989798651 -0.809235886 229 -0.889300152 -0.989798651 230 -0.580778643 -0.889300152 231 0.005358535 -0.580778643 232 0.227858946 0.005358535 233 0.306854861 0.227858946 234 0.047365797 0.306854861 235 0.155062270 0.047365797 236 -0.168843098 0.155062270 237 -0.006688865 -0.168843098 238 0.066133540 -0.006688865 239 1.359033188 0.066133540 240 1.283898524 1.359033188 241 1.258936591 1.283898524 242 0.613134341 1.258936591 243 0.585747728 0.613134341 244 0.430372830 0.585747728 245 0.488380440 0.430372830 246 0.865036623 0.488380440 247 0.586317600 0.865036623 248 0.335169139 0.586317600 249 0.080919731 0.335169139 250 -0.368430368 0.080919731 251 -0.990166748 -0.368430368 252 -0.966601129 -0.990166748 253 -0.484321766 -0.966601129 254 -0.264536356 -0.484321766 255 -0.753973978 -0.264536356 256 -0.806886458 -0.753973978 257 -0.733875342 -0.806886458 258 -0.968038517 -0.733875342 259 -1.297763254 -0.968038517 260 -1.382669325 -1.297763254 261 -1.481995410 -1.382669325 262 -1.133688266 -1.481995410 263 -0.911590774 -1.133688266 264 -1.227349541 -0.911590774 265 -1.050177374 -1.227349541 266 -0.950138175 -1.050177374 267 -1.000906290 -0.950138175 268 -1.137736820 -1.000906290 269 -1.332466672 -1.137736820 270 -0.784585475 -1.332466672 271 -0.424261737 -0.784585475 272 -0.063442687 -0.424261737 273 -0.126174111 -0.063442687 274 -0.305232930 -0.126174111 275 -0.386814616 -0.305232930 276 -0.375651131 -0.386814616 277 -0.674711301 -0.375651131 278 -0.957081491 -0.674711301 279 -0.926340455 -0.957081491 280 -0.356766199 -0.926340455 281 -0.319198381 -0.356766199 282 -0.485914388 -0.319198381 283 -0.606016622 -0.485914388 284 -0.806959383 -0.606016622 285 -0.932638128 -0.806959383 286 -0.910700965 -0.932638128 287 -0.943220260 -0.910700965 288 -1.235668552 -0.943220260 289 -1.338371592 -1.235668552 290 -1.629368690 -1.338371592 291 -1.477784843 -1.629368690 292 -1.351689645 -1.477784843 293 -1.296184899 -1.351689645 294 -1.153095226 -1.296184899 295 -1.087624192 -1.153095226 296 -1.217019288 -1.087624192 297 -1.260373459 -1.217019288 298 -1.172655244 -1.260373459 299 -1.149267945 -1.172655244 300 -1.381097772 -1.149267945 301 -1.605131673 -1.381097772 302 -1.528795824 -1.605131673 303 -1.873275502 -1.528795824 304 -2.026103477 -1.873275502 305 -2.113811995 -2.026103477 306 -2.086765409 -2.113811995 307 -2.101218029 -2.086765409 308 -2.268636343 -2.101218029 309 -2.249363135 -2.268636343 310 -2.411726607 -2.249363135 311 -1.871199077 -2.411726607 312 -2.092034366 -1.871199077 313 -2.347265150 -2.092034366 314 -2.550364482 -2.347265150 315 -2.645757419 -2.550364482 316 -2.511380737 -2.645757419 317 -2.669992470 -2.511380737 318 -2.802924992 -2.669992470 319 -2.945034019 -2.802924992 320 -3.206514477 -2.945034019 321 -3.170761507 -3.206514477 322 -2.939118519 -3.170761507 323 -3.172513631 -2.939118519 324 -3.297929885 -3.172513631 325 -3.096235396 -3.297929885 326 -3.281703489 -3.096235396 327 -3.230037420 -3.281703489 328 -3.135503360 -3.230037420 329 -2.863756187 -3.135503360 330 -3.160000964 -2.863756187 331 -3.287348680 -3.160000964 332 -3.240873709 -3.287348680 333 -3.453319231 -3.240873709 334 -3.596641159 -3.453319231 335 -3.455355051 -3.596641159 336 -3.375141704 -3.455355051 337 -3.482455375 -3.375141704 338 -3.575277068 -3.482455375 339 -3.480038465 -3.575277068 340 -3.475100457 -3.480038465 341 -3.294161483 -3.475100457 342 -2.635370056 -3.294161483 343 -2.145494350 -2.635370056 344 -2.712545696 -2.145494350 345 -2.574232417 -2.712545696 346 -3.022943197 -2.574232417 347 -3.137445418 -3.022943197 348 -2.880439940 -3.137445418 349 -2.731559518 -2.880439940 350 -2.951093926 -2.731559518 351 -2.305207603 -2.951093926 352 -1.472042613 -2.305207603 353 -1.072676718 -1.472042613 354 -1.548666823 -1.072676718 355 -1.137721289 -1.548666823 356 -1.649438105 -1.137721289 357 -2.040938725 -1.649438105 358 -2.292960073 -2.040938725 359 -2.084803269 -2.292960073 360 -1.746373839 -2.084803269 361 -1.982531293 -1.746373839 362 -2.280411729 -1.982531293 363 -2.146359156 -2.280411729 364 -1.613909103 -2.146359156 365 -1.491769680 -1.613909103 366 -1.560680036 -1.491769680 367 -1.523710417 -1.560680036 368 -1.741688219 -1.523710417 369 -1.637111588 -1.741688219 370 -1.497345718 -1.637111588 371 -1.689407589 -1.497345718 372 -1.966128349 -1.689407589 373 -1.689626770 -1.966128349 374 -2.084019782 -1.689626770 375 -1.872570424 -2.084019782 376 -1.785733187 -1.872570424 377 -1.432114075 -1.785733187 378 -1.427306874 -1.432114075 379 -1.531242298 -1.427306874 380 -1.328063563 -1.531242298 381 -1.274320939 -1.328063563 382 -0.905443645 -1.274320939 383 -0.744436718 -0.905443645 384 -0.786097918 -0.744436718 385 -0.818573483 -0.786097918 386 -0.636162256 -0.818573483 387 -0.960252415 -0.636162256 388 -1.114927447 -0.960252415 389 -1.278052043 -1.114927447 390 -1.295115154 -1.278052043 391 -1.411735494 -1.295115154 392 -1.235699697 -1.411735494 393 -0.897584818 -1.235699697 394 -0.904269212 -0.897584818 395 -0.896772153 -0.904269212 396 -1.124488034 -0.896772153 397 -1.053218043 -1.124488034 398 -1.194947956 -1.053218043 399 -1.142411923 -1.194947956 400 -1.225477111 -1.142411923 401 -1.298328105 -1.225477111 402 -1.122212868 -1.298328105 403 -0.993779503 -1.122212868 404 -1.091805169 -0.993779503 405 -1.232633302 -1.091805169 406 -1.200052036 -1.232633302 407 -1.322092168 -1.200052036 408 -1.195462169 -1.322092168 409 -0.985313783 -1.195462169 410 -1.068605664 -0.985313783 411 -0.743028069 -1.068605664 412 -0.677678033 -0.743028069 413 -0.793680952 -0.677678033 414 -0.746508800 -0.793680952 415 -0.805142671 -0.746508800 416 -0.804874209 -0.805142671 417 -0.841901947 -0.804874209 418 -0.779023487 -0.841901947 419 -0.738739764 -0.779023487 420 -0.356713897 -0.738739764 421 -0.336940220 -0.356713897 422 -0.403784571 -0.336940220 423 -0.515292718 -0.403784571 424 -0.546812446 -0.515292718 425 -0.372270528 -0.546812446 426 -0.234750066 -0.372270528 427 -0.249870851 -0.234750066 428 -0.343896455 -0.249870851 429 -0.308033578 -0.343896455 430 -0.218273705 -0.308033578 431 -0.158447668 -0.218273705 432 -0.126697034 -0.158447668 433 -0.179410519 -0.126697034 434 -0.281313227 -0.179410519 435 0.016307076 -0.281313227 436 -0.069247826 0.016307076 437 0.254857500 -0.069247826 438 0.152501676 0.254857500 439 0.198951738 0.152501676 440 0.428665471 0.198951738 441 0.834327107 0.428665471 442 0.843631267 0.834327107 443 0.724285935 0.843631267 444 0.501659874 0.724285935 445 0.477966840 0.501659874 446 0.549132336 0.477966840 447 0.878252241 0.549132336 448 0.871362794 0.878252241 449 0.907263999 0.871362794 450 0.612224452 0.907263999 451 0.275666519 0.612224452 452 0.426787405 0.275666519 453 0.611519713 0.426787405 454 0.564094832 0.611519713 455 0.838164773 0.564094832 456 0.987107075 0.838164773 457 1.213827899 0.987107075 458 1.486350198 1.213827899 459 1.540194239 1.486350198 460 1.464804854 1.540194239 461 1.950533686 1.464804854 462 2.043865359 1.950533686 463 1.592240601 2.043865359 464 1.568427553 1.592240601 465 1.447972511 1.568427553 466 1.569065643 1.447972511 467 1.810016229 1.569065643 468 2.112452560 1.810016229 469 2.221467970 2.112452560 470 2.262989001 2.221467970 471 2.227515073 2.262989001 472 2.164061733 2.227515073 473 2.214780332 2.164061733 474 2.006347567 2.214780332 475 1.944639871 2.006347567 476 1.445106726 1.944639871 477 1.215619391 1.445106726 478 1.486340115 1.215619391 479 1.443999721 1.486340115 480 1.589280178 1.443999721 481 1.739751310 1.589280178 482 1.571635541 1.739751310 483 1.579761180 1.571635541 484 1.557372565 1.579761180 485 1.720370667 1.557372565 486 1.467538051 1.720370667 487 1.302821189 1.467538051 488 1.291249621 1.302821189 489 0.968954202 1.291249621 490 1.073429670 0.968954202 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7rr0b1292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/8rr0b1292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/9k0he1292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/rcomp/tmp/10k0he1292971268.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/116jyj1292971268.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12r1w71292971268.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13y2t11292971268.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/14qcb41292971268.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/15cc9a1292971268.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/16iwrv1292971269.tab") + } > > try(system("convert tmp/1vz2k1292971268.ps tmp/1vz2k1292971268.png",intern=TRUE)) character(0) > try(system("convert tmp/2682n1292971268.ps tmp/2682n1292971268.png",intern=TRUE)) character(0) > try(system("convert tmp/3682n1292971268.ps tmp/3682n1292971268.png",intern=TRUE)) character(0) > try(system("convert tmp/4682n1292971268.ps tmp/4682n1292971268.png",intern=TRUE)) character(0) > try(system("convert tmp/5682n1292971268.ps tmp/5682n1292971268.png",intern=TRUE)) character(0) > try(system("convert tmp/6hz181292971268.ps tmp/6hz181292971268.png",intern=TRUE)) character(0) > try(system("convert tmp/7rr0b1292971268.ps tmp/7rr0b1292971268.png",intern=TRUE)) character(0) > try(system("convert tmp/8rr0b1292971268.ps tmp/8rr0b1292971268.png",intern=TRUE)) character(0) > try(system("convert tmp/9k0he1292971268.ps tmp/9k0he1292971268.png",intern=TRUE)) character(0) > try(system("convert tmp/10k0he1292971268.ps tmp/10k0he1292971268.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 13.703 2.390 37.805