R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(210907 + ,1 + ,1 + ,2 + ,2 + ,3 + ,4 + ,3 + ,120982 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,176508 + ,2 + ,2 + ,2 + ,4 + ,2 + ,2 + ,1 + ,179321 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,123185 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,52746 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,385534 + ,1 + ,1 + ,2 + ,2 + ,2 + ,3 + ,2 + ,33170 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,101645 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,149061 + ,1 + ,1 + ,2 + ,2 + ,3 + ,4 + ,1 + ,165446 + ,2 + ,2 + ,2 + ,3 + ,3 + ,3 + ,2 + ,237213 + ,1 + ,1 + ,1 + ,2 + ,3 + ,4 + ,1 + ,173326 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,133131 + ,2 + ,1 + ,1 + ,1 + ,1 + ,4 + ,2 + ,258873 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,180083 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,324799 + ,1 + ,1 + ,1 + ,1 + ,1 + ,3 + ,1 + ,230964 + ,4 + ,4 + ,4 + ,4 + ,3 + ,2 + ,4 + ,236785 + ,2 + ,1 + ,1 + ,1 + ,2 + ,4 + ,2 + ,135473 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,202925 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,215147 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,1 + ,344297 + ,1 + ,1 + ,1 + ,4 + ,1 + ,4 + ,1 + ,153935 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,132943 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,174724 + ,1 + ,1 + ,1 + ,1 + ,1 + ,3 + ,1 + ,174415 + ,1 + ,1 + ,1 + ,3 + ,3 + ,3 + ,2 + ,225548 + ,4 + ,4 + ,4 + ,4 + ,3 + ,3 + ,4 + ,223632 + ,1 + ,1 + ,1 + ,2 + ,2 + ,3 + ,2 + ,124817 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,2 + ,221698 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,210767 + ,4 + ,3 + ,3 + ,2 + ,2 + ,2 + ,3 + ,170266 + ,1 + ,1 + ,1 + ,2 + ,1 + ,4 + ,2 + ,260561 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,84853 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,294424 + ,1 + ,1 + ,1 + ,1 + ,2 + ,1 + ,2 + ,101011 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,215641 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,325107 + ,2 + ,1 + ,2 + ,2 + ,2 + ,3 + ,3 + ,7176 + ,2 + ,1 + ,2 + ,2 + ,3 + ,3 + ,2 + ,167542 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,106408 + ,3 + ,4 + ,4 + ,2 + ,4 + ,2 + ,4 + ,96560 + ,2 + ,3 + ,3 + ,3 + ,3 + ,2 + ,4 + ,265769 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,269651 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,149112 + ,2 + ,1 + ,2 + ,2 + ,2 + ,3 + ,3 + ,175824 + ,1 + ,1 + ,1 + ,2 + ,1 + ,1 + ,1 + ,152871 + ,1 + ,1 + ,1 + ,3 + ,1 + ,4 + ,1 + ,111665 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,1 + ,116408 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,362301 + ,2 + ,3 + ,3 + ,3 + ,2 + ,3 + ,2 + ,78800 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,183167 + ,2 + ,3 + ,3 + ,3 + ,2 + ,3 + ,3 + ,277965 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,150629 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,168809 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,24188 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,329267 + ,2 + ,3 + ,2 + ,2 + ,2 + ,1 + ,1 + ,65029 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,101097 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,218946 + ,3 + ,3 + ,3 + ,4 + ,4 + ,2 + ,3 + ,244052 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,1 + ,341570 + ,2 + ,2 + ,2 + ,2 + ,2 + ,3 + ,2 + ,103597 + ,2 + ,1 + ,2 + ,2 + ,3 + ,4 + ,2 + ,233328 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,256462 + ,2 + ,2 + ,2 + ,3 + ,3 + ,2 + ,2 + ,206161 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,311473 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,235800 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,177939 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,207176 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,196553 + ,2 + ,1 + ,1 + ,2 + ,2 + ,4 + ,1 + ,174184 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,143246 + ,1 + ,1 + ,2 + ,1 + ,1 + ,4 + ,1 + ,187559 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,187681 + ,2 + ,2 + ,2 + ,3 + ,3 + ,2 + ,2 + ,119016 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,182192 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,73566 + ,1 + ,1 + ,1 + ,3 + ,3 + ,4 + ,1 + ,194979 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,167488 + ,1 + ,1 + ,2 + ,1 + ,2 + ,3 + ,1 + ,143756 + ,3 + ,1 + ,2 + ,3 + ,4 + ,3 + ,4 + ,275541 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,243199 + ,1 + ,2 + ,2 + ,3 + ,2 + ,3 + ,3 + ,182999 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,135649 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,152299 + ,1 + ,1 + ,2 + ,2 + ,2 + ,3 + ,1 + ,120221 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,346485 + ,1 + ,1 + ,2 + ,3 + ,4 + ,3 + ,3 + ,145790 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,193339 + ,1 + ,1 + ,1 + ,3 + ,1 + ,2 + ,2 + ,80953 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,122774 + ,1 + ,1 + ,2 + ,1 + ,3 + ,3 + ,1 + ,130585 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,3 + ,112611 + ,1 + ,1 + ,1 + ,2 + ,1 + ,4 + ,3 + ,286468 + ,2 + ,2 + ,2 + ,2 + ,2 + ,3 + ,1 + ,241066 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,148446 + ,2 + ,2 + ,2 + ,4 + ,3 + ,3 + ,3 + ,204713 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,182079 + ,2 + ,1 + ,1 + ,3 + ,4 + ,4 + ,2 + ,140344 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,1 + ,220516 + ,1 + ,1 + ,1 + ,2 + ,2 + ,3 + ,2 + ,243060 + ,1 + ,1 + ,2 + ,3 + ,1 + ,3 + ,1 + ,162765 + ,1 + ,1 + ,2 + ,3 + ,2 + ,3 + ,2 + ,182613 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,232138 + ,2 + ,1 + ,2 + ,2 + ,3 + ,3 + ,2 + ,265318 + ,3 + ,3 + ,3 + ,4 + ,2 + ,2 + ,1 + ,85574 + ,1 + ,2 + ,1 + ,2 + ,4 + ,3 + ,3 + ,310839 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,225060 + ,2 + ,2 + ,2 + ,3 + ,1 + ,3 + ,3 + ,232317 + ,3 + ,3 + ,3 + ,2 + ,3 + ,1 + ,4 + ,144966 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,43287 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,155754 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,164709 + ,3 + ,2 + ,2 + ,2 + ,2 + ,3 + ,3 + ,201940 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,235454 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,220801 + ,2 + ,1 + ,2 + ,2 + ,1 + ,4 + ,3 + ,99466 + ,1 + ,2 + ,2 + ,1 + ,2 + ,3 + ,2 + ,92661 + ,2 + ,1 + ,3 + ,3 + ,2 + ,2 + ,3 + ,133328 + ,1 + ,1 + ,2 + ,3 + ,3 + ,4 + ,3 + ,61361 + ,2 + ,2 + ,2 + ,3 + ,3 + ,3 + ,1 + ,125930 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,100750 + ,1 + ,1 + ,1 + ,2 + ,4 + ,3 + ,1 + ,224549 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,82316 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,102010 + ,2 + ,1 + ,2 + ,1 + ,2 + ,4 + ,3 + ,101523 + ,2 + ,1 + ,1 + ,2 + ,3 + ,4 + ,3 + ,243511 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,1 + ,22938 + ,4 + ,2 + ,2 + ,3 + ,3 + ,2 + ,3 + ,41566 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,152474 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,1 + ,61857 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,99923 + ,3 + ,1 + ,2 + ,2 + ,2 + ,3 + ,4 + ,132487 + ,3 + ,3 + ,3 + ,2 + ,3 + ,2 + ,3 + ,317394 + ,1 + ,1 + ,2 + ,2 + ,2 + ,4 + ,3 + ,21054 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,209641 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,22648 + ,2 + ,1 + ,1 + ,2 + ,2 + ,3 + ,4 + ,31414 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,46698 + ,2 + ,2 + ,3 + ,2 + ,3 + ,2 + ,3 + ,131698 + ,1 + ,1 + ,1 + ,2 + ,3 + ,4 + ,1 + ,91735 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,244749 + ,2 + ,2 + ,2 + ,3 + ,2 + ,3 + ,3 + ,184510 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,79863 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,128423 + ,3 + ,1 + ,1 + ,2 + ,1 + ,3 + ,1 + ,97839 + ,4 + ,4 + ,4 + ,4 + ,4 + ,3 + ,4 + ,38214 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,151101 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,272458 + ,1 + ,1 + ,2 + ,1 + ,1 + ,4 + ,2 + ,172494 + ,3 + ,1 + ,1 + ,2 + ,2 + ,1 + ,2 + ,108043 + ,1 + ,2 + ,2 + ,2 + ,3 + ,4 + ,1 + ,328107 + ,2 + ,3 + ,3 + ,2 + ,4 + ,1 + ,3 + ,250579 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,351067 + ,2 + ,2 + ,3 + ,3 + ,3 + ,2 + ,4 + ,158015 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,98866 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,85439 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,229242 + ,3 + ,2 + ,3 + ,4 + ,4 + ,3 + ,2 + ,351619 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,84207 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,120445 + ,2 + ,1 + ,1 + ,1 + ,2 + ,4 + ,2 + ,324598 + ,3 + ,1 + ,2 + ,3 + ,1 + ,3 + ,3 + ,131069 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,204271 + ,3 + ,2 + ,3 + ,3 + ,3 + ,2 + ,4 + ,165543 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,141722 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,116048 + ,2 + ,1 + ,2 + ,3 + ,3 + ,3 + ,3 + ,250047 + ,1 + ,1 + ,1 + ,2 + ,1 + ,3 + ,1 + ,299775 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,195838 + ,3 + ,2 + ,2 + ,3 + ,4 + ,3 + ,1 + ,173260 + ,1 + ,1 + ,1 + ,1 + ,1 + ,3 + ,3 + ,254488 + ,2 + ,3 + ,4 + ,3 + ,3 + ,2 + ,1 + ,104389 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,136084 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,199476 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,92499 + ,1 + ,1 + ,2 + ,2 + ,2 + ,3 + ,3 + ,224330 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,1 + ,135781 + ,1 + ,1 + ,2 + ,4 + ,1 + ,4 + ,1 + ,74408 + ,1 + ,1 + ,2 + ,2 + ,2 + ,3 + ,1 + ,81240 + ,3 + ,4 + ,3 + ,3 + ,4 + ,1 + ,4 + ,14688 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,181633 + ,1 + ,1 + ,2 + ,1 + ,2 + ,4 + ,1 + ,271856 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,2 + ,7199 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,46660 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,17547 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,133368 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,95227 + ,2 + ,2 + ,2 + ,2 + ,2 + ,4 + ,1 + ,152601 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,98146 + ,1 + ,1 + ,2 + ,3 + ,1 + ,3 + ,2 + ,79619 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,59194 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,4 + ,139942 + ,3 + ,1 + ,2 + ,3 + ,3 + ,2 + ,2 + ,118612 + ,1 + ,1 + ,2 + ,2 + ,1 + ,4 + ,1 + ,72880 + ,2 + ,2 + ,2 + ,2 + ,2 + ,4 + ,2 + ,65475 + ,1 + ,1 + ,2 + ,2 + ,4 + ,3 + ,2 + ,99643 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,71965 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,77272 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,49289 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,135131 + ,2 + ,3 + ,3 + ,3 + ,3 + ,2 + ,3 + ,108446 + ,2 + ,2 + ,2 + ,2 + ,2 + ,3 + ,2 + ,89746 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,44296 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,77648 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,181528 + ,2 + ,1 + ,1 + ,2 + ,2 + ,4 + ,1 + ,134019 + ,1 + ,1 + ,1 + ,1 + ,1 + ,3 + ,1 + ,124064 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,92630 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,121848 + ,2 + ,3 + ,3 + ,1 + ,1 + ,3 + ,1 + ,52915 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,81872 + ,2 + ,1 + ,2 + ,2 + ,1 + ,4 + ,2 + ,58981 + ,1 + ,1 + ,2 + ,1 + ,2 + ,4 + ,1 + ,53515 + ,1 + ,1 + ,3 + ,2 + ,2 + ,4 + ,4 + ,60812 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,56375 + ,1 + ,1 + ,2 + ,1 + ,1 + ,4 + ,1 + ,65490 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,80949 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,76302 + ,1 + ,1 + ,2 + ,3 + ,3 + ,3 + ,4 + ,104011 + ,1 + ,2 + ,2 + ,3 + ,3 + ,3 + ,1 + ,98104 + ,1 + ,1 + ,1 + ,2 + ,4 + ,3 + ,3 + ,67989 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,30989 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,135458 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,73504 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,63123 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,2 + ,61254 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,74914 + ,1 + ,1 + ,1 + ,1 + ,3 + ,4 + ,3 + ,31774 + ,1 + ,1 + ,2 + ,2 + ,2 + ,4 + ,2 + ,81437 + ,0 + ,1 + ,2 + ,2 + ,2 + ,4 + ,3 + ,87186 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,50090 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,65745 + ,1 + ,1 + ,2 + ,2 + ,3 + ,3 + ,2 + ,56653 + ,1 + ,1 + ,2 + ,1 + ,3 + ,3 + ,2 + ,158399 + ,2 + ,1 + ,2 + ,2 + ,3 + ,3 + ,1 + ,46455 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,73624 + ,1 + ,1 + ,2 + ,1 + ,1 + ,3 + ,1 + ,38395 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,91899 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,139526 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,52164 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,51567 + ,3 + ,2 + ,3 + ,2 + ,3 + ,3 + ,1 + ,70551 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,84856 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,102538 + ,1 + ,1 + ,2 + ,1 + ,1 + ,3 + ,2 + ,86678 + ,1 + ,1 + ,2 + ,2 + ,3 + ,3 + ,1 + ,85709 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,34662 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,150580 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,3 + ,99611 + ,2 + ,1 + ,2 + ,2 + ,1 + ,3 + ,3 + ,19349 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,99373 + ,3 + ,1 + ,3 + ,2 + ,4 + ,3 + ,2 + ,86230 + ,3 + ,3 + ,2 + ,2 + ,1 + ,4 + ,3 + ,30837 + ,1 + ,1 + ,2 + ,3 + ,3 + ,4 + ,4 + ,31706 + ,2 + ,1 + ,1 + ,2 + ,1 + ,4 + ,3 + ,89806 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,62088 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,40151 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,27634 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,76990 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,37460 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,54157 + ,1 + ,1 + ,1 + ,1 + ,2 + ,3 + ,3 + ,49862 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,84337 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,64175 + ,2 + ,2 + ,2 + ,2 + ,3 + ,3 + ,2 + ,59382 + ,2 + ,2 + ,2 + ,2 + ,3 + ,4 + ,2 + ,119308 + ,2 + ,1 + ,1 + ,1 + ,3 + ,4 + ,1 + ,76702 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,103425 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,70344 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,43410 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,104838 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,62215 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,69304 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,53117 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,19764 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,86680 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,84105 + ,3 + ,2 + ,3 + ,2 + ,1 + ,3 + ,2 + ,77945 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,89113 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,91005 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,40248 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,64187 + ,1 + ,2 + ,2 + ,2 + ,3 + ,3 + ,2 + ,50857 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,56613 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,62792 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,72535 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA) + ,dim=c(8 + ,289) + ,dimnames=list(c('Time' + ,'Q1' + ,'Q2' + ,'Q3' + ,'Q4' + ,'Q5' + ,'Q6' + ,'Q7') + ,1:289)) > y <- array(NA,dim=c(8,289),dimnames=list(c('Time','Q1','Q2','Q3','Q4','Q5','Q6','Q7'),1:289)) > 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' > 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, 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 Time Q1 Q2 Q3 Q4 Q5 Q6 Q7 1 210907 1 1 2 2 3 4 3 2 120982 NA NA NA NA NA NA NA 3 176508 2 2 2 4 2 2 1 4 179321 NA NA NA NA NA NA NA 5 123185 NA NA NA NA NA NA NA 6 52746 NA NA NA NA NA NA NA 7 385534 1 1 2 2 2 3 2 8 33170 NA NA NA NA NA NA NA 9 101645 NA NA NA NA NA NA NA 10 149061 1 1 2 2 3 4 1 11 165446 2 2 2 3 3 3 2 12 237213 1 1 1 2 3 4 1 13 173326 NA NA NA NA NA NA NA 14 133131 2 1 1 1 1 4 2 15 258873 NA NA NA NA NA NA NA 16 180083 NA NA NA NA NA NA NA 17 324799 1 1 1 1 1 3 1 18 230964 4 4 4 4 3 2 4 19 236785 2 1 1 1 2 4 2 20 135473 NA NA NA NA NA NA NA 21 202925 NA NA NA NA NA NA NA 22 215147 1 1 1 1 1 4 1 23 344297 1 1 1 4 1 4 1 24 153935 NA NA NA NA NA NA NA 25 132943 NA NA NA NA NA NA NA 26 174724 1 1 1 1 1 3 1 27 174415 1 1 1 3 3 3 2 28 225548 4 4 4 4 3 3 4 29 223632 1 1 1 2 2 3 2 30 124817 1 1 1 1 1 4 2 31 221698 NA NA NA NA NA NA NA 32 210767 4 3 3 2 2 2 3 33 170266 1 1 1 2 1 4 2 34 260561 NA NA NA NA NA NA NA 35 84853 NA NA NA NA NA NA NA 36 294424 1 1 1 1 2 1 2 37 101011 NA NA NA NA NA NA NA 38 215641 NA NA NA NA NA NA NA 39 325107 2 1 2 2 2 3 3 40 7176 2 1 2 2 3 3 2 41 167542 NA NA NA NA NA NA NA 42 106408 3 4 4 2 4 2 4 43 96560 2 3 3 3 3 2 4 44 265769 0 0 0 0 0 0 0 45 269651 NA NA NA NA NA NA NA 46 149112 2 1 2 2 2 3 3 47 175824 1 1 1 2 1 1 1 48 152871 1 1 1 3 1 4 1 49 111665 1 1 1 1 2 4 1 50 116408 NA NA NA NA NA NA NA 51 362301 2 3 3 3 2 3 2 52 78800 NA NA NA NA NA NA NA 53 183167 2 3 3 3 2 3 3 54 277965 NA NA NA NA NA NA NA 55 150629 NA NA NA NA NA NA NA 56 168809 NA NA NA NA NA NA NA 57 24188 NA NA NA NA NA NA NA 58 329267 2 3 2 2 2 1 1 59 65029 NA NA NA NA NA NA NA 60 101097 NA NA NA NA NA NA NA 61 218946 3 3 3 4 4 2 3 62 244052 1 1 1 2 2 4 1 63 341570 2 2 2 2 2 3 2 64 103597 2 1 2 2 3 4 2 65 233328 NA NA NA NA NA NA NA 66 256462 2 2 2 3 3 2 2 67 206161 NA NA NA NA NA NA NA 68 311473 NA NA NA NA NA NA NA 69 235800 NA NA NA NA NA NA NA 70 177939 NA NA NA NA NA NA NA 71 207176 NA NA NA NA NA NA NA 72 196553 2 1 1 2 2 4 1 73 174184 NA NA NA NA NA NA NA 74 143246 1 1 2 1 1 4 1 75 187559 NA NA NA NA NA NA NA 76 187681 2 2 2 3 3 2 2 77 119016 NA NA NA NA NA NA NA 78 182192 NA NA NA NA NA NA NA 79 73566 1 1 1 3 3 4 1 80 194979 NA NA NA NA NA NA NA 81 167488 1 1 2 1 2 3 1 82 143756 3 1 2 3 4 3 4 83 275541 NA NA NA NA NA NA NA 84 243199 1 2 2 3 2 3 3 85 182999 NA NA NA NA NA NA NA 86 135649 NA NA NA NA NA NA NA 87 152299 1 1 2 2 2 3 1 88 120221 NA NA NA NA NA NA NA 89 346485 1 1 2 3 4 3 3 90 145790 NA NA NA NA NA NA NA 91 193339 1 1 1 3 1 2 2 92 80953 NA NA NA NA NA NA NA 93 122774 1 1 2 1 3 3 1 94 130585 1 1 1 2 2 4 3 95 112611 1 1 1 2 1 4 3 96 286468 2 2 2 2 2 3 1 97 241066 NA NA NA NA NA NA NA 98 148446 2 2 2 4 3 3 3 99 204713 NA NA NA NA NA NA NA 100 182079 2 1 1 3 4 4 2 101 140344 1 1 1 1 2 4 1 102 220516 1 1 1 2 2 3 2 103 243060 1 1 2 3 1 3 1 104 162765 1 1 2 3 2 3 2 105 182613 NA NA NA NA NA NA NA 106 232138 2 1 2 2 3 3 2 107 265318 3 3 3 4 2 2 1 108 85574 1 2 1 2 4 3 3 109 310839 NA NA NA NA NA NA NA 110 225060 2 2 2 3 1 3 3 111 232317 3 3 3 2 3 1 4 112 144966 NA NA NA NA NA NA NA 113 43287 NA NA NA NA NA NA NA 114 155754 NA NA NA NA NA NA NA 115 164709 3 2 2 2 2 3 3 116 201940 NA NA NA NA NA NA NA 117 235454 NA NA NA NA NA NA NA 118 220801 2 1 2 2 1 4 3 119 99466 1 2 2 1 2 3 2 120 92661 2 1 3 3 2 2 3 121 133328 1 1 2 3 3 4 3 122 61361 2 2 2 3 3 3 1 123 125930 NA NA NA NA NA NA NA 124 100750 1 1 1 2 4 3 1 125 224549 NA NA NA NA NA NA NA 126 82316 NA NA NA NA NA NA NA 127 102010 2 1 2 1 2 4 3 128 101523 2 1 1 2 3 4 3 129 243511 1 1 1 1 1 4 1 130 22938 4 2 2 3 3 2 3 131 41566 NA NA NA NA NA NA NA 132 152474 1 1 1 1 2 4 1 133 61857 NA NA NA NA NA NA NA 134 99923 3 1 2 2 2 3 4 135 132487 3 3 3 2 3 2 3 136 317394 1 1 2 2 2 4 3 137 21054 NA NA NA NA NA NA NA 138 209641 NA NA NA NA NA NA NA 139 22648 2 1 1 2 2 3 4 140 31414 NA NA NA NA NA NA NA 141 46698 2 2 3 2 3 2 3 142 131698 1 1 1 2 3 4 1 143 91735 NA NA NA NA NA NA NA 144 244749 2 2 2 3 2 3 3 145 184510 NA NA NA NA NA NA NA 146 79863 NA NA NA NA NA NA NA 147 128423 3 1 1 2 1 3 1 148 97839 4 4 4 4 4 3 4 149 38214 NA NA NA NA NA NA NA 150 151101 NA NA NA NA NA NA NA 151 272458 1 1 2 1 1 4 2 152 172494 3 1 1 2 2 1 2 153 108043 1 2 2 2 3 4 1 154 328107 2 3 3 2 4 1 3 155 250579 NA NA NA NA NA NA NA 156 351067 2 2 3 3 3 2 4 157 158015 NA NA NA NA NA NA NA 158 98866 NA NA NA NA NA NA NA 159 85439 NA NA NA NA NA NA NA 160 229242 3 2 3 4 4 3 2 161 351619 NA NA NA NA NA NA NA 162 84207 NA NA NA NA NA NA NA 163 120445 2 1 1 1 2 4 2 164 324598 3 1 2 3 1 3 3 165 131069 1 1 1 2 2 4 2 166 204271 3 2 3 3 3 2 4 167 165543 NA NA NA NA NA NA NA 168 141722 NA NA NA NA NA NA NA 169 116048 2 1 2 3 3 3 3 170 250047 1 1 1 2 1 3 1 171 299775 NA NA NA NA NA NA NA 172 195838 3 2 2 3 4 3 1 173 173260 1 1 1 1 1 3 3 174 254488 2 3 4 3 3 2 1 175 104389 NA NA NA NA NA NA NA 176 136084 NA NA NA NA NA NA NA 177 199476 NA NA NA NA NA NA NA 178 92499 1 1 2 2 2 3 3 179 224330 1 1 1 1 2 4 1 180 135781 1 1 2 4 1 4 1 181 74408 1 1 2 2 2 3 1 182 81240 3 4 3 3 4 1 4 183 14688 NA NA NA NA NA NA NA 184 181633 1 1 2 1 2 4 1 185 271856 1 1 1 1 2 4 2 186 7199 NA NA NA NA NA NA NA 187 46660 NA NA NA NA NA NA NA 188 17547 NA NA NA NA NA NA NA 189 133368 NA NA NA NA NA NA NA 190 95227 2 2 2 2 2 4 1 191 152601 NA NA NA NA NA NA NA 192 98146 1 1 2 3 1 3 2 193 79619 NA NA NA NA NA NA NA 194 59194 1 1 1 1 1 4 4 195 139942 3 1 2 3 3 2 2 196 118612 1 1 2 2 1 4 1 197 72880 2 2 2 2 2 4 2 198 65475 1 1 2 2 4 3 2 199 99643 NA NA NA NA NA NA NA 200 71965 NA NA NA NA NA NA NA 201 77272 NA NA NA NA NA NA NA 202 49289 NA NA NA NA NA NA NA 203 135131 2 3 3 3 3 2 3 204 108446 2 2 2 2 2 3 2 205 89746 NA NA NA NA NA NA NA 206 44296 NA NA NA NA NA NA NA 207 77648 NA NA NA NA NA NA NA 208 181528 2 1 1 2 2 4 1 209 134019 1 1 1 1 1 3 1 210 124064 NA NA NA NA NA NA NA 211 92630 NA NA NA NA NA NA NA 212 121848 2 3 3 1 1 3 1 213 52915 NA NA NA NA NA NA NA 214 81872 2 1 2 2 1 4 2 215 58981 1 1 2 1 2 4 1 216 53515 1 1 3 2 2 4 4 217 60812 NA NA NA NA NA NA NA 218 56375 1 1 2 1 1 4 1 219 65490 NA NA NA NA NA NA NA 220 80949 NA NA NA NA NA NA NA 221 76302 1 1 2 3 3 3 4 222 104011 1 2 2 3 3 3 1 223 98104 1 1 1 2 4 3 3 224 67989 NA NA NA NA NA NA NA 225 30989 NA NA NA NA NA NA NA 226 135458 NA NA NA NA NA NA NA 227 73504 NA NA NA NA NA NA NA 228 63123 1 1 1 1 1 4 2 229 61254 NA NA NA NA NA NA NA 230 74914 1 1 1 1 3 4 3 231 31774 1 1 2 2 2 4 2 232 81437 0 1 2 2 2 4 3 233 87186 NA NA NA NA NA NA NA 234 50090 NA NA NA NA NA NA NA 235 65745 1 1 2 2 3 3 2 236 56653 1 1 2 1 3 3 2 237 158399 2 1 2 2 3 3 1 238 46455 NA NA NA NA NA NA NA 239 73624 1 1 2 1 1 3 1 240 38395 NA NA NA NA NA NA NA 241 91899 NA NA NA NA NA NA NA 242 139526 NA NA NA NA NA NA NA 243 52164 NA NA NA NA NA NA NA 244 51567 3 2 3 2 3 3 1 245 70551 NA NA NA NA NA NA NA 246 84856 NA NA NA NA NA NA NA 247 102538 1 1 2 1 1 3 2 248 86678 1 1 2 2 3 3 1 249 85709 NA NA NA NA NA NA NA 250 34662 NA NA NA NA NA NA NA 251 150580 1 1 1 2 2 4 3 252 99611 2 1 2 2 1 3 3 253 19349 NA NA NA NA NA NA NA 254 99373 3 1 3 2 4 3 2 255 86230 3 3 2 2 1 4 3 256 30837 1 1 2 3 3 4 4 257 31706 2 1 1 2 1 4 3 258 89806 NA NA NA NA NA NA NA 259 62088 NA NA NA NA NA NA NA 260 40151 NA NA NA NA NA NA NA 261 27634 NA NA NA NA NA NA NA 262 76990 NA NA NA NA NA NA NA 263 37460 NA NA NA NA NA NA NA 264 54157 1 1 1 1 2 3 3 265 49862 NA NA NA NA NA NA NA 266 84337 NA NA NA NA NA NA NA 267 64175 2 2 2 2 3 3 2 268 59382 2 2 2 2 3 4 2 269 119308 2 1 1 1 3 4 1 270 76702 NA NA NA NA NA NA NA 271 103425 NA NA NA NA NA NA NA 272 70344 NA NA NA NA NA NA NA 273 43410 NA NA NA NA NA NA NA 274 104838 NA NA NA NA NA NA NA 275 62215 NA NA NA NA NA NA NA 276 69304 NA NA NA NA NA NA NA 277 53117 NA NA NA NA NA NA NA 278 19764 NA NA NA NA NA NA NA 279 86680 NA NA NA NA NA NA NA 280 84105 3 2 3 2 1 3 2 281 77945 NA NA NA NA NA NA NA 282 89113 NA NA NA NA NA NA NA 283 91005 NA NA NA NA NA NA NA 284 40248 NA NA NA NA NA NA NA 285 64187 1 2 2 2 3 3 2 286 50857 NA NA NA NA NA NA NA 287 56613 NA NA NA NA NA NA NA 288 62792 NA NA NA NA NA NA NA 289 72535 NA NA NA NA NA NA NA > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Q1 Q2 Q3 Q4 Q5 263754 -6604 14991 -13540 22773 -17537 Q6 Q7 -24895 -11074 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -149995 -59081 -9182 48738 226836 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 263754 41573 6.344 2.94e-09 *** Q1 -6604 10865 -0.608 0.54430 Q2 14991 13670 1.097 0.27471 Q3 -13540 13416 -1.009 0.31460 Q4 22773 9254 2.461 0.01509 * Q5 -17537 7882 -2.225 0.02770 * Q6 -24895 8722 -2.854 0.00497 ** Q7 -11074 7357 -1.505 0.13452 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 80780 on 139 degrees of freedom (142 observations deleted due to missingness) Multiple R-squared: 0.1435, Adjusted R-squared: 0.1004 F-statistic: 3.327 on 7 and 139 DF, p-value: 0.00263 > 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 + } Error in if (gqarr[mypoint - kp3 + 1, 2] < 0.01) numsignificant1 <- numsignificant1 + : missing value where TRUE/FALSE needed Execution halted