R version 2.8.0 (2008-10-20) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale 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(8.155 + ,14.296 + ,133 + ,19.573 + ,278.295 + ,469.137 + ,94 + ,13.804 + ,7.014 + ,12.242 + ,117 + ,17.224 + ,9.682 + ,16.380 + ,100 + ,14.634 + ,6.196 + ,10.178 + ,111 + ,16.378 + ,21.303 + ,37.161 + ,122 + ,17.943 + ,15.031 + ,26.949 + ,109 + ,16.013 + ,12.637 + ,21.502 + ,126 + ,18.559 + ,9.353 + ,17.279 + ,99 + ,14.479 + ,5.606 + ,9.808 + ,107 + ,15.726 + ,4.581 + ,8.312 + ,138 + ,20.202 + ,9.655 + ,17.541 + ,113 + ,16.608 + ,14.851 + ,26.173 + ,123 + ,18.126 + ,11.332 + ,20.266 + ,122 + ,17.923 + ,4.509 + ,8.062 + ,117 + ,17.127 + ,8.043 + ,14.191 + ,107 + ,15.687 + ,14.626 + ,24.423 + ,118 + ,17.389 + ,5.377 + ,9.038 + ,102 + ,15.016 + ,10.169 + ,18.004 + ,114 + ,16.792 + ,8.452 + ,14.671 + ,116 + ,17.073 + ,4.503 + ,7.974 + ,113 + ,16.639 + ,10.763 + ,19.526 + ,131 + ,19.183 + ,19.207 + ,33.194 + ,114 + ,16.816 + ,10.194 + ,17.682 + ,118 + ,17.296 + ,5.282 + ,8.888 + ,112 + ,16.470 + ,7.120 + ,12.152 + ,111 + ,16.367 + ,10.374 + ,18.721 + ,100 + ,14.764 + ,7.099 + ,12.301 + ,113 + ,16.530 + ,11.843 + ,20.903 + ,118 + ,17.395 + ,11.191 + ,18.386 + ,116 + ,17.027 + ,6.202 + ,10.699 + ,104 + ,15.346 + ,8.259 + ,14.528 + ,130 + ,19.150 + ,11.673 + ,20.291 + ,114 + ,16.731 + ,9.642 + ,16.187 + ,113 + ,16.555 + ,23.073 + ,38.960 + ,105 + ,15.444 + ,20.102 + ,33.375 + ,111 + ,16.238 + ,47.035 + ,79.201 + ,106 + ,15.607 + ,12.181 + ,21.094 + ,104 + ,15.305 + ,9.293 + ,15.860 + ,108 + ,15.866 + ,9.341 + ,16.324 + ,114 + ,16.778 + ,4.389 + ,7.853 + ,109 + ,16.071 + ,11.293 + ,19.809 + ,111 + ,16.360 + ,13.594 + ,23.482 + ,105 + ,15.369 + ,6.969 + ,12.423 + ,95 + ,14.026 + ,1.306 + ,2.361 + ,94 + ,13.780 + ,11.855 + ,20.511 + ,102 + ,14.993 + ,9.266 + ,16.678 + ,108 + ,15.892 + ,4.983 + ,8.898 + ,100 + ,14.749 + ,21.305 + ,35.758 + ,108 + ,15.835 + ,6.105 + ,10.808 + ,110 + ,16.087 + ,16.081 + ,26.210 + ,108 + ,15.796 + ,4.980 + ,8.498 + ,103 + ,15.190 + ,8.225 + ,14.018 + ,103 + ,15.153 + ,10.576 + ,19.049 + ,98 + ,14.458 + ,5.576 + ,9.482 + ,105 + ,15.360 + ,10.143 + ,18.020 + ,106 + ,15.584 + ,8.884 + ,15.039 + ,104 + ,15.207 + ,8.985 + ,15.716 + ,105 + ,15.432 + ,5.634 + ,9.520 + ,106 + ,15.552 + ,4.485 + ,8.295 + ,101 + ,14.851 + ,19.285 + ,33.230 + ,106 + ,15.557 + ,6.678 + ,11.507 + ,108 + ,15.858 + ,7.101 + ,12.655 + ,112 + ,16.401 + ,7.739 + ,14.079 + ,95 + ,13.941 + ,5.742 + ,10.515 + ,98 + ,14.427 + ,6.192 + ,10.914 + ,95 + ,13.912 + ,24.554 + ,39.966 + ,103 + ,15.194 + ,4.386 + ,7.409 + ,103 + ,15.146 + ,5.697 + ,10.231 + ,120 + ,17.576 + ,13.625 + ,23.331 + ,109 + ,15.960 + ,55.457 + ,98.343 + ,75 + ,10.980 + ,16.181 + ,29.883 + ,107 + ,15.669 + ,11.927 + ,20.703 + ,99 + ,14.554 + ,79.009 + ,147.395 + ,76 + ,11.237 + ,21.129 + ,42.622 + ,81 + ,11.836 + ,19.132 + ,34.427 + ,86 + ,12.651 + ,27.088 + ,48.595 + ,87 + ,12.749 + ,13.241 + ,21.569 + ,97 + ,14.227 + ,41.602 + ,78.928 + ,88 + ,12.897 + ,26.141 + ,44.082 + ,93 + ,13.602 + ,10.678 + ,18.780 + ,82 + ,12.059 + ,42.356 + ,82.653 + ,65 + ,9.496 + ,23.658 + ,45.001 + ,69 + ,10.126 + ,11.936 + ,23.931 + ,55 + ,8.100 + ,59.836 + ,114.766 + ,72 + ,10.518 + ,39.267 + ,76.654 + ,110 + ,16.101 + ,13.664 + ,24.127 + ,113 + ,16.628 + ,25.699 + ,48.788 + ,102 + ,14.940 + ,19.016 + ,38.602 + ,114 + ,16.719 + ,6.849 + ,12.199 + ,118 + ,17.379 + ,16.942 + ,29.446 + ,116 + ,17.066 + ,12.988 + ,23.582 + ,120 + ,17.652 + ,1.130 + ,2.070 + ,112 + ,16.435 + ,22.943 + ,39.619 + ,123 + ,18.074 + ,2.835 + ,4.857 + ,99 + ,14.546 + ,4.732 + ,8.221 + ,116 + ,17.085 + ,5.016 + ,8.905 + ,119 + ,17.517 + ,20.042 + ,34.423 + ,127 + ,18.613 + ,20.786 + ,35.229 + ,111 + ,16.257 + ,3.666 + ,6.524 + ,114 + ,16.746 + ,5.381 + ,10.142 + ,120 + ,17.593 + ,6.272 + ,11.075 + ,125 + ,18.415 + ,5.235 + ,8.942 + ,120 + ,17.629 + ,6.097 + ,13.092 + ,119 + ,17.446 + ,4.931 + ,8.738 + ,127 + ,18.693 + ,7.093 + ,12.162 + ,110 + ,16.149 + ,2.664 + ,4.715 + ,123 + ,18.125 + ,10.084 + ,17.488 + ,117 + ,17.217 + ,7.252 + ,12.872 + ,102 + ,15.044 + ,10.439 + ,18.505 + ,135 + ,19.792 + ,8.663 + ,15.081 + ,116 + ,17.012 + ,7.147 + ,12.484 + ,119 + ,17.443 + ,12.040 + ,24.279 + ,123 + ,18.024 + ,2.479 + ,4.379 + ,122 + ,17.955 + ,6.110 + ,10.861 + ,116 + ,17.069 + ,9.125 + ,17.998 + ,135 + ,19.777 + ,17.724 + ,31.024 + ,107 + ,15.756 + ,6.361 + ,11.094 + ,126 + ,18.557 + ,8.481 + ,14.773 + ,120 + ,17.659 + ,21.670 + ,38.260 + ,104 + ,15.314 + ,8.455 + ,14.956 + ,125 + ,18.294 + ,6.431 + ,13.515 + ,119 + ,17.543 + ,15.695 + ,29.267 + ,111 + ,16.322 + ,12.215 + ,21.626 + ,129 + ,19.020 + ,16.985 + ,28.075 + ,112 + ,16.459 + ,5.509 + ,9.557 + ,112 + ,16.453 + ,3.349 + ,5.920 + ,104 + ,15.257 + ,5.100 + ,9.365 + ,123 + ,18.042 + ,5.098 + ,9.256 + ,128 + ,18.799 + ,6.498 + ,11.617 + ,129 + ,18.907 + ,4.213 + ,7.528 + ,125 + ,18.364 + ,13.679 + ,22.935 + ,108 + ,15.869 + ,3.339 + ,5.811 + ,102 + ,14.936 + ,2.865 + ,5.231 + ,111 + ,16.293 + ,7.810 + ,13.811 + ,118 + ,17.378 + ,10.848 + ,19.479 + ,131 + ,19.293 + ,3.560 + ,6.291 + ,119 + ,17.547 + ,5.120 + ,9.427 + ,122 + ,17.955 + ,4.994 + ,9.280 + ,121 + ,17.780 + ,6.921 + ,12.608 + ,142 + ,20.906 + ,4.308 + ,7.582 + ,105 + ,15.361 + ,10.372 + ,18.543 + ,125 + ,18.361 + ,8.788 + ,14.967 + ,112 + ,16.471 + ,55.767 + ,92.323 + ,118 + ,17.289 + ,4.019 + ,7.043 + ,111 + ,16.299 + ,7.486 + ,13.762 + ,131 + ,19.288 + ,5.889 + ,10.899 + ,138 + ,20.228 + ,8.479 + ,15.222 + ,120 + ,17.668 + ,13.187 + ,22.208 + ,107 + ,15.654 + ,9.933 + ,20.892 + ,120 + ,17.612 + ,5.821 + ,10.161 + ,113 + ,16.607 + ,19.645 + ,32.145 + ,110 + ,16.143 + ,7.991 + ,13.977 + ,115 + ,16.870 + ,4.779 + ,8.015 + ,106 + ,15.506 + ,3.602 + ,6.608 + ,117 + ,17.222 + ,21.505 + ,37.620 + ,120 + ,17.653 + ,5.378 + ,9.548 + ,113 + ,16.574 + ,3.641 + ,6.855 + ,119 + ,17.439 + ,5.885 + ,11.099 + ,126 + ,18.449 + ,5.154 + ,9.577 + ,112 + ,16.437 + ,7.978 + ,14.341 + ,114 + ,16.700 + ,6.874 + ,12.540 + ,126 + ,18.500 + ,1.784 + ,3.128 + ,107 + ,15.673 + ,2.514 + ,4.679 + ,111 + ,16.298 + ,3.423 + ,6.112 + ,114 + ,16.739 + ,7.271 + ,12.703 + ,107 + ,15.671 + ,3.965 + ,7.326 + ,123 + ,18.075 + ,7.124 + ,14.018 + ,139 + ,20.448 + ,3.645 + ,6.586 + ,117 + ,17.116 + ,15.202 + ,24.828 + ,115 + ,16.931 + ,4.408 + ,7.993 + ,107 + ,15.695 + ,16.341 + ,30.235 + ,103 + ,15.135 + ,4.294 + ,7.697 + ,110 + ,16.126 + ,3.175 + ,5.900 + ,105 + ,15.462 + ,5.761 + ,10.246 + ,104 + ,15.242 + ,11.727 + ,21.498 + ,124 + ,18.159 + ,13.005 + ,22.829 + ,97 + ,14.226 + ,5.386 + ,9.778 + ,113 + ,16.536 + ,3.269 + ,6.124 + ,121 + ,17.751 + ,14.968 + ,29.423 + ,116 + ,17.095 + ,18.943 + ,32.680 + ,114 + ,16.706 + ,8.405 + ,14.748 + ,104 + ,15.282 + ,12.020 + ,18.389 + ,97 + ,14.314 + ,70.511 + ,117.027 + ,108 + ,15.823 + ,6.158 + ,10.864 + ,102 + ,14.918 + ,7.398 + ,13.645 + ,108 + ,15.843 + ,21.308 + ,34.079 + ,133 + ,19.477 + ,12.316 + ,22.029 + ,106 + ,15.514 + ,11.570 + ,19.647 + ,101 + ,14.889 + ,11.978 + ,22.000 + ,100 + ,14.699 + ,1.546 + ,2.755 + ,107 + ,15.663 + ,9.130 + ,15.881 + ,90 + ,13.183 + ,5.226 + ,9.205 + ,87 + ,12.722 + ,4.837 + ,8.333 + ,92 + ,13.448 + ,6.837 + ,12.039 + ,91 + ,13.358 + ,1.839 + ,3.331 + ,81 + ,11.891 + ,4.550 + ,8.172 + ,83 + ,12.135 + ,20.470 + ,34.865 + ,95 + ,14.030 + ,4.273 + ,7.783 + ,88 + ,12.938 + ,560 + ,990 + ,78 + ,11.519 + ,11.147 + ,19.772 + ,90 + ,13.197 + ,2.042 + ,3.656 + ,80 + ,11.769 + ,10.107 + ,17.732 + ,92 + ,13.462 + ,6.693 + ,12.000 + ,91 + ,13.413 + ,7.696 + ,14.110 + ,102 + ,14.966 + ,5.520 + ,9.560 + ,98 + ,14.371 + ,6.305 + ,11.327 + ,101 + ,14.797 + ,15.090 + ,26.311 + ,99 + ,14.607 + ,43.413 + ,73.859 + ,103 + ,15.168 + ,7.131 + ,12.637 + ,99 + ,14.599 + ,3.092 + ,5.439 + ,99 + ,14.556 + ,18.743 + ,32.444 + ,91 + ,13.417 + ,1.138 + ,2.048 + ,92 + ,13.573 + ,20.382 + ,35.897 + ,103 + ,15.112 + ,17.304 + ,30.928 + ,99 + ,14.517 + ,13.142 + ,23.853 + ,102 + ,15.022 + ,9.196 + ,15.460 + ,102 + ,14.965 + ,7.585 + ,12.210 + ,108 + ,15.926 + ,6.365 + ,11.277 + ,96 + ,14.100 + ,7.570 + ,13.533 + ,92 + ,13.539 + ,11.344 + ,18.229 + ,96 + ,14.117 + ,45.037 + ,69.145 + ,100 + ,14.636 + ,5.037 + ,8.951 + ,98 + ,14.354 + ,5.425 + ,9.897 + ,97 + ,14.217 + ,5.923 + ,10.633 + ,97 + ,14.182 + ,15.461 + ,26.633 + ,101 + ,14.775 + ,5.195 + ,9.375 + ,94 + ,13.745 + ,4.642 + ,8.463 + ,97 + ,14.254 + ,6.105 + ,10.768 + ,93 + ,13.704 + ,32.793 + ,56.407 + ,102 + ,14.965 + ,6.098 + ,11.017 + ,88 + ,12.996 + ,5.197 + ,9.098 + ,96 + ,14.140 + ,4.161 + ,8.207 + ,93 + ,13.727 + ,6.140 + ,11.007 + ,92 + ,13.545 + ,4.092 + ,7.480 + ,94 + ,13.881 + ,3.597 + ,6.597 + ,94 + ,13.777 + ,2.885 + ,5.133 + ,95 + ,13.903 + ,10.799 + ,19.415 + ,100 + ,14.672 + ,4.972 + ,9.047 + ,94 + ,13.876 + ,7.354 + ,13.394 + ,93 + ,13.654 + ,2.746 + ,4.959 + ,83 + ,12.205 + ,6.520 + ,10.219 + ,95 + ,13.968 + ,13.284 + ,21.467 + ,130 + ,19.027 + ,6.979 + ,11.001 + ,98 + ,14.375 + ,6.922 + ,11.790 + ,93 + ,13.589 + ,47.352 + ,78.030 + ,111 + ,16.324 + ,10.535 + ,17.715 + ,112 + ,16.442 + ,11.341 + ,19.128 + ,117 + ,17.130 + ,18.588 + ,31.605 + ,105 + ,15.498 + ,10.325 + ,17.485 + ,115 + ,16.868 + ,9.841 + ,16.910 + ,112 + ,16.509 + ,9.930 + ,17.334 + ,114 + ,16.678 + ,21.171 + ,36.046 + ,109 + ,15.940 + ,5.516 + ,9.507 + ,109 + ,16.023 + ,14.675 + ,24.619 + ,116 + ,16.998 + ,8.212 + ,14.362 + ,107 + ,15.734 + ,8.027 + ,13.931 + ,114 + ,16.803 + ,25.978 + ,43.569 + ,107 + ,15.730 + ,13.565 + ,23.457 + ,98 + ,14.362 + ,6.887 + ,11.809 + ,115 + ,16.899 + ,10.218 + ,17.647 + ,107 + ,15.704 + ,5.733 + ,10.320 + ,118 + ,17.368 + ,14.076 + ,23.359 + ,109 + ,15.985 + ,6.197 + ,11.081 + ,110 + ,16.216 + ,11.413 + ,20.450 + ,97 + ,14.246 + ,7.861 + ,13.610 + ,104 + ,15.291 + ,12.047 + ,19.735 + ,101 + ,14.832 + ,3.494 + ,6.171 + ,106 + ,15.632 + ,12.806 + ,22.381 + ,99 + ,14.592 + ,3.800 + ,6.591 + ,99 + ,14.513 + ,7.337 + ,12.216 + ,102 + ,15.027 + ,10.494 + ,18.957 + ,108 + ,15.910 + ,5.660 + ,10.262 + ,139 + ,20.393 + ,16.514 + ,28.696 + ,108 + ,15.806 + ,9.992 + ,17.345 + ,120 + ,17.691 + ,18.510 + ,32.459 + ,110 + ,16.126 + ,6.947 + ,12.202 + ,113 + ,16.593 + ,145.933 + ,236.196 + ,105 + ,15.391 + ,4.574 + ,7.946 + ,99 + ,14.576 + ,11.405 + ,20.564 + ,115 + ,16.838 + ,5.257 + ,9.405 + ,118 + ,17.358 + ,6.127 + ,10.685 + ,116 + ,17.048 + ,12.857 + ,22.564 + ,121 + ,17.845 + ,3.343 + ,5.923 + ,108 + ,15.888 + ,6.343 + ,11.021 + ,113 + ,16.618 + ,6.317 + ,11.403 + ,108 + ,15.841 + ,7.342 + ,13.217 + ,114 + ,16.701 + ,4.605 + ,8.253 + ,152 + ,22.269 + ,4.506 + ,7.766 + ,103 + ,15.095 + ,4.008 + ,6.950 + ,109 + ,15.959 + ,4.658 + ,8.018 + ,105 + ,15.367 + ,8.334 + ,14.671 + ,102 + ,14.979 + ,8.154 + ,13.820 + ,109 + ,15.964 + ,1.143 + ,2.063 + ,109 + ,16.059 + ,3.608 + ,6.320 + ,101 + ,14.833 + ,4.462 + ,8.150 + ,108 + ,15.813 + ,3.623 + ,6.464 + ,110 + ,16.226 + ,3.699 + ,6.453 + ,107 + ,15.782 + ,17.103 + ,28.935 + ,109 + ,15.996 + ,13.826 + ,24.532 + ,87 + ,12.736 + ,3.366 + ,6.110 + ,111 + ,16.249 + ,3.932 + ,6.887 + ,108 + ,15.901 + ,4.406 + ,7.877 + ,111 + ,16.343 + ,25.953 + ,45.829 + ,111 + ,16.337 + ,8.746 + ,15.408 + ,107 + ,15.747 + ,21.624 + ,38.382 + ,101 + ,14.822 + ,10.162 + ,18.275 + ,105 + ,15.374 + ,41.500 + ,70.233 + ,102 + ,15.004 + ,9.663 + ,17.068 + ,103 + ,15.127 + ,15.252 + ,26.976 + ,105 + ,15.430 + ,16.211 + ,27.363 + ,100 + ,14.706 + ,7.783 + ,13.492 + ,90 + ,13.274 + ,6.593 + ,11.562 + ,86 + ,12.680 + ,1.932 + ,3.391 + ,94 + ,13.775 + ,3.561 + ,6.369 + ,96 + ,14.067 + ,3.317 + ,5.747 + ,101 + ,14.888 + ,1.905 + ,3.278 + ,100 + ,14.654 + ,6.355 + ,11.071 + ,95 + ,13.995 + ,5.976 + ,10.720 + ,85 + ,12.529 + ,8.174 + ,14.142 + ,83 + ,12.240 + ,117.696 + ,201.571 + ,78 + ,11.433 + ,20.776 + ,35.784 + ,78 + ,11.482 + ,17.180 + ,29.879 + ,85 + ,12.512 + ,6.020 + ,11.034 + ,70 + ,10.226 + ,12.990 + ,22.269 + ,87 + ,12.840 + ,9.788 + ,16.775 + ,84 + ,12.384 + ,6.853 + ,12.048 + ,114 + ,16.795 + ,5.112 + ,8.963 + ,106 + ,15.607 + ,12.900 + ,22.421 + ,84 + ,12.274 + ,6.021 + ,10.192 + ,115 + ,16.905 + ,9.198 + ,16.423 + ,102 + ,14.971 + ,6.033 + ,10.677 + ,102 + ,14.946 + ,11.794 + ,20.155 + ,79 + ,11.629 + ,11.416 + ,20.021 + ,74 + ,10.941 + ,9.873 + ,16.890 + ,83 + ,12.123 + ,12.235 + ,20.728 + ,86 + ,12.630 + ,3.810 + ,6.728 + ,84 + ,12.331 + ,2.843 + ,5.017 + ,95 + ,13.905 + ,4.920 + ,9.720 + ,105 + ,15.369 + ,2.187 + ,4.080 + ,97 + ,14.254 + ,53.206 + ,91.174 + ,89 + ,13.011 + ,10.872 + ,18.714 + ,78 + ,11.447 + ,4.402 + ,7.746 + ,102 + ,15.029 + ,3.940 + ,6.556 + ,83 + ,12.215 + ,13.356 + ,22.622 + ,93 + ,13.643 + ,9.899 + ,17.597 + ,82 + ,12.015 + ,30.531 + ,53.467 + ,84 + ,12.336 + ,11.700 + ,20.620 + ,102 + ,14.933 + ,5.871 + ,10.207 + ,99 + ,14.506 + ,7.185 + ,12.387 + ,109 + ,16.084 + ,44.601 + ,77.562 + ,83 + ,12.134 + ,4.475 + ,7.994 + ,102 + ,14.971 + ,10.543 + ,18.085 + ,94 + ,13.874 + ,4.321 + ,8.042 + ,110 + ,16.202 + ,14.135 + ,25.785 + ,96 + ,14.079 + ,6.825 + ,11.632 + ,82 + ,12.102 + ,3.782 + ,6.755 + ,87 + ,12.798 + ,19.290 + ,32.591 + ,87 + ,12.713 + ,5.622 + ,9.898 + ,79 + ,11.559 + ,5.501 + ,9.588 + ,82 + ,12.085 + ,4.216 + ,7.455 + ,93 + ,13.627 + ,2.085 + ,3.604 + ,81 + ,11.893 + ,7.508 + ,13.377 + ,114 + ,16.791 + ,3.083 + ,5.550 + ,101 + ,14.790 + ,2.187 + ,4.130 + ,83 + ,12.129 + ,2.914 + ,5.148 + ,75 + ,11.016 + ,11.033 + ,18.768 + ,85 + ,12.556 + ,2.657 + ,4.688 + ,84 + ,12.363 + ,8.308 + ,14.640 + ,102 + ,14.943 + ,4.482 + ,7.550 + ,88 + ,12.947 + ,4.325 + ,7.726 + ,88 + ,12.937 + ,3.056 + ,5.449 + ,91 + ,13.430 + ,5.351 + ,9.826 + ,90 + ,13.216 + ,7.756 + ,13.315 + ,93 + ,13.721 + ,1.918 + ,3.423 + ,95 + ,13.997 + ,2.933 + ,5.324 + ,92 + ,13.548 + ,9.489 + ,16.883 + ,82 + ,12.087 + ,2.888 + ,5.046 + ,91 + ,13.368 + ,40.378 + ,68.018 + ,94 + ,13.788 + ,7.744 + ,13.213 + ,92 + ,13.543 + ,2.261 + ,4.007 + ,104 + ,15.341 + ,1.583 + ,2.870 + ,104 + ,15.319 + ,2.380 + ,4.238 + ,96 + ,14.168 + ,3.280 + ,5.742 + ,86 + ,12.691 + ,2.551 + ,4.519 + ,96 + ,14.045 + ,2.116 + ,3.618 + ,93 + ,13.695 + ,2.532 + ,4.593 + ,100 + ,14.646 + ,12.499 + ,20.263 + ,93 + ,13.722 + ,2.896 + ,5.151 + ,95 + ,14.021 + ,2.091 + ,3.736 + ,103 + ,15.122 + ,3.169 + ,5.669 + ,118 + ,17.390 + ,1.487 + ,2.582 + ,97 + ,14.321 + ,1.288 + ,2.397 + ,102 + ,15.047 + ,1.946 + ,3.650 + ,105 + ,15.422 + ,3.374 + ,6.142 + ,99 + ,14.491 + ,7.220 + ,12.870 + ,101 + ,14.792 + ,16.349 + ,27.483 + ,94 + ,13.826 + ,4.996 + ,8.737 + ,105 + ,15.358 + ,6.435 + ,11.160 + ,95 + ,13.886 + ,4.905 + ,8.439 + ,97 + ,14.201 + ,6.779 + ,11.798 + ,91 + ,13.363 + ,7.304 + ,12.932 + ,106 + ,15.569 + ,11.935 + ,20.934 + ,123 + ,18.061 + ,3.173 + ,5.383 + ,88 + ,12.898 + ,3.677 + ,6.616 + ,108 + ,15.890 + ,7.851 + ,13.069 + ,113 + ,16.556 + ,14.750 + ,24.992 + ,89 + ,13.136 + ,9.395 + ,16.150 + ,94 + ,13.751 + ,12.687 + ,21.833 + ,86 + ,12.615 + ,22.124 + ,37.702 + ,82 + ,12.102 + ,4.886 + ,8.555 + ,107 + ,15.728 + ,112.237 + ,189.504 + ,85 + ,12.439 + ,5.420 + ,9.716 + ,127 + ,18.611 + ,13.897 + ,23.651 + ,96 + ,14.067 + ,13.110 + ,22.591 + ,79 + ,11.556 + ,36.631 + ,61.447 + ,81 + ,11.932 + ,8.813 + ,15.654 + ,97 + ,14.234 + ,7.325 + ,13.080 + ,107 + ,15.725 + ,4.334 + ,7.718 + ,97 + ,14.213 + ,9.665 + ,16.877 + ,92 + ,13.573 + ,2.352 + ,4.156 + ,99 + ,14.498 + ,2.307 + ,4.133 + ,97 + ,14.318 + ,8.020 + ,14.467 + ,75 + ,10.952 + ,9.416 + ,16.747 + ,99 + ,14.597 + ,4.524 + ,8.022 + ,107 + ,15.685 + ,1.873 + ,3.351 + ,90 + ,13.158 + ,3.194 + ,5.662 + ,91 + ,13.312 + ,6.869 + ,11.910 + ,92 + ,13.571 + ,2.099 + ,3.785 + ,118 + ,17.363 + ,5.376 + ,9.574 + ,90 + ,13.211 + ,5.404 + ,9.773 + ,88 + ,12.967 + ,6.601 + ,10.510 + ,96 + ,14.036 + ,3.926 + ,6.687 + ,91 + ,13.386 + ,1.688 + ,3.010 + ,101 + ,14.849 + ,6.668 + ,11.717 + ,102 + ,15.051 + ,2.837 + ,5.380 + ,104 + ,15.237 + ,1.452 + ,2.516 + ,90 + ,13.255 + ,31.389 + ,54.334 + ,82 + ,12.019 + ,3.784 + ,6.820 + ,89 + ,13.004 + ,5.429 + ,9.290 + ,96 + ,14.109 + ,2.686 + ,5.338 + ,83 + ,12.250 + ,2.601 + ,5.462 + ,79 + ,11.634 + ,1.996 + ,3.952 + ,82 + ,12.084 + ,2.930 + ,5.597 + ,82 + ,12.070 + ,9.787 + ,18.360 + ,98 + ,14.437 + ,5.359 + ,10.519 + ,84 + ,12.291 + ,2.611 + ,5.221 + ,93 + ,13.664 + ,5.028 + ,10.255 + ,97 + ,14.230 + ,4.939 + ,9.216 + ,93 + ,13.589 + ,1.565 + ,2.833 + ,104 + ,15.294 + ,3.166 + ,5.639 + ,102 + ,15.010 + ,1.601 + ,2.867 + ,110 + ,16.195 + ,1.574 + ,2.859 + ,111 + ,16.338 + ,1.988 + ,3.577 + ,110 + ,16.218 + ,1.732 + ,3.056 + ,104 + ,15.233 + ,1.585 + ,2.962 + ,105 + ,15.403 + ,8.232 + ,14.596 + ,105 + ,15.455 + ,1.738 + ,3.047 + ,104 + ,15.335 + ,2.003 + ,3.506 + ,105 + ,15.366 + ,2.944 + ,5.164 + ,109 + ,15.965 + ,3.847 + ,6.563 + ,91 + ,13.437 + ,8.635 + ,14.252 + ,104 + ,15.286 + ,1.424 + ,2.581 + ,94 + ,13.821 + ,4.251 + ,7.592 + ,99 + ,14.599 + ,23.050 + ,41.757 + ,92 + ,13.570 + ,10.052 + ,17.937 + ,105 + ,15.437 + ,35.964 + ,64.194 + ,89 + ,13.023 + ,4.536 + ,7.922 + ,103 + ,15.160 + ,5.056 + ,8.745 + ,103 + ,15.070 + ,5.611 + ,9.899 + ,101 + ,14.853 + ,43.970 + ,71.063 + ,112 + ,16.460 + ,6.788 + ,11.916 + ,105 + ,15.395 + ,17.047 + ,31.147 + ,93 + ,13.713 + ,8.488 + ,14.554 + ,96 + ,14.071 + ,7.928 + ,13.905 + ,104 + ,15.311 + ,3.879 + ,6.696 + ,99 + ,14.577 + ,5.280 + ,9.762 + ,97 + ,14.269 + ,23.525 + ,38.627 + ,103 + ,15.102 + ,9.981 + ,17.113 + ,107 + ,15.785 + ,11.471 + ,20.183 + ,103 + ,15.061 + ,3.918 + ,6.970 + ,101 + ,14.859 + ,6.876 + ,12.535 + ,89 + ,13.140 + ,8.429 + ,14.644 + ,97 + ,14.210 + ,10.793 + ,19.233 + ,91 + ,13.345 + ,7.693 + ,13.880 + ,90 + ,13.157 + ,6.390 + ,11.564 + ,95 + ,13.934 + ,16.513 + ,30.021 + ,89 + ,13.091 + ,6.516 + ,12.108 + ,89 + ,13.021 + ,17.937 + ,32.282 + ,95 + ,13.927 + ,13.589 + ,24.009 + ,96 + ,14.072 + ,6.857 + ,12.637 + ,95 + ,13.908 + ,9.005 + ,16.201 + ,96 + ,14.170 + ,7.427 + ,13.650 + ,95 + ,13.993 + ,9.015 + ,15.969 + ,97 + ,14.200 + ,6.216 + ,11.121 + ,105 + ,15.470 + ,17.455 + ,30.278 + ,98 + ,14.327 + ,6.014 + ,10.272 + ,101 + ,14.778 + ,3.978 + ,6.855 + ,98 + ,14.356 + ,56 + ,84 + ,117 + ,17.213 + ,5.409 + ,9.377 + ,100 + ,14.699 + ,4.685 + ,8.152 + ,104 + ,15.302 + ,14.213 + ,24.829 + ,101 + ,14.801 + ,20.890 + ,36.558 + ,86 + ,12.655 + ,8.987 + ,16.110 + ,96 + ,14.140 + ,17.813 + ,29.832 + ,101 + ,14.853 + ,2.298 + ,4.234 + ,93 + ,13.651 + ,4.065 + ,7.050 + ,102 + ,15.047 + ,14.443 + ,26.738 + ,119 + ,17.495 + ,2.249 + ,4.873 + ,122 + ,17.969 + ,8.045 + ,15.263 + ,97 + ,14.211 + ,883 + ,1.575 + ,95 + ,13.964 + ,3.955 + ,7.543 + ,118 + ,17.334 + ,7.905 + ,14.481 + ,88 + ,12.997 + ,1.576 + ,3.025 + ,85 + ,12.454 + ,1.073 + ,2.093 + ,103 + ,15.123 + ,2.606 + ,4.803 + ,82 + ,12.043 + ,2.755 + ,4.810 + ,87 + ,12.748 + ,1.310 + ,2.325 + ,85 + ,12.544 + ,2.487 + ,4.817 + ,87 + ,12.844 + ,4.285 + ,7.372 + ,89 + ,13.091 + ,6.276 + ,10.706 + ,86 + ,12.697 + ,1.672 + ,2.948 + ,86 + ,12.667 + ,2.807 + ,5.076 + ,83 + ,12.150 + ,2.521 + ,4.349 + ,82 + ,11.993 + ,1.869 + ,3.207 + ,89 + ,13.102 + ,9.606 + ,17.100 + ,88 + ,12.948 + ,2.772 + ,5.069 + ,88 + ,12.984 + ,1.371 + ,2.316 + ,83 + ,12.213 + ,1.469 + ,2.585 + ,90 + ,13.150 + ,4.678 + ,8.187 + ,88 + ,12.989 + ,3.278 + ,5.489 + ,85 + ,12.427 + ,778 + ,1.369 + ,86 + ,12.563 + ,893 + ,1.547 + ,87 + ,12.718 + ,2.216 + ,4.259 + ,93 + ,13.676 + ,2.623 + ,4.681 + ,92 + ,13.554 + ,5.567 + ,10.016 + ,95 + ,13.897 + ,3.660 + ,6.676 + ,90 + ,13.240 + ,2.829 + ,5.072 + ,84 + ,12.327 + ,3.294 + ,5.774 + ,90 + ,13.250 + ,1.321 + ,2.392 + ,89 + ,13.019 + ,1.630 + ,2.974 + ,83 + ,12.255 + ,2.758 + ,5.071 + ,91 + ,13.297 + ,2.641 + ,5.429 + ,111 + ,16.308 + ,3.259 + ,5.447 + ,87 + ,12.801 + ,4.097 + ,7.964 + ,108 + ,15.891 + ,1.474 + ,2.699 + ,92 + ,13.533 + ,2.167 + ,4.264 + ,97 + ,14.205 + ,1.080 + ,2.037 + ,86 + ,12.631 + ,1.647 + ,3.212 + ,107 + ,15.668 + ,1.907 + ,3.758 + ,91 + ,13.332 + ,5.890 + ,11.230 + ,91 + ,13.326 + ,4.022 + ,6.986 + ,90 + ,13.196 + ,4.852 + ,8.494 + ,84 + ,12.410 + ,1.822 + ,3.163 + ,81 + ,11.949 + ,8.816 + ,15.197 + ,93 + ,13.641 + ,7.872 + ,13.164 + ,88 + ,12.950 + ,2.531 + ,4.451 + ,84 + ,12.384 + ,3.595 + ,6.722 + ,96 + ,14.121 + ,3.261 + ,5.348 + ,84 + ,12.326 + ,2.702 + ,4.953 + ,90 + ,13.197 + ,2.613 + ,4.482 + ,90 + ,13.232 + ,1.747 + ,3.131 + ,90 + ,13.190 + ,6.816 + ,12.115 + ,85 + ,12.550 + ,2.631 + ,4.741 + ,87 + ,12.772 + ,1.586 + ,2.831 + ,76 + ,11.231 + ,4.831 + ,8.586 + ,101 + ,14.830 + ,14.237 + ,24.686 + ,89 + ,13.061 + ,3.447 + ,6.316 + ,110 + ,16.146 + ,8.178 + ,14.680 + ,106 + ,15.568 + ,3.671 + ,6.804 + ,104 + ,15.311 + ,4.194 + ,7.536 + ,103 + ,15.095 + ,5.132 + ,9.460 + ,92 + ,13.511 + ,12.580 + ,22.252 + ,110 + ,16.091 + ,3.578 + ,6.503 + ,102 + ,14.943 + ,10.006 + ,18.109 + ,94 + ,13.848 + ,4.502 + ,8.428 + ,113 + ,16.629 + ,6.730 + ,12.135 + ,93 + ,13.637 + ,64.771 + ,107.796 + ,100 + ,14.749 + ,2.453 + ,4.387 + ,93 + ,13.602 + ,6.369 + ,11.462 + ,104 + ,15.350 + ,15.968 + ,27.083 + ,89 + ,13.096 + ,4.247 + ,7.766 + ,100 + ,14.720 + ,2.579 + ,4.601 + ,84 + ,12.411 + ,7.872 + ,13.512 + ,80 + ,11.793 + ,1.640 + ,2.913 + ,84 + ,12.307 + ,6.160 + ,10.853 + ,91 + ,13.376 + ,4.989 + ,8.554 + ,88 + ,12.870 + ,3.493 + ,5.743 + ,79 + ,11.667 + ,9.952 + ,17.728 + ,98 + ,14.428) + ,dim=c(4 + ,589) + ,dimnames=list(c('Totaal' + ,'inwoners' + ,'welvaart' + ,'gem_ink ') + ,1:589)) > y <- array(NA,dim=c(4,589),dimnames=list(c('Totaal','inwoners','welvaart','gem_ink '),1:589)) > 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 = '1' > #'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 Totaal inwoners welvaart gem_ink\r 1 8.155 14.296 133 19.573 2 278.295 469.137 94 13.804 3 7.014 12.242 117 17.224 4 9.682 16.380 100 14.634 5 6.196 10.178 111 16.378 6 21.303 37.161 122 17.943 7 15.031 26.949 109 16.013 8 12.637 21.502 126 18.559 9 9.353 17.279 99 14.479 10 5.606 9.808 107 15.726 11 4.581 8.312 138 20.202 12 9.655 17.541 113 16.608 13 14.851 26.173 123 18.126 14 11.332 20.266 122 17.923 15 4.509 8.062 117 17.127 16 8.043 14.191 107 15.687 17 14.626 24.423 118 17.389 18 5.377 9.038 102 15.016 19 10.169 18.004 114 16.792 20 8.452 14.671 116 17.073 21 4.503 7.974 113 16.639 22 10.763 19.526 131 19.183 23 19.207 33.194 114 16.816 24 10.194 17.682 118 17.296 25 5.282 8.888 112 16.470 26 7.120 12.152 111 16.367 27 10.374 18.721 100 14.764 28 7.099 12.301 113 16.530 29 11.843 20.903 118 17.395 30 11.191 18.386 116 17.027 31 6.202 10.699 104 15.346 32 8.259 14.528 130 19.150 33 11.673 20.291 114 16.731 34 9.642 16.187 113 16.555 35 23.073 38.960 105 15.444 36 20.102 33.375 111 16.238 37 47.035 79.201 106 15.607 38 12.181 21.094 104 15.305 39 9.293 15.860 108 15.866 40 9.341 16.324 114 16.778 41 4.389 7.853 109 16.071 42 11.293 19.809 111 16.360 43 13.594 23.482 105 15.369 44 6.969 12.423 95 14.026 45 1.306 2.361 94 13.780 46 11.855 20.511 102 14.993 47 9.266 16.678 108 15.892 48 4.983 8.898 100 14.749 49 21.305 35.758 108 15.835 50 6.105 10.808 110 16.087 51 16.081 26.210 108 15.796 52 4.980 8.498 103 15.190 53 8.225 14.018 103 15.153 54 10.576 19.049 98 14.458 55 5.576 9.482 105 15.360 56 10.143 18.020 106 15.584 57 8.884 15.039 104 15.207 58 8.985 15.716 105 15.432 59 5.634 9.520 106 15.552 60 4.485 8.295 101 14.851 61 19.285 33.230 106 15.557 62 6.678 11.507 108 15.858 63 7.101 12.655 112 16.401 64 7.739 14.079 95 13.941 65 5.742 10.515 98 14.427 66 6.192 10.914 95 13.912 67 24.554 39.966 103 15.194 68 4.386 7.409 103 15.146 69 5.697 10.231 120 17.576 70 13.625 23.331 109 15.960 71 55.457 98.343 75 10.980 72 16.181 29.883 107 15.669 73 11.927 20.703 99 14.554 74 79.009 147.395 76 11.237 75 21.129 42.622 81 11.836 76 19.132 34.427 86 12.651 77 27.088 48.595 87 12.749 78 13.241 21.569 97 14.227 79 41.602 78.928 88 12.897 80 26.141 44.082 93 13.602 81 10.678 18.780 82 12.059 82 42.356 82.653 65 9.496 83 23.658 45.001 69 10.126 84 11.936 23.931 55 8.100 85 59.836 114.766 72 10.518 86 39.267 76.654 110 16.101 87 13.664 24.127 113 16.628 88 25.699 48.788 102 14.940 89 19.016 38.602 114 16.719 90 6.849 12.199 118 17.379 91 16.942 29.446 116 17.066 92 12.988 23.582 120 17.652 93 1.130 2.070 112 16.435 94 22.943 39.619 123 18.074 95 2.835 4.857 99 14.546 96 4.732 8.221 116 17.085 97 5.016 8.905 119 17.517 98 20.042 34.423 127 18.613 99 20.786 35.229 111 16.257 100 3.666 6.524 114 16.746 101 5.381 10.142 120 17.593 102 6.272 11.075 125 18.415 103 5.235 8.942 120 17.629 104 6.097 13.092 119 17.446 105 4.931 8.738 127 18.693 106 7.093 12.162 110 16.149 107 2.664 4.715 123 18.125 108 10.084 17.488 117 17.217 109 7.252 12.872 102 15.044 110 10.439 18.505 135 19.792 111 8.663 15.081 116 17.012 112 7.147 12.484 119 17.443 113 12.040 24.279 123 18.024 114 2.479 4.379 122 17.955 115 6.110 10.861 116 17.069 116 9.125 17.998 135 19.777 117 17.724 31.024 107 15.756 118 6.361 11.094 126 18.557 119 8.481 14.773 120 17.659 120 21.670 38.260 104 15.314 121 8.455 14.956 125 18.294 122 6.431 13.515 119 17.543 123 15.695 29.267 111 16.322 124 12.215 21.626 129 19.020 125 16.985 28.075 112 16.459 126 5.509 9.557 112 16.453 127 3.349 5.920 104 15.257 128 5.100 9.365 123 18.042 129 5.098 9.256 128 18.799 130 6.498 11.617 129 18.907 131 4.213 7.528 125 18.364 132 13.679 22.935 108 15.869 133 3.339 5.811 102 14.936 134 2.865 5.231 111 16.293 135 7.810 13.811 118 17.378 136 10.848 19.479 131 19.293 137 3.560 6.291 119 17.547 138 5.120 9.427 122 17.955 139 4.994 9.280 121 17.780 140 6.921 12.608 142 20.906 141 4.308 7.582 105 15.361 142 10.372 18.543 125 18.361 143 8.788 14.967 112 16.471 144 55.767 92.323 118 17.289 145 4.019 7.043 111 16.299 146 7.486 13.762 131 19.288 147 5.889 10.899 138 20.228 148 8.479 15.222 120 17.668 149 13.187 22.208 107 15.654 150 9.933 20.892 120 17.612 151 5.821 10.161 113 16.607 152 19.645 32.145 110 16.143 153 7.991 13.977 115 16.870 154 4.779 8.015 106 15.506 155 3.602 6.608 117 17.222 156 21.505 37.620 120 17.653 157 5.378 9.548 113 16.574 158 3.641 6.855 119 17.439 159 5.885 11.099 126 18.449 160 5.154 9.577 112 16.437 161 7.978 14.341 114 16.700 162 6.874 12.540 126 18.500 163 1.784 3.128 107 15.673 164 2.514 4.679 111 16.298 165 3.423 6.112 114 16.739 166 7.271 12.703 107 15.671 167 3.965 7.326 123 18.075 168 7.124 14.018 139 20.448 169 3.645 6.586 117 17.116 170 15.202 24.828 115 16.931 171 4.408 7.993 107 15.695 172 16.341 30.235 103 15.135 173 4.294 7.697 110 16.126 174 3.175 5.900 105 15.462 175 5.761 10.246 104 15.242 176 11.727 21.498 124 18.159 177 13.005 22.829 97 14.226 178 5.386 9.778 113 16.536 179 3.269 6.124 121 17.751 180 14.968 29.423 116 17.095 181 18.943 32.680 114 16.706 182 8.405 14.748 104 15.282 183 12.020 18.389 97 14.314 184 70.511 117.027 108 15.823 185 6.158 10.864 102 14.918 186 7.398 13.645 108 15.843 187 21.308 34.079 133 19.477 188 12.316 22.029 106 15.514 189 11.570 19.647 101 14.889 190 11.978 22.000 100 14.699 191 1.546 2.755 107 15.663 192 9.130 15.881 90 13.183 193 5.226 9.205 87 12.722 194 4.837 8.333 92 13.448 195 6.837 12.039 91 13.358 196 1.839 3.331 81 11.891 197 4.550 8.172 83 12.135 198 20.470 34.865 95 14.030 199 4.273 7.783 88 12.938 200 560.000 990.000 78 11.519 201 11.147 19.772 90 13.197 202 2.042 3.656 80 11.769 203 10.107 17.732 92 13.462 204 6.693 12.000 91 13.413 205 7.696 14.110 102 14.966 206 5.520 9.560 98 14.371 207 6.305 11.327 101 14.797 208 15.090 26.311 99 14.607 209 43.413 73.859 103 15.168 210 7.131 12.637 99 14.599 211 3.092 5.439 99 14.556 212 18.743 32.444 91 13.417 213 1.138 2.048 92 13.573 214 20.382 35.897 103 15.112 215 17.304 30.928 99 14.517 216 13.142 23.853 102 15.022 217 9.196 15.460 102 14.965 218 7.585 12.210 108 15.926 219 6.365 11.277 96 14.100 220 7.570 13.533 92 13.539 221 11.344 18.229 96 14.117 222 45.037 69.145 100 14.636 223 5.037 8.951 98 14.354 224 5.425 9.897 97 14.217 225 5.923 10.633 97 14.182 226 15.461 26.633 101 14.775 227 5.195 9.375 94 13.745 228 4.642 8.463 97 14.254 229 6.105 10.768 93 13.704 230 32.793 56.407 102 14.965 231 6.098 11.017 88 12.996 232 5.197 9.098 96 14.140 233 4.161 8.207 93 13.727 234 6.140 11.007 92 13.545 235 4.092 7.480 94 13.881 236 3.597 6.597 94 13.777 237 2.885 5.133 95 13.903 238 10.799 19.415 100 14.672 239 4.972 9.047 94 13.876 240 7.354 13.394 93 13.654 241 2.746 4.959 83 12.205 242 6.520 10.219 95 13.968 243 13.284 21.467 130 19.027 244 6.979 11.001 98 14.375 245 6.922 11.790 93 13.589 246 47.352 78.030 111 16.324 247 10.535 17.715 112 16.442 248 11.341 19.128 117 17.130 249 18.588 31.605 105 15.498 250 10.325 17.485 115 16.868 251 9.841 16.910 112 16.509 252 9.930 17.334 114 16.678 253 21.171 36.046 109 15.940 254 5.516 9.507 109 16.023 255 14.675 24.619 116 16.998 256 8.212 14.362 107 15.734 257 8.027 13.931 114 16.803 258 25.978 43.569 107 15.730 259 13.565 23.457 98 14.362 260 6.887 11.809 115 16.899 261 10.218 17.647 107 15.704 262 5.733 10.320 118 17.368 263 14.076 23.359 109 15.985 264 6.197 11.081 110 16.216 265 11.413 20.450 97 14.246 266 7.861 13.610 104 15.291 267 12.047 19.735 101 14.832 268 3.494 6.171 106 15.632 269 12.806 22.381 99 14.592 270 3.800 6.591 99 14.513 271 7.337 12.216 102 15.027 272 10.494 18.957 108 15.910 273 5.660 10.262 139 20.393 274 16.514 28.696 108 15.806 275 9.992 17.345 120 17.691 276 18.510 32.459 110 16.126 277 6.947 12.202 113 16.593 278 145.933 236.196 105 15.391 279 4.574 7.946 99 14.576 280 11.405 20.564 115 16.838 281 5.257 9.405 118 17.358 282 6.127 10.685 116 17.048 283 12.857 22.564 121 17.845 284 3.343 5.923 108 15.888 285 6.343 11.021 113 16.618 286 6.317 11.403 108 15.841 287 7.342 13.217 114 16.701 288 4.605 8.253 152 22.269 289 4.506 7.766 103 15.095 290 4.008 6.950 109 15.959 291 4.658 8.018 105 15.367 292 8.334 14.671 102 14.979 293 8.154 13.820 109 15.964 294 1.143 2.063 109 16.059 295 3.608 6.320 101 14.833 296 4.462 8.150 108 15.813 297 3.623 6.464 110 16.226 298 3.699 6.453 107 15.782 299 17.103 28.935 109 15.996 300 13.826 24.532 87 12.736 301 3.366 6.110 111 16.249 302 3.932 6.887 108 15.901 303 4.406 7.877 111 16.343 304 25.953 45.829 111 16.337 305 8.746 15.408 107 15.747 306 21.624 38.382 101 14.822 307 10.162 18.275 105 15.374 308 41.500 70.233 102 15.004 309 9.663 17.068 103 15.127 310 15.252 26.976 105 15.430 311 16.211 27.363 100 14.706 312 7.783 13.492 90 13.274 313 6.593 11.562 86 12.680 314 1.932 3.391 94 13.775 315 3.561 6.369 96 14.067 316 3.317 5.747 101 14.888 317 1.905 3.278 100 14.654 318 6.355 11.071 95 13.995 319 5.976 10.720 85 12.529 320 8.174 14.142 83 12.240 321 117.696 201.571 78 11.433 322 20.776 35.784 78 11.482 323 17.180 29.879 85 12.512 324 6.020 11.034 70 10.226 325 12.990 22.269 87 12.840 326 9.788 16.775 84 12.384 327 6.853 12.048 114 16.795 328 5.112 8.963 106 15.607 329 12.900 22.421 84 12.274 330 6.021 10.192 115 16.905 331 9.198 16.423 102 14.971 332 6.033 10.677 102 14.946 333 11.794 20.155 79 11.629 334 11.416 20.021 74 10.941 335 9.873 16.890 83 12.123 336 12.235 20.728 86 12.630 337 3.810 6.728 84 12.331 338 2.843 5.017 95 13.905 339 4.920 9.720 105 15.369 340 2.187 4.080 97 14.254 341 53.206 91.174 89 13.011 342 10.872 18.714 78 11.447 343 4.402 7.746 102 15.029 344 3.940 6.556 83 12.215 345 13.356 22.622 93 13.643 346 9.899 17.597 82 12.015 347 30.531 53.467 84 12.336 348 11.700 20.620 102 14.933 349 5.871 10.207 99 14.506 350 7.185 12.387 109 16.084 351 44.601 77.562 83 12.134 352 4.475 7.994 102 14.971 353 10.543 18.085 94 13.874 354 4.321 8.042 110 16.202 355 14.135 25.785 96 14.079 356 6.825 11.632 82 12.102 357 3.782 6.755 87 12.798 358 19.290 32.591 87 12.713 359 5.622 9.898 79 11.559 360 5.501 9.588 82 12.085 361 4.216 7.455 93 13.627 362 2.085 3.604 81 11.893 363 7.508 13.377 114 16.791 364 3.083 5.550 101 14.790 365 2.187 4.130 83 12.129 366 2.914 5.148 75 11.016 367 11.033 18.768 85 12.556 368 2.657 4.688 84 12.363 369 8.308 14.640 102 14.943 370 4.482 7.550 88 12.947 371 4.325 7.726 88 12.937 372 3.056 5.449 91 13.430 373 5.351 9.826 90 13.216 374 7.756 13.315 93 13.721 375 1.918 3.423 95 13.997 376 2.933 5.324 92 13.548 377 9.489 16.883 82 12.087 378 2.888 5.046 91 13.368 379 40.378 68.018 94 13.788 380 7.744 13.213 92 13.543 381 2.261 4.007 104 15.341 382 1.583 2.870 104 15.319 383 2.380 4.238 96 14.168 384 3.280 5.742 86 12.691 385 2.551 4.519 96 14.045 386 2.116 3.618 93 13.695 387 2.532 4.593 100 14.646 388 12.499 20.263 93 13.722 389 2.896 5.151 95 14.021 390 2.091 3.736 103 15.122 391 3.169 5.669 118 17.390 392 1.487 2.582 97 14.321 393 1.288 2.397 102 15.047 394 1.946 3.650 105 15.422 395 3.374 6.142 99 14.491 396 7.220 12.870 101 14.792 397 16.349 27.483 94 13.826 398 4.996 8.737 105 15.358 399 6.435 11.160 95 13.886 400 4.905 8.439 97 14.201 401 6.779 11.798 91 13.363 402 7.304 12.932 106 15.569 403 11.935 20.934 123 18.061 404 3.173 5.383 88 12.898 405 3.677 6.616 108 15.890 406 7.851 13.069 113 16.556 407 14.750 24.992 89 13.136 408 9.395 16.150 94 13.751 409 12.687 21.833 86 12.615 410 22.124 37.702 82 12.102 411 4.886 8.555 107 15.728 412 112.237 189.504 85 12.439 413 5.420 9.716 127 18.611 414 13.897 23.651 96 14.067 415 13.110 22.591 79 11.556 416 36.631 61.447 81 11.932 417 8.813 15.654 97 14.234 418 7.325 13.080 107 15.725 419 4.334 7.718 97 14.213 420 9.665 16.877 92 13.573 421 2.352 4.156 99 14.498 422 2.307 4.133 97 14.318 423 8.020 14.467 75 10.952 424 9.416 16.747 99 14.597 425 4.524 8.022 107 15.685 426 1.873 3.351 90 13.158 427 3.194 5.662 91 13.312 428 6.869 11.910 92 13.571 429 2.099 3.785 118 17.363 430 5.376 9.574 90 13.211 431 5.404 9.773 88 12.967 432 6.601 10.510 96 14.036 433 3.926 6.687 91 13.386 434 1.688 3.010 101 14.849 435 6.668 11.717 102 15.051 436 2.837 5.380 104 15.237 437 1.452 2.516 90 13.255 438 31.389 54.334 82 12.019 439 3.784 6.820 89 13.004 440 5.429 9.290 96 14.109 441 2.686 5.338 83 12.250 442 2.601 5.462 79 11.634 443 1.996 3.952 82 12.084 444 2.930 5.597 82 12.070 445 9.787 18.360 98 14.437 446 5.359 10.519 84 12.291 447 2.611 5.221 93 13.664 448 5.028 10.255 97 14.230 449 4.939 9.216 93 13.589 450 1.565 2.833 104 15.294 451 3.166 5.639 102 15.010 452 1.601 2.867 110 16.195 453 1.574 2.859 111 16.338 454 1.988 3.577 110 16.218 455 1.732 3.056 104 15.233 456 1.585 2.962 105 15.403 457 8.232 14.596 105 15.455 458 1.738 3.047 104 15.335 459 2.003 3.506 105 15.366 460 2.944 5.164 109 15.965 461 3.847 6.563 91 13.437 462 8.635 14.252 104 15.286 463 1.424 2.581 94 13.821 464 4.251 7.592 99 14.599 465 23.050 41.757 92 13.570 466 10.052 17.937 105 15.437 467 35.964 64.194 89 13.023 468 4.536 7.922 103 15.160 469 5.056 8.745 103 15.070 470 5.611 9.899 101 14.853 471 43.970 71.063 112 16.460 472 6.788 11.916 105 15.395 473 17.047 31.147 93 13.713 474 8.488 14.554 96 14.071 475 7.928 13.905 104 15.311 476 3.879 6.696 99 14.577 477 5.280 9.762 97 14.269 478 23.525 38.627 103 15.102 479 9.981 17.113 107 15.785 480 11.471 20.183 103 15.061 481 3.918 6.970 101 14.859 482 6.876 12.535 89 13.140 483 8.429 14.644 97 14.210 484 10.793 19.233 91 13.345 485 7.693 13.880 90 13.157 486 6.390 11.564 95 13.934 487 16.513 30.021 89 13.091 488 6.516 12.108 89 13.021 489 17.937 32.282 95 13.927 490 13.589 24.009 96 14.072 491 6.857 12.637 95 13.908 492 9.005 16.201 96 14.170 493 7.427 13.650 95 13.993 494 9.015 15.969 97 14.200 495 6.216 11.121 105 15.470 496 17.455 30.278 98 14.327 497 6.014 10.272 101 14.778 498 3.978 6.855 98 14.356 499 56.000 84.000 117 17.213 500 5.409 9.377 100 14.699 501 4.685 8.152 104 15.302 502 14.213 24.829 101 14.801 503 20.890 36.558 86 12.655 504 8.987 16.110 96 14.140 505 17.813 29.832 101 14.853 506 2.298 4.234 93 13.651 507 4.065 7.050 102 15.047 508 14.443 26.738 119 17.495 509 2.249 4.873 122 17.969 510 8.045 15.263 97 14.211 511 883.000 1.575 95 13.964 512 3.955 7.543 118 17.334 513 7.905 14.481 88 12.997 514 1.576 3.025 85 12.454 515 1.073 2.093 103 15.123 516 2.606 4.803 82 12.043 517 2.755 4.810 87 12.748 518 1.310 2.325 85 12.544 519 2.487 4.817 87 12.844 520 4.285 7.372 89 13.091 521 6.276 10.706 86 12.697 522 1.672 2.948 86 12.667 523 2.807 5.076 83 12.150 524 2.521 4.349 82 11.993 525 1.869 3.207 89 13.102 526 9.606 17.100 88 12.948 527 2.772 5.069 88 12.984 528 1.371 2.316 83 12.213 529 1.469 2.585 90 13.150 530 4.678 8.187 88 12.989 531 3.278 5.489 85 12.427 532 778.000 1.369 86 12.563 533 893.000 1.547 87 12.718 534 2.216 4.259 93 13.676 535 2.623 4.681 92 13.554 536 5.567 10.016 95 13.897 537 3.660 6.676 90 13.240 538 2.829 5.072 84 12.327 539 3.294 5.774 90 13.250 540 1.321 2.392 89 13.019 541 1.630 2.974 83 12.255 542 2.758 5.071 91 13.297 543 2.641 5.429 111 16.308 544 3.259 5.447 87 12.801 545 4.097 7.964 108 15.891 546 1.474 2.699 92 13.533 547 2.167 4.264 97 14.205 548 1.080 2.037 86 12.631 549 1.647 3.212 107 15.668 550 1.907 3.758 91 13.332 551 5.890 11.230 91 13.326 552 4.022 6.986 90 13.196 553 4.852 8.494 84 12.410 554 1.822 3.163 81 11.949 555 8.816 15.197 93 13.641 556 7.872 13.164 88 12.950 557 2.531 4.451 84 12.384 558 3.595 6.722 96 14.121 559 3.261 5.348 84 12.326 560 2.702 4.953 90 13.197 561 2.613 4.482 90 13.232 562 1.747 3.131 90 13.190 563 6.816 12.115 85 12.550 564 2.631 4.741 87 12.772 565 1.586 2.831 76 11.231 566 4.831 8.586 101 14.830 567 14.237 24.686 89 13.061 568 3.447 6.316 110 16.146 569 8.178 14.680 106 15.568 570 3.671 6.804 104 15.311 571 4.194 7.536 103 15.095 572 5.132 9.460 92 13.511 573 12.580 22.252 110 16.091 574 3.578 6.503 102 14.943 575 10.006 18.109 94 13.848 576 4.502 8.428 113 16.629 577 6.730 12.135 93 13.637 578 64.771 107.796 100 14.749 579 2.453 4.387 93 13.602 580 6.369 11.462 104 15.350 581 15.968 27.083 89 13.096 582 4.247 7.766 100 14.720 583 2.579 4.601 84 12.411 584 7.872 13.512 80 11.793 585 1.640 2.913 84 12.307 586 6.160 10.853 91 13.376 587 4.989 8.554 88 12.870 588 3.493 5.743 79 11.667 589 9.952 17.728 98 14.428 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) inwoners welvaart `gem_ink\r` 36.0532 0.5335 14.1896 -98.6667 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -23.1764 -8.5935 -4.2577 -0.1482 876.4713 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 36.05322 19.25165 1.873 0.0616 . inwoners 0.53350 0.05096 10.469 <2e-16 *** welvaart 14.18958 8.57765 1.654 0.0986 . `gem_ink\r` -98.66670 58.40897 -1.689 0.0917 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 60.57 on 585 degrees of freedom Multiple R-squared: 0.1691, Adjusted R-squared: 0.1649 F-statistic: 39.69 on 3 and 585 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,] 7.532592e-07 1.506518e-06 9.999992e-01 [2,] 5.890917e-09 1.178183e-08 1.000000e+00 [3,] 5.024446e-11 1.004889e-10 1.000000e+00 [4,] 3.176764e-13 6.353528e-13 1.000000e+00 [5,] 3.507894e-15 7.015788e-15 1.000000e+00 [6,] 3.329559e-17 6.659119e-17 1.000000e+00 [7,] 3.272509e-19 6.545019e-19 1.000000e+00 [8,] 2.396383e-21 4.792766e-21 1.000000e+00 [9,] 1.844039e-23 3.688077e-23 1.000000e+00 [10,] 1.091026e-25 2.182053e-25 1.000000e+00 [11,] 1.393604e-27 2.787208e-27 1.000000e+00 [12,] 1.151508e-29 2.303015e-29 1.000000e+00 [13,] 7.972160e-32 1.594432e-31 1.000000e+00 [14,] 4.505685e-34 9.011371e-34 1.000000e+00 [15,] 2.474656e-36 4.949311e-36 1.000000e+00 [16,] 1.594931e-38 3.189861e-38 1.000000e+00 [17,] 1.138653e-40 2.277306e-40 1.000000e+00 [18,] 7.077629e-43 1.415526e-42 1.000000e+00 [19,] 6.514484e-45 1.302897e-44 1.000000e+00 [20,] 3.607511e-47 7.215022e-47 1.000000e+00 [21,] 7.714890e-49 1.542978e-48 1.000000e+00 [22,] 5.811815e-51 1.162363e-50 1.000000e+00 [23,] 3.922319e-53 7.844638e-53 1.000000e+00 [24,] 1.508931e-54 3.017862e-54 1.000000e+00 [25,] 8.693993e-57 1.738799e-56 1.000000e+00 [26,] 4.620291e-59 9.240581e-59 1.000000e+00 [27,] 2.409734e-61 4.819468e-61 1.000000e+00 [28,] 2.722274e-63 5.444549e-63 1.000000e+00 [29,] 1.823665e-65 3.647329e-65 1.000000e+00 [30,] 5.189291e-67 1.037858e-66 1.000000e+00 [31,] 3.946859e-69 7.893718e-69 1.000000e+00 [32,] 2.174696e-71 4.349392e-71 1.000000e+00 [33,] 1.256644e-73 2.513288e-73 1.000000e+00 [34,] 6.559737e-76 1.311947e-75 1.000000e+00 [35,] 3.340734e-78 6.681467e-78 1.000000e+00 [36,] 1.955718e-80 3.911436e-80 1.000000e+00 [37,] 1.013899e-82 2.027798e-82 1.000000e+00 [38,] 5.903572e-85 1.180714e-84 1.000000e+00 [39,] 3.088687e-87 6.177374e-87 1.000000e+00 [40,] 1.525423e-89 3.050846e-89 1.000000e+00 [41,] 1.371916e-91 2.743832e-91 1.000000e+00 [42,] 6.574605e-94 1.314921e-93 1.000000e+00 [43,] 6.119673e-96 1.223935e-95 1.000000e+00 [44,] 2.933361e-98 5.866723e-98 1.000000e+00 [45,] 2.624454e-99 5.248909e-99 1.000000e+00 [46,] 1.627924e-101 3.255849e-101 1.000000e+00 [47,] 8.840465e-104 1.768093e-103 1.000000e+00 [48,] 1.353018e-105 2.706036e-105 1.000000e+00 [49,] 7.196931e-108 1.439386e-107 1.000000e+00 [50,] 5.379970e-110 1.075994e-109 1.000000e+00 [51,] 2.817941e-112 5.635882e-112 1.000000e+00 [52,] 1.413753e-114 2.827507e-114 1.000000e+00 [53,] 8.518469e-117 1.703694e-116 1.000000e+00 [54,] 4.966743e-119 9.933486e-119 1.000000e+00 [55,] 2.933461e-121 5.866923e-121 1.000000e+00 [56,] 1.458059e-123 2.916118e-123 1.000000e+00 [57,] 8.235522e-126 1.647104e-125 1.000000e+00 [58,] 9.771557e-128 1.954311e-127 1.000000e+00 [59,] 5.899401e-130 1.179880e-129 1.000000e+00 [60,] 2.843758e-132 5.687516e-132 1.000000e+00 [61,] 9.617521e-132 1.923504e-131 1.000000e+00 [62,] 6.420417e-134 1.284083e-133 1.000000e+00 [63,] 3.355075e-136 6.710151e-136 1.000000e+00 [64,] 1.677662e-138 3.355324e-138 1.000000e+00 [65,] 3.592298e-131 7.184596e-131 1.000000e+00 [66,] 7.834344e-132 1.566869e-131 1.000000e+00 [67,] 5.619149e-134 1.123830e-133 1.000000e+00 [68,] 1.430916e-113 2.861831e-113 1.000000e+00 [69,] 9.446323e-114 1.889265e-113 1.000000e+00 [70,] 1.482647e-115 2.965293e-115 1.000000e+00 [71,] 2.408910e-117 4.817819e-117 1.000000e+00 [72,] 9.170315e-119 1.834063e-118 1.000000e+00 [73,] 8.304581e-118 1.660916e-117 1.000000e+00 [74,] 2.545238e-119 5.090476e-119 1.000000e+00 [75,] 7.755839e-121 1.551168e-120 1.000000e+00 [76,] 4.592259e-120 9.184518e-120 1.000000e+00 [77,] 8.350797e-122 1.670159e-121 1.000000e+00 [78,] 2.474397e-123 4.948795e-123 1.000000e+00 [79,] 4.370288e-121 8.740577e-121 1.000000e+00 [80,] 2.161561e-119 4.323122e-119 1.000000e+00 [81,] 4.845713e-121 9.691426e-121 1.000000e+00 [82,] 2.765536e-122 5.531073e-122 1.000000e+00 [83,] 9.592969e-123 1.918594e-122 1.000000e+00 [84,] 2.255732e-124 4.511465e-124 1.000000e+00 [85,] 5.075163e-126 1.015033e-125 1.000000e+00 [86,] 1.458328e-127 2.916656e-127 1.000000e+00 [87,] 3.220671e-129 6.441341e-129 1.000000e+00 [88,] 7.555226e-131 1.511045e-130 1.000000e+00 [89,] 2.061833e-132 4.123666e-132 1.000000e+00 [90,] 4.409165e-134 8.818329e-134 1.000000e+00 [91,] 9.699023e-136 1.939805e-135 1.000000e+00 [92,] 2.134464e-137 4.268928e-137 1.000000e+00 [93,] 5.053205e-139 1.010641e-138 1.000000e+00 [94,] 1.048607e-140 2.097215e-140 1.000000e+00 [95,] 2.319212e-142 4.638423e-142 1.000000e+00 [96,] 5.599683e-144 1.119937e-143 1.000000e+00 [97,] 1.144726e-145 2.289453e-145 1.000000e+00 [98,] 4.076143e-147 8.152286e-147 1.000000e+00 [99,] 9.591187e-149 1.918237e-148 1.000000e+00 [100,] 2.056822e-150 4.113644e-150 1.000000e+00 [101,] 4.443286e-152 8.886573e-152 1.000000e+00 [102,] 8.949297e-154 1.789859e-153 1.000000e+00 [103,] 1.832139e-155 3.664277e-155 1.000000e+00 [104,] 4.782063e-157 9.564127e-157 1.000000e+00 [105,] 9.464955e-159 1.892991e-158 1.000000e+00 [106,] 1.854228e-160 3.708456e-160 1.000000e+00 [107,] 1.309545e-161 2.619091e-161 1.000000e+00 [108,] 2.634743e-163 5.269486e-163 1.000000e+00 [109,] 5.176985e-165 1.035397e-164 1.000000e+00 [110,] 2.480412e-166 4.960824e-166 1.000000e+00 [111,] 4.820209e-168 9.640419e-168 1.000000e+00 [112,] 1.024893e-169 2.049786e-169 1.000000e+00 [113,] 2.008503e-171 4.017005e-171 1.000000e+00 [114,] 3.872325e-173 7.744651e-173 1.000000e+00 [115,] 7.392153e-175 1.478431e-174 1.000000e+00 [116,] 3.138779e-176 6.277557e-176 1.000000e+00 [117,] 8.185365e-178 1.637073e-177 1.000000e+00 [118,] 2.131468e-179 4.262935e-179 1.000000e+00 [119,] 5.333622e-181 1.066724e-180 1.000000e+00 [120,] 1.046753e-182 2.093506e-182 1.000000e+00 [121,] 2.346516e-184 4.693033e-184 1.000000e+00 [122,] 4.478952e-186 8.957905e-186 1.000000e+00 [123,] 9.010262e-188 1.802052e-187 1.000000e+00 [124,] 1.731615e-189 3.463230e-189 1.000000e+00 [125,] 3.258078e-191 6.516157e-191 1.000000e+00 [126,] 7.483462e-193 1.496692e-192 1.000000e+00 [127,] 1.887909e-194 3.775817e-194 1.000000e+00 [128,] 3.573356e-196 7.146713e-196 1.000000e+00 [129,] 6.619545e-198 1.323909e-197 1.000000e+00 [130,] 1.727083e-199 3.454166e-199 1.000000e+00 [131,] 3.157772e-201 6.315543e-201 1.000000e+00 [132,] 6.065895e-203 1.213179e-202 1.000000e+00 [133,] 1.125120e-204 2.250240e-204 1.000000e+00 [134,] 3.550548e-206 7.101095e-206 1.000000e+00 [135,] 7.945464e-208 1.589093e-207 1.000000e+00 [136,] 1.524889e-209 3.049779e-209 1.000000e+00 [137,] 2.858382e-211 5.716764e-211 1.000000e+00 [138,] 1.995230e-212 3.990460e-212 1.000000e+00 [139,] 3.728681e-214 7.457362e-214 1.000000e+00 [140,] 8.988442e-216 1.797688e-215 1.000000e+00 [141,] 2.009533e-217 4.019066e-217 1.000000e+00 [142,] 3.694260e-219 7.388521e-219 1.000000e+00 [143,] 9.143334e-221 1.828667e-220 1.000000e+00 [144,] 6.659445e-222 1.331889e-221 1.000000e+00 [145,] 1.196161e-223 2.392322e-223 1.000000e+00 [146,] 3.995911e-225 7.991822e-225 1.000000e+00 [147,] 7.071704e-227 1.414341e-226 1.000000e+00 [148,] 1.665598e-228 3.331195e-228 1.000000e+00 [149,] 2.882595e-230 5.765191e-230 1.000000e+00 [150,] 5.280786e-232 1.056157e-231 1.000000e+00 [151,] 9.289824e-234 1.857965e-233 1.000000e+00 [152,] 1.584229e-235 3.168459e-235 1.000000e+00 [153,] 2.914244e-237 5.828488e-237 1.000000e+00 [154,] 4.930314e-239 9.860628e-239 1.000000e+00 [155,] 8.333601e-241 1.666720e-240 1.000000e+00 [156,] 1.509395e-242 3.018791e-242 1.000000e+00 [157,] 3.034022e-244 6.068045e-244 1.000000e+00 [158,] 5.206966e-246 1.041393e-245 1.000000e+00 [159,] 8.831528e-248 1.766306e-247 1.000000e+00 [160,] 1.652853e-249 3.305706e-249 1.000000e+00 [161,] 2.789217e-251 5.578433e-251 1.000000e+00 [162,] 1.177172e-252 2.354344e-252 1.000000e+00 [163,] 1.982322e-254 3.964644e-254 1.000000e+00 [164,] 4.447241e-256 8.894481e-256 1.000000e+00 [165,] 7.717720e-258 1.543544e-257 1.000000e+00 [166,] 1.441092e-259 2.882185e-259 1.000000e+00 [167,] 2.471231e-261 4.942461e-261 1.000000e+00 [168,] 4.141893e-263 8.283786e-263 1.000000e+00 [169,] 7.517069e-265 1.503414e-264 1.000000e+00 [170,] 1.422476e-266 2.844952e-266 1.000000e+00 [171,] 2.653411e-268 5.306823e-268 1.000000e+00 [172,] 4.320318e-270 8.640636e-270 1.000000e+00 [173,] 6.867068e-272 1.373414e-271 1.000000e+00 [174,] 5.316945e-273 1.063389e-272 1.000000e+00 [175,] 8.691874e-275 1.738375e-274 1.000000e+00 [176,] 1.497243e-276 2.994486e-276 1.000000e+00 [177,] 1.121238e-277 2.242477e-277 1.000000e+00 [178,] 2.183588e-278 4.367176e-278 1.000000e+00 [179,] 4.190117e-280 8.380234e-280 1.000000e+00 [180,] 6.678780e-282 1.335756e-281 1.000000e+00 [181,] 1.813852e-283 3.627704e-283 1.000000e+00 [182,] 2.908619e-285 5.817238e-285 1.000000e+00 [183,] 5.661574e-287 1.132315e-286 1.000000e+00 [184,] 8.955902e-289 1.791180e-288 1.000000e+00 [185,] 1.637205e-290 3.274411e-290 1.000000e+00 [186,] 3.810996e-292 7.621993e-292 1.000000e+00 [187,] 9.737303e-294 1.947461e-293 1.000000e+00 [188,] 2.376261e-295 4.752522e-295 1.000000e+00 [189,] 4.823541e-297 9.647083e-297 1.000000e+00 [190,] 1.301821e-298 2.603643e-298 1.000000e+00 [191,] 3.115047e-300 6.230094e-300 1.000000e+00 [192,] 6.191167e-302 1.238233e-301 1.000000e+00 [193,] 1.176348e-303 2.352695e-303 1.000000e+00 [194,] 2.774333e-294 5.548666e-294 1.000000e+00 [195,] 5.121365e-296 1.024273e-295 1.000000e+00 [196,] 1.022060e-297 2.044120e-297 1.000000e+00 [197,] 1.924537e-299 3.849073e-299 1.000000e+00 [198,] 3.423319e-301 6.846639e-301 1.000000e+00 [199,] 5.935661e-303 1.187132e-302 1.000000e+00 [200,] 1.080726e-304 2.161451e-304 1.000000e+00 [201,] 1.875877e-306 3.751755e-306 1.000000e+00 [202,] 3.337427e-308 6.674855e-308 1.000000e+00 [203,] 1.120583e-309 2.241167e-309 1.000000e+00 [204,] 1.924300e-311 3.848599e-311 1.000000e+00 [205,] 3.341554e-313 6.683109e-313 1.000000e+00 [206,] 6.386088e-315 1.277218e-314 1.000000e+00 [207,] 1.112356e-316 2.224712e-316 1.000000e+00 [208,] 1.874925e-318 3.749850e-318 1.000000e+00 [209,] 3.140281e-320 6.280562e-320 1.000000e+00 [210,] 5.335909e-322 1.067182e-321 1.000000e+00 [211,] 9.881313e-324 1.976263e-323 1.000000e+00 [212,] 0.000000e+00 0.000000e+00 1.000000e+00 [213,] 0.000000e+00 0.000000e+00 1.000000e+00 [214,] 0.000000e+00 0.000000e+00 1.000000e+00 [215,] 0.000000e+00 0.000000e+00 1.000000e+00 [216,] 0.000000e+00 0.000000e+00 1.000000e+00 [217,] 0.000000e+00 0.000000e+00 1.000000e+00 [218,] 0.000000e+00 0.000000e+00 1.000000e+00 [219,] 0.000000e+00 0.000000e+00 1.000000e+00 [220,] 0.000000e+00 0.000000e+00 1.000000e+00 [221,] 0.000000e+00 0.000000e+00 1.000000e+00 [222,] 0.000000e+00 0.000000e+00 1.000000e+00 [223,] 0.000000e+00 0.000000e+00 1.000000e+00 [224,] 0.000000e+00 0.000000e+00 1.000000e+00 [225,] 0.000000e+00 0.000000e+00 1.000000e+00 [226,] 0.000000e+00 0.000000e+00 1.000000e+00 [227,] 0.000000e+00 0.000000e+00 1.000000e+00 [228,] 0.000000e+00 0.000000e+00 1.000000e+00 [229,] 0.000000e+00 0.000000e+00 1.000000e+00 [230,] 0.000000e+00 0.000000e+00 1.000000e+00 [231,] 0.000000e+00 0.000000e+00 1.000000e+00 [232,] 0.000000e+00 0.000000e+00 1.000000e+00 [233,] 0.000000e+00 0.000000e+00 1.000000e+00 [234,] 0.000000e+00 0.000000e+00 1.000000e+00 [235,] 0.000000e+00 0.000000e+00 1.000000e+00 [236,] 0.000000e+00 0.000000e+00 1.000000e+00 [237,] 0.000000e+00 0.000000e+00 1.000000e+00 [238,] 0.000000e+00 0.000000e+00 1.000000e+00 [239,] 0.000000e+00 0.000000e+00 1.000000e+00 [240,] 0.000000e+00 0.000000e+00 1.000000e+00 [241,] 0.000000e+00 0.000000e+00 1.000000e+00 [242,] 0.000000e+00 0.000000e+00 1.000000e+00 [243,] 0.000000e+00 0.000000e+00 1.000000e+00 [244,] 0.000000e+00 0.000000e+00 1.000000e+00 [245,] 0.000000e+00 0.000000e+00 1.000000e+00 [246,] 0.000000e+00 0.000000e+00 1.000000e+00 [247,] 0.000000e+00 0.000000e+00 1.000000e+00 [248,] 0.000000e+00 0.000000e+00 1.000000e+00 [249,] 0.000000e+00 0.000000e+00 1.000000e+00 [250,] 0.000000e+00 0.000000e+00 1.000000e+00 [251,] 0.000000e+00 0.000000e+00 1.000000e+00 [252,] 0.000000e+00 0.000000e+00 1.000000e+00 [253,] 0.000000e+00 0.000000e+00 1.000000e+00 [254,] 0.000000e+00 0.000000e+00 1.000000e+00 [255,] 0.000000e+00 0.000000e+00 1.000000e+00 [256,] 0.000000e+00 0.000000e+00 1.000000e+00 [257,] 0.000000e+00 0.000000e+00 1.000000e+00 [258,] 0.000000e+00 0.000000e+00 1.000000e+00 [259,] 0.000000e+00 0.000000e+00 1.000000e+00 [260,] 0.000000e+00 0.000000e+00 1.000000e+00 [261,] 0.000000e+00 0.000000e+00 1.000000e+00 [262,] 0.000000e+00 0.000000e+00 1.000000e+00 [263,] 0.000000e+00 0.000000e+00 1.000000e+00 [264,] 0.000000e+00 0.000000e+00 1.000000e+00 [265,] 0.000000e+00 0.000000e+00 1.000000e+00 [266,] 0.000000e+00 0.000000e+00 1.000000e+00 [267,] 0.000000e+00 0.000000e+00 1.000000e+00 [268,] 0.000000e+00 0.000000e+00 1.000000e+00 [269,] 0.000000e+00 0.000000e+00 1.000000e+00 [270,] 0.000000e+00 0.000000e+00 1.000000e+00 [271,] 0.000000e+00 0.000000e+00 1.000000e+00 [272,] 0.000000e+00 0.000000e+00 1.000000e+00 [273,] 0.000000e+00 0.000000e+00 1.000000e+00 [274,] 0.000000e+00 0.000000e+00 1.000000e+00 [275,] 0.000000e+00 0.000000e+00 1.000000e+00 [276,] 0.000000e+00 0.000000e+00 1.000000e+00 [277,] 0.000000e+00 0.000000e+00 1.000000e+00 [278,] 0.000000e+00 0.000000e+00 1.000000e+00 [279,] 0.000000e+00 0.000000e+00 1.000000e+00 [280,] 0.000000e+00 0.000000e+00 1.000000e+00 [281,] 0.000000e+00 0.000000e+00 1.000000e+00 [282,] 0.000000e+00 0.000000e+00 1.000000e+00 [283,] 0.000000e+00 0.000000e+00 1.000000e+00 [284,] 0.000000e+00 0.000000e+00 1.000000e+00 [285,] 0.000000e+00 0.000000e+00 1.000000e+00 [286,] 0.000000e+00 0.000000e+00 1.000000e+00 [287,] 0.000000e+00 0.000000e+00 1.000000e+00 [288,] 0.000000e+00 0.000000e+00 1.000000e+00 [289,] 0.000000e+00 0.000000e+00 1.000000e+00 [290,] 0.000000e+00 0.000000e+00 1.000000e+00 [291,] 0.000000e+00 0.000000e+00 1.000000e+00 [292,] 0.000000e+00 0.000000e+00 1.000000e+00 [293,] 0.000000e+00 0.000000e+00 1.000000e+00 [294,] 0.000000e+00 0.000000e+00 1.000000e+00 [295,] 0.000000e+00 0.000000e+00 1.000000e+00 [296,] 0.000000e+00 0.000000e+00 1.000000e+00 [297,] 0.000000e+00 0.000000e+00 1.000000e+00 [298,] 0.000000e+00 0.000000e+00 1.000000e+00 [299,] 0.000000e+00 0.000000e+00 1.000000e+00 [300,] 0.000000e+00 0.000000e+00 1.000000e+00 [301,] 0.000000e+00 0.000000e+00 1.000000e+00 [302,] 0.000000e+00 0.000000e+00 1.000000e+00 [303,] 0.000000e+00 0.000000e+00 1.000000e+00 [304,] 0.000000e+00 0.000000e+00 1.000000e+00 [305,] 0.000000e+00 0.000000e+00 1.000000e+00 [306,] 0.000000e+00 0.000000e+00 1.000000e+00 [307,] 0.000000e+00 0.000000e+00 1.000000e+00 [308,] 0.000000e+00 0.000000e+00 1.000000e+00 [309,] 0.000000e+00 0.000000e+00 1.000000e+00 [310,] 0.000000e+00 0.000000e+00 1.000000e+00 [311,] 0.000000e+00 0.000000e+00 1.000000e+00 [312,] 0.000000e+00 0.000000e+00 1.000000e+00 [313,] 0.000000e+00 0.000000e+00 1.000000e+00 [314,] 0.000000e+00 0.000000e+00 1.000000e+00 [315,] 0.000000e+00 0.000000e+00 1.000000e+00 [316,] 0.000000e+00 0.000000e+00 1.000000e+00 [317,] 0.000000e+00 0.000000e+00 1.000000e+00 [318,] 0.000000e+00 0.000000e+00 1.000000e+00 [319,] 0.000000e+00 0.000000e+00 1.000000e+00 [320,] 0.000000e+00 0.000000e+00 1.000000e+00 [321,] 0.000000e+00 0.000000e+00 1.000000e+00 [322,] 0.000000e+00 0.000000e+00 1.000000e+00 [323,] 0.000000e+00 0.000000e+00 1.000000e+00 [324,] 0.000000e+00 0.000000e+00 1.000000e+00 [325,] 0.000000e+00 0.000000e+00 1.000000e+00 [326,] 0.000000e+00 0.000000e+00 1.000000e+00 [327,] 0.000000e+00 0.000000e+00 1.000000e+00 [328,] 0.000000e+00 0.000000e+00 1.000000e+00 [329,] 0.000000e+00 0.000000e+00 1.000000e+00 [330,] 0.000000e+00 0.000000e+00 1.000000e+00 [331,] 0.000000e+00 0.000000e+00 1.000000e+00 [332,] 0.000000e+00 0.000000e+00 1.000000e+00 [333,] 0.000000e+00 0.000000e+00 1.000000e+00 [334,] 0.000000e+00 0.000000e+00 1.000000e+00 [335,] 0.000000e+00 0.000000e+00 1.000000e+00 [336,] 0.000000e+00 0.000000e+00 1.000000e+00 [337,] 0.000000e+00 0.000000e+00 1.000000e+00 [338,] 0.000000e+00 0.000000e+00 1.000000e+00 [339,] 0.000000e+00 0.000000e+00 1.000000e+00 [340,] 0.000000e+00 0.000000e+00 1.000000e+00 [341,] 0.000000e+00 0.000000e+00 1.000000e+00 [342,] 0.000000e+00 0.000000e+00 1.000000e+00 [343,] 0.000000e+00 0.000000e+00 1.000000e+00 [344,] 0.000000e+00 0.000000e+00 1.000000e+00 [345,] 0.000000e+00 0.000000e+00 1.000000e+00 [346,] 0.000000e+00 0.000000e+00 1.000000e+00 [347,] 0.000000e+00 0.000000e+00 1.000000e+00 [348,] 0.000000e+00 0.000000e+00 1.000000e+00 [349,] 0.000000e+00 0.000000e+00 1.000000e+00 [350,] 0.000000e+00 0.000000e+00 1.000000e+00 [351,] 0.000000e+00 0.000000e+00 1.000000e+00 [352,] 0.000000e+00 0.000000e+00 1.000000e+00 [353,] 0.000000e+00 0.000000e+00 1.000000e+00 [354,] 0.000000e+00 0.000000e+00 1.000000e+00 [355,] 0.000000e+00 0.000000e+00 1.000000e+00 [356,] 0.000000e+00 0.000000e+00 1.000000e+00 [357,] 0.000000e+00 0.000000e+00 1.000000e+00 [358,] 0.000000e+00 0.000000e+00 1.000000e+00 [359,] 0.000000e+00 0.000000e+00 1.000000e+00 [360,] 0.000000e+00 0.000000e+00 1.000000e+00 [361,] 0.000000e+00 0.000000e+00 1.000000e+00 [362,] 0.000000e+00 0.000000e+00 1.000000e+00 [363,] 0.000000e+00 0.000000e+00 1.000000e+00 [364,] 0.000000e+00 0.000000e+00 1.000000e+00 [365,] 0.000000e+00 0.000000e+00 1.000000e+00 [366,] 0.000000e+00 0.000000e+00 1.000000e+00 [367,] 0.000000e+00 0.000000e+00 1.000000e+00 [368,] 0.000000e+00 0.000000e+00 1.000000e+00 [369,] 0.000000e+00 0.000000e+00 1.000000e+00 [370,] 0.000000e+00 0.000000e+00 1.000000e+00 [371,] 0.000000e+00 0.000000e+00 1.000000e+00 [372,] 0.000000e+00 0.000000e+00 1.000000e+00 [373,] 0.000000e+00 0.000000e+00 1.000000e+00 [374,] 0.000000e+00 0.000000e+00 1.000000e+00 [375,] 0.000000e+00 0.000000e+00 1.000000e+00 [376,] 0.000000e+00 0.000000e+00 1.000000e+00 [377,] 0.000000e+00 0.000000e+00 1.000000e+00 [378,] 0.000000e+00 0.000000e+00 1.000000e+00 [379,] 0.000000e+00 0.000000e+00 1.000000e+00 [380,] 0.000000e+00 0.000000e+00 1.000000e+00 [381,] 0.000000e+00 0.000000e+00 1.000000e+00 [382,] 0.000000e+00 0.000000e+00 1.000000e+00 [383,] 0.000000e+00 0.000000e+00 1.000000e+00 [384,] 0.000000e+00 0.000000e+00 1.000000e+00 [385,] 0.000000e+00 0.000000e+00 1.000000e+00 [386,] 0.000000e+00 0.000000e+00 1.000000e+00 [387,] 0.000000e+00 0.000000e+00 1.000000e+00 [388,] 0.000000e+00 0.000000e+00 1.000000e+00 [389,] 0.000000e+00 0.000000e+00 1.000000e+00 [390,] 0.000000e+00 0.000000e+00 1.000000e+00 [391,] 0.000000e+00 0.000000e+00 1.000000e+00 [392,] 0.000000e+00 0.000000e+00 1.000000e+00 [393,] 0.000000e+00 0.000000e+00 1.000000e+00 [394,] 0.000000e+00 0.000000e+00 1.000000e+00 [395,] 0.000000e+00 0.000000e+00 1.000000e+00 [396,] 0.000000e+00 0.000000e+00 1.000000e+00 [397,] 0.000000e+00 0.000000e+00 1.000000e+00 [398,] 0.000000e+00 0.000000e+00 1.000000e+00 [399,] 0.000000e+00 0.000000e+00 1.000000e+00 [400,] 0.000000e+00 0.000000e+00 1.000000e+00 [401,] 0.000000e+00 0.000000e+00 1.000000e+00 [402,] 0.000000e+00 0.000000e+00 1.000000e+00 [403,] 0.000000e+00 0.000000e+00 1.000000e+00 [404,] 0.000000e+00 0.000000e+00 1.000000e+00 [405,] 0.000000e+00 0.000000e+00 1.000000e+00 [406,] 0.000000e+00 0.000000e+00 1.000000e+00 [407,] 0.000000e+00 0.000000e+00 1.000000e+00 [408,] 0.000000e+00 0.000000e+00 1.000000e+00 [409,] 0.000000e+00 0.000000e+00 1.000000e+00 [410,] 0.000000e+00 0.000000e+00 1.000000e+00 [411,] 0.000000e+00 0.000000e+00 1.000000e+00 [412,] 0.000000e+00 0.000000e+00 1.000000e+00 [413,] 0.000000e+00 0.000000e+00 1.000000e+00 [414,] 0.000000e+00 0.000000e+00 1.000000e+00 [415,] 0.000000e+00 0.000000e+00 1.000000e+00 [416,] 0.000000e+00 0.000000e+00 1.000000e+00 [417,] 0.000000e+00 0.000000e+00 1.000000e+00 [418,] 0.000000e+00 0.000000e+00 1.000000e+00 [419,] 0.000000e+00 0.000000e+00 1.000000e+00 [420,] 0.000000e+00 0.000000e+00 1.000000e+00 [421,] 0.000000e+00 0.000000e+00 1.000000e+00 [422,] 0.000000e+00 0.000000e+00 1.000000e+00 [423,] 0.000000e+00 0.000000e+00 1.000000e+00 [424,] 0.000000e+00 0.000000e+00 1.000000e+00 [425,] 0.000000e+00 0.000000e+00 1.000000e+00 [426,] 0.000000e+00 0.000000e+00 1.000000e+00 [427,] 0.000000e+00 0.000000e+00 1.000000e+00 [428,] 0.000000e+00 0.000000e+00 1.000000e+00 [429,] 0.000000e+00 0.000000e+00 1.000000e+00 [430,] 0.000000e+00 0.000000e+00 1.000000e+00 [431,] 0.000000e+00 0.000000e+00 1.000000e+00 [432,] 0.000000e+00 0.000000e+00 1.000000e+00 [433,] 0.000000e+00 0.000000e+00 1.000000e+00 [434,] 0.000000e+00 0.000000e+00 1.000000e+00 [435,] 0.000000e+00 0.000000e+00 1.000000e+00 [436,] 0.000000e+00 0.000000e+00 1.000000e+00 [437,] 0.000000e+00 0.000000e+00 1.000000e+00 [438,] 0.000000e+00 0.000000e+00 1.000000e+00 [439,] 0.000000e+00 0.000000e+00 1.000000e+00 [440,] 0.000000e+00 0.000000e+00 1.000000e+00 [441,] 0.000000e+00 0.000000e+00 1.000000e+00 [442,] 0.000000e+00 0.000000e+00 1.000000e+00 [443,] 0.000000e+00 0.000000e+00 1.000000e+00 [444,] 0.000000e+00 0.000000e+00 1.000000e+00 [445,] 0.000000e+00 0.000000e+00 1.000000e+00 [446,] 0.000000e+00 0.000000e+00 1.000000e+00 [447,] 0.000000e+00 0.000000e+00 1.000000e+00 [448,] 0.000000e+00 0.000000e+00 1.000000e+00 [449,] 0.000000e+00 0.000000e+00 1.000000e+00 [450,] 0.000000e+00 0.000000e+00 1.000000e+00 [451,] 0.000000e+00 0.000000e+00 1.000000e+00 [452,] 0.000000e+00 0.000000e+00 1.000000e+00 [453,] 0.000000e+00 0.000000e+00 1.000000e+00 [454,] 0.000000e+00 0.000000e+00 1.000000e+00 [455,] 0.000000e+00 0.000000e+00 1.000000e+00 [456,] 0.000000e+00 0.000000e+00 1.000000e+00 [457,] 0.000000e+00 0.000000e+00 1.000000e+00 [458,] 0.000000e+00 0.000000e+00 1.000000e+00 [459,] 0.000000e+00 0.000000e+00 1.000000e+00 [460,] 0.000000e+00 0.000000e+00 1.000000e+00 [461,] 0.000000e+00 0.000000e+00 1.000000e+00 [462,] 0.000000e+00 0.000000e+00 1.000000e+00 [463,] 0.000000e+00 0.000000e+00 1.000000e+00 [464,] 0.000000e+00 0.000000e+00 1.000000e+00 [465,] 0.000000e+00 0.000000e+00 1.000000e+00 [466,] 0.000000e+00 0.000000e+00 1.000000e+00 [467,] 0.000000e+00 0.000000e+00 1.000000e+00 [468,] 0.000000e+00 0.000000e+00 1.000000e+00 [469,] 0.000000e+00 0.000000e+00 1.000000e+00 [470,] 0.000000e+00 0.000000e+00 1.000000e+00 [471,] 0.000000e+00 0.000000e+00 1.000000e+00 [472,] 0.000000e+00 0.000000e+00 1.000000e+00 [473,] 0.000000e+00 0.000000e+00 1.000000e+00 [474,] 0.000000e+00 0.000000e+00 1.000000e+00 [475,] 0.000000e+00 0.000000e+00 1.000000e+00 [476,] 0.000000e+00 0.000000e+00 1.000000e+00 [477,] 0.000000e+00 0.000000e+00 1.000000e+00 [478,] 0.000000e+00 0.000000e+00 1.000000e+00 [479,] 0.000000e+00 0.000000e+00 1.000000e+00 [480,] 0.000000e+00 0.000000e+00 1.000000e+00 [481,] 0.000000e+00 0.000000e+00 1.000000e+00 [482,] 0.000000e+00 0.000000e+00 1.000000e+00 [483,] 0.000000e+00 0.000000e+00 1.000000e+00 [484,] 0.000000e+00 0.000000e+00 1.000000e+00 [485,] 0.000000e+00 0.000000e+00 1.000000e+00 [486,] 0.000000e+00 0.000000e+00 1.000000e+00 [487,] 0.000000e+00 0.000000e+00 1.000000e+00 [488,] 0.000000e+00 0.000000e+00 1.000000e+00 [489,] 0.000000e+00 0.000000e+00 1.000000e+00 [490,] 0.000000e+00 0.000000e+00 1.000000e+00 [491,] 0.000000e+00 0.000000e+00 1.000000e+00 [492,] 0.000000e+00 0.000000e+00 1.000000e+00 [493,] 0.000000e+00 0.000000e+00 1.000000e+00 [494,] 0.000000e+00 0.000000e+00 1.000000e+00 [495,] 0.000000e+00 0.000000e+00 1.000000e+00 [496,] 0.000000e+00 0.000000e+00 1.000000e+00 [497,] 0.000000e+00 0.000000e+00 1.000000e+00 [498,] 0.000000e+00 0.000000e+00 1.000000e+00 [499,] 0.000000e+00 0.000000e+00 1.000000e+00 [500,] 0.000000e+00 0.000000e+00 1.000000e+00 [501,] 0.000000e+00 0.000000e+00 1.000000e+00 [502,] 0.000000e+00 0.000000e+00 1.000000e+00 [503,] 0.000000e+00 0.000000e+00 1.000000e+00 [504,] 0.000000e+00 0.000000e+00 1.000000e+00 [505,] 2.214401e-67 4.428802e-67 1.000000e+00 [506,] 7.465581e-68 1.493116e-67 1.000000e+00 [507,] 2.473496e-68 4.946991e-68 1.000000e+00 [508,] 1.290618e-68 2.581236e-68 1.000000e+00 [509,] 3.785678e-69 7.571355e-69 1.000000e+00 [510,] 1.499277e-69 2.998555e-69 1.000000e+00 [511,] 7.503248e-70 1.500650e-69 1.000000e+00 [512,] 2.179510e-70 4.359021e-70 1.000000e+00 [513,] 6.586168e-71 1.317234e-70 1.000000e+00 [514,] 1.913064e-71 3.826127e-71 1.000000e+00 [515,] 5.654130e-72 1.130826e-71 1.000000e+00 [516,] 1.575255e-72 3.150510e-72 1.000000e+00 [517,] 1.211956e-72 2.423912e-72 1.000000e+00 [518,] 1.517707e-72 3.035414e-72 1.000000e+00 [519,] 4.146659e-73 8.293318e-73 1.000000e+00 [520,] 1.160395e-73 2.320789e-73 1.000000e+00 [521,] 3.255026e-74 6.510053e-74 1.000000e+00 [522,] 9.959524e-75 1.991905e-74 1.000000e+00 [523,] 1.414577e-74 2.829154e-74 1.000000e+00 [524,] 4.080413e-75 8.160825e-75 1.000000e+00 [525,] 8.479807e-75 1.695961e-74 1.000000e+00 [526,] 6.030996e-25 1.206199e-24 1.000000e+00 [527,] 1.000000e+00 2.707954e-121 1.353977e-121 [528,] 1.000000e+00 9.199342e-119 4.599671e-119 [529,] 1.000000e+00 2.454403e-116 1.227202e-116 [530,] 1.000000e+00 6.898294e-114 3.449147e-114 [531,] 1.000000e+00 2.160068e-111 1.080034e-111 [532,] 1.000000e+00 6.992105e-109 3.496053e-109 [533,] 1.000000e+00 1.932026e-106 9.660132e-107 [534,] 1.000000e+00 4.849570e-104 2.424785e-104 [535,] 1.000000e+00 1.524778e-101 7.623892e-102 [536,] 1.000000e+00 4.827657e-99 2.413828e-99 [537,] 1.000000e+00 1.507211e-96 7.536054e-97 [538,] 1.000000e+00 2.406110e-94 1.203055e-94 [539,] 1.000000e+00 5.985939e-92 2.992969e-92 [540,] 1.000000e+00 1.257432e-89 6.287161e-90 [541,] 1.000000e+00 3.856379e-87 1.928189e-87 [542,] 1.000000e+00 1.036964e-84 5.184821e-85 [543,] 1.000000e+00 1.488756e-82 7.443781e-83 [544,] 1.000000e+00 4.437535e-80 2.218767e-80 [545,] 1.000000e+00 2.399356e-79 1.199678e-79 [546,] 1.000000e+00 7.093273e-77 3.546637e-77 [547,] 1.000000e+00 2.293084e-74 1.146542e-74 [548,] 1.000000e+00 6.391683e-72 3.195841e-72 [549,] 1.000000e+00 2.018376e-69 1.009188e-69 [550,] 1.000000e+00 3.414601e-67 1.707301e-67 [551,] 1.000000e+00 9.434516e-65 4.717258e-65 [552,] 1.000000e+00 2.845469e-62 1.422735e-62 [553,] 1.000000e+00 3.937164e-60 1.968582e-60 [554,] 1.000000e+00 1.192682e-57 5.963411e-58 [555,] 1.000000e+00 1.763345e-55 8.816724e-56 [556,] 1.000000e+00 3.550326e-53 1.775163e-53 [557,] 1.000000e+00 7.130529e-51 3.565264e-51 [558,] 1.000000e+00 2.005242e-48 1.002621e-48 [559,] 1.000000e+00 5.640158e-46 2.820079e-46 [560,] 1.000000e+00 1.307685e-43 6.538427e-44 [561,] 1.000000e+00 1.186298e-41 5.931492e-42 [562,] 1.000000e+00 1.502606e-39 7.513030e-40 [563,] 1.000000e+00 4.142390e-37 2.071195e-37 [564,] 1.000000e+00 9.328551e-35 4.664275e-35 [565,] 1.000000e+00 1.767472e-32 8.837361e-33 [566,] 1.000000e+00 3.432838e-30 1.716419e-30 [567,] 1.000000e+00 5.332257e-28 2.666128e-28 [568,] 1.000000e+00 1.249985e-25 6.249926e-26 [569,] 1.000000e+00 1.223713e-24 6.118566e-25 [570,] 1.000000e+00 2.827164e-22 1.413582e-22 [571,] 1.000000e+00 3.358738e-20 1.679369e-20 [572,] 1.000000e+00 9.960217e-19 4.980108e-19 [573,] 1.000000e+00 4.918272e-16 2.459136e-16 [574,] 1.000000e+00 1.840632e-13 9.203161e-14 [575,] 1.000000e+00 1.425178e-11 7.125889e-12 [576,] 1.000000e+00 7.480532e-09 3.740266e-09 > postscript(file="/var/www/html/freestat/rcomp/tmp/1fedx1293645852.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/freestat/rcomp/tmp/28oc01293645852.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/freestat/rcomp/tmp/38oc01293645852.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/freestat/rcomp/tmp/48oc01293645852.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/freestat/rcomp/tmp/58oc01293645852.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 = 589 Frequency = 1 1 2 3 4 5 8.464328e+00 2.013082e+01 3.684305e+00 -1.017926e+01 5.632891e+00 6 7 8 9 10 4.672424e+00 -2.113670e+00 8.380918e+00 -1.209164e+01 -2.332091e+00 11 12 13 14 15 -8.037254e-01 -5.221019e-01 7.948979e+00 1.741616e+00 -6.161324e+00 16 17 18 19 20 -6.081434e+00 6.888138e+00 -1.255763e+00 3.709982e+00 3.117333e+00 21 22 23 24 25 2.488585e+00 -1.818748e+00 7.012078e+00 -3.123524e+00 2.948687e-01 26 27 28 29 30 4.418424e+00 2.090478e+00 -7.978551e+00 6.575067e+00 -6.642968e-01 31 32 33 34 35 2.863945e+00 9.277274e+00 -2.024807e+00 -5.042074e+00 1.373882e-01 36 37 38 39 40 -6.650105e+00 4.523806e+00 -7.481487e-01 -2.250088e+00 2.396932e+00 41 42 43 44 45 3.154763e+00 3.815728e+00 -8.484061e+00 1.773405e-01 -1.019999e+01 46 47 48 49 50 -3.167972e+00 -1.481587e-01 4.600735e-01 -3.912389e+00 -9.316647e+00 51 52 53 54 55 -7.890508e+00 1.613757e+00 -1.736845e+00 3.046343e-01 -9.921026e+00 56 57 58 59 60 -1.997308e+00 -1.048413e+01 -2.733878e+00 -5.128870e+00 -3.841798e+00 61 62 63 64 65 -3.633882e+00 -3.332085e+00 -6.703837e+00 -8.322809e+00 -3.035122e+00 66 67 68 69 70 -1.104261e+01 4.794165e+00 -2.740594e+00 -4.397875e+00 -6.818793e+00 71 72 73 74 75 -1.392042e+01 -8.091156e+00 -3.944353e+00 -5.370021e+00 -1.919989e+01 76 77 78 79 80 -7.359364e+00 -1.148227e+01 -6.977222e+00 -1.273791e+01 -1.099634e+01 81 82 83 84 85 -9.118026e+00 -2.317637e+01 -1.638522e+01 -1.811097e+01 -2.131841e+01 86 87 88 89 90 -9.902323e+00 1.946584e+00 -9.639159e+00 -5.634773e+00 4.646007e+00 91 92 93 94 95 3.034166e+00 3.269001e+00 -3.673045e+00 3.736835e+00 -5.371805e+00 96 97 98 99 100 4.022425e+00 3.996791e+00 3.095573e-02 -5.080552e+00 -1.207077e+00 101 102 103 104 105 -2.989060e+00 7.560322e+00 1.057145e+00 -4.161320e+00 6.516305e+00 106 107 108 109 110 -2.933674e+00 7.111210e+00 3.264884e+00 1.336456e+00 1.731668e+00 111 112 113 114 115 -2.909071e+00 -3.082950e+00 -3.915570e+00 4.521706e+00 2.413311e+00 116 117 118 119 120 -7.918462e-01 1.427120e+00 7.460279e+00 4.152292e+00 4.707465e-01 121 122 123 124 125 -4.265872e+00 5.517678e+00 -5.774741e-01 1.080938e+01 6.762213e-01 126 127 128 129 130 -1.512378e+00 -6.220782e+00 -1.122912e+00 2.676043e+00 -7.171300e-01 131 132 133 134 135 2.361654e+00 -1.342618e+00 -9.465486e+00 -3.445541e+00 4.648334e+00 136 137 138 139 140 9.144664e+00 6.895368e+00 4.469585e+00 1.344915e+00 1.194739e+01 141 142 143 144 145 -1.007670e+01 2.348123e+00 6.563733e-01 1.937644e+00 -2.666247e+00 146 147 148 149 150 8.339364e+00 1.689438e+00 4.798750e+00 -8.470525e+00 -2.297545e+00 151 152 153 154 155 -5.175196e-01 -1.634656e+00 -2.813178e+00 -9.719617e+00 3.080725e+00 156 157 158 159 160 4.395359e+00 -3.889484e+00 -3.980531e+00 -3.674392e+00 -3.456716e+00 161 162 163 164 165 -5.604135e+00 1.577833e+00 -7.819628e+00 -3.008714e+00 -1.920941e+00 166 167 168 169 170 -7.638249e+00 2.085900e+00 8.777539e+00 -7.323208e+00 4.627454e+00 171 172 173 174 175 -5.620451e+00 -4.048657e+00 -5.619919e+00 -3.470164e-01 -7.596715e+00 176 177 178 179 180 -3.614473e+00 -7.984102e+00 -7.753524e+00 -1.557685e+00 3.933771e+00 181 182 183 184 185 -3.831038e+00 -3.407875e+00 2.082319e+00 7.523901e-01 -1.111827e+01 186 187 188 189 190 -5.232714e+00 1.591044e+00 -8.869788e+00 9.362148e-01 -4.468212e+00 191 192 193 194 195 -8.845299e+00 -1.173465e+01 -1.499360e+01 -1.423325e+01 -8.900840e+00 196 197 198 199 200 -1.210137e+01 -1.627754e+01 2.100143e+00 -8.065534e+00 2.553390e+01 201 202 203 204 205 -1.041218e+01 -9.919521e+00 -1.259631e+01 -3.597365e+00 -6.576023e+00 206 207 208 209 210 -8.272962e+00 -8.967380e+00 1.456099e+00 3.005828e+00 2.879939e-03 211 212 213 214 215 -4.438637e+00 -2.059624e+00 -2.245853e+00 -5.297683e+00 -7.673085e+00 216 217 218 219 220 -8.026034e-01 -5.894918e+00 3.909198e+00 -6.703507e+00 -5.295797e+00 221 222 223 224 225 -3.756083e+00 -2.777192e+00 -1.010839e+01 -9.552847e+00 -1.290084e+01 226 227 228 229 230 -1.014784e+01 -1.350631e+01 -5.920137e+00 -3.195235e+00 -4.143248e+00 231 232 233 234 235 -2.243213e+00 -2.762337e+00 -1.503601e+00 -4.786170e+00 -1.796502e-01 236 237 238 239 240 -1.046490e+01 -1.215343e+01 -6.932109e+00 -6.289822e-01 -8.280547e+00 241 242 243 244 245 -9.460726e+00 -4.818488e+00 -1.535704e+00 -7.188073e+00 -1.427014e+01 246 247 248 249 250 5.261674e+00 -1.924026e+00 -4.937063e+00 4.904301e+00 -2.548039e+00 251 252 253 254 255 4.422112e+00 -7.419575e+00 -8.029612e+00 -1.336652e+00 -3.366953e+00 256 257 258 259 260 -1.366328e+00 4.826271e+00 4.229960e-01 -8.530048e+00 1.007894e-01 261 262 263 264 265 -4.072885e+00 3.447124e+00 -3.916064e+00 3.357711e+00 -6.333565e+00 266 267 268 269 270 -2.456749e+00 -4.257735e+00 2.411166e+00 -2.112360e-01 -8.587899e+00 271 272 273 274 275 9.409951e-02 1.639989e+00 3.890707e+00 -7.797128e+00 7.448458e+00 276 277 278 279 280 -4.614510e+00 -1.861732e+00 1.254214e+01 -2.320794e+00 -6.070694e+00 281 282 283 284 285 2.472612e+00 4.522071e-01 8.534202e+00 -7.280054e-01 6.310018e-01 286 287 288 289 290 -5.314934e+00 -5.541811e+00 4.541734e+00 -7.843056e+00 -7.795154e+00 291 292 293 294 295 -9.367311e+00 -4.954651e+00 -6.820983e+00 1.813743e+00 -5.441131e+00 296 297 298 299 300 -8.197118e+00 4.233559e+00 3.076145e+00 -2.778540e+00 -1.318926e+01 301 302 303 304 305 -7.754824e+00 6.293652e-01 1.617146e+00 2.324657e+00 -1.077044e-01 306 307 308 309 310 -5.615624e+00 -8.644780e+00 1.035547e+00 -4.491362e+00 -2.671451e+00 311 312 313 314 315 -2.405719e+00 -2.828444e+00 -4.838494e+00 -1.061683e+01 -1.014508e+01 316 317 318 319 320 2.338463e-04 -8.992978e+00 -2.774032e+00 -5.715379e+00 -5.478546e+00 321 322 323 324 325 -4.626538e+00 -8.264079e+00 -6.410088e+00 -2.022465e+01 -2.556611e+00 326 327 328 329 330 -5.250831e+00 3.867523e+00 7.295919e-02 -1.600432e+01 6.894633e-01 331 332 333 334 335 -5.814681e+00 -8.380842e+00 -8.593540e+00 -5.834853e+00 -1.678961e+01 336 337 338 339 340 -9.019913e+00 -1.109807e+01 -1.193621e+01 -9.815999e+00 -6.036795e+00 341 342 343 344 345 -1.060875e+01 -1.251452e+01 -2.588106e-01 -8.132063e+00 -8.287043e+00 346 347 348 349 350 -1.360723e+01 -8.819109e+00 -9.301126e+00 -9.136712e+00 4.814530e+00 351 352 353 354 355 -1.334495e+01 -6.040788e+00 -7.711183e-02 1.721691e+00 -8.745563e+00 356 357 358 359 360 -4.914882e+00 -7.631854e+00 -1.429410e+01 -1.620007e+01 -6.825736e+00 361 362 363 364 365 -1.091408e+01 -1.180369e+01 3.418831e+00 -9.798002e+00 -1.707612e+01 366 367 368 369 370 -1.319164e+01 -2.288007e+00 -8.005386e+00 -8.516113e+00 -6.844228e+00 371 372 373 374 375 -8.081791e+00 -2.062056e+00 -9.027292e+00 -1.225732e+00 -2.933471e+00 376 377 378 379 380 -4.665275e+00 -6.532304e+00 -8.132389e+00 -5.366831e+00 -4.556410e+00 381 382 383 384 385 1.999811e+00 -2.422642e-01 -2.238471e-01 -3.961175e+00 -1.233876e+01 386 387 388 389 390 -4.257692e+00 -9.856868e+00 -9.084125e-02 -5.093623e-01 -5.444040e+00 391 392 393 394 395 5.535112e+00 6.730607e-01 1.256895e+00 -4.322303e+00 -1.094502e+01 396 397 398 399 400 -9.368908e+00 -4.020970e+00 -1.030090e+01 -1.349618e+01 -1.087367e+01 401 402 403 404 405 -8.336932e+00 -3.601847e+00 1.414664e+00 -1.183180e+01 -5.663892e-01 406 407 408 409 410 -5.070947e+00 -1.423151e+00 -1.232879e+01 -1.063743e+01 -3.524293e+00 411 412 413 414 415 -2.186279e+00 -3.716105e+00 -1.607130e+00 -9.029069e+00 -1.577982e+01 416 417 418 419 420 -4.269073e+00 -7.558887e+00 -2.457378e+00 -9.876012e+00 -1.630162e+00 421 422 423 424 425 -1.021682e+01 3.695981e-01 -1.937202e+01 -1.021490e-01 -6.506590e+00 426 427 428 429 430 -1.477353e+01 -1.368036e+01 -1.973588e+00 2.806230e+00 -9.361183e+00 431 432 433 434 435 -5.134870e+00 -1.237298e+01 -6.193866e+00 -4.016571e+00 2.059318e+00 436 437 438 439 440 -8.418025e+00 -5.178387e+00 -1.132184e+01 -1.571835e+01 -5.691437e+00 441 442 443 444 445 -5.282922e+00 -9.454454e+00 -7.422582e+00 -8.747528e+00 -2.188783e+00 446 447 448 449 450 -1.551824e+01 -7.676564e+00 -8.858174e+00 -1.487991e+01 -2.707192e+00 451 452 453 454 455 -2.245388e+00 1.071900e+00 9.689288e-01 3.349447e+00 -8.677832e+00 456 457 458 459 460 -6.190921e+00 -6.200216e-01 1.396973e+00 -9.713814e+00 -7.314318e+00 461 462 463 464 465 -1.174711e+00 -2.518591e+00 -6.154023e+00 -1.855996e-01 -1.814706e+00 466 467 468 469 470 -2.358454e+00 -1.227286e+01 -1.482947e+00 -1.028202e+01 -3.374203e+00 471 472 473 474 475 4.825680e+00 -6.554237e+00 -2.237483e+00 -9.190129e+00 -5.737986e-01 476 477 478 479 480 -2.250249e+00 -4.495156e+00 -4.597812e+00 3.967008e+00 -1.085722e+01 481 482 483 484 485 -2.912574e+00 -2.256643e+00 -9.772051e+00 -1.006552e+01 -1.466945e+01 486 487 488 489 490 -9.020717e+00 -6.783137e+00 -1.413017e+01 -9.217490e+00 -9.034729e+00 491 492 493 494 495 -1.169150e+01 2.161951e-01 -3.275268e+00 -1.087961e+01 6.979003e-01 496 497 498 499 500 -1.173240e+01 -1.057020e+01 -9.851838e+00 1.330189e+01 -4.302809e+00 501 502 503 504 505 -1.635559e+00 -7.868065e+00 -6.343591e+00 -2.713257e+00 -1.806510e+00 506 507 508 509 510 -8.745664e+00 1.551508e+00 1.739173e+00 5.409489e+00 -1.038762e+01 511 512 513 514 515 8.758784e+02 -2.040067e-01 -2.185599e+00 -1.341008e+01 -5.486828e+00 516 517 518 519 520 -1.131193e+01 -1.255453e+01 -4.422624e+00 -3.354258e+00 -6.927837e+00 521 522 523 524 525 -3.021482e+00 -6.446570e+00 -1.488881e+01 -1.608805e+01 -6.036465e+00 526 527 528 529 530 -6.716511e+00 -3.579940e+00 -8.636345e+00 -1.555820e+01 -2.844068e+00 531 532 533 534 535 -1.568663e+01 7.604625e+02 8.764713e+02 -6.374334e+00 -4.040232e+00 536 537 538 539 540 -1.266852e+01 -6.669758e+00 -1.159025e+01 -5.567872e+00 -1.433900e+01 541 542 543 544 545 -4.584389e+00 -1.528106e+01 -2.295174e+00 -7.161032e+00 -7.668840e-01 546 547 548 549 550 -6.203831e+00 -1.098963e+01 -1.010455e+01 -8.494776e+00 -1.197824e+01 551 552 553 554 555 -1.257357e+01 -1.081448e+01 -3.203562e+00 -6.306076e+00 -9.063120e+00 556 557 558 559 560 -6.153311e+00 -5.932945e+00 -4.971402e+00 -1.140417e+01 -1.095120e+01 561 562 563 564 565 -7.335587e+00 -1.162483e+01 -3.547614e+00 -1.027371e+01 -6.259756e+00 566 567 568 569 570 -5.723048e+00 -9.172901e+00 -3.756818e+00 -3.759076e+00 -1.042397e+00 571 572 573 574 575 -8.032350e+00 -8.323509e+00 -8.552383e+00 -8.905003e+00 -3.192250e+00 576 577 578 579 580 1.258708e+00 -9.910202e+00 7.485737e+00 -1.350696e+01 3.018549e+00 581 582 583 584 585 -5.267373e+00 -2.533336e+00 -3.300970e+00 -6.979721e+00 -1.360075e+01 586 587 588 589 -7.169105e+00 -1.447020e+01 -5.456367e+00 -2.574610e+00 > postscript(file="/var/www/html/freestat/rcomp/tmp/60xbl1293645852.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 = 589 Frequency = 1 lag(myerror, k = 1) myerror 0 8.464328e+00 NA 1 2.013082e+01 8.464328e+00 2 3.684305e+00 2.013082e+01 3 -1.017926e+01 3.684305e+00 4 5.632891e+00 -1.017926e+01 5 4.672424e+00 5.632891e+00 6 -2.113670e+00 4.672424e+00 7 8.380918e+00 -2.113670e+00 8 -1.209164e+01 8.380918e+00 9 -2.332091e+00 -1.209164e+01 10 -8.037254e-01 -2.332091e+00 11 -5.221019e-01 -8.037254e-01 12 7.948979e+00 -5.221019e-01 13 1.741616e+00 7.948979e+00 14 -6.161324e+00 1.741616e+00 15 -6.081434e+00 -6.161324e+00 16 6.888138e+00 -6.081434e+00 17 -1.255763e+00 6.888138e+00 18 3.709982e+00 -1.255763e+00 19 3.117333e+00 3.709982e+00 20 2.488585e+00 3.117333e+00 21 -1.818748e+00 2.488585e+00 22 7.012078e+00 -1.818748e+00 23 -3.123524e+00 7.012078e+00 24 2.948687e-01 -3.123524e+00 25 4.418424e+00 2.948687e-01 26 2.090478e+00 4.418424e+00 27 -7.978551e+00 2.090478e+00 28 6.575067e+00 -7.978551e+00 29 -6.642968e-01 6.575067e+00 30 2.863945e+00 -6.642968e-01 31 9.277274e+00 2.863945e+00 32 -2.024807e+00 9.277274e+00 33 -5.042074e+00 -2.024807e+00 34 1.373882e-01 -5.042074e+00 35 -6.650105e+00 1.373882e-01 36 4.523806e+00 -6.650105e+00 37 -7.481487e-01 4.523806e+00 38 -2.250088e+00 -7.481487e-01 39 2.396932e+00 -2.250088e+00 40 3.154763e+00 2.396932e+00 41 3.815728e+00 3.154763e+00 42 -8.484061e+00 3.815728e+00 43 1.773405e-01 -8.484061e+00 44 -1.019999e+01 1.773405e-01 45 -3.167972e+00 -1.019999e+01 46 -1.481587e-01 -3.167972e+00 47 4.600735e-01 -1.481587e-01 48 -3.912389e+00 4.600735e-01 49 -9.316647e+00 -3.912389e+00 50 -7.890508e+00 -9.316647e+00 51 1.613757e+00 -7.890508e+00 52 -1.736845e+00 1.613757e+00 53 3.046343e-01 -1.736845e+00 54 -9.921026e+00 3.046343e-01 55 -1.997308e+00 -9.921026e+00 56 -1.048413e+01 -1.997308e+00 57 -2.733878e+00 -1.048413e+01 58 -5.128870e+00 -2.733878e+00 59 -3.841798e+00 -5.128870e+00 60 -3.633882e+00 -3.841798e+00 61 -3.332085e+00 -3.633882e+00 62 -6.703837e+00 -3.332085e+00 63 -8.322809e+00 -6.703837e+00 64 -3.035122e+00 -8.322809e+00 65 -1.104261e+01 -3.035122e+00 66 4.794165e+00 -1.104261e+01 67 -2.740594e+00 4.794165e+00 68 -4.397875e+00 -2.740594e+00 69 -6.818793e+00 -4.397875e+00 70 -1.392042e+01 -6.818793e+00 71 -8.091156e+00 -1.392042e+01 72 -3.944353e+00 -8.091156e+00 73 -5.370021e+00 -3.944353e+00 74 -1.919989e+01 -5.370021e+00 75 -7.359364e+00 -1.919989e+01 76 -1.148227e+01 -7.359364e+00 77 -6.977222e+00 -1.148227e+01 78 -1.273791e+01 -6.977222e+00 79 -1.099634e+01 -1.273791e+01 80 -9.118026e+00 -1.099634e+01 81 -2.317637e+01 -9.118026e+00 82 -1.638522e+01 -2.317637e+01 83 -1.811097e+01 -1.638522e+01 84 -2.131841e+01 -1.811097e+01 85 -9.902323e+00 -2.131841e+01 86 1.946584e+00 -9.902323e+00 87 -9.639159e+00 1.946584e+00 88 -5.634773e+00 -9.639159e+00 89 4.646007e+00 -5.634773e+00 90 3.034166e+00 4.646007e+00 91 3.269001e+00 3.034166e+00 92 -3.673045e+00 3.269001e+00 93 3.736835e+00 -3.673045e+00 94 -5.371805e+00 3.736835e+00 95 4.022425e+00 -5.371805e+00 96 3.996791e+00 4.022425e+00 97 3.095573e-02 3.996791e+00 98 -5.080552e+00 3.095573e-02 99 -1.207077e+00 -5.080552e+00 100 -2.989060e+00 -1.207077e+00 101 7.560322e+00 -2.989060e+00 102 1.057145e+00 7.560322e+00 103 -4.161320e+00 1.057145e+00 104 6.516305e+00 -4.161320e+00 105 -2.933674e+00 6.516305e+00 106 7.111210e+00 -2.933674e+00 107 3.264884e+00 7.111210e+00 108 1.336456e+00 3.264884e+00 109 1.731668e+00 1.336456e+00 110 -2.909071e+00 1.731668e+00 111 -3.082950e+00 -2.909071e+00 112 -3.915570e+00 -3.082950e+00 113 4.521706e+00 -3.915570e+00 114 2.413311e+00 4.521706e+00 115 -7.918462e-01 2.413311e+00 116 1.427120e+00 -7.918462e-01 117 7.460279e+00 1.427120e+00 118 4.152292e+00 7.460279e+00 119 4.707465e-01 4.152292e+00 120 -4.265872e+00 4.707465e-01 121 5.517678e+00 -4.265872e+00 122 -5.774741e-01 5.517678e+00 123 1.080938e+01 -5.774741e-01 124 6.762213e-01 1.080938e+01 125 -1.512378e+00 6.762213e-01 126 -6.220782e+00 -1.512378e+00 127 -1.122912e+00 -6.220782e+00 128 2.676043e+00 -1.122912e+00 129 -7.171300e-01 2.676043e+00 130 2.361654e+00 -7.171300e-01 131 -1.342618e+00 2.361654e+00 132 -9.465486e+00 -1.342618e+00 133 -3.445541e+00 -9.465486e+00 134 4.648334e+00 -3.445541e+00 135 9.144664e+00 4.648334e+00 136 6.895368e+00 9.144664e+00 137 4.469585e+00 6.895368e+00 138 1.344915e+00 4.469585e+00 139 1.194739e+01 1.344915e+00 140 -1.007670e+01 1.194739e+01 141 2.348123e+00 -1.007670e+01 142 6.563733e-01 2.348123e+00 143 1.937644e+00 6.563733e-01 144 -2.666247e+00 1.937644e+00 145 8.339364e+00 -2.666247e+00 146 1.689438e+00 8.339364e+00 147 4.798750e+00 1.689438e+00 148 -8.470525e+00 4.798750e+00 149 -2.297545e+00 -8.470525e+00 150 -5.175196e-01 -2.297545e+00 151 -1.634656e+00 -5.175196e-01 152 -2.813178e+00 -1.634656e+00 153 -9.719617e+00 -2.813178e+00 154 3.080725e+00 -9.719617e+00 155 4.395359e+00 3.080725e+00 156 -3.889484e+00 4.395359e+00 157 -3.980531e+00 -3.889484e+00 158 -3.674392e+00 -3.980531e+00 159 -3.456716e+00 -3.674392e+00 160 -5.604135e+00 -3.456716e+00 161 1.577833e+00 -5.604135e+00 162 -7.819628e+00 1.577833e+00 163 -3.008714e+00 -7.819628e+00 164 -1.920941e+00 -3.008714e+00 165 -7.638249e+00 -1.920941e+00 166 2.085900e+00 -7.638249e+00 167 8.777539e+00 2.085900e+00 168 -7.323208e+00 8.777539e+00 169 4.627454e+00 -7.323208e+00 170 -5.620451e+00 4.627454e+00 171 -4.048657e+00 -5.620451e+00 172 -5.619919e+00 -4.048657e+00 173 -3.470164e-01 -5.619919e+00 174 -7.596715e+00 -3.470164e-01 175 -3.614473e+00 -7.596715e+00 176 -7.984102e+00 -3.614473e+00 177 -7.753524e+00 -7.984102e+00 178 -1.557685e+00 -7.753524e+00 179 3.933771e+00 -1.557685e+00 180 -3.831038e+00 3.933771e+00 181 -3.407875e+00 -3.831038e+00 182 2.082319e+00 -3.407875e+00 183 7.523901e-01 2.082319e+00 184 -1.111827e+01 7.523901e-01 185 -5.232714e+00 -1.111827e+01 186 1.591044e+00 -5.232714e+00 187 -8.869788e+00 1.591044e+00 188 9.362148e-01 -8.869788e+00 189 -4.468212e+00 9.362148e-01 190 -8.845299e+00 -4.468212e+00 191 -1.173465e+01 -8.845299e+00 192 -1.499360e+01 -1.173465e+01 193 -1.423325e+01 -1.499360e+01 194 -8.900840e+00 -1.423325e+01 195 -1.210137e+01 -8.900840e+00 196 -1.627754e+01 -1.210137e+01 197 2.100143e+00 -1.627754e+01 198 -8.065534e+00 2.100143e+00 199 2.553390e+01 -8.065534e+00 200 -1.041218e+01 2.553390e+01 201 -9.919521e+00 -1.041218e+01 202 -1.259631e+01 -9.919521e+00 203 -3.597365e+00 -1.259631e+01 204 -6.576023e+00 -3.597365e+00 205 -8.272962e+00 -6.576023e+00 206 -8.967380e+00 -8.272962e+00 207 1.456099e+00 -8.967380e+00 208 3.005828e+00 1.456099e+00 209 2.879939e-03 3.005828e+00 210 -4.438637e+00 2.879939e-03 211 -2.059624e+00 -4.438637e+00 212 -2.245853e+00 -2.059624e+00 213 -5.297683e+00 -2.245853e+00 214 -7.673085e+00 -5.297683e+00 215 -8.026034e-01 -7.673085e+00 216 -5.894918e+00 -8.026034e-01 217 3.909198e+00 -5.894918e+00 218 -6.703507e+00 3.909198e+00 219 -5.295797e+00 -6.703507e+00 220 -3.756083e+00 -5.295797e+00 221 -2.777192e+00 -3.756083e+00 222 -1.010839e+01 -2.777192e+00 223 -9.552847e+00 -1.010839e+01 224 -1.290084e+01 -9.552847e+00 225 -1.014784e+01 -1.290084e+01 226 -1.350631e+01 -1.014784e+01 227 -5.920137e+00 -1.350631e+01 228 -3.195235e+00 -5.920137e+00 229 -4.143248e+00 -3.195235e+00 230 -2.243213e+00 -4.143248e+00 231 -2.762337e+00 -2.243213e+00 232 -1.503601e+00 -2.762337e+00 233 -4.786170e+00 -1.503601e+00 234 -1.796502e-01 -4.786170e+00 235 -1.046490e+01 -1.796502e-01 236 -1.215343e+01 -1.046490e+01 237 -6.932109e+00 -1.215343e+01 238 -6.289822e-01 -6.932109e+00 239 -8.280547e+00 -6.289822e-01 240 -9.460726e+00 -8.280547e+00 241 -4.818488e+00 -9.460726e+00 242 -1.535704e+00 -4.818488e+00 243 -7.188073e+00 -1.535704e+00 244 -1.427014e+01 -7.188073e+00 245 5.261674e+00 -1.427014e+01 246 -1.924026e+00 5.261674e+00 247 -4.937063e+00 -1.924026e+00 248 4.904301e+00 -4.937063e+00 249 -2.548039e+00 4.904301e+00 250 4.422112e+00 -2.548039e+00 251 -7.419575e+00 4.422112e+00 252 -8.029612e+00 -7.419575e+00 253 -1.336652e+00 -8.029612e+00 254 -3.366953e+00 -1.336652e+00 255 -1.366328e+00 -3.366953e+00 256 4.826271e+00 -1.366328e+00 257 4.229960e-01 4.826271e+00 258 -8.530048e+00 4.229960e-01 259 1.007894e-01 -8.530048e+00 260 -4.072885e+00 1.007894e-01 261 3.447124e+00 -4.072885e+00 262 -3.916064e+00 3.447124e+00 263 3.357711e+00 -3.916064e+00 264 -6.333565e+00 3.357711e+00 265 -2.456749e+00 -6.333565e+00 266 -4.257735e+00 -2.456749e+00 267 2.411166e+00 -4.257735e+00 268 -2.112360e-01 2.411166e+00 269 -8.587899e+00 -2.112360e-01 270 9.409951e-02 -8.587899e+00 271 1.639989e+00 9.409951e-02 272 3.890707e+00 1.639989e+00 273 -7.797128e+00 3.890707e+00 274 7.448458e+00 -7.797128e+00 275 -4.614510e+00 7.448458e+00 276 -1.861732e+00 -4.614510e+00 277 1.254214e+01 -1.861732e+00 278 -2.320794e+00 1.254214e+01 279 -6.070694e+00 -2.320794e+00 280 2.472612e+00 -6.070694e+00 281 4.522071e-01 2.472612e+00 282 8.534202e+00 4.522071e-01 283 -7.280054e-01 8.534202e+00 284 6.310018e-01 -7.280054e-01 285 -5.314934e+00 6.310018e-01 286 -5.541811e+00 -5.314934e+00 287 4.541734e+00 -5.541811e+00 288 -7.843056e+00 4.541734e+00 289 -7.795154e+00 -7.843056e+00 290 -9.367311e+00 -7.795154e+00 291 -4.954651e+00 -9.367311e+00 292 -6.820983e+00 -4.954651e+00 293 1.813743e+00 -6.820983e+00 294 -5.441131e+00 1.813743e+00 295 -8.197118e+00 -5.441131e+00 296 4.233559e+00 -8.197118e+00 297 3.076145e+00 4.233559e+00 298 -2.778540e+00 3.076145e+00 299 -1.318926e+01 -2.778540e+00 300 -7.754824e+00 -1.318926e+01 301 6.293652e-01 -7.754824e+00 302 1.617146e+00 6.293652e-01 303 2.324657e+00 1.617146e+00 304 -1.077044e-01 2.324657e+00 305 -5.615624e+00 -1.077044e-01 306 -8.644780e+00 -5.615624e+00 307 1.035547e+00 -8.644780e+00 308 -4.491362e+00 1.035547e+00 309 -2.671451e+00 -4.491362e+00 310 -2.405719e+00 -2.671451e+00 311 -2.828444e+00 -2.405719e+00 312 -4.838494e+00 -2.828444e+00 313 -1.061683e+01 -4.838494e+00 314 -1.014508e+01 -1.061683e+01 315 2.338463e-04 -1.014508e+01 316 -8.992978e+00 2.338463e-04 317 -2.774032e+00 -8.992978e+00 318 -5.715379e+00 -2.774032e+00 319 -5.478546e+00 -5.715379e+00 320 -4.626538e+00 -5.478546e+00 321 -8.264079e+00 -4.626538e+00 322 -6.410088e+00 -8.264079e+00 323 -2.022465e+01 -6.410088e+00 324 -2.556611e+00 -2.022465e+01 325 -5.250831e+00 -2.556611e+00 326 3.867523e+00 -5.250831e+00 327 7.295919e-02 3.867523e+00 328 -1.600432e+01 7.295919e-02 329 6.894633e-01 -1.600432e+01 330 -5.814681e+00 6.894633e-01 331 -8.380842e+00 -5.814681e+00 332 -8.593540e+00 -8.380842e+00 333 -5.834853e+00 -8.593540e+00 334 -1.678961e+01 -5.834853e+00 335 -9.019913e+00 -1.678961e+01 336 -1.109807e+01 -9.019913e+00 337 -1.193621e+01 -1.109807e+01 338 -9.815999e+00 -1.193621e+01 339 -6.036795e+00 -9.815999e+00 340 -1.060875e+01 -6.036795e+00 341 -1.251452e+01 -1.060875e+01 342 -2.588106e-01 -1.251452e+01 343 -8.132063e+00 -2.588106e-01 344 -8.287043e+00 -8.132063e+00 345 -1.360723e+01 -8.287043e+00 346 -8.819109e+00 -1.360723e+01 347 -9.301126e+00 -8.819109e+00 348 -9.136712e+00 -9.301126e+00 349 4.814530e+00 -9.136712e+00 350 -1.334495e+01 4.814530e+00 351 -6.040788e+00 -1.334495e+01 352 -7.711183e-02 -6.040788e+00 353 1.721691e+00 -7.711183e-02 354 -8.745563e+00 1.721691e+00 355 -4.914882e+00 -8.745563e+00 356 -7.631854e+00 -4.914882e+00 357 -1.429410e+01 -7.631854e+00 358 -1.620007e+01 -1.429410e+01 359 -6.825736e+00 -1.620007e+01 360 -1.091408e+01 -6.825736e+00 361 -1.180369e+01 -1.091408e+01 362 3.418831e+00 -1.180369e+01 363 -9.798002e+00 3.418831e+00 364 -1.707612e+01 -9.798002e+00 365 -1.319164e+01 -1.707612e+01 366 -2.288007e+00 -1.319164e+01 367 -8.005386e+00 -2.288007e+00 368 -8.516113e+00 -8.005386e+00 369 -6.844228e+00 -8.516113e+00 370 -8.081791e+00 -6.844228e+00 371 -2.062056e+00 -8.081791e+00 372 -9.027292e+00 -2.062056e+00 373 -1.225732e+00 -9.027292e+00 374 -2.933471e+00 -1.225732e+00 375 -4.665275e+00 -2.933471e+00 376 -6.532304e+00 -4.665275e+00 377 -8.132389e+00 -6.532304e+00 378 -5.366831e+00 -8.132389e+00 379 -4.556410e+00 -5.366831e+00 380 1.999811e+00 -4.556410e+00 381 -2.422642e-01 1.999811e+00 382 -2.238471e-01 -2.422642e-01 383 -3.961175e+00 -2.238471e-01 384 -1.233876e+01 -3.961175e+00 385 -4.257692e+00 -1.233876e+01 386 -9.856868e+00 -4.257692e+00 387 -9.084125e-02 -9.856868e+00 388 -5.093623e-01 -9.084125e-02 389 -5.444040e+00 -5.093623e-01 390 5.535112e+00 -5.444040e+00 391 6.730607e-01 5.535112e+00 392 1.256895e+00 6.730607e-01 393 -4.322303e+00 1.256895e+00 394 -1.094502e+01 -4.322303e+00 395 -9.368908e+00 -1.094502e+01 396 -4.020970e+00 -9.368908e+00 397 -1.030090e+01 -4.020970e+00 398 -1.349618e+01 -1.030090e+01 399 -1.087367e+01 -1.349618e+01 400 -8.336932e+00 -1.087367e+01 401 -3.601847e+00 -8.336932e+00 402 1.414664e+00 -3.601847e+00 403 -1.183180e+01 1.414664e+00 404 -5.663892e-01 -1.183180e+01 405 -5.070947e+00 -5.663892e-01 406 -1.423151e+00 -5.070947e+00 407 -1.232879e+01 -1.423151e+00 408 -1.063743e+01 -1.232879e+01 409 -3.524293e+00 -1.063743e+01 410 -2.186279e+00 -3.524293e+00 411 -3.716105e+00 -2.186279e+00 412 -1.607130e+00 -3.716105e+00 413 -9.029069e+00 -1.607130e+00 414 -1.577982e+01 -9.029069e+00 415 -4.269073e+00 -1.577982e+01 416 -7.558887e+00 -4.269073e+00 417 -2.457378e+00 -7.558887e+00 418 -9.876012e+00 -2.457378e+00 419 -1.630162e+00 -9.876012e+00 420 -1.021682e+01 -1.630162e+00 421 3.695981e-01 -1.021682e+01 422 -1.937202e+01 3.695981e-01 423 -1.021490e-01 -1.937202e+01 424 -6.506590e+00 -1.021490e-01 425 -1.477353e+01 -6.506590e+00 426 -1.368036e+01 -1.477353e+01 427 -1.973588e+00 -1.368036e+01 428 2.806230e+00 -1.973588e+00 429 -9.361183e+00 2.806230e+00 430 -5.134870e+00 -9.361183e+00 431 -1.237298e+01 -5.134870e+00 432 -6.193866e+00 -1.237298e+01 433 -4.016571e+00 -6.193866e+00 434 2.059318e+00 -4.016571e+00 435 -8.418025e+00 2.059318e+00 436 -5.178387e+00 -8.418025e+00 437 -1.132184e+01 -5.178387e+00 438 -1.571835e+01 -1.132184e+01 439 -5.691437e+00 -1.571835e+01 440 -5.282922e+00 -5.691437e+00 441 -9.454454e+00 -5.282922e+00 442 -7.422582e+00 -9.454454e+00 443 -8.747528e+00 -7.422582e+00 444 -2.188783e+00 -8.747528e+00 445 -1.551824e+01 -2.188783e+00 446 -7.676564e+00 -1.551824e+01 447 -8.858174e+00 -7.676564e+00 448 -1.487991e+01 -8.858174e+00 449 -2.707192e+00 -1.487991e+01 450 -2.245388e+00 -2.707192e+00 451 1.071900e+00 -2.245388e+00 452 9.689288e-01 1.071900e+00 453 3.349447e+00 9.689288e-01 454 -8.677832e+00 3.349447e+00 455 -6.190921e+00 -8.677832e+00 456 -6.200216e-01 -6.190921e+00 457 1.396973e+00 -6.200216e-01 458 -9.713814e+00 1.396973e+00 459 -7.314318e+00 -9.713814e+00 460 -1.174711e+00 -7.314318e+00 461 -2.518591e+00 -1.174711e+00 462 -6.154023e+00 -2.518591e+00 463 -1.855996e-01 -6.154023e+00 464 -1.814706e+00 -1.855996e-01 465 -2.358454e+00 -1.814706e+00 466 -1.227286e+01 -2.358454e+00 467 -1.482947e+00 -1.227286e+01 468 -1.028202e+01 -1.482947e+00 469 -3.374203e+00 -1.028202e+01 470 4.825680e+00 -3.374203e+00 471 -6.554237e+00 4.825680e+00 472 -2.237483e+00 -6.554237e+00 473 -9.190129e+00 -2.237483e+00 474 -5.737986e-01 -9.190129e+00 475 -2.250249e+00 -5.737986e-01 476 -4.495156e+00 -2.250249e+00 477 -4.597812e+00 -4.495156e+00 478 3.967008e+00 -4.597812e+00 479 -1.085722e+01 3.967008e+00 480 -2.912574e+00 -1.085722e+01 481 -2.256643e+00 -2.912574e+00 482 -9.772051e+00 -2.256643e+00 483 -1.006552e+01 -9.772051e+00 484 -1.466945e+01 -1.006552e+01 485 -9.020717e+00 -1.466945e+01 486 -6.783137e+00 -9.020717e+00 487 -1.413017e+01 -6.783137e+00 488 -9.217490e+00 -1.413017e+01 489 -9.034729e+00 -9.217490e+00 490 -1.169150e+01 -9.034729e+00 491 2.161951e-01 -1.169150e+01 492 -3.275268e+00 2.161951e-01 493 -1.087961e+01 -3.275268e+00 494 6.979003e-01 -1.087961e+01 495 -1.173240e+01 6.979003e-01 496 -1.057020e+01 -1.173240e+01 497 -9.851838e+00 -1.057020e+01 498 1.330189e+01 -9.851838e+00 499 -4.302809e+00 1.330189e+01 500 -1.635559e+00 -4.302809e+00 501 -7.868065e+00 -1.635559e+00 502 -6.343591e+00 -7.868065e+00 503 -2.713257e+00 -6.343591e+00 504 -1.806510e+00 -2.713257e+00 505 -8.745664e+00 -1.806510e+00 506 1.551508e+00 -8.745664e+00 507 1.739173e+00 1.551508e+00 508 5.409489e+00 1.739173e+00 509 -1.038762e+01 5.409489e+00 510 8.758784e+02 -1.038762e+01 511 -2.040067e-01 8.758784e+02 512 -2.185599e+00 -2.040067e-01 513 -1.341008e+01 -2.185599e+00 514 -5.486828e+00 -1.341008e+01 515 -1.131193e+01 -5.486828e+00 516 -1.255453e+01 -1.131193e+01 517 -4.422624e+00 -1.255453e+01 518 -3.354258e+00 -4.422624e+00 519 -6.927837e+00 -3.354258e+00 520 -3.021482e+00 -6.927837e+00 521 -6.446570e+00 -3.021482e+00 522 -1.488881e+01 -6.446570e+00 523 -1.608805e+01 -1.488881e+01 524 -6.036465e+00 -1.608805e+01 525 -6.716511e+00 -6.036465e+00 526 -3.579940e+00 -6.716511e+00 527 -8.636345e+00 -3.579940e+00 528 -1.555820e+01 -8.636345e+00 529 -2.844068e+00 -1.555820e+01 530 -1.568663e+01 -2.844068e+00 531 7.604625e+02 -1.568663e+01 532 8.764713e+02 7.604625e+02 533 -6.374334e+00 8.764713e+02 534 -4.040232e+00 -6.374334e+00 535 -1.266852e+01 -4.040232e+00 536 -6.669758e+00 -1.266852e+01 537 -1.159025e+01 -6.669758e+00 538 -5.567872e+00 -1.159025e+01 539 -1.433900e+01 -5.567872e+00 540 -4.584389e+00 -1.433900e+01 541 -1.528106e+01 -4.584389e+00 542 -2.295174e+00 -1.528106e+01 543 -7.161032e+00 -2.295174e+00 544 -7.668840e-01 -7.161032e+00 545 -6.203831e+00 -7.668840e-01 546 -1.098963e+01 -6.203831e+00 547 -1.010455e+01 -1.098963e+01 548 -8.494776e+00 -1.010455e+01 549 -1.197824e+01 -8.494776e+00 550 -1.257357e+01 -1.197824e+01 551 -1.081448e+01 -1.257357e+01 552 -3.203562e+00 -1.081448e+01 553 -6.306076e+00 -3.203562e+00 554 -9.063120e+00 -6.306076e+00 555 -6.153311e+00 -9.063120e+00 556 -5.932945e+00 -6.153311e+00 557 -4.971402e+00 -5.932945e+00 558 -1.140417e+01 -4.971402e+00 559 -1.095120e+01 -1.140417e+01 560 -7.335587e+00 -1.095120e+01 561 -1.162483e+01 -7.335587e+00 562 -3.547614e+00 -1.162483e+01 563 -1.027371e+01 -3.547614e+00 564 -6.259756e+00 -1.027371e+01 565 -5.723048e+00 -6.259756e+00 566 -9.172901e+00 -5.723048e+00 567 -3.756818e+00 -9.172901e+00 568 -3.759076e+00 -3.756818e+00 569 -1.042397e+00 -3.759076e+00 570 -8.032350e+00 -1.042397e+00 571 -8.323509e+00 -8.032350e+00 572 -8.552383e+00 -8.323509e+00 573 -8.905003e+00 -8.552383e+00 574 -3.192250e+00 -8.905003e+00 575 1.258708e+00 -3.192250e+00 576 -9.910202e+00 1.258708e+00 577 7.485737e+00 -9.910202e+00 578 -1.350696e+01 7.485737e+00 579 3.018549e+00 -1.350696e+01 580 -5.267373e+00 3.018549e+00 581 -2.533336e+00 -5.267373e+00 582 -3.300970e+00 -2.533336e+00 583 -6.979721e+00 -3.300970e+00 584 -1.360075e+01 -6.979721e+00 585 -7.169105e+00 -1.360075e+01 586 -1.447020e+01 -7.169105e+00 587 -5.456367e+00 -1.447020e+01 588 -2.574610e+00 -5.456367e+00 589 NA -2.574610e+00 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2.013082e+01 8.464328e+00 [2,] 3.684305e+00 2.013082e+01 [3,] -1.017926e+01 3.684305e+00 [4,] 5.632891e+00 -1.017926e+01 [5,] 4.672424e+00 5.632891e+00 [6,] -2.113670e+00 4.672424e+00 [7,] 8.380918e+00 -2.113670e+00 [8,] -1.209164e+01 8.380918e+00 [9,] -2.332091e+00 -1.209164e+01 [10,] -8.037254e-01 -2.332091e+00 [11,] -5.221019e-01 -8.037254e-01 [12,] 7.948979e+00 -5.221019e-01 [13,] 1.741616e+00 7.948979e+00 [14,] -6.161324e+00 1.741616e+00 [15,] -6.081434e+00 -6.161324e+00 [16,] 6.888138e+00 -6.081434e+00 [17,] -1.255763e+00 6.888138e+00 [18,] 3.709982e+00 -1.255763e+00 [19,] 3.117333e+00 3.709982e+00 [20,] 2.488585e+00 3.117333e+00 [21,] -1.818748e+00 2.488585e+00 [22,] 7.012078e+00 -1.818748e+00 [23,] -3.123524e+00 7.012078e+00 [24,] 2.948687e-01 -3.123524e+00 [25,] 4.418424e+00 2.948687e-01 [26,] 2.090478e+00 4.418424e+00 [27,] -7.978551e+00 2.090478e+00 [28,] 6.575067e+00 -7.978551e+00 [29,] -6.642968e-01 6.575067e+00 [30,] 2.863945e+00 -6.642968e-01 [31,] 9.277274e+00 2.863945e+00 [32,] -2.024807e+00 9.277274e+00 [33,] -5.042074e+00 -2.024807e+00 [34,] 1.373882e-01 -5.042074e+00 [35,] -6.650105e+00 1.373882e-01 [36,] 4.523806e+00 -6.650105e+00 [37,] -7.481487e-01 4.523806e+00 [38,] -2.250088e+00 -7.481487e-01 [39,] 2.396932e+00 -2.250088e+00 [40,] 3.154763e+00 2.396932e+00 [41,] 3.815728e+00 3.154763e+00 [42,] -8.484061e+00 3.815728e+00 [43,] 1.773405e-01 -8.484061e+00 [44,] -1.019999e+01 1.773405e-01 [45,] -3.167972e+00 -1.019999e+01 [46,] -1.481587e-01 -3.167972e+00 [47,] 4.600735e-01 -1.481587e-01 [48,] -3.912389e+00 4.600735e-01 [49,] -9.316647e+00 -3.912389e+00 [50,] -7.890508e+00 -9.316647e+00 [51,] 1.613757e+00 -7.890508e+00 [52,] -1.736845e+00 1.613757e+00 [53,] 3.046343e-01 -1.736845e+00 [54,] -9.921026e+00 3.046343e-01 [55,] -1.997308e+00 -9.921026e+00 [56,] -1.048413e+01 -1.997308e+00 [57,] -2.733878e+00 -1.048413e+01 [58,] -5.128870e+00 -2.733878e+00 [59,] -3.841798e+00 -5.128870e+00 [60,] -3.633882e+00 -3.841798e+00 [61,] -3.332085e+00 -3.633882e+00 [62,] -6.703837e+00 -3.332085e+00 [63,] -8.322809e+00 -6.703837e+00 [64,] -3.035122e+00 -8.322809e+00 [65,] -1.104261e+01 -3.035122e+00 [66,] 4.794165e+00 -1.104261e+01 [67,] -2.740594e+00 4.794165e+00 [68,] -4.397875e+00 -2.740594e+00 [69,] -6.818793e+00 -4.397875e+00 [70,] -1.392042e+01 -6.818793e+00 [71,] -8.091156e+00 -1.392042e+01 [72,] -3.944353e+00 -8.091156e+00 [73,] -5.370021e+00 -3.944353e+00 [74,] -1.919989e+01 -5.370021e+00 [75,] -7.359364e+00 -1.919989e+01 [76,] -1.148227e+01 -7.359364e+00 [77,] -6.977222e+00 -1.148227e+01 [78,] -1.273791e+01 -6.977222e+00 [79,] -1.099634e+01 -1.273791e+01 [80,] -9.118026e+00 -1.099634e+01 [81,] -2.317637e+01 -9.118026e+00 [82,] -1.638522e+01 -2.317637e+01 [83,] -1.811097e+01 -1.638522e+01 [84,] -2.131841e+01 -1.811097e+01 [85,] -9.902323e+00 -2.131841e+01 [86,] 1.946584e+00 -9.902323e+00 [87,] -9.639159e+00 1.946584e+00 [88,] -5.634773e+00 -9.639159e+00 [89,] 4.646007e+00 -5.634773e+00 [90,] 3.034166e+00 4.646007e+00 [91,] 3.269001e+00 3.034166e+00 [92,] -3.673045e+00 3.269001e+00 [93,] 3.736835e+00 -3.673045e+00 [94,] -5.371805e+00 3.736835e+00 [95,] 4.022425e+00 -5.371805e+00 [96,] 3.996791e+00 4.022425e+00 [97,] 3.095573e-02 3.996791e+00 [98,] -5.080552e+00 3.095573e-02 [99,] -1.207077e+00 -5.080552e+00 [100,] -2.989060e+00 -1.207077e+00 [101,] 7.560322e+00 -2.989060e+00 [102,] 1.057145e+00 7.560322e+00 [103,] -4.161320e+00 1.057145e+00 [104,] 6.516305e+00 -4.161320e+00 [105,] -2.933674e+00 6.516305e+00 [106,] 7.111210e+00 -2.933674e+00 [107,] 3.264884e+00 7.111210e+00 [108,] 1.336456e+00 3.264884e+00 [109,] 1.731668e+00 1.336456e+00 [110,] -2.909071e+00 1.731668e+00 [111,] -3.082950e+00 -2.909071e+00 [112,] -3.915570e+00 -3.082950e+00 [113,] 4.521706e+00 -3.915570e+00 [114,] 2.413311e+00 4.521706e+00 [115,] -7.918462e-01 2.413311e+00 [116,] 1.427120e+00 -7.918462e-01 [117,] 7.460279e+00 1.427120e+00 [118,] 4.152292e+00 7.460279e+00 [119,] 4.707465e-01 4.152292e+00 [120,] -4.265872e+00 4.707465e-01 [121,] 5.517678e+00 -4.265872e+00 [122,] -5.774741e-01 5.517678e+00 [123,] 1.080938e+01 -5.774741e-01 [124,] 6.762213e-01 1.080938e+01 [125,] -1.512378e+00 6.762213e-01 [126,] -6.220782e+00 -1.512378e+00 [127,] -1.122912e+00 -6.220782e+00 [128,] 2.676043e+00 -1.122912e+00 [129,] -7.171300e-01 2.676043e+00 [130,] 2.361654e+00 -7.171300e-01 [131,] -1.342618e+00 2.361654e+00 [132,] -9.465486e+00 -1.342618e+00 [133,] -3.445541e+00 -9.465486e+00 [134,] 4.648334e+00 -3.445541e+00 [135,] 9.144664e+00 4.648334e+00 [136,] 6.895368e+00 9.144664e+00 [137,] 4.469585e+00 6.895368e+00 [138,] 1.344915e+00 4.469585e+00 [139,] 1.194739e+01 1.344915e+00 [140,] -1.007670e+01 1.194739e+01 [141,] 2.348123e+00 -1.007670e+01 [142,] 6.563733e-01 2.348123e+00 [143,] 1.937644e+00 6.563733e-01 [144,] -2.666247e+00 1.937644e+00 [145,] 8.339364e+00 -2.666247e+00 [146,] 1.689438e+00 8.339364e+00 [147,] 4.798750e+00 1.689438e+00 [148,] -8.470525e+00 4.798750e+00 [149,] -2.297545e+00 -8.470525e+00 [150,] -5.175196e-01 -2.297545e+00 [151,] -1.634656e+00 -5.175196e-01 [152,] -2.813178e+00 -1.634656e+00 [153,] -9.719617e+00 -2.813178e+00 [154,] 3.080725e+00 -9.719617e+00 [155,] 4.395359e+00 3.080725e+00 [156,] -3.889484e+00 4.395359e+00 [157,] -3.980531e+00 -3.889484e+00 [158,] -3.674392e+00 -3.980531e+00 [159,] -3.456716e+00 -3.674392e+00 [160,] -5.604135e+00 -3.456716e+00 [161,] 1.577833e+00 -5.604135e+00 [162,] -7.819628e+00 1.577833e+00 [163,] -3.008714e+00 -7.819628e+00 [164,] -1.920941e+00 -3.008714e+00 [165,] -7.638249e+00 -1.920941e+00 [166,] 2.085900e+00 -7.638249e+00 [167,] 8.777539e+00 2.085900e+00 [168,] -7.323208e+00 8.777539e+00 [169,] 4.627454e+00 -7.323208e+00 [170,] -5.620451e+00 4.627454e+00 [171,] -4.048657e+00 -5.620451e+00 [172,] -5.619919e+00 -4.048657e+00 [173,] -3.470164e-01 -5.619919e+00 [174,] -7.596715e+00 -3.470164e-01 [175,] -3.614473e+00 -7.596715e+00 [176,] -7.984102e+00 -3.614473e+00 [177,] -7.753524e+00 -7.984102e+00 [178,] -1.557685e+00 -7.753524e+00 [179,] 3.933771e+00 -1.557685e+00 [180,] -3.831038e+00 3.933771e+00 [181,] -3.407875e+00 -3.831038e+00 [182,] 2.082319e+00 -3.407875e+00 [183,] 7.523901e-01 2.082319e+00 [184,] -1.111827e+01 7.523901e-01 [185,] -5.232714e+00 -1.111827e+01 [186,] 1.591044e+00 -5.232714e+00 [187,] -8.869788e+00 1.591044e+00 [188,] 9.362148e-01 -8.869788e+00 [189,] -4.468212e+00 9.362148e-01 [190,] -8.845299e+00 -4.468212e+00 [191,] -1.173465e+01 -8.845299e+00 [192,] -1.499360e+01 -1.173465e+01 [193,] -1.423325e+01 -1.499360e+01 [194,] -8.900840e+00 -1.423325e+01 [195,] -1.210137e+01 -8.900840e+00 [196,] -1.627754e+01 -1.210137e+01 [197,] 2.100143e+00 -1.627754e+01 [198,] -8.065534e+00 2.100143e+00 [199,] 2.553390e+01 -8.065534e+00 [200,] -1.041218e+01 2.553390e+01 [201,] -9.919521e+00 -1.041218e+01 [202,] -1.259631e+01 -9.919521e+00 [203,] -3.597365e+00 -1.259631e+01 [204,] -6.576023e+00 -3.597365e+00 [205,] -8.272962e+00 -6.576023e+00 [206,] -8.967380e+00 -8.272962e+00 [207,] 1.456099e+00 -8.967380e+00 [208,] 3.005828e+00 1.456099e+00 [209,] 2.879939e-03 3.005828e+00 [210,] -4.438637e+00 2.879939e-03 [211,] -2.059624e+00 -4.438637e+00 [212,] -2.245853e+00 -2.059624e+00 [213,] -5.297683e+00 -2.245853e+00 [214,] -7.673085e+00 -5.297683e+00 [215,] -8.026034e-01 -7.673085e+00 [216,] -5.894918e+00 -8.026034e-01 [217,] 3.909198e+00 -5.894918e+00 [218,] -6.703507e+00 3.909198e+00 [219,] -5.295797e+00 -6.703507e+00 [220,] -3.756083e+00 -5.295797e+00 [221,] -2.777192e+00 -3.756083e+00 [222,] -1.010839e+01 -2.777192e+00 [223,] -9.552847e+00 -1.010839e+01 [224,] -1.290084e+01 -9.552847e+00 [225,] -1.014784e+01 -1.290084e+01 [226,] -1.350631e+01 -1.014784e+01 [227,] -5.920137e+00 -1.350631e+01 [228,] -3.195235e+00 -5.920137e+00 [229,] -4.143248e+00 -3.195235e+00 [230,] -2.243213e+00 -4.143248e+00 [231,] -2.762337e+00 -2.243213e+00 [232,] -1.503601e+00 -2.762337e+00 [233,] -4.786170e+00 -1.503601e+00 [234,] -1.796502e-01 -4.786170e+00 [235,] -1.046490e+01 -1.796502e-01 [236,] -1.215343e+01 -1.046490e+01 [237,] -6.932109e+00 -1.215343e+01 [238,] -6.289822e-01 -6.932109e+00 [239,] -8.280547e+00 -6.289822e-01 [240,] -9.460726e+00 -8.280547e+00 [241,] -4.818488e+00 -9.460726e+00 [242,] -1.535704e+00 -4.818488e+00 [243,] -7.188073e+00 -1.535704e+00 [244,] -1.427014e+01 -7.188073e+00 [245,] 5.261674e+00 -1.427014e+01 [246,] -1.924026e+00 5.261674e+00 [247,] -4.937063e+00 -1.924026e+00 [248,] 4.904301e+00 -4.937063e+00 [249,] -2.548039e+00 4.904301e+00 [250,] 4.422112e+00 -2.548039e+00 [251,] -7.419575e+00 4.422112e+00 [252,] -8.029612e+00 -7.419575e+00 [253,] -1.336652e+00 -8.029612e+00 [254,] -3.366953e+00 -1.336652e+00 [255,] -1.366328e+00 -3.366953e+00 [256,] 4.826271e+00 -1.366328e+00 [257,] 4.229960e-01 4.826271e+00 [258,] -8.530048e+00 4.229960e-01 [259,] 1.007894e-01 -8.530048e+00 [260,] -4.072885e+00 1.007894e-01 [261,] 3.447124e+00 -4.072885e+00 [262,] -3.916064e+00 3.447124e+00 [263,] 3.357711e+00 -3.916064e+00 [264,] -6.333565e+00 3.357711e+00 [265,] -2.456749e+00 -6.333565e+00 [266,] -4.257735e+00 -2.456749e+00 [267,] 2.411166e+00 -4.257735e+00 [268,] -2.112360e-01 2.411166e+00 [269,] -8.587899e+00 -2.112360e-01 [270,] 9.409951e-02 -8.587899e+00 [271,] 1.639989e+00 9.409951e-02 [272,] 3.890707e+00 1.639989e+00 [273,] -7.797128e+00 3.890707e+00 [274,] 7.448458e+00 -7.797128e+00 [275,] -4.614510e+00 7.448458e+00 [276,] -1.861732e+00 -4.614510e+00 [277,] 1.254214e+01 -1.861732e+00 [278,] -2.320794e+00 1.254214e+01 [279,] -6.070694e+00 -2.320794e+00 [280,] 2.472612e+00 -6.070694e+00 [281,] 4.522071e-01 2.472612e+00 [282,] 8.534202e+00 4.522071e-01 [283,] -7.280054e-01 8.534202e+00 [284,] 6.310018e-01 -7.280054e-01 [285,] -5.314934e+00 6.310018e-01 [286,] -5.541811e+00 -5.314934e+00 [287,] 4.541734e+00 -5.541811e+00 [288,] -7.843056e+00 4.541734e+00 [289,] -7.795154e+00 -7.843056e+00 [290,] -9.367311e+00 -7.795154e+00 [291,] -4.954651e+00 -9.367311e+00 [292,] -6.820983e+00 -4.954651e+00 [293,] 1.813743e+00 -6.820983e+00 [294,] -5.441131e+00 1.813743e+00 [295,] -8.197118e+00 -5.441131e+00 [296,] 4.233559e+00 -8.197118e+00 [297,] 3.076145e+00 4.233559e+00 [298,] -2.778540e+00 3.076145e+00 [299,] -1.318926e+01 -2.778540e+00 [300,] -7.754824e+00 -1.318926e+01 [301,] 6.293652e-01 -7.754824e+00 [302,] 1.617146e+00 6.293652e-01 [303,] 2.324657e+00 1.617146e+00 [304,] -1.077044e-01 2.324657e+00 [305,] -5.615624e+00 -1.077044e-01 [306,] -8.644780e+00 -5.615624e+00 [307,] 1.035547e+00 -8.644780e+00 [308,] -4.491362e+00 1.035547e+00 [309,] -2.671451e+00 -4.491362e+00 [310,] -2.405719e+00 -2.671451e+00 [311,] -2.828444e+00 -2.405719e+00 [312,] -4.838494e+00 -2.828444e+00 [313,] -1.061683e+01 -4.838494e+00 [314,] -1.014508e+01 -1.061683e+01 [315,] 2.338463e-04 -1.014508e+01 [316,] -8.992978e+00 2.338463e-04 [317,] -2.774032e+00 -8.992978e+00 [318,] -5.715379e+00 -2.774032e+00 [319,] -5.478546e+00 -5.715379e+00 [320,] -4.626538e+00 -5.478546e+00 [321,] -8.264079e+00 -4.626538e+00 [322,] -6.410088e+00 -8.264079e+00 [323,] -2.022465e+01 -6.410088e+00 [324,] -2.556611e+00 -2.022465e+01 [325,] -5.250831e+00 -2.556611e+00 [326,] 3.867523e+00 -5.250831e+00 [327,] 7.295919e-02 3.867523e+00 [328,] -1.600432e+01 7.295919e-02 [329,] 6.894633e-01 -1.600432e+01 [330,] -5.814681e+00 6.894633e-01 [331,] -8.380842e+00 -5.814681e+00 [332,] -8.593540e+00 -8.380842e+00 [333,] -5.834853e+00 -8.593540e+00 [334,] -1.678961e+01 -5.834853e+00 [335,] -9.019913e+00 -1.678961e+01 [336,] -1.109807e+01 -9.019913e+00 [337,] -1.193621e+01 -1.109807e+01 [338,] -9.815999e+00 -1.193621e+01 [339,] -6.036795e+00 -9.815999e+00 [340,] -1.060875e+01 -6.036795e+00 [341,] -1.251452e+01 -1.060875e+01 [342,] -2.588106e-01 -1.251452e+01 [343,] -8.132063e+00 -2.588106e-01 [344,] -8.287043e+00 -8.132063e+00 [345,] -1.360723e+01 -8.287043e+00 [346,] -8.819109e+00 -1.360723e+01 [347,] -9.301126e+00 -8.819109e+00 [348,] -9.136712e+00 -9.301126e+00 [349,] 4.814530e+00 -9.136712e+00 [350,] -1.334495e+01 4.814530e+00 [351,] -6.040788e+00 -1.334495e+01 [352,] -7.711183e-02 -6.040788e+00 [353,] 1.721691e+00 -7.711183e-02 [354,] -8.745563e+00 1.721691e+00 [355,] -4.914882e+00 -8.745563e+00 [356,] -7.631854e+00 -4.914882e+00 [357,] -1.429410e+01 -7.631854e+00 [358,] -1.620007e+01 -1.429410e+01 [359,] -6.825736e+00 -1.620007e+01 [360,] -1.091408e+01 -6.825736e+00 [361,] -1.180369e+01 -1.091408e+01 [362,] 3.418831e+00 -1.180369e+01 [363,] -9.798002e+00 3.418831e+00 [364,] -1.707612e+01 -9.798002e+00 [365,] -1.319164e+01 -1.707612e+01 [366,] -2.288007e+00 -1.319164e+01 [367,] -8.005386e+00 -2.288007e+00 [368,] -8.516113e+00 -8.005386e+00 [369,] -6.844228e+00 -8.516113e+00 [370,] -8.081791e+00 -6.844228e+00 [371,] -2.062056e+00 -8.081791e+00 [372,] -9.027292e+00 -2.062056e+00 [373,] -1.225732e+00 -9.027292e+00 [374,] -2.933471e+00 -1.225732e+00 [375,] -4.665275e+00 -2.933471e+00 [376,] -6.532304e+00 -4.665275e+00 [377,] -8.132389e+00 -6.532304e+00 [378,] -5.366831e+00 -8.132389e+00 [379,] -4.556410e+00 -5.366831e+00 [380,] 1.999811e+00 -4.556410e+00 [381,] -2.422642e-01 1.999811e+00 [382,] -2.238471e-01 -2.422642e-01 [383,] -3.961175e+00 -2.238471e-01 [384,] -1.233876e+01 -3.961175e+00 [385,] -4.257692e+00 -1.233876e+01 [386,] -9.856868e+00 -4.257692e+00 [387,] -9.084125e-02 -9.856868e+00 [388,] -5.093623e-01 -9.084125e-02 [389,] -5.444040e+00 -5.093623e-01 [390,] 5.535112e+00 -5.444040e+00 [391,] 6.730607e-01 5.535112e+00 [392,] 1.256895e+00 6.730607e-01 [393,] -4.322303e+00 1.256895e+00 [394,] -1.094502e+01 -4.322303e+00 [395,] -9.368908e+00 -1.094502e+01 [396,] -4.020970e+00 -9.368908e+00 [397,] -1.030090e+01 -4.020970e+00 [398,] -1.349618e+01 -1.030090e+01 [399,] -1.087367e+01 -1.349618e+01 [400,] -8.336932e+00 -1.087367e+01 [401,] -3.601847e+00 -8.336932e+00 [402,] 1.414664e+00 -3.601847e+00 [403,] -1.183180e+01 1.414664e+00 [404,] -5.663892e-01 -1.183180e+01 [405,] -5.070947e+00 -5.663892e-01 [406,] -1.423151e+00 -5.070947e+00 [407,] -1.232879e+01 -1.423151e+00 [408,] -1.063743e+01 -1.232879e+01 [409,] -3.524293e+00 -1.063743e+01 [410,] -2.186279e+00 -3.524293e+00 [411,] -3.716105e+00 -2.186279e+00 [412,] -1.607130e+00 -3.716105e+00 [413,] -9.029069e+00 -1.607130e+00 [414,] -1.577982e+01 -9.029069e+00 [415,] -4.269073e+00 -1.577982e+01 [416,] -7.558887e+00 -4.269073e+00 [417,] -2.457378e+00 -7.558887e+00 [418,] -9.876012e+00 -2.457378e+00 [419,] -1.630162e+00 -9.876012e+00 [420,] -1.021682e+01 -1.630162e+00 [421,] 3.695981e-01 -1.021682e+01 [422,] -1.937202e+01 3.695981e-01 [423,] -1.021490e-01 -1.937202e+01 [424,] -6.506590e+00 -1.021490e-01 [425,] -1.477353e+01 -6.506590e+00 [426,] -1.368036e+01 -1.477353e+01 [427,] -1.973588e+00 -1.368036e+01 [428,] 2.806230e+00 -1.973588e+00 [429,] -9.361183e+00 2.806230e+00 [430,] -5.134870e+00 -9.361183e+00 [431,] -1.237298e+01 -5.134870e+00 [432,] -6.193866e+00 -1.237298e+01 [433,] -4.016571e+00 -6.193866e+00 [434,] 2.059318e+00 -4.016571e+00 [435,] -8.418025e+00 2.059318e+00 [436,] -5.178387e+00 -8.418025e+00 [437,] -1.132184e+01 -5.178387e+00 [438,] -1.571835e+01 -1.132184e+01 [439,] -5.691437e+00 -1.571835e+01 [440,] -5.282922e+00 -5.691437e+00 [441,] -9.454454e+00 -5.282922e+00 [442,] -7.422582e+00 -9.454454e+00 [443,] -8.747528e+00 -7.422582e+00 [444,] -2.188783e+00 -8.747528e+00 [445,] -1.551824e+01 -2.188783e+00 [446,] -7.676564e+00 -1.551824e+01 [447,] -8.858174e+00 -7.676564e+00 [448,] -1.487991e+01 -8.858174e+00 [449,] -2.707192e+00 -1.487991e+01 [450,] -2.245388e+00 -2.707192e+00 [451,] 1.071900e+00 -2.245388e+00 [452,] 9.689288e-01 1.071900e+00 [453,] 3.349447e+00 9.689288e-01 [454,] -8.677832e+00 3.349447e+00 [455,] -6.190921e+00 -8.677832e+00 [456,] -6.200216e-01 -6.190921e+00 [457,] 1.396973e+00 -6.200216e-01 [458,] -9.713814e+00 1.396973e+00 [459,] -7.314318e+00 -9.713814e+00 [460,] -1.174711e+00 -7.314318e+00 [461,] -2.518591e+00 -1.174711e+00 [462,] -6.154023e+00 -2.518591e+00 [463,] -1.855996e-01 -6.154023e+00 [464,] -1.814706e+00 -1.855996e-01 [465,] -2.358454e+00 -1.814706e+00 [466,] -1.227286e+01 -2.358454e+00 [467,] -1.482947e+00 -1.227286e+01 [468,] -1.028202e+01 -1.482947e+00 [469,] -3.374203e+00 -1.028202e+01 [470,] 4.825680e+00 -3.374203e+00 [471,] -6.554237e+00 4.825680e+00 [472,] -2.237483e+00 -6.554237e+00 [473,] -9.190129e+00 -2.237483e+00 [474,] -5.737986e-01 -9.190129e+00 [475,] -2.250249e+00 -5.737986e-01 [476,] -4.495156e+00 -2.250249e+00 [477,] -4.597812e+00 -4.495156e+00 [478,] 3.967008e+00 -4.597812e+00 [479,] -1.085722e+01 3.967008e+00 [480,] -2.912574e+00 -1.085722e+01 [481,] -2.256643e+00 -2.912574e+00 [482,] -9.772051e+00 -2.256643e+00 [483,] -1.006552e+01 -9.772051e+00 [484,] -1.466945e+01 -1.006552e+01 [485,] -9.020717e+00 -1.466945e+01 [486,] -6.783137e+00 -9.020717e+00 [487,] -1.413017e+01 -6.783137e+00 [488,] -9.217490e+00 -1.413017e+01 [489,] -9.034729e+00 -9.217490e+00 [490,] -1.169150e+01 -9.034729e+00 [491,] 2.161951e-01 -1.169150e+01 [492,] -3.275268e+00 2.161951e-01 [493,] -1.087961e+01 -3.275268e+00 [494,] 6.979003e-01 -1.087961e+01 [495,] -1.173240e+01 6.979003e-01 [496,] -1.057020e+01 -1.173240e+01 [497,] -9.851838e+00 -1.057020e+01 [498,] 1.330189e+01 -9.851838e+00 [499,] -4.302809e+00 1.330189e+01 [500,] -1.635559e+00 -4.302809e+00 [501,] -7.868065e+00 -1.635559e+00 [502,] -6.343591e+00 -7.868065e+00 [503,] -2.713257e+00 -6.343591e+00 [504,] -1.806510e+00 -2.713257e+00 [505,] -8.745664e+00 -1.806510e+00 [506,] 1.551508e+00 -8.745664e+00 [507,] 1.739173e+00 1.551508e+00 [508,] 5.409489e+00 1.739173e+00 [509,] -1.038762e+01 5.409489e+00 [510,] 8.758784e+02 -1.038762e+01 [511,] -2.040067e-01 8.758784e+02 [512,] -2.185599e+00 -2.040067e-01 [513,] -1.341008e+01 -2.185599e+00 [514,] -5.486828e+00 -1.341008e+01 [515,] -1.131193e+01 -5.486828e+00 [516,] -1.255453e+01 -1.131193e+01 [517,] -4.422624e+00 -1.255453e+01 [518,] -3.354258e+00 -4.422624e+00 [519,] -6.927837e+00 -3.354258e+00 [520,] -3.021482e+00 -6.927837e+00 [521,] -6.446570e+00 -3.021482e+00 [522,] -1.488881e+01 -6.446570e+00 [523,] -1.608805e+01 -1.488881e+01 [524,] -6.036465e+00 -1.608805e+01 [525,] -6.716511e+00 -6.036465e+00 [526,] -3.579940e+00 -6.716511e+00 [527,] -8.636345e+00 -3.579940e+00 [528,] -1.555820e+01 -8.636345e+00 [529,] -2.844068e+00 -1.555820e+01 [530,] -1.568663e+01 -2.844068e+00 [531,] 7.604625e+02 -1.568663e+01 [532,] 8.764713e+02 7.604625e+02 [533,] -6.374334e+00 8.764713e+02 [534,] -4.040232e+00 -6.374334e+00 [535,] -1.266852e+01 -4.040232e+00 [536,] -6.669758e+00 -1.266852e+01 [537,] -1.159025e+01 -6.669758e+00 [538,] -5.567872e+00 -1.159025e+01 [539,] -1.433900e+01 -5.567872e+00 [540,] -4.584389e+00 -1.433900e+01 [541,] -1.528106e+01 -4.584389e+00 [542,] -2.295174e+00 -1.528106e+01 [543,] -7.161032e+00 -2.295174e+00 [544,] -7.668840e-01 -7.161032e+00 [545,] -6.203831e+00 -7.668840e-01 [546,] -1.098963e+01 -6.203831e+00 [547,] -1.010455e+01 -1.098963e+01 [548,] -8.494776e+00 -1.010455e+01 [549,] -1.197824e+01 -8.494776e+00 [550,] -1.257357e+01 -1.197824e+01 [551,] -1.081448e+01 -1.257357e+01 [552,] -3.203562e+00 -1.081448e+01 [553,] -6.306076e+00 -3.203562e+00 [554,] -9.063120e+00 -6.306076e+00 [555,] -6.153311e+00 -9.063120e+00 [556,] -5.932945e+00 -6.153311e+00 [557,] -4.971402e+00 -5.932945e+00 [558,] -1.140417e+01 -4.971402e+00 [559,] -1.095120e+01 -1.140417e+01 [560,] -7.335587e+00 -1.095120e+01 [561,] -1.162483e+01 -7.335587e+00 [562,] -3.547614e+00 -1.162483e+01 [563,] -1.027371e+01 -3.547614e+00 [564,] -6.259756e+00 -1.027371e+01 [565,] -5.723048e+00 -6.259756e+00 [566,] -9.172901e+00 -5.723048e+00 [567,] -3.756818e+00 -9.172901e+00 [568,] -3.759076e+00 -3.756818e+00 [569,] -1.042397e+00 -3.759076e+00 [570,] -8.032350e+00 -1.042397e+00 [571,] -8.323509e+00 -8.032350e+00 [572,] -8.552383e+00 -8.323509e+00 [573,] -8.905003e+00 -8.552383e+00 [574,] -3.192250e+00 -8.905003e+00 [575,] 1.258708e+00 -3.192250e+00 [576,] -9.910202e+00 1.258708e+00 [577,] 7.485737e+00 -9.910202e+00 [578,] -1.350696e+01 7.485737e+00 [579,] 3.018549e+00 -1.350696e+01 [580,] -5.267373e+00 3.018549e+00 [581,] -2.533336e+00 -5.267373e+00 [582,] -3.300970e+00 -2.533336e+00 [583,] -6.979721e+00 -3.300970e+00 [584,] -1.360075e+01 -6.979721e+00 [585,] -7.169105e+00 -1.360075e+01 [586,] -1.447020e+01 -7.169105e+00 [587,] -5.456367e+00 -1.447020e+01 [588,] -2.574610e+00 -5.456367e+00 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2.013082e+01 8.464328e+00 2 3.684305e+00 2.013082e+01 3 -1.017926e+01 3.684305e+00 4 5.632891e+00 -1.017926e+01 5 4.672424e+00 5.632891e+00 6 -2.113670e+00 4.672424e+00 7 8.380918e+00 -2.113670e+00 8 -1.209164e+01 8.380918e+00 9 -2.332091e+00 -1.209164e+01 10 -8.037254e-01 -2.332091e+00 11 -5.221019e-01 -8.037254e-01 12 7.948979e+00 -5.221019e-01 13 1.741616e+00 7.948979e+00 14 -6.161324e+00 1.741616e+00 15 -6.081434e+00 -6.161324e+00 16 6.888138e+00 -6.081434e+00 17 -1.255763e+00 6.888138e+00 18 3.709982e+00 -1.255763e+00 19 3.117333e+00 3.709982e+00 20 2.488585e+00 3.117333e+00 21 -1.818748e+00 2.488585e+00 22 7.012078e+00 -1.818748e+00 23 -3.123524e+00 7.012078e+00 24 2.948687e-01 -3.123524e+00 25 4.418424e+00 2.948687e-01 26 2.090478e+00 4.418424e+00 27 -7.978551e+00 2.090478e+00 28 6.575067e+00 -7.978551e+00 29 -6.642968e-01 6.575067e+00 30 2.863945e+00 -6.642968e-01 31 9.277274e+00 2.863945e+00 32 -2.024807e+00 9.277274e+00 33 -5.042074e+00 -2.024807e+00 34 1.373882e-01 -5.042074e+00 35 -6.650105e+00 1.373882e-01 36 4.523806e+00 -6.650105e+00 37 -7.481487e-01 4.523806e+00 38 -2.250088e+00 -7.481487e-01 39 2.396932e+00 -2.250088e+00 40 3.154763e+00 2.396932e+00 41 3.815728e+00 3.154763e+00 42 -8.484061e+00 3.815728e+00 43 1.773405e-01 -8.484061e+00 44 -1.019999e+01 1.773405e-01 45 -3.167972e+00 -1.019999e+01 46 -1.481587e-01 -3.167972e+00 47 4.600735e-01 -1.481587e-01 48 -3.912389e+00 4.600735e-01 49 -9.316647e+00 -3.912389e+00 50 -7.890508e+00 -9.316647e+00 51 1.613757e+00 -7.890508e+00 52 -1.736845e+00 1.613757e+00 53 3.046343e-01 -1.736845e+00 54 -9.921026e+00 3.046343e-01 55 -1.997308e+00 -9.921026e+00 56 -1.048413e+01 -1.997308e+00 57 -2.733878e+00 -1.048413e+01 58 -5.128870e+00 -2.733878e+00 59 -3.841798e+00 -5.128870e+00 60 -3.633882e+00 -3.841798e+00 61 -3.332085e+00 -3.633882e+00 62 -6.703837e+00 -3.332085e+00 63 -8.322809e+00 -6.703837e+00 64 -3.035122e+00 -8.322809e+00 65 -1.104261e+01 -3.035122e+00 66 4.794165e+00 -1.104261e+01 67 -2.740594e+00 4.794165e+00 68 -4.397875e+00 -2.740594e+00 69 -6.818793e+00 -4.397875e+00 70 -1.392042e+01 -6.818793e+00 71 -8.091156e+00 -1.392042e+01 72 -3.944353e+00 -8.091156e+00 73 -5.370021e+00 -3.944353e+00 74 -1.919989e+01 -5.370021e+00 75 -7.359364e+00 -1.919989e+01 76 -1.148227e+01 -7.359364e+00 77 -6.977222e+00 -1.148227e+01 78 -1.273791e+01 -6.977222e+00 79 -1.099634e+01 -1.273791e+01 80 -9.118026e+00 -1.099634e+01 81 -2.317637e+01 -9.118026e+00 82 -1.638522e+01 -2.317637e+01 83 -1.811097e+01 -1.638522e+01 84 -2.131841e+01 -1.811097e+01 85 -9.902323e+00 -2.131841e+01 86 1.946584e+00 -9.902323e+00 87 -9.639159e+00 1.946584e+00 88 -5.634773e+00 -9.639159e+00 89 4.646007e+00 -5.634773e+00 90 3.034166e+00 4.646007e+00 91 3.269001e+00 3.034166e+00 92 -3.673045e+00 3.269001e+00 93 3.736835e+00 -3.673045e+00 94 -5.371805e+00 3.736835e+00 95 4.022425e+00 -5.371805e+00 96 3.996791e+00 4.022425e+00 97 3.095573e-02 3.996791e+00 98 -5.080552e+00 3.095573e-02 99 -1.207077e+00 -5.080552e+00 100 -2.989060e+00 -1.207077e+00 101 7.560322e+00 -2.989060e+00 102 1.057145e+00 7.560322e+00 103 -4.161320e+00 1.057145e+00 104 6.516305e+00 -4.161320e+00 105 -2.933674e+00 6.516305e+00 106 7.111210e+00 -2.933674e+00 107 3.264884e+00 7.111210e+00 108 1.336456e+00 3.264884e+00 109 1.731668e+00 1.336456e+00 110 -2.909071e+00 1.731668e+00 111 -3.082950e+00 -2.909071e+00 112 -3.915570e+00 -3.082950e+00 113 4.521706e+00 -3.915570e+00 114 2.413311e+00 4.521706e+00 115 -7.918462e-01 2.413311e+00 116 1.427120e+00 -7.918462e-01 117 7.460279e+00 1.427120e+00 118 4.152292e+00 7.460279e+00 119 4.707465e-01 4.152292e+00 120 -4.265872e+00 4.707465e-01 121 5.517678e+00 -4.265872e+00 122 -5.774741e-01 5.517678e+00 123 1.080938e+01 -5.774741e-01 124 6.762213e-01 1.080938e+01 125 -1.512378e+00 6.762213e-01 126 -6.220782e+00 -1.512378e+00 127 -1.122912e+00 -6.220782e+00 128 2.676043e+00 -1.122912e+00 129 -7.171300e-01 2.676043e+00 130 2.361654e+00 -7.171300e-01 131 -1.342618e+00 2.361654e+00 132 -9.465486e+00 -1.342618e+00 133 -3.445541e+00 -9.465486e+00 134 4.648334e+00 -3.445541e+00 135 9.144664e+00 4.648334e+00 136 6.895368e+00 9.144664e+00 137 4.469585e+00 6.895368e+00 138 1.344915e+00 4.469585e+00 139 1.194739e+01 1.344915e+00 140 -1.007670e+01 1.194739e+01 141 2.348123e+00 -1.007670e+01 142 6.563733e-01 2.348123e+00 143 1.937644e+00 6.563733e-01 144 -2.666247e+00 1.937644e+00 145 8.339364e+00 -2.666247e+00 146 1.689438e+00 8.339364e+00 147 4.798750e+00 1.689438e+00 148 -8.470525e+00 4.798750e+00 149 -2.297545e+00 -8.470525e+00 150 -5.175196e-01 -2.297545e+00 151 -1.634656e+00 -5.175196e-01 152 -2.813178e+00 -1.634656e+00 153 -9.719617e+00 -2.813178e+00 154 3.080725e+00 -9.719617e+00 155 4.395359e+00 3.080725e+00 156 -3.889484e+00 4.395359e+00 157 -3.980531e+00 -3.889484e+00 158 -3.674392e+00 -3.980531e+00 159 -3.456716e+00 -3.674392e+00 160 -5.604135e+00 -3.456716e+00 161 1.577833e+00 -5.604135e+00 162 -7.819628e+00 1.577833e+00 163 -3.008714e+00 -7.819628e+00 164 -1.920941e+00 -3.008714e+00 165 -7.638249e+00 -1.920941e+00 166 2.085900e+00 -7.638249e+00 167 8.777539e+00 2.085900e+00 168 -7.323208e+00 8.777539e+00 169 4.627454e+00 -7.323208e+00 170 -5.620451e+00 4.627454e+00 171 -4.048657e+00 -5.620451e+00 172 -5.619919e+00 -4.048657e+00 173 -3.470164e-01 -5.619919e+00 174 -7.596715e+00 -3.470164e-01 175 -3.614473e+00 -7.596715e+00 176 -7.984102e+00 -3.614473e+00 177 -7.753524e+00 -7.984102e+00 178 -1.557685e+00 -7.753524e+00 179 3.933771e+00 -1.557685e+00 180 -3.831038e+00 3.933771e+00 181 -3.407875e+00 -3.831038e+00 182 2.082319e+00 -3.407875e+00 183 7.523901e-01 2.082319e+00 184 -1.111827e+01 7.523901e-01 185 -5.232714e+00 -1.111827e+01 186 1.591044e+00 -5.232714e+00 187 -8.869788e+00 1.591044e+00 188 9.362148e-01 -8.869788e+00 189 -4.468212e+00 9.362148e-01 190 -8.845299e+00 -4.468212e+00 191 -1.173465e+01 -8.845299e+00 192 -1.499360e+01 -1.173465e+01 193 -1.423325e+01 -1.499360e+01 194 -8.900840e+00 -1.423325e+01 195 -1.210137e+01 -8.900840e+00 196 -1.627754e+01 -1.210137e+01 197 2.100143e+00 -1.627754e+01 198 -8.065534e+00 2.100143e+00 199 2.553390e+01 -8.065534e+00 200 -1.041218e+01 2.553390e+01 201 -9.919521e+00 -1.041218e+01 202 -1.259631e+01 -9.919521e+00 203 -3.597365e+00 -1.259631e+01 204 -6.576023e+00 -3.597365e+00 205 -8.272962e+00 -6.576023e+00 206 -8.967380e+00 -8.272962e+00 207 1.456099e+00 -8.967380e+00 208 3.005828e+00 1.456099e+00 209 2.879939e-03 3.005828e+00 210 -4.438637e+00 2.879939e-03 211 -2.059624e+00 -4.438637e+00 212 -2.245853e+00 -2.059624e+00 213 -5.297683e+00 -2.245853e+00 214 -7.673085e+00 -5.297683e+00 215 -8.026034e-01 -7.673085e+00 216 -5.894918e+00 -8.026034e-01 217 3.909198e+00 -5.894918e+00 218 -6.703507e+00 3.909198e+00 219 -5.295797e+00 -6.703507e+00 220 -3.756083e+00 -5.295797e+00 221 -2.777192e+00 -3.756083e+00 222 -1.010839e+01 -2.777192e+00 223 -9.552847e+00 -1.010839e+01 224 -1.290084e+01 -9.552847e+00 225 -1.014784e+01 -1.290084e+01 226 -1.350631e+01 -1.014784e+01 227 -5.920137e+00 -1.350631e+01 228 -3.195235e+00 -5.920137e+00 229 -4.143248e+00 -3.195235e+00 230 -2.243213e+00 -4.143248e+00 231 -2.762337e+00 -2.243213e+00 232 -1.503601e+00 -2.762337e+00 233 -4.786170e+00 -1.503601e+00 234 -1.796502e-01 -4.786170e+00 235 -1.046490e+01 -1.796502e-01 236 -1.215343e+01 -1.046490e+01 237 -6.932109e+00 -1.215343e+01 238 -6.289822e-01 -6.932109e+00 239 -8.280547e+00 -6.289822e-01 240 -9.460726e+00 -8.280547e+00 241 -4.818488e+00 -9.460726e+00 242 -1.535704e+00 -4.818488e+00 243 -7.188073e+00 -1.535704e+00 244 -1.427014e+01 -7.188073e+00 245 5.261674e+00 -1.427014e+01 246 -1.924026e+00 5.261674e+00 247 -4.937063e+00 -1.924026e+00 248 4.904301e+00 -4.937063e+00 249 -2.548039e+00 4.904301e+00 250 4.422112e+00 -2.548039e+00 251 -7.419575e+00 4.422112e+00 252 -8.029612e+00 -7.419575e+00 253 -1.336652e+00 -8.029612e+00 254 -3.366953e+00 -1.336652e+00 255 -1.366328e+00 -3.366953e+00 256 4.826271e+00 -1.366328e+00 257 4.229960e-01 4.826271e+00 258 -8.530048e+00 4.229960e-01 259 1.007894e-01 -8.530048e+00 260 -4.072885e+00 1.007894e-01 261 3.447124e+00 -4.072885e+00 262 -3.916064e+00 3.447124e+00 263 3.357711e+00 -3.916064e+00 264 -6.333565e+00 3.357711e+00 265 -2.456749e+00 -6.333565e+00 266 -4.257735e+00 -2.456749e+00 267 2.411166e+00 -4.257735e+00 268 -2.112360e-01 2.411166e+00 269 -8.587899e+00 -2.112360e-01 270 9.409951e-02 -8.587899e+00 271 1.639989e+00 9.409951e-02 272 3.890707e+00 1.639989e+00 273 -7.797128e+00 3.890707e+00 274 7.448458e+00 -7.797128e+00 275 -4.614510e+00 7.448458e+00 276 -1.861732e+00 -4.614510e+00 277 1.254214e+01 -1.861732e+00 278 -2.320794e+00 1.254214e+01 279 -6.070694e+00 -2.320794e+00 280 2.472612e+00 -6.070694e+00 281 4.522071e-01 2.472612e+00 282 8.534202e+00 4.522071e-01 283 -7.280054e-01 8.534202e+00 284 6.310018e-01 -7.280054e-01 285 -5.314934e+00 6.310018e-01 286 -5.541811e+00 -5.314934e+00 287 4.541734e+00 -5.541811e+00 288 -7.843056e+00 4.541734e+00 289 -7.795154e+00 -7.843056e+00 290 -9.367311e+00 -7.795154e+00 291 -4.954651e+00 -9.367311e+00 292 -6.820983e+00 -4.954651e+00 293 1.813743e+00 -6.820983e+00 294 -5.441131e+00 1.813743e+00 295 -8.197118e+00 -5.441131e+00 296 4.233559e+00 -8.197118e+00 297 3.076145e+00 4.233559e+00 298 -2.778540e+00 3.076145e+00 299 -1.318926e+01 -2.778540e+00 300 -7.754824e+00 -1.318926e+01 301 6.293652e-01 -7.754824e+00 302 1.617146e+00 6.293652e-01 303 2.324657e+00 1.617146e+00 304 -1.077044e-01 2.324657e+00 305 -5.615624e+00 -1.077044e-01 306 -8.644780e+00 -5.615624e+00 307 1.035547e+00 -8.644780e+00 308 -4.491362e+00 1.035547e+00 309 -2.671451e+00 -4.491362e+00 310 -2.405719e+00 -2.671451e+00 311 -2.828444e+00 -2.405719e+00 312 -4.838494e+00 -2.828444e+00 313 -1.061683e+01 -4.838494e+00 314 -1.014508e+01 -1.061683e+01 315 2.338463e-04 -1.014508e+01 316 -8.992978e+00 2.338463e-04 317 -2.774032e+00 -8.992978e+00 318 -5.715379e+00 -2.774032e+00 319 -5.478546e+00 -5.715379e+00 320 -4.626538e+00 -5.478546e+00 321 -8.264079e+00 -4.626538e+00 322 -6.410088e+00 -8.264079e+00 323 -2.022465e+01 -6.410088e+00 324 -2.556611e+00 -2.022465e+01 325 -5.250831e+00 -2.556611e+00 326 3.867523e+00 -5.250831e+00 327 7.295919e-02 3.867523e+00 328 -1.600432e+01 7.295919e-02 329 6.894633e-01 -1.600432e+01 330 -5.814681e+00 6.894633e-01 331 -8.380842e+00 -5.814681e+00 332 -8.593540e+00 -8.380842e+00 333 -5.834853e+00 -8.593540e+00 334 -1.678961e+01 -5.834853e+00 335 -9.019913e+00 -1.678961e+01 336 -1.109807e+01 -9.019913e+00 337 -1.193621e+01 -1.109807e+01 338 -9.815999e+00 -1.193621e+01 339 -6.036795e+00 -9.815999e+00 340 -1.060875e+01 -6.036795e+00 341 -1.251452e+01 -1.060875e+01 342 -2.588106e-01 -1.251452e+01 343 -8.132063e+00 -2.588106e-01 344 -8.287043e+00 -8.132063e+00 345 -1.360723e+01 -8.287043e+00 346 -8.819109e+00 -1.360723e+01 347 -9.301126e+00 -8.819109e+00 348 -9.136712e+00 -9.301126e+00 349 4.814530e+00 -9.136712e+00 350 -1.334495e+01 4.814530e+00 351 -6.040788e+00 -1.334495e+01 352 -7.711183e-02 -6.040788e+00 353 1.721691e+00 -7.711183e-02 354 -8.745563e+00 1.721691e+00 355 -4.914882e+00 -8.745563e+00 356 -7.631854e+00 -4.914882e+00 357 -1.429410e+01 -7.631854e+00 358 -1.620007e+01 -1.429410e+01 359 -6.825736e+00 -1.620007e+01 360 -1.091408e+01 -6.825736e+00 361 -1.180369e+01 -1.091408e+01 362 3.418831e+00 -1.180369e+01 363 -9.798002e+00 3.418831e+00 364 -1.707612e+01 -9.798002e+00 365 -1.319164e+01 -1.707612e+01 366 -2.288007e+00 -1.319164e+01 367 -8.005386e+00 -2.288007e+00 368 -8.516113e+00 -8.005386e+00 369 -6.844228e+00 -8.516113e+00 370 -8.081791e+00 -6.844228e+00 371 -2.062056e+00 -8.081791e+00 372 -9.027292e+00 -2.062056e+00 373 -1.225732e+00 -9.027292e+00 374 -2.933471e+00 -1.225732e+00 375 -4.665275e+00 -2.933471e+00 376 -6.532304e+00 -4.665275e+00 377 -8.132389e+00 -6.532304e+00 378 -5.366831e+00 -8.132389e+00 379 -4.556410e+00 -5.366831e+00 380 1.999811e+00 -4.556410e+00 381 -2.422642e-01 1.999811e+00 382 -2.238471e-01 -2.422642e-01 383 -3.961175e+00 -2.238471e-01 384 -1.233876e+01 -3.961175e+00 385 -4.257692e+00 -1.233876e+01 386 -9.856868e+00 -4.257692e+00 387 -9.084125e-02 -9.856868e+00 388 -5.093623e-01 -9.084125e-02 389 -5.444040e+00 -5.093623e-01 390 5.535112e+00 -5.444040e+00 391 6.730607e-01 5.535112e+00 392 1.256895e+00 6.730607e-01 393 -4.322303e+00 1.256895e+00 394 -1.094502e+01 -4.322303e+00 395 -9.368908e+00 -1.094502e+01 396 -4.020970e+00 -9.368908e+00 397 -1.030090e+01 -4.020970e+00 398 -1.349618e+01 -1.030090e+01 399 -1.087367e+01 -1.349618e+01 400 -8.336932e+00 -1.087367e+01 401 -3.601847e+00 -8.336932e+00 402 1.414664e+00 -3.601847e+00 403 -1.183180e+01 1.414664e+00 404 -5.663892e-01 -1.183180e+01 405 -5.070947e+00 -5.663892e-01 406 -1.423151e+00 -5.070947e+00 407 -1.232879e+01 -1.423151e+00 408 -1.063743e+01 -1.232879e+01 409 -3.524293e+00 -1.063743e+01 410 -2.186279e+00 -3.524293e+00 411 -3.716105e+00 -2.186279e+00 412 -1.607130e+00 -3.716105e+00 413 -9.029069e+00 -1.607130e+00 414 -1.577982e+01 -9.029069e+00 415 -4.269073e+00 -1.577982e+01 416 -7.558887e+00 -4.269073e+00 417 -2.457378e+00 -7.558887e+00 418 -9.876012e+00 -2.457378e+00 419 -1.630162e+00 -9.876012e+00 420 -1.021682e+01 -1.630162e+00 421 3.695981e-01 -1.021682e+01 422 -1.937202e+01 3.695981e-01 423 -1.021490e-01 -1.937202e+01 424 -6.506590e+00 -1.021490e-01 425 -1.477353e+01 -6.506590e+00 426 -1.368036e+01 -1.477353e+01 427 -1.973588e+00 -1.368036e+01 428 2.806230e+00 -1.973588e+00 429 -9.361183e+00 2.806230e+00 430 -5.134870e+00 -9.361183e+00 431 -1.237298e+01 -5.134870e+00 432 -6.193866e+00 -1.237298e+01 433 -4.016571e+00 -6.193866e+00 434 2.059318e+00 -4.016571e+00 435 -8.418025e+00 2.059318e+00 436 -5.178387e+00 -8.418025e+00 437 -1.132184e+01 -5.178387e+00 438 -1.571835e+01 -1.132184e+01 439 -5.691437e+00 -1.571835e+01 440 -5.282922e+00 -5.691437e+00 441 -9.454454e+00 -5.282922e+00 442 -7.422582e+00 -9.454454e+00 443 -8.747528e+00 -7.422582e+00 444 -2.188783e+00 -8.747528e+00 445 -1.551824e+01 -2.188783e+00 446 -7.676564e+00 -1.551824e+01 447 -8.858174e+00 -7.676564e+00 448 -1.487991e+01 -8.858174e+00 449 -2.707192e+00 -1.487991e+01 450 -2.245388e+00 -2.707192e+00 451 1.071900e+00 -2.245388e+00 452 9.689288e-01 1.071900e+00 453 3.349447e+00 9.689288e-01 454 -8.677832e+00 3.349447e+00 455 -6.190921e+00 -8.677832e+00 456 -6.200216e-01 -6.190921e+00 457 1.396973e+00 -6.200216e-01 458 -9.713814e+00 1.396973e+00 459 -7.314318e+00 -9.713814e+00 460 -1.174711e+00 -7.314318e+00 461 -2.518591e+00 -1.174711e+00 462 -6.154023e+00 -2.518591e+00 463 -1.855996e-01 -6.154023e+00 464 -1.814706e+00 -1.855996e-01 465 -2.358454e+00 -1.814706e+00 466 -1.227286e+01 -2.358454e+00 467 -1.482947e+00 -1.227286e+01 468 -1.028202e+01 -1.482947e+00 469 -3.374203e+00 -1.028202e+01 470 4.825680e+00 -3.374203e+00 471 -6.554237e+00 4.825680e+00 472 -2.237483e+00 -6.554237e+00 473 -9.190129e+00 -2.237483e+00 474 -5.737986e-01 -9.190129e+00 475 -2.250249e+00 -5.737986e-01 476 -4.495156e+00 -2.250249e+00 477 -4.597812e+00 -4.495156e+00 478 3.967008e+00 -4.597812e+00 479 -1.085722e+01 3.967008e+00 480 -2.912574e+00 -1.085722e+01 481 -2.256643e+00 -2.912574e+00 482 -9.772051e+00 -2.256643e+00 483 -1.006552e+01 -9.772051e+00 484 -1.466945e+01 -1.006552e+01 485 -9.020717e+00 -1.466945e+01 486 -6.783137e+00 -9.020717e+00 487 -1.413017e+01 -6.783137e+00 488 -9.217490e+00 -1.413017e+01 489 -9.034729e+00 -9.217490e+00 490 -1.169150e+01 -9.034729e+00 491 2.161951e-01 -1.169150e+01 492 -3.275268e+00 2.161951e-01 493 -1.087961e+01 -3.275268e+00 494 6.979003e-01 -1.087961e+01 495 -1.173240e+01 6.979003e-01 496 -1.057020e+01 -1.173240e+01 497 -9.851838e+00 -1.057020e+01 498 1.330189e+01 -9.851838e+00 499 -4.302809e+00 1.330189e+01 500 -1.635559e+00 -4.302809e+00 501 -7.868065e+00 -1.635559e+00 502 -6.343591e+00 -7.868065e+00 503 -2.713257e+00 -6.343591e+00 504 -1.806510e+00 -2.713257e+00 505 -8.745664e+00 -1.806510e+00 506 1.551508e+00 -8.745664e+00 507 1.739173e+00 1.551508e+00 508 5.409489e+00 1.739173e+00 509 -1.038762e+01 5.409489e+00 510 8.758784e+02 -1.038762e+01 511 -2.040067e-01 8.758784e+02 512 -2.185599e+00 -2.040067e-01 513 -1.341008e+01 -2.185599e+00 514 -5.486828e+00 -1.341008e+01 515 -1.131193e+01 -5.486828e+00 516 -1.255453e+01 -1.131193e+01 517 -4.422624e+00 -1.255453e+01 518 -3.354258e+00 -4.422624e+00 519 -6.927837e+00 -3.354258e+00 520 -3.021482e+00 -6.927837e+00 521 -6.446570e+00 -3.021482e+00 522 -1.488881e+01 -6.446570e+00 523 -1.608805e+01 -1.488881e+01 524 -6.036465e+00 -1.608805e+01 525 -6.716511e+00 -6.036465e+00 526 -3.579940e+00 -6.716511e+00 527 -8.636345e+00 -3.579940e+00 528 -1.555820e+01 -8.636345e+00 529 -2.844068e+00 -1.555820e+01 530 -1.568663e+01 -2.844068e+00 531 7.604625e+02 -1.568663e+01 532 8.764713e+02 7.604625e+02 533 -6.374334e+00 8.764713e+02 534 -4.040232e+00 -6.374334e+00 535 -1.266852e+01 -4.040232e+00 536 -6.669758e+00 -1.266852e+01 537 -1.159025e+01 -6.669758e+00 538 -5.567872e+00 -1.159025e+01 539 -1.433900e+01 -5.567872e+00 540 -4.584389e+00 -1.433900e+01 541 -1.528106e+01 -4.584389e+00 542 -2.295174e+00 -1.528106e+01 543 -7.161032e+00 -2.295174e+00 544 -7.668840e-01 -7.161032e+00 545 -6.203831e+00 -7.668840e-01 546 -1.098963e+01 -6.203831e+00 547 -1.010455e+01 -1.098963e+01 548 -8.494776e+00 -1.010455e+01 549 -1.197824e+01 -8.494776e+00 550 -1.257357e+01 -1.197824e+01 551 -1.081448e+01 -1.257357e+01 552 -3.203562e+00 -1.081448e+01 553 -6.306076e+00 -3.203562e+00 554 -9.063120e+00 -6.306076e+00 555 -6.153311e+00 -9.063120e+00 556 -5.932945e+00 -6.153311e+00 557 -4.971402e+00 -5.932945e+00 558 -1.140417e+01 -4.971402e+00 559 -1.095120e+01 -1.140417e+01 560 -7.335587e+00 -1.095120e+01 561 -1.162483e+01 -7.335587e+00 562 -3.547614e+00 -1.162483e+01 563 -1.027371e+01 -3.547614e+00 564 -6.259756e+00 -1.027371e+01 565 -5.723048e+00 -6.259756e+00 566 -9.172901e+00 -5.723048e+00 567 -3.756818e+00 -9.172901e+00 568 -3.759076e+00 -3.756818e+00 569 -1.042397e+00 -3.759076e+00 570 -8.032350e+00 -1.042397e+00 571 -8.323509e+00 -8.032350e+00 572 -8.552383e+00 -8.323509e+00 573 -8.905003e+00 -8.552383e+00 574 -3.192250e+00 -8.905003e+00 575 1.258708e+00 -3.192250e+00 576 -9.910202e+00 1.258708e+00 577 7.485737e+00 -9.910202e+00 578 -1.350696e+01 7.485737e+00 579 3.018549e+00 -1.350696e+01 580 -5.267373e+00 3.018549e+00 581 -2.533336e+00 -5.267373e+00 582 -3.300970e+00 -2.533336e+00 583 -6.979721e+00 -3.300970e+00 584 -1.360075e+01 -6.979721e+00 585 -7.169105e+00 -1.360075e+01 586 -1.447020e+01 -7.169105e+00 587 -5.456367e+00 -1.447020e+01 588 -2.574610e+00 -5.456367e+00 > 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/freestat/rcomp/tmp/7boso1293645852.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/freestat/rcomp/tmp/8boso1293645852.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/freestat/rcomp/tmp/9boso1293645852.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/freestat/rcomp/tmp/104fsr1293645852.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/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/11pgqw1293645852.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/freestat/rcomp/tmp/12ty721293645852.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/freestat/rcomp/tmp/137qnb1293645852.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/freestat/rcomp/tmp/14ar3h1293645852.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/freestat/rcomp/tmp/15wrkn1293645852.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/freestat/rcomp/tmp/16hs0t1293645852.tab") + } > > try(system("convert tmp/1fedx1293645852.ps tmp/1fedx1293645852.png",intern=TRUE)) character(0) > try(system("convert tmp/28oc01293645852.ps tmp/28oc01293645852.png",intern=TRUE)) character(0) > try(system("convert tmp/38oc01293645852.ps tmp/38oc01293645852.png",intern=TRUE)) character(0) > try(system("convert tmp/48oc01293645852.ps tmp/48oc01293645852.png",intern=TRUE)) character(0) > try(system("convert tmp/58oc01293645852.ps tmp/58oc01293645852.png",intern=TRUE)) character(0) > try(system("convert tmp/60xbl1293645852.ps tmp/60xbl1293645852.png",intern=TRUE)) character(0) > try(system("convert tmp/7boso1293645852.ps tmp/7boso1293645852.png",intern=TRUE)) character(0) > try(system("convert tmp/8boso1293645852.ps tmp/8boso1293645852.png",intern=TRUE)) character(0) > try(system("convert tmp/9boso1293645852.ps tmp/9boso1293645852.png",intern=TRUE)) character(0) > try(system("convert tmp/104fsr1293645852.ps tmp/104fsr1293645852.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 20.424 3.493 20.928