Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
2441 1 0 0 0 3406 1 0 0 0 4029 0 1 0 0 1924 0 1 0 0 2319 0 1 0 0 2156 0 1 0 0 2117 0 1 0 0 2189 0 1 0 0 2625 0 0 1 0 1959 0 0 1 0 3096 0 0 1 0 1997 0 0 1 0 1813 0 0 0 1 2648 0 0 0 1 2648 0 0 0 1 5782 0 0 0 1 2218 0 0 0 1 25289 1 0 0 0 25389 1 0 0 0 25196 0 1 0 0 24153 0 1 0 0 26079 0 1 0 0 26540 0 1 0 0 25099 0 0 1 0 27402 0 0 1 0 26206 0 0 1 0 26526 0 0 1 0 27154 0 0 1 0 24464 0 0 1 0 24967 0 0 0 1 24713 0 0 0 1 26733 0 0 0 1 24578 0 0 0 1 26092 0 0 0 1 24486 0 1 0 0 24630 0 1 0 0 26233 0 1 0 0 26492 0 1 0 0 24973 0 1 0 0 27437 0 1 0 0 23829 0 1 0 0 27158 0 0 1 0 25670 0 0 1 0 23530 0 0 1 0 24474 0 0 1 0 26668 0 0 1 0 26060 0 0 0 1 24856 0 0 0 1 24067 0 0 0 1 25545 0 0 0 1 24213 0 0 0 1 23703 1 0 0 0 23566 1 0 0 0 22876 0 1 0 0 22744 0 1 0 0 27615 0 1 0 0 24421 0 1 0 0 24728 0 1 0 0 25732 0 1 0 0 24204 0 1 0 0 23869 0 0 1 0 24120 0 0 1 0 22474 0 0 1 0 26406 0 0 1 0 22440 0 0 0 1 20387 0 0 0 0 21609 0 0 0 0 24905 0 0 0 0 21584 0 0 0 0 20920 0 0 0 0 5042 1 0 0 0 4353 1 0 0 0 7996 0 1 0 0 3998 0 1 0 0 4697 0 1 0 0 7837 0 1 0 0 3512 0 0 1 0 3503 0 0 1 0 3572 0 0 1 0 3918 0 0 1 0 4767 0 0 1 0 5833 0 0 0 1 4154 0 0 0 1 3894 0 0 0 1 4133 0 0 0 1 4273 0 0 0 1 5574 0 0 0 1 5029 1 0 0 0 5279 0 1 0 0 4876 0 1 0 0 3850 0 1 0 0 4109 0 1 0 0 4137 0 1 0 0 3725 0 1 0 0 5675 0 0 1 0 3405 0 0 1 0 3568 0 0 1 0 3408 0 0 1 0 7203 0 0 1 0 5392 0 0 0 1 4053 0 0 0 1 7863 0 0 0 1 3716 0 0 0 1 4027 0 0 0 1 3608 1 0 0 0 3333 1 0 0 0 3014 0 1 0 0 5014 0 1 0 0 4328 0 1 0 0 2956 0 1 0 0 6535 0 1 0 0 3153 0 1 0 0 3081 0 1 0 0 2996 0 0 1 0 3150 1 0 0 0 3673 1 0 0 0 2870 0 1 0 0 3230 0 1 0 0 3821 0 1 0 0 3178 0 1 0 0 2988 0 1 0 0 2347 0 0 1 0 2891 0 0 1 0 4775 0 0 1 0 4758 0 0 1 0 2962 0 0 1 0 2687 0 0 0 1 2825 0 0 0 1 4201 0 0 0 1 2545 0 0 0 1 2626 0 0 0 1 3556 1 0 0 0 6069 1 0 0 0 2795 0 1 0 0 2763 0 1 0 0 3024 0 1 0 0 2622 0 1 0 0 3800 0 0 1 0 5217 0 0 1 0 3163 0 0 1 0 3765 0 0 1 0 2991 0 0 1 0 4856 0 0 0 1 5752 0 0 0 1 3351 0 0 0 1 3392 0 0 0 1 3145 0 0 0 1 3820 0 0 0 1 4790 1 0 0 0 2729 1 0 0 0 3025 0 1 0 0 2428 0 1 0 0 2981 0 1 0 0 3051 0 1 0 0 6330 0 1 0 0 3006 0 0 1 0 3301 0 0 1 0 5265 0 0 1 0 3975 0 0 1 0 2643 0 0 1 0 3130 0 0 0 1 3832 0 0 0 1 3819 0 0 0 1 3037 0 0 0 1 4272 0 0 0 1 10589 1 0 0 0 8945 1 0 0 0 7764 0 1 0 0 8704 0 1 0 0 7546 0 1 0 0 7694 0 1 0 0 10499 0 1 0 0 7614 0 0 1 0 8248 0 0 1 0 8158 0 0 1 0 8174 0 0 1 0 8097 0 0 1 0 9154 0 0 0 1 10287 0 0 0 1 7972 0 0 0 1 7518 0 0 0 1 9492 0 0 0 1 8317 1 0 0 0 8158 0 1 0 0 9174 0 1 0 0 8262 0 1 0 0 10533 0 1 0 0 10434 0 1 0 0 8047 0 1 0 0 7831 0 0 1 0 8062 0 0 1 0 8834 0 0 1 0 8957 0 0 1 0 8753 0 0 1 0 7663 0 0 1 0 8290 0 0 0 1 8435 0 0 0 1 10802 0 0 0 1 9391 0 0 0 1 10280 1 0 0 0 8461 1 0 0 0 9152 1 0 0 0 8380 0 1 0 0 8171 0 1 0 0 8386 0 1 0 0 8212 0 1 0 0 9103 0 1 0 0 8461 0 1 0 0 8443 0 0 1 0 9253 0 0 1 0 8220 0 0 1 0 10435 0 0 1 0 8627 0 0 1 0 8196 0 0 0 1 9431 0 0 0 1 7917 0 0 0 0 8186 0 0 0 0 4350 0 0 0 0 9341 1 0 0 0 9545 0 0 0 1 10624 1 0 0 0 10665 1 0 0 0 11698 0 1 0 0 9516 0 1 0 0 8815 0 1 0 0 8389 0 1 0 0 10475 0 1 0 0 10170 0 0 1 0 9192 0 0 1 0 9198 0 0 1 0 8764 0 0 1 0 9996 0 0 1 0 9219 0 0 0 1 10801 0 0 0 1 8631 0 0 0 1 11110 0 0 0 1 8101 0 0 0 1 9696 1 0 0 0 10542 0 1 0 0 10069 0 1 0 0 11789 0 1 0 0 9416 0 1 0 0 9543 0 1 0 0 8919 0 1 0 0 8958 0 0 1 0 8933 0 0 1 0 11251 0 0 1 0 9589 0 0 1 0 8870 0 0 0 1 9108 0 0 0 1 9544 0 0 0 1 9611 0 0 0 1 11798 0 0 0 1 11269 1 0 0 0 10411 1 0 0 0 9690 0 1 0 0 9625 0 1 0 0 9522 0 1 0 0 10330 0 1 0 0 10803 0 1 0 0 9946 0 1 0 0 9782 0 0 1 0 11660 0 0 1 0 9960 0 0 1 0 10286 0 0 1 0 10790 0 0 1 0 10188 0 0 1 0 9465 0 0 0 0 7791 0 0 0 0 7793 0 0 0 0 8175 0 0 0 0 10328 0 0 0 0 4510 1 0 0 0 3589 0 1 0 0 4039 0 1 0 0 7656 0 1 0 0 4662 0 1 0 0 5001 0 0 1 0 7089 0 0 1 0 4103 0 0 1 0 4314 0 0 1 0 7187 0 0 1 0 5954 0 0 1 0 3597 0 0 1 0 3647 0 0 0 1 8287 0 0 0 1 4192 0 0 0 1 4046 0 0 0 1 5195 0 0 0 1 7626 1 0 0 0 5232 1 0 0 0 5251 1 0 0 0 5043 0 1 0 0 5842 0 1 0 0 4879 0 1 0 0 5429 0 1 0 0 4772 0 1 0 0 6159 0 0 1 0 3761 0 0 1 0 8832 0 0 1 0 4337 0 0 1 0 3979 0 0 1 0 4886 0 0 1 0 6057 0 0 0 1 4922 0 0 0 0 4650 0 0 0 0 4938 0 0 0 0 6610 0 0 0 0 6041 0 0 0 0 11379 1 0 0 0 10702 1 0 0 0 7455 0 1 0 0 8425 0 1 0 0 7679 0 1 0 0 8312 0 1 0 0 7238 0 0 1 0 9412 0 0 1 0 7698 0 0 1 0 7776 0 0 1 0 7870 0 0 1 0 8122 0 0 0 1 9138 0 0 0 1 10187 0 0 0 1 8315 0 0 0 1 8424 0 0 0 1 7731 0 0 0 1 8079 1 0 0 0 7926 0 1 0 0 9975 0 1 0 0 8397 0 1 0 0 8572 0 1 0 0 8157 0 1 0 0 7856 0 1 0 0 9835 0 1 0 0 9524 0 0 1 0 7750 0 0 1 0 8221 0 0 1 0 8998 0 0 1 0 9875 0 0 0 1 8015 0 0 0 1 7749 0 0 0 1 8174 0 0 0 1 8504 0 0 0 1 11129 1 0 0 0 12615 1 0 0 0 12219 1 0 0 0 10828 0 1 0 0 11463 0 1 0 0 12524 0 1 0 0 10638 0 1 0 0 11085 0 0 1 0 10831 0 0 1 0 12022 0 0 1 0 10544 0 0 1 0 11569 0 0 1 0 10889 0 0 1 0 11064 0 0 0 1 11221 0 0 0 1 10339 0 0 0 0 10652 0 0 0 0 11155 0 0 0 0 3597 1 0 0 0 2768 1 0 0 0 2812 0 1 0 0 3781 0 1 0 0 7789 0 1 0 0 2886 0 1 0 0 2283 0 0 1 0 2389 0 0 1 0 3784 0 0 1 0 2990 0 0 1 0 3615 0 0 1 0 2767 0 0 0 1 2673 0 0 0 1 3068 0 0 0 1 2894 0 0 0 1 2621 0 0 0 1 6440 0 0 0 1 3082 1 0 0 0 5532 0 1 0 0 2421 0 1 0 0 3653 0 1 0 0 2656 0 1 0 0 3059 0 1 0 0 3341 0 1 0 0 2387 0 0 1 0 2469 0 0 1 0 2758 0 0 1 0 2254 0 0 1 0 2305 0 0 1 0 7075 0 0 0 1 2260 0 0 0 1 2988 0 0 0 1 2091 0 0 0 1 2169 0 0 0 1 15711 1 0 0 0 14409 0 0 1 0 17306 1 0 0 0 17157 1 0 0 0 17611 0 1 0 0 20394 0 1 0 0 18757 0 1 0 0 20250 0 1 0 0 17622 0 1 0 0 17270 0 0 1 0 18330 0 0 1 0 17580 0 0 1 0 18128 0 0 1 0 17261 0 0 0 1 17287 0 0 0 1 17433 0 0 0 1 17518 0 0 0 1 16890 0 0 0 1 18728 0 0 0 1 16953 1 0 0 0 17970 0 1 0 0 16920 0 1 0 0 19400 0 1 0 0 15769 0 1 0 0 17431 0 1 0 0 16058 0 0 1 0 15312 0 0 1 0 16214 0 0 1 0 15962 0 0 1 0 15852 0 0 0 1 15634 0 0 0 1 17699 0 0 0 1 16100 0 0 0 1 16252 0 0 0 1 17874 0 0 0 1 14058 1 0 0 0 14466 1 0 0 0 14531 1 0 0 0 14102 0 1 0 0 14014 0 1 0 0 16871 0 1 0 0 14903 0 1 0 0 16411 0 1 0 0 14687 0 0 1 0 14363 0 0 1 0 16062 0 0 1 0 15361 0 0 1 0 16134 0 0 1 0 14256 0 0 0 1 15863 0 0 0 0 14196 0 0 0 0 14120 0 0 0 0 14825 0 0 0 0 16946 0 0 0 0 23867 1 0 0 0 24107 1 0 0 0 24041 0 1 0 0 24415 0 1 0 0 24496 0 1 0 0 24022 0 1 0 0 24367 0 1 0 0 23869 0 0 1 0 24495 0 0 1 0 23818 0 0 1 0 24081 0 0 1 0 24132 0 0 0 1 23651 0 0 0 1 23622 0 0 0 1 23726 0 0 0 1 23942 0 0 0 1 24573 0 0 0 1 23085 1 0 0 0 22612 1 0 0 0 22960 0 1 0 0 22921 0 1 0 0 23510 0 1 0 0 22729 0 1 0 0 23047 0 1 0 0 22850 0 0 1 0 23426 0 0 1 0 22812 0 0 1 0 22446 0 0 1 0 23567 0 0 1 0 23185 0 0 0 1 22777 0 0 0 1 23508 0 0 0 1 23193 0 0 0 1 23006 0 0 0 1 22332 1 0 0 0 22347 1 0 0 0 23061 1 0 0 0 22887 0 1 0 0 22890 0 1 0 0 22701 0 1 0 0 22467 0 1 0 0 22357 0 1 0 0 22443 0 1 0 0 22824 0 0 1 0 22906 0 0 1 0 23059 0 0 1 0 23055 0 0 1 0 22564 0 0 0 1 18570 0 0 0 0 20329 0 0 0 0 19279 0 0 0 0 19541 0 0 0 0 19517 0 0 0 0 6519 1 0 0 0 7169 1 0 0 0 8107 1 0 0 0 10668 0 1 0 0 6650 0 1 0 0 5726 0 1 0 0 5224 0 0 1 0 6297 0 0 1 0 5011 0 0 1 0 5075 0 0 1 0 5434 0 0 1 0 11758 0 0 0 1 4531 0 0 0 1 5373 0 0 0 1 6343 0 0 0 1 20051 0 0 0 1 5482 0 0 0 1 5066 1 0 0 0 5040 0 1 0 0 5100 0 1 0 0 4679 0 1 0 0 10940 0 1 0 0 6035 0 1 0 0 5364 0 0 1 0 4424 0 0 1 0 4486 0 0 1 0 4962 0 0 1 0 10445 0 0 1 0 5973 0 0 1 0 5415 0 0 0 1 17792 0 0 0 1 5581 0 0 0 1 4997 0 0 0 1 6893 0 0 0 1 10181 1 0 0 0 7007 1 0 0 0 6621 0 1 0 0 7309 0 1 0 0 6114 0 1 0 0 5521 0 1 0 0 5263 0 1 0 0 5400 0 0 1 0 6141 0 0 1 0 5736 0 0 1 0 16104 0 0 1 0 10810 0 0 0 1 5057 0 0 0 1 5732 0 0 0 1 4000 0 0 0 0 4000 0 0 0 0 4200 0 0 0 0 3551 1 0 0 0 4025 1 0 0 0 5591 0 1 0 0 3868 0 1 0 0 3566 0 1 0 0 4525 0 1 0 0 3752 0 0 1 0 3182 0 0 1 0 6152 0 0 1 0 3548 0 0 1 0 6876 0 0 1 0 3199 0 0 0 1 3386 0 0 0 1 3411 0 0 0 1 6892 0 0 0 1 4920 0 0 0 1 3193 0 0 0 1 3054 1 0 0 0 3262 0 1 0 0 3509 0 1 0 0 3471 0 1 0 0 3101 0 1 0 0 5956 0 1 0 0 6232 0 0 1 0 5456 0 0 1 0 3186 0 0 1 0 3751 0 0 1 0 2973 0 0 1 0 5548 0 0 0 1 3219 0 0 0 1 6595 0 0 0 1 4886 0 0 0 1 3082 0 0 0 1 3516 0 0 0 1 3807 1 0 0 0 3607 1 0 0 0 3163 1 0 0 0 4981 0 1 0 0 3276 0 1 0 0 3278 0 1 0 0 3850 0 1 0 0 3439 0 1 0 0 5545 0 0 1 0 4749 0 0 1 0 3656 0 0 1 0 3520 0 0 1 0 4392 0 0 1 0 3057 0 0 1 0 6542 0 0 0 1 2785 0 0 0 0 3057 0 0 0 0 4379 0 0 0 0 23934 1 0 0 0 23625 1 0 0 0 26185 0 1 0 0 23777 0 1 0 0 24586 0 1 0 0 25439 0 1 0 0 24037 0 1 0 0 25403 0 0 1 0 25133 0 0 1 0 24023 0 0 1 0 23901 0 0 1 0 24892 0 0 1 0 24560 0 0 0 1 24226 0 0 0 1 24885 0 0 0 1 25466 0 0 0 1 24903 0 0 0 1 23761 1 0 0 0 23868 0 1 0 0 26118 0 1 0 0 25120 0 1 0 0 26119 0 1 0 0 25440 0 1 0 0 24206 0 0 1 0 25312 0 0 1 0 24499 0 0 1 0 24330 0 0 1 0 24217 0 0 1 0 25047 0 0 0 1 25817 0 0 0 1 25466 0 0 0 1 25410 0 0 0 1 26246 0 0 0 1 25718 0 0 0 1 26543 1 0 0 0 26723 1 0 0 0 26700 0 1 0 0 24743 0 1 0 0 25520 0 1 0 0 25298 0 1 0 0 26382 0 1 0 0 25719 0 0 1 0 24547 0 0 1 0 25640 0 0 1 0 26103 0 0 1 0 24911 0 0 1 0 25199 0 0 1 0 25308 0 0 0 1 17000 0 0 0 0 12000 0 0 0 0 20000 0 0 0 0 3915 1 0 0 0 3229 1 0 0 0 6671 0 1 0 0 5937 0 1 0 0 3639 0 1 0 0 4274 0 1 0 0 3781 0 0 1 0 5612 0 0 1 0 4498 0 0 1 0 3520 0 0 1 0 6323 0 0 1 0 3622 0 0 1 0 4085 0 0 0 1 3978 0 0 0 1 3788 0 0 0 1 3973 0 0 0 1 3268 0 0 0 1 6852 1 0 0 0 6237 1 0 0 0 9194 0 1 0 0 9177 0 1 0 0 6170 0 1 0 0 6295 0 1 0 0 5878 0 1 0 0 7172 0 1 0 0 5741 0 0 1 0 7093 0 0 1 0 5774 0 0 1 0 5690 0 0 1 0 7717 0 0 1 0 6511 0 0 0 1 6989 0 0 0 0 9006 0 0 0 0 6052 0 0 0 0 5094 0 0 0 0 6198 0 0 0 0 8845 0 0 0 0 6219 1 0 0 0 5984 1 0 0 0 7303 0 1 0 0 6887 0 1 0 0 8083 0 1 0 0 18978 0 1 0 0 25222 0 1 0 0 6093 0 0 1 0 7206 0 0 1 0 8070 0 0 1 0 16129 0 0 1 0 7646 0 0 1 0 5415 0 0 0 1 10480 0 0 0 1 5998 0 0 0 1 6289 0 0 0 1 6146 0 0 0 1 12308 1 0 0 0 7128 1 0 0 0 7653 0 1 0 0 10130 0 1 0 0 8741 0 1 0 0 7719 0 1 0 0 8167 0 1 0 0 7786 0 0 1 0 8091 0 0 1 0 8089 0 0 1 0 7219 0 0 1 0 7373 0 0 1 0 20659 0 0 0 1 7158 0 0 0 1 7503 0 0 0 1 7654 0 0 0 1 7747 0 0 0 1 8631 1 0 0 0 6867 1 0 0 0 17989 1 0 0 0 6798 0 1 0 0 7485 0 1 0 0 7085 0 1 0 0 7235 0 1 0 0 8534 0 0 1 0 7785 0 0 1 0 7017 0 0 1 0 5951 0 0 1 0 6709 0 0 1 0 6999 0 0 0 1 6175 0 0 0 1 5927 0 0 0 0 6703 0 0 0 0 6118 0 0 0 0 4565 1 0 0 0 6178 1 0 0 0 5389 0 1 0 0 5079 0 1 0 0 7382 0 1 0 0 4305 0 1 0 0 4216 0 1 0 0 7849 0 0 1 0 4248 0 0 1 0 4238 0 0 1 0 4746 0 0 1 0 4227 0 0 1 0 4946 0 0 0 1 4234 0 0 0 1 4379 0 0 0 1 5464 0 0 0 1 4240 0 0 0 1 5465 1 0 0 0 5634 1 0 0 0 3984 0 1 0 0 6957 0 1 0 0 4492 0 1 0 0 3863 0 1 0 0 3845 0 1 0 0 3768 0 0 1 0 6071 0 0 1 0 3794 0 0 1 0 4078 0 0 1 0 3927 0 0 1 0 3931 0 0 0 1 5368 0 0 0 1 5142 0 0 0 1 5165 0 0 0 1 4432 0 0 0 1 5082 1 0 0 0 6087 1 0 0 0 4434 1 0 0 0 4360 0 1 0 0 5634 0 1 0 0 7836 0 1 0 0 5394 0 1 0 0 4327 0 0 1 0 4142 0 0 1 0 5251 0 0 1 0 4951 0 0 1 0 4565 0 0 1 0 4463 0 0 1 0 5922 0 0 1 0 4581 0 0 0 1 7566 1 0 0 0 5500 1 0 0 0 5745 0 1 0 0 6924 0 1 0 0 5354 0 1 0 0 5563 0 1 0 0 5369 0 1 0 0 5658 0 0 1 0 5215 0 0 1 0 5824 0 0 1 0 6667 0 0 1 0 7795 0 0 1 0 5490 0 0 0 1 5232 0 0 0 1 7739 0 0 0 1 5404 0 0 0 1 6045 0 0 0 1 6012 1 0 0 0 6287 0 1 0 0 5185 0 1 0 0 8080 0 1 0 0 7229 0 1 0 0 5602 0 1 0 0 5329 0 1 0 0 5401 0 0 1 0 8283 0 0 1 0 6359 0 0 1 0 5457 0 0 1 0 5654 0 0 1 0 6391 0 0 0 1 5765 0 0 0 1 6707 0 0 0 1 8214 0 0 0 1 5621 0 0 0 1 6387 1 0 0 0 8299 1 0 0 0 6526 0 1 0 0 5514 0 1 0 0 6659 0 1 0 0 6023 0 1 0 0 5701 0 1 0 0 6628 0 1 0 0 5845 0 1 0 0 5778 0 0 1 0 5668 0 0 1 0 5982 0 0 1 0 8294 0 0 1 0 5970 0 0 1 0 7440 0 0 0 1 5385 0 0 0 0 6226 0 0 0 0 6905 0 0 0 0 7566 0 0 0 0 6033 0 0 0 0 3338 1 0 0 0 2778 1 0 0 0 2876 0 1 0 0 3059 0 1 0 0 2827 0 1 0 0 3819 0 1 0 0 3319 0 1 0 0 5529 0 0 1 0 2791 0 0 1 0 6521 0 0 1 0 2959 0 0 1 0 4378 0 0 1 0 6042 0 0 0 1 3715 0 0 0 1 6219 0 0 0 1 2890 0 0 0 1 3134 0 0 0 1 3544 1 0 0 0 3915 1 0 0 0 3139 0 1 0 0 2989 0 1 0 0 2856 0 1 0 0 5619 0 1 0 0 3955 0 1 0 0 3027 0 1 0 0 3760 0 0 1 0 6323 0 0 1 0 3362 0 0 1 0 6263 0 0 1 0 5720 0 0 0 1 3035 0 0 0 1 6509 0 0 0 1 3123 0 0 0 1 3332 0 0 0 1 3298 1 0 0 0 4579 1 0 0 0 2963 1 0 0 0 5861 0 1 0 0 4549 0 1 0 0 6211 0 1 0 0 2942 0 1 0 0 3181 0 1 0 0 5019 0 1 0 0 6590 0 0 1 0 4528 0 0 1 0 3744 0 0 1 0 3096 0 0 1 0 2893 0 0 0 1 3946 0 0 0 0 2838 0 0 0 0 2804 0 0 0 0
Names of X columns:
Toeschouwers juli-aug sep-nov dec-feb ma-eind
Sample Range:
(leave blank to include all observations)
From:
To:
Column Number of Endogenous Series
(?)
Fixed Seasonal Effects
Do not include Seasonal Dummies
Do not include Seasonal Dummies
Include Seasonal Dummies
Type of Equation
No Linear Trend
No Linear Trend
Linear Trend
First Differences
Seasonal Differences (s)
First and Seasonal Differences (s)
Degree of Predetermination (lagged endogenous variables)
Degree of Seasonal Predetermination
Seasonality
12
1
2
3
4
5
6
7
8
9
10
11
12
Chart options
R Code
library(lattice) library(lmtest) n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test mywarning <- '' par1 <- as.numeric(par1) if(is.na(par1)) { par1 <- 1 mywarning = 'Warning: you did not specify the column number of the endogenous series! The first column was selected by default.' } if (par4=='') par4 <- 0 par4 <- as.numeric(par4) if (par5=='') par5 <- 0 par5 <- as.numeric(par5) x <- na.omit(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'){ (n <- n -1) x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B)',colnames(x),sep=''))) for (i in 1:n) { for (j in 1:k) { x2[i,j] <- x[i+1,j] - x[i,j] } } x <- x2 } if (par3 == 'Seasonal Differences (s=12)'){ (n <- n - 12) x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B12)',colnames(x),sep=''))) for (i in 1:n) { for (j in 1:k) { x2[i,j] <- x[i+12,j] - x[i,j] } } x <- x2 } if (par3 == 'First and Seasonal Differences (s=12)'){ (n <- n -1) x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B)',colnames(x),sep=''))) for (i in 1:n) { for (j in 1:k) { x2[i,j] <- x[i+1,j] - x[i,j] } } x <- x2 (n <- n - 12) x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B12)',colnames(x),sep=''))) for (i in 1:n) { for (j in 1:k) { x2[i,j] <- x[i+12,j] - x[i,j] } } x <- x2 } if(par4 > 0) { x2 <- array(0, dim=c(n-par4,par4), dimnames=list(1:(n-par4), paste(colnames(x)[par1],'(t-',1:par4,')',sep=''))) for (i in 1:(n-par4)) { for (j in 1:par4) { x2[i,j] <- x[i+par4-j,par1] } } x <- cbind(x[(par4+1):n,], x2) n <- n - par4 } if(par5 > 0) { x2 <- array(0, dim=c(n-par5*12,par5), dimnames=list(1:(n-par5*12), paste(colnames(x)[par1],'(t-',1:par5,'s)',sep=''))) for (i in 1:(n-par5*12)) { for (j in 1:par5) { x2[i,j] <- x[i+par5*12-j*12,par1] } } x <- cbind(x[(par5*12+1):n,], x2) n <- n - par5*12 } 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[n,])) if (par3 == 'Linear Trend'){ x <- cbind(x, c(1:n)) colnames(x)[k+1] <- 't' } x (k <- length(x[n,])) head(x) df <- as.data.frame(x) (mylm <- lm(df)) (mysum <- summary(mylm)) 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 } bitmap(file='test0.png') 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() bitmap(file='test1.png') plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') grid() dev.off() bitmap(file='test2.png') hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') grid() dev.off() bitmap(file='test3.png') densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') dev.off() bitmap(file='test4.png') qqnorm(mysum$resid, main='Residual Normal Q-Q Plot') qqline(mysum$resid) grid() dev.off() (myerror <- as.ts(mysum$resid)) bitmap(file='test5.png') dum <- cbind(lag(myerror,k=1),myerror) dum dum1 <- dum[2:length(myerror),] dum1 z <- as.data.frame(dum1) z 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() bitmap(file='test6.png') acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') grid() dev.off() bitmap(file='test7.png') pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') grid() dev.off() bitmap(file='test8.png') opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) plot(mylm, las = 1, sub='Residual Diagnostics') par(opar) dev.off() if (n > n25) { bitmap(file='test9.png') plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') grid() dev.off() } load(file='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, signif(mysum$coefficients[i,1],6), 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.row.start(a) a<-table.element(a, mywarning) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable1.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<br />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,formatC(signif(mysum$coefficients[i,1],5),format='g',flag='+')) a<-table.element(a,formatC(signif(mysum$coefficients[i,2],5),format='g',flag=' ')) a<-table.element(a,formatC(signif(mysum$coefficients[i,3],4),format='e',flag='+')) a<-table.element(a,formatC(signif(mysum$coefficients[i,4],4),format='g',flag=' ')) a<-table.element(a,formatC(signif(mysum$coefficients[i,4]/2,4),format='g',flag=' ')) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable2.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,formatC(signif(sqrt(mysum$r.squared),6),format='g',flag=' ')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'R-squared',1,TRUE) a<-table.element(a,formatC(signif(mysum$r.squared,6),format='g',flag=' ')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Adjusted R-squared',1,TRUE) a<-table.element(a,formatC(signif(mysum$adj.r.squared,6),format='g',flag=' ')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'F-TEST (value)',1,TRUE) a<-table.element(a,formatC(signif(mysum$fstatistic[1],6),format='g',flag=' ')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) a<-table.element(a, signif(mysum$fstatistic[2],6)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) a<-table.element(a, signif(mysum$fstatistic[3],6)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'p-value',1,TRUE) a<-table.element(a,formatC(signif(1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]),6),format='g',flag=' ')) 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,formatC(signif(mysum$sigma,6),format='g',flag=' ')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Sum Squared Residuals',1,TRUE) a<-table.element(a,formatC(signif(sum(myerror*myerror),6),format='g',flag=' ')) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable3.tab') if(n < 200) { 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<br />Forecast', 1, TRUE) a<-table.element(a, 'Residuals<br />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,formatC(signif(x[i],6),format='g',flag=' ')) a<-table.element(a,formatC(signif(x[i]-mysum$resid[i],6),format='g',flag=' ')) a<-table.element(a,formatC(signif(mysum$resid[i],6),format='g',flag=' ')) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable4.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,formatC(signif(gqarr[mypoint-kp3+1,1],6),format='g',flag=' ')) a<-table.element(a,formatC(signif(gqarr[mypoint-kp3+1,2],6),format='g',flag=' ')) a<-table.element(a,formatC(signif(gqarr[mypoint-kp3+1,3],6),format='g',flag=' ')) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable5.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,signif(numsignificant1,6)) a<-table.element(a,formatC(signif(numsignificant1/numgqtests,6),format='g',flag=' ')) 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,signif(numsignificant5,6)) a<-table.element(a,signif(numsignificant5/numgqtests,6)) 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,signif(numsignificant10,6)) a<-table.element(a,signif(numsignificant10/numgqtests,6)) 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='mytable6.tab') } }
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation