R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(17080 + ,9 + ,57 + ,0 + ,0 + ,17240 + ,8 + ,68 + ,0 + ,0 + ,17200 + ,7 + ,55 + ,0 + ,0 + ,15580 + ,9 + ,53 + ,0 + ,1 + ,18950 + ,9 + ,57 + ,0 + ,1 + ,23360 + ,8 + ,61 + ,0 + ,0 + ,19190 + ,6 + ,58 + ,0 + ,0 + ,14150 + ,6 + ,56 + ,0 + ,0 + ,19870 + ,8 + ,59 + ,0 + ,0 + ,19420 + ,9 + ,60 + ,0 + ,0 + ,16020 + ,6 + ,53 + ,0 + ,0 + ,17350 + ,8 + ,54 + ,0 + ,1 + ,21930 + ,8 + ,59 + ,0 + ,0 + ,21180 + ,8 + ,61 + ,0 + ,1 + ,22580 + ,8 + ,58 + ,0 + ,1 + ,19320 + ,7 + ,53 + ,0 + ,1 + ,14720 + ,5 + ,50 + ,0 + ,1 + ,18780 + ,6 + ,53 + ,0 + ,0 + ,23520 + ,9 + ,59 + ,0 + ,1 + ,26040 + ,9 + ,62 + ,0 + ,1 + ,14000 + ,5 + ,49 + ,0 + ,0 + ,12560 + ,5 + ,53 + ,0 + ,0 + ,25780 + ,7 + ,63 + ,0 + ,1 + ,29880 + ,9 + ,65 + ,0 + ,0 + ,14040 + ,3 + ,52 + ,0 + ,1 + ,23480 + ,9 + ,60 + ,0 + ,1 + ,17550 + ,5 + ,52 + ,0 + ,1 + ,29800 + ,8 + ,60 + ,0 + ,0 + ,21000 + ,9 + ,60 + ,0 + ,0 + ,12820 + ,5 + ,49 + ,0 + ,0 + ,30000 + ,9 + ,66 + ,0 + ,1 + ,26730 + ,8 + ,60 + ,0 + ,0 + ,20930 + ,7 + ,58 + ,0 + ,0 + ,16120 + ,5 + ,52 + ,0 + ,0 + ,21750 + ,8 + ,59 + ,0 + ,0 + ,27250 + ,9 + ,59 + ,0 + ,1 + ,20710 + ,8 + ,55 + ,0 + ,1 + ,15470 + ,8 + ,57 + ,0 + ,1 + ,20040 + ,8 + ,57 + ,0 + ,1 + ,31350 + ,9 + ,60 + ,0 + ,0 + ,24200 + ,8 + ,59 + ,0 + ,1 + ,17760 + ,5 + ,51 + ,0 + ,1 + ,19310 + ,8 + ,57 + ,0 + ,0 + ,13430 + ,5 + ,50 + ,0 + ,0 + ,20760 + ,9 + ,57 + ,0 + ,0 + ,16240 + ,7 + ,54 + ,0 + ,1 + ,13440 + ,8 + ,53 + ,0 + ,0 + ,16500 + ,6 + ,55 + ,0 + ,1 + ,27320 + ,8 + ,61 + ,0 + ,1 + ,20170 + ,5 + ,55 + ,0 + ,1 + ,27970 + ,9 + ,62 + ,0 + ,0 + ,35560 + ,9 + ,62 + ,0 + ,1 + ,17030 + ,8 + ,55 + ,0 + ,1 + ,16340 + ,6 + ,54 + ,0 + ,1 + ,25700 + ,9 + ,57 + ,0 + ,1 + ,30160 + ,9 + ,63 + ,0 + ,0 + ,24190 + ,7 + ,60 + ,0 + ,0 + ,15690 + ,4 + ,50 + ,0 + ,0 + ,16980 + ,8 + ,58 + ,0 + ,0 + ,21230 + ,8 + ,60 + ,0 + ,1 + ,24810 + ,8 + ,60 + ,0 + ,0 + ,14810 + ,6 + ,51 + ,0 + ,0 + ,15770 + ,4 + ,49 + ,0 + ,0 + ,19400 + ,8 + ,59 + ,0 + ,1 + ,17470 + ,6 + ,58 + ,0 + ,1 + ,20690 + ,9 + ,58 + ,0 + ,1 + ,16310 + ,7 + ,56 + ,0 + ,0 + ,15360 + ,5 + ,52 + ,0 + ,0 + ,25600 + ,9 + ,61 + ,0 + ,0 + ,19620 + ,8 + ,57 + ,0 + ,1 + ,25310 + ,8 + ,58 + ,0 + ,0 + ,27150 + ,9 + ,60 + ,0 + ,1 + ,24570 + ,9 + ,59 + ,0 + ,1 + ,20900 + ,9 + ,60 + ,0 + ,1 + ,17890 + ,7 + ,56 + ,0 + ,1 + ,18580 + ,5 + ,53 + ,0 + ,1 + ,14520 + ,5 + ,51 + ,0 + ,1 + ,38420 + ,9 + ,69 + ,0 + ,1 + ,17190 + ,6 + ,53 + ,0 + ,0 + ,21110 + ,7 + ,57 + ,0 + ,0 + ,16950 + ,6 + ,53 + ,0 + ,0 + ,22110 + ,8 + ,63 + ,0 + ,1 + ,17940 + ,8 + ,55 + ,0 + ,1 + ,19170 + ,7 + ,58 + ,0 + ,0 + ,21440 + ,8 + ,63 + ,0 + ,0 + ,12530 + ,7 + ,52 + ,0 + ,1 + ,26590 + ,9 + ,62 + ,0 + ,1 + ,15800 + ,5 + ,53 + ,0 + ,1 + ,21260 + ,9 + ,62 + ,0 + ,1 + ,30290 + ,9 + ,62 + ,0 + ,0 + ,29640 + ,9 + ,65 + ,0 + ,1 + ,16110 + ,7 + ,58 + ,0 + ,1 + ,22150 + ,8 + ,60 + ,0 + ,0 + ,23880 + ,8 + ,60 + ,0 + ,0 + ,21960 + ,9 + ,61 + ,0 + ,1 + ,17510 + ,9 + ,58 + ,0 + ,1 + ,21650 + ,9 + ,62 + ,0 + ,1 + ,16820 + ,7 + ,55 + ,0 + ,1 + ,15230 + ,8 + ,55 + ,0 + ,1 + ,12920 + ,8 + ,52 + ,0 + ,0 + ,16490 + ,7 + ,54 + ,0 + ,1 + ,25880 + ,9 + ,63 + ,0 + ,1 + ,25740 + ,9 + ,61 + ,0 + ,0 + ,19790 + ,6 + ,56 + ,0 + ,1 + ,23540 + ,8 + ,59 + ,0 + ,1 + ,17180 + ,6 + ,55 + ,0 + ,1 + ,17420 + ,7 + ,59 + ,0 + ,0 + ,16030 + ,7 + ,51 + ,0 + ,0 + ,26390 + ,8 + ,60 + ,0 + ,0 + ,18290 + ,7 + ,54 + ,0 + ,0 + ,20840 + ,7 + ,58 + ,0 + ,1 + ,22200 + ,7 + ,58 + ,0 + ,1 + ,14730 + ,7 + ,53 + ,0 + ,0 + ,23410 + ,8 + ,61 + ,0 + ,0 + ,16980 + ,7 + ,55 + ,0 + ,0 + ,11960 + ,5 + ,47 + ,0 + ,0 + ,18720 + ,8 + ,57 + ,0 + ,0 + ,22190 + ,7 + ,55 + ,0 + ,1 + ,24200 + ,9 + ,57 + ,0 + ,1 + ,18270 + ,7 + ,55 + ,0 + ,0 + ,14610 + ,7 + ,54 + ,0 + ,0 + ,13380 + ,6 + ,53 + ,0 + ,1 + ,20900 + ,8 + ,57 + ,0 + ,1 + ,16970 + ,8 + ,59 + ,0 + ,0 + ,15620 + ,8 + ,55 + ,0 + ,1 + ,20400 + ,9 + ,56 + ,0 + ,0 + ,16090 + ,7 + ,52 + ,0 + ,0 + ,24580 + ,8 + ,61 + ,0 + ,0 + ,26500 + ,9 + ,64 + ,0 + ,1 + ,14290 + ,8 + ,58 + ,0 + ,1 + ,16750 + ,8 + ,53 + ,0 + ,1 + ,19470 + ,9 + ,57 + ,0 + ,0 + ,20690 + ,8 + ,54 + ,0 + ,1 + ,15720 + ,6 + ,52 + ,0 + ,1 + ,13480 + ,6 + ,53 + ,0 + ,1 + ,22880 + ,8 + ,62 + ,0 + ,0 + ,17730 + ,9 + ,59 + ,0 + ,1 + ,19050 + ,7 + ,58 + ,0 + ,1 + ,24630 + ,9 + ,61 + ,0 + ,0 + ,14310 + ,6 + ,51 + ,0 + ,1 + ,26310 + ,9 + ,62 + ,0 + ,0 + ,31140 + ,9 + ,65 + ,0 + ,1 + ,21350 + ,9 + ,59 + ,0 + ,1 + ,15270 + ,6 + ,53 + ,0 + ,1 + ,22930 + ,8 + ,58 + ,0 + ,0 + ,30420 + ,9 + ,66 + ,0 + ,0 + ,29270 + ,8 + ,64 + ,0 + ,1 + ,26650 + ,8 + ,64 + ,0 + ,0 + ,23010 + ,9 + ,59 + ,0 + ,1 + ,24600 + ,9 + ,64 + ,0 + ,1 + ,25920 + ,9 + ,61 + ,0 + ,0 + ,17500 + ,7 + ,55 + ,0 + ,0 + ,17590 + ,8 + ,53 + ,0 + ,1 + ,15360 + ,6 + ,48 + ,0 + ,1 + ,22590 + ,9 + ,59 + ,0 + ,0 + ,20480 + ,9 + ,65 + ,0 + ,0 + ,25710 + ,9 + ,61 + ,0 + ,1 + ,20460 + ,7 + ,56 + ,0 + ,1 + ,17800 + ,8 + ,59 + ,0 + ,0 + ,15520 + ,5 + ,54 + ,0 + ,0 + ,19530 + ,8 + ,58 + ,0 + ,0 + ,28930 + ,9 + ,65 + ,0 + ,1 + ,17130 + ,6 + ,51 + ,0 + ,1 + ,28510 + ,9 + ,60 + ,0 + ,0 + ,16240 + ,6 + ,52 + ,0 + ,1 + ,26310 + ,8 + ,59 + ,0 + ,1 + ,18190 + ,5 + ,53 + ,0 + ,1 + ,16580 + ,7 + ,53 + ,0 + ,1 + ,21580 + ,7 + ,54 + ,0 + ,1 + ,17890 + ,4 + ,52 + ,0 + ,1 + ,30040 + ,9 + ,64 + ,0 + ,0 + ,25030 + ,8 + ,63 + ,0 + ,1 + ,19330 + ,9 + ,58 + ,0 + ,0 + ,20910 + ,9 + ,59 + ,0 + ,0 + ,23160 + ,9 + ,60 + ,0 + ,0 + ,17040 + ,5 + ,51 + ,0 + ,0 + ,16060 + ,9 + ,58 + ,0 + ,0 + ,11650 + ,7 + ,47 + ,0 + ,1 + ,21020 + ,6 + ,56 + ,0 + ,0 + ,23200 + ,9 + ,57 + ,0 + ,0 + ,22300 + ,9 + ,61 + ,0 + ,1 + ,17160 + ,9 + ,56 + ,0 + ,1 + ,17900 + ,7 + ,54 + ,0 + ,1 + ,11460 + ,5 + ,50 + ,0 + ,0 + ,21870 + ,8 + ,62 + ,0 + ,0 + ,27170 + ,9 + ,62 + ,0 + ,1 + ,17960 + ,7 + ,55 + ,0 + ,1 + ,19530 + ,9 + ,58 + ,1 + ,1 + ,13350 + ,8 + ,57 + ,0 + ,0 + ,21190 + ,9 + ,57 + ,0 + ,1 + ,16660 + ,6 + ,52 + ,0 + ,1 + ,18260 + ,6 + ,53 + ,0 + ,1 + ,27090 + ,8 + ,63 + ,0 + ,0 + ,28710 + ,9 + ,65 + ,0 + ,1 + ,10920 + ,5 + ,50 + ,0 + ,0 + ,22620 + ,6 + ,58 + ,0 + ,1 + ,21040 + ,6 + ,57 + ,0 + ,1 + ,21660 + ,9 + ,58 + ,0 + ,0 + ,16900 + ,7 + ,54 + ,0 + ,0 + ,29730 + ,9 + ,60 + ,0 + ,1 + ,21450 + ,8 + ,60 + ,0 + ,0 + ,19710 + ,5 + ,58 + ,0 + ,1 + ,20950 + ,7 + ,57 + ,0 + ,0 + ,16970 + ,6 + ,55 + ,0 + ,0 + ,24550 + ,9 + ,60 + ,0 + ,0 + ,19200 + ,7 + ,57 + ,0 + ,1 + ,21640 + ,9 + ,60 + ,0 + ,1 + ,21300 + ,9 + ,59 + ,0 + ,0 + ,29930 + ,8 + ,63 + ,0 + ,0 + ,25290 + ,9 + ,59 + ,0 + ,0 + ,17260 + ,7 + ,53 + ,0 + ,0 + ,24420 + ,9 + ,62 + ,0 + ,0 + ,11020 + ,4 + ,48 + ,0 + ,0 + ,20560 + ,9 + ,63 + ,0 + ,0 + ,18080 + ,5 + ,56 + ,0 + ,1 + ,23050 + ,8 + ,65 + ,0 + ,0 + ,19690 + ,9 + ,59 + ,0 + ,0 + ,15560 + ,8 + ,59 + ,0 + ,0 + ,10720 + ,3 + ,46 + ,0 + ,0 + ,20420 + ,9 + ,62 + ,0 + ,1 + ,15120 + ,8 + ,53 + ,0 + ,0 + ,14230 + ,6 + ,50 + ,0 + ,1 + ,36810 + ,9 + ,68 + ,0 + ,1 + ,19910 + ,8 + ,60 + ,0 + ,1 + ,18970 + ,8 + ,56 + ,0 + ,1 + ,13700 + ,7 + ,55 + ,0 + ,0 + ,13380 + ,6 + ,52 + ,0 + ,0 + ,20160 + ,8 + ,56 + ,0 + ,1 + ,26390 + ,9 + ,63 + ,0 + ,0 + ,13890 + ,4 + ,48 + ,0 + ,0 + ,16120 + ,7 + ,57 + ,0 + ,1 + ,21350 + ,8 + ,59 + ,0 + ,0 + ,26810 + ,8 + ,61 + ,0 + ,1 + ,32230 + ,9 + ,65 + ,0 + ,0 + ,17960 + ,6 + ,55 + ,0 + ,0 + ,20100 + ,8 + ,55 + ,0 + ,1 + ,15230 + ,6 + ,51 + ,0 + ,0 + ,17440 + ,8 + ,53 + ,0 + ,1 + ,24850 + ,9 + ,64 + ,0 + ,0 + ,23350 + ,8 + ,59 + ,0 + ,0 + ,14150 + ,7 + ,54 + ,0 + ,0 + ,20760 + ,9 + ,61 + ,0 + ,1 + ,24350 + ,8 + ,60 + ,0 + ,1 + ,17280 + ,7 + ,57 + ,0 + ,0 + ,28500 + ,9 + ,63 + ,0 + ,0 + ,18440 + ,8 + ,57 + ,0 + ,0 + ,17540 + ,9 + ,62 + ,0 + ,0 + ,13430 + ,6 + ,52 + ,0 + ,0 + ,23030 + ,8 + ,57 + ,0 + ,1 + ,22460 + ,9 + ,64 + ,0 + ,1 + ,24760 + ,8 + ,63 + ,0 + ,0 + ,32390 + ,9 + ,65 + ,0 + ,1 + ,24570 + ,9 + ,62 + ,0 + ,1 + ,23820 + ,8 + ,62 + ,0 + ,0 + ,16400 + ,7 + ,55 + ,0 + ,0 + ,15890 + ,5 + ,51 + ,0 + ,0 + ,20560 + ,7 + ,54 + ,0 + ,1 + ,22260 + ,8 + ,57 + ,0 + ,1 + ,18860 + ,9 + ,56 + ,0 + ,0 + ,28330 + ,9 + ,62 + ,0 + ,1 + ,17150 + ,6 + ,53 + ,0 + ,1 + ,26310 + ,8 + ,59 + ,0 + ,1 + ,25500 + ,7 + ,56 + ,0 + ,1 + ,19120 + ,9 + ,59 + ,0 + ,0 + ,18770 + ,7 + ,53 + ,0 + ,0 + ,19350 + ,7 + ,53 + ,0 + ,0 + ,15390 + ,5 + ,50 + ,0 + ,0 + ,28030 + ,9 + ,60 + ,0 + ,1 + ,29230 + ,9 + ,64 + ,0 + ,1 + ,23580 + ,8 + ,61 + ,0 + ,0 + ,20940 + ,8 + ,58 + ,0 + ,1 + ,18550 + ,9 + ,60 + ,0 + ,1 + ,15350 + ,6 + ,55 + ,0 + ,0 + ,21350 + ,7 + ,56 + ,0 + ,1 + ,19300 + ,5 + ,51 + ,0 + ,1 + ,21820 + ,9 + ,60 + ,0 + ,0 + ,13590 + ,5 + ,51 + ,0 + ,1 + ,20020 + ,7 + ,58 + ,0 + ,0 + ,16990 + ,6 + ,54 + ,0 + ,1 + ,25000 + ,8 + ,57 + ,0 + ,1 + ,23660 + ,7 + ,58 + ,0 + ,0 + ,20690 + ,8 + ,60 + ,0 + ,0 + ,14180 + ,4 + ,49 + ,0 + ,0 + ,23330 + ,8 + ,57 + ,0 + ,0 + ,15140 + ,5 + ,52 + ,0 + ,1 + ,17580 + ,8 + ,52 + ,0 + ,0 + ,25350 + ,7 + ,60 + ,0 + ,1 + ,25640 + ,7 + ,58 + ,0 + ,0 + ,24870 + ,9 + ,64 + ,0 + ,0 + ,15910 + ,9 + ,57 + ,0 + ,0 + ,16240 + ,8 + ,53 + ,0 + ,1 + ,27980 + ,9 + ,62 + ,0 + ,1 + ,16910 + ,6 + ,53 + ,0 + ,1 + ,19990 + ,8 + ,57 + ,0 + ,0 + ,18690 + ,9 + ,57 + ,0 + ,1 + ,10040 + ,4 + ,48 + ,0 + ,1 + ,14270 + ,6 + ,50 + ,0 + ,1 + ,18260 + ,7 + ,51 + ,0 + ,1 + ,26880 + ,9 + ,60 + ,0 + ,0 + ,16570 + ,8 + ,56 + ,0 + ,1 + ,16720 + ,6 + ,54 + ,0 + ,0 + ,20150 + ,8 + ,58 + ,0 + ,0 + ,23710 + ,7 + ,56 + ,0 + ,0 + ,21150 + ,5 + ,50 + ,0 + ,1 + ,23280 + ,8 + ,60 + ,0 + ,0 + ,14950 + ,7 + ,57 + ,0 + ,0 + ,28840 + ,11 + ,69 + ,0 + ,1 + ,23280 + ,10 + ,64 + ,0 + ,1 + ,33810 + ,14 + ,63 + ,0 + ,1 + ,21700 + ,11 + ,58 + ,0 + ,0 + ,34700 + ,11 + ,67 + ,0 + ,1 + ,30580 + ,12 + ,61 + ,0 + ,0 + ,18110 + ,10 + ,57 + ,0 + ,1 + ,25240 + ,11 + ,64 + ,0 + ,1 + ,26420 + ,10 + ,61 + ,0 + ,0 + ,37410 + ,14 + ,69 + ,0 + ,1 + ,43360 + ,13 + ,70 + ,0 + ,1 + ,48420 + ,14 + ,72 + ,0 + ,1 + ,45500 + ,12 + ,71 + ,0 + ,1 + ,28410 + ,12 + ,63 + ,0 + ,0 + ,31660 + ,10 + ,62 + ,0 + ,0 + ,38160 + ,13 + ,64 + ,0 + ,0 + ,25610 + ,10 + ,62 + ,0 + ,1 + ,36540 + ,11 + ,65 + ,0 + ,0 + ,24810 + ,10 + ,61 + ,0 + ,1 + ,26650 + ,11 + ,63 + ,0 + ,0 + ,32030 + ,10 + ,66 + ,0 + ,1 + ,35490 + ,13 + ,68 + ,0 + ,1 + ,22360 + ,14 + ,66 + ,1 + ,0 + ,32220 + ,11 + ,72 + ,0 + ,1 + ,31110 + ,10 + ,66 + ,0 + ,1 + ,34900 + ,11 + ,67 + ,0 + ,0 + ,31470 + ,13 + ,64 + ,0 + ,0 + ,25200 + ,10 + ,61 + ,0 + ,0 + ,22920 + ,10 + ,63 + ,0 + ,1 + ,28890 + ,12 + ,64 + ,0 + ,0 + ,22460 + ,10 + ,61 + ,0 + ,1 + ,19370 + ,10 + ,62 + ,0 + ,1 + ,26460 + ,10 + ,60 + ,0 + ,1 + ,29570 + ,11 + ,65 + ,0 + ,1 + ,40070 + ,11 + ,67 + ,0 + ,1 + ,23860 + ,11 + ,62 + ,0 + ,0 + ,32510 + ,10 + ,66 + ,0 + ,1 + ,27620 + ,11 + ,60 + ,0 + ,0 + ,30110 + ,11 + ,64 + ,0 + ,0 + ,43050 + ,13 + ,69 + ,0 + ,1 + ,39060 + ,13 + ,67 + ,0 + ,1 + ,35830 + ,11 + ,67 + ,0 + ,1 + ,32360 + ,11 + ,66 + ,0 + ,0 + ,34360 + ,14 + ,63 + ,0 + ,1 + ,30580 + ,11 + ,61 + ,0 + ,1 + ,30070 + ,10 + ,62 + ,0 + ,1 + ,34890 + ,10 + ,67 + ,0 + ,1 + ,28640 + ,10 + ,60 + ,0 + ,0 + ,34280 + ,14 + ,64 + ,1 + ,0 + ,28190 + ,13 + ,62 + ,0 + ,0 + ,22500 + ,10 + ,58 + ,0 + ,0 + ,46830 + ,14 + ,69 + ,0 + ,1 + ,23520 + ,10 + ,62 + ,0 + ,1 + ,31080 + ,11 + ,65 + ,0 + ,1 + ,39940 + ,13 + ,67 + ,0 + ,1 + ,43930 + ,12 + ,69 + ,0 + ,1 + ,32080 + ,13 + ,61 + ,1 + ,0 + ,25920 + ,10 + ,65 + ,0 + ,1 + ,31930 + ,13 + ,70 + ,0 + ,1 + ,16940 + ,11 + ,60 + ,1 + ,1 + ,39570 + ,14 + ,72 + ,1 + ,1 + ,23460 + ,11 + ,59 + ,0 + ,0 + ,47890 + ,13 + ,69 + ,1 + ,1 + ,35150 + ,11 + ,68 + ,0 + ,1 + ,27540 + ,11 + ,66 + ,0 + ,0 + ,27200 + ,10 + ,66 + ,0 + ,1 + ,24630 + ,11 + ,65 + ,0 + ,1 + ,26330 + ,11 + ,62 + ,0 + ,0 + ,30480 + ,10 + ,66 + ,0 + ,0 + ,31110 + ,11 + ,68 + ,0 + ,1 + ,37450 + ,13 + ,68 + ,0 + ,0 + ,23840 + ,12 + ,64 + ,1 + ,0 + ,20940 + ,10 + ,59 + ,0 + ,1 + ,31830 + ,10 + ,66 + ,0 + ,0 + ,30740 + ,14 + ,65 + ,1 + ,0 + ,39770 + ,11 + ,71 + ,0 + ,1 + ,33540 + ,10 + ,63 + ,0 + ,1 + ,34110 + ,11 + ,64 + ,0 + ,0 + ,23870 + ,10 + ,66 + ,1 + ,0 + ,31710 + ,11 + ,63 + ,0 + ,0 + ,38870 + ,13 + ,68 + ,0 + ,1 + ,26460 + ,13 + ,62 + ,0 + ,0 + ,25040 + ,10 + ,60 + ,0 + ,0 + ,35870 + ,11 + ,65 + ,0 + ,1 + ,38450 + ,11 + ,69 + ,0 + ,1 + ,29710 + ,12 + ,65 + ,0 + ,1 + ,28910 + ,10 + ,61 + ,0 + ,0 + ,18230 + ,10 + ,57 + ,0 + ,0 + ,24170 + ,11 + ,63 + ,0 + ,1 + ,21750 + ,10 + ,58 + ,0 + ,0 + ,27350 + ,11 + ,63 + ,0 + ,0 + ,42730 + ,14 + ,73 + ,0 + ,1 + ,29760 + ,13 + ,66 + ,0 + ,1 + ,38350 + ,12 + ,70 + ,1 + ,0 + ,40650 + ,11 + ,67 + ,0 + ,1 + ,23180 + ,11 + ,59 + ,0 + ,0 + ,35960 + ,11 + ,68 + ,0 + ,1 + ,33950 + ,14 + ,67 + ,0 + ,0 + ,27510 + ,12 + ,63 + ,0 + ,0 + ,26730 + ,10 + ,65 + ,0 + ,0 + ,25560 + ,12 + ,62 + ,0 + ,0 + ,25420 + ,11 + ,62 + ,0 + ,0 + ,26080 + ,10 + ,66 + ,0 + ,1 + ,23540 + ,11 + ,62 + ,0 + ,0 + ,25990 + ,13 + ,63 + ,1 + ,0 + ,14580 + ,10 + ,57 + ,0 + ,0 + ,37950 + ,10 + ,69 + ,0 + ,1 + ,24910 + ,11 + ,59 + ,0 + ,0 + ,30600 + ,13 + ,62 + ,0 + ,0 + ,25450 + ,10 + ,65 + ,0 + ,1 + ,29930 + ,11 + ,67 + ,0 + ,1 + ,33050 + ,10 + ,65 + ,0 + ,0 + ,47560 + ,13 + ,68 + ,1 + ,1 + ,37740 + ,11 + ,67 + ,0 + ,0 + ,28550 + ,10 + ,65 + ,0 + ,1 + ,29880 + ,11 + ,70 + ,0 + ,1 + ,24980 + ,11 + ,60 + ,0 + ,1 + ,31690 + ,14 + ,64 + ,0 + ,0 + ,28870 + ,11 + ,63 + ,0 + ,1 + ,27040 + ,13 + ,61 + ,0 + ,0 + ,35150 + ,11 + ,64 + ,0 + ,0 + ,34250 + ,11 + ,66 + ,0 + ,1 + ,22870 + ,10 + ,61 + ,0 + ,0 + ,24340 + ,13 + ,65 + ,0 + ,0 + ,23650 + ,10 + ,64 + ,0 + ,0 + ,30860 + ,13 + ,68 + ,1 + ,0 + ,26960 + ,10 + ,66 + ,0 + ,1 + ,28680 + ,12 + ,62 + ,0 + ,0 + ,28130 + ,10 + ,62 + ,0 + ,0 + ,43090 + ,14 + ,69 + ,1 + ,1 + ,32550 + ,12 + ,66 + ,0 + ,0 + ,34130 + ,10 + ,66 + ,1 + ,0 + ,45930 + ,11 + ,69 + ,0 + ,1 + ,41110 + ,14 + ,71 + ,0 + ,1 + ,19160 + ,12 + ,61 + ,0 + ,1 + ,18580 + ,10 + ,58 + ,0 + ,1 + ,29750 + ,10 + ,63 + ,1 + ,0 + ,33500 + ,10 + ,69 + ,0 + ,1 + ,29010 + ,10 + ,60 + ,0 + ,1 + ,22410 + ,12 + ,64 + ,0 + ,1 + ,42250 + ,13 + ,74 + ,0 + ,1 + ,32230 + ,11 + ,65 + ,0 + ,0 + ,52240 + ,12 + ,70 + ,0 + ,1 + ,40730 + ,11 + ,67 + ,0 + ,1 + ,40800 + ,12 + ,65 + ,0 + ,1 + ,26060 + ,11 + ,65 + ,0 + ,0 + ,31690 + ,11 + ,63 + ,1 + ,0 + ,44110 + ,12 + ,68 + ,0 + ,1 + ,37910 + ,12 + ,69 + ,0 + ,1 + ,30890 + ,13 + ,68 + ,0 + ,1 + ,24650 + ,11 + ,60 + ,0 + ,1 + ,33430 + ,12 + ,68 + ,1 + ,1 + ,32000 + ,10 + ,65 + ,0 + ,1 + ,29130 + ,12 + ,64 + ,0 + ,1 + ,48770 + ,13 + ,73 + ,0 + ,1 + ,23580 + ,10 + ,59 + ,0 + ,0 + ,32790 + ,12 + ,71 + ,0 + ,1 + ,25810 + ,10 + ,66 + ,0 + ,1 + ,23470 + ,12 + ,62 + ,0 + ,0 + ,26910 + ,10 + ,67 + ,0 + ,0 + ,28270 + ,11 + ,63 + ,0 + ,0 + ,18730 + ,10 + ,53 + ,0 + ,1 + ,37510 + ,12 + ,72 + ,1 + ,1 + ,25380 + ,14 + ,71 + ,0 + ,0 + ,27580 + ,10 + ,66 + ,0 + ,1 + ,30500 + ,10 + ,60 + ,0 + ,0 + ,30790 + ,12 + ,60 + ,0 + ,0 + ,22010 + ,10 + ,61 + ,0 + ,1 + ,18580 + ,10 + ,59 + ,0 + ,1 + ,22160 + ,13 + ,68 + ,1 + ,0 + ,34030 + ,12 + ,62 + ,0 + ,0 + ,35010 + ,12 + ,65 + ,0 + ,0 + ,25780 + ,11 + ,63 + ,0 + ,0 + ,30780 + ,13 + ,66 + ,1 + ,0 + ,31860 + ,12 + ,67 + ,1 + ,0 + ,16650 + ,10 + ,57 + ,0 + ,1 + ,20810 + ,11 + ,63 + ,0 + ,0 + ,29740 + ,11 + ,62 + ,0 + ,0 + ,32970 + ,13 + ,65 + ,1 + ,0 + ,40730 + ,12 + ,69 + ,0 + ,1 + ,44480 + ,13 + ,69 + ,0 + ,1 + ,39840 + ,13 + ,71 + ,0 + ,1 + ,22500 + ,10 + ,58 + ,0 + ,0 + ,27520 + ,12 + ,64 + ,0 + ,0 + ,23040 + ,12 + ,67 + ,1 + ,1 + ,36800 + ,14 + ,67 + ,0 + ,1 + ,31020 + ,11 + ,64 + ,1 + ,0 + ,28620 + ,10 + ,61 + ,0 + ,0 + ,26770 + ,13 + ,67 + ,1 + ,0 + ,30230 + ,11 + ,68 + ,0 + ,0 + ,36810 + ,11 + ,68 + ,0 + ,0 + ,32550 + ,13 + ,67 + ,0 + ,0 + ,36920 + ,12 + ,67 + ,0 + ,1 + ,23560 + ,10 + ,61 + ,0 + ,0 + ,45910 + ,10 + ,70 + ,0 + ,1 + ,30820 + ,12 + ,64 + ,0 + ,0 + ,32970 + ,13 + ,65 + ,1 + ,0 + ,32580 + ,11 + ,63 + ,0 + ,0 + ,22160 + ,10 + ,61 + ,0 + ,1 + ,32470 + ,11 + ,66 + ,0 + ,1 + ,43240 + ,11 + ,68 + ,0 + ,1 + ,23620 + ,11 + ,61 + ,0 + ,0 + ,25630 + ,11 + ,63 + ,0 + ,0 + ,32060 + ,11 + ,64 + ,0 + ,1 + ,35850 + ,14 + ,70 + ,0 + ,1 + ,47200 + ,12 + ,72 + ,0 + ,1 + ,33310 + ,13 + ,66 + ,0 + ,0 + ,50830 + ,13 + ,74 + ,0 + ,1 + ,34980 + ,10 + ,68 + ,1 + ,1 + ,24170 + ,12 + ,61 + ,0 + ,0 + ,23640 + ,10 + ,61 + ,0 + ,1 + ,23410 + ,10 + ,61 + ,0 + ,1 + ,27590 + ,12 + ,62 + ,1 + ,0 + ,29530 + ,11 + ,67 + ,1 + ,0 + ,32310 + ,12 + ,63 + ,0 + ,1 + ,30780 + ,11 + ,68 + ,0 + ,1 + ,33690 + ,11 + ,71 + ,0 + ,1 + ,35290 + ,12 + ,71 + ,0 + ,1 + ,28660 + ,12 + ,62 + ,0 + ,0 + ,28910 + ,14 + ,62 + ,0 + ,0 + ,30220 + ,11 + ,62 + ,0 + ,0 + ,31270 + ,10 + ,62 + ,0 + ,1 + ,28660 + ,11 + ,61 + ,0 + ,0 + ,26050 + ,12 + ,63 + ,0 + ,0 + ,30560 + ,13 + ,63 + ,0 + ,0 + ,25690 + ,12 + ,63 + ,0 + ,0 + ,25010 + ,11 + ,62 + ,0 + ,0 + ,33200 + ,11 + ,66 + ,0 + ,1 + ,21230 + ,11 + ,65 + ,0 + ,1 + ,37800 + ,14 + ,70 + ,0 + ,1 + ,38470 + ,11 + ,66 + ,0 + ,1 + ,37850 + ,13 + ,63 + ,1 + ,0 + ,39240 + ,12 + ,68 + ,0 + ,1 + ,21320 + ,10 + ,59 + ,0 + ,1 + ,27520 + ,12 + ,69 + ,0 + ,1 + ,24490 + ,13 + ,63 + ,0 + ,0 + ,34560 + ,10 + ,63 + ,0 + ,1 + ,30730 + ,10 + ,66 + ,0 + ,0 + ,26880 + ,10 + ,62 + ,0 + ,0 + ,33290 + ,10 + ,68 + ,0 + ,1 + ,42710 + ,14 + ,73 + ,0 + ,1 + ,35300 + ,12 + ,64 + ,0 + ,1 + ,29280 + ,11 + ,66 + ,0 + ,1 + ,26890 + ,11 + ,62 + ,0 + ,0 + ,23320 + ,12 + ,57 + ,0 + ,1 + ,29340 + ,14 + ,64 + ,0 + ,0 + ,22760 + ,14 + ,66 + ,1 + ,1 + ,31100 + ,10 + ,65 + ,0 + ,1 + ,28940 + ,11 + ,67 + ,0 + ,1 + ,46370 + ,11 + ,72 + ,1 + ,1 + ,24350 + ,10 + ,65 + ,0 + ,0 + ,28380 + ,10 + ,63 + ,0 + ,0 + ,30350 + ,12 + ,62 + ,0 + ,0 + ,48310 + ,12 + ,71 + ,0 + ,1 + ,28120 + ,11 + ,61 + ,0 + ,1 + ,27140 + ,12 + ,66 + ,0 + ,0 + ,30860 + ,10 + ,62 + ,0 + ,0 + ,35190 + ,12 + ,66 + ,0 + ,0 + ,42320 + ,13 + ,71 + ,0 + ,1 + ,27700 + ,10 + ,62 + ,0 + ,1 + ,33410 + ,12 + ,66 + ,0 + ,0 + ,30900 + ,10 + ,65 + ,0 + ,1 + ,25310 + ,13 + ,61 + ,0 + ,1 + ,28220 + ,12 + ,70 + ,0 + ,1 + ,30380 + ,10 + ,65 + ,1 + ,0 + ,29350 + ,12 + ,66 + ,0 + ,1 + ,25680 + ,10 + ,64 + ,0 + ,0 + ,23870 + ,11 + ,61 + ,0 + ,1 + ,24990 + ,12 + ,65 + ,0 + ,1 + ,41300 + ,11 + ,67 + ,0 + ,1 + ,30010 + ,12 + ,64 + ,0 + ,0 + ,31320 + ,10 + ,60 + ,0 + ,0 + ,35770 + ,13 + ,64 + ,0 + ,0 + ,32220 + ,12 + ,61 + ,0 + ,0 + ,32800 + ,11 + ,66 + ,0 + ,1 + ,26590 + ,11 + ,64 + ,0 + ,1 + ,28220 + ,11 + ,62 + ,0 + ,0 + ,21400 + ,11 + ,61 + ,0 + ,0 + ,42030 + ,12 + ,71 + ,0 + ,1 + ,29970 + ,14 + ,65 + ,0 + ,0 + ,31200 + ,11 + ,61 + ,1 + ,0 + ,25620 + ,11 + ,63 + ,0 + ,0 + ,30820 + ,12 + ,65 + ,0 + ,0 + ,38060 + ,14 + ,68 + ,0 + ,1 + ,33390 + ,11 + ,69 + ,1 + ,1 + ,31520 + ,13 + ,62 + ,1 + ,0 + ,24580 + ,11 + ,60 + ,0 + ,0 + ,23910 + ,10 + ,60 + ,0 + ,1 + ,31410 + ,13 + ,61 + ,0 + ,0 + ,25790 + ,12 + ,63 + ,0 + ,0 + ,31040 + ,11 + ,68 + ,1 + ,0 + ,40450 + ,13 + ,69 + ,1 + ,1 + ,47630 + ,14 + ,68 + ,1 + ,1 + ,21000 + ,10 + ,58 + ,0 + ,1 + ,30690 + ,11 + ,65 + ,1 + ,0 + ,27850 + ,11 + ,69 + ,0 + ,1 + ,42840 + ,15 + ,70 + ,0 + ,1 + ,45060 + ,15 + ,71 + ,1 + ,1 + ,29060 + ,18 + ,66 + ,0 + ,0 + ,51020 + ,19 + ,72 + ,0 + ,1 + ,35190 + ,19 + ,66 + ,1 + ,0 + ,36880 + ,16 + ,68 + ,1 + ,1 + ,44290 + ,17 + ,70 + ,0 + ,1 + ,42790 + ,15 + ,68 + ,0 + ,1 + ,45000 + ,15 + ,70 + ,0 + ,1 + ,26350 + ,15 + ,64 + ,0 + ,0 + ,26790 + ,15 + ,66 + ,1 + ,0 + ,21980 + ,15 + ,62 + ,1 + ,0 + ,33450 + ,19 + ,66 + ,1 + ,0 + ,30820 + ,18 + ,65 + ,0 + ,0 + ,33870 + ,16 + ,67 + ,0 + ,0 + ,30820 + ,17 + ,67 + ,1 + ,1 + ,29030 + ,16 + ,63 + ,1 + ,0 + ,30040 + ,15 + ,64 + ,1 + ,0 + ,57930 + ,15 + ,69 + ,0 + ,1 + ,39850 + ,15 + ,71 + ,0 + ,1 + ,42200 + ,18 + ,68 + ,0 + ,1 + ,47240 + ,17 + ,71 + ,0 + ,1 + ,37310 + ,15 + ,67 + ,0 + ,1 + ,34060 + ,17 + ,69 + ,1 + ,1 + ,35000 + ,17 + ,62 + ,0 + ,0 + ,36740 + ,16 + ,68 + ,0 + ,0 + ,56330 + ,17 + ,73 + ,0 + ,1 + ,31220 + ,15 + ,64 + ,1 + ,0 + ,33300 + ,15 + ,69 + ,1 + ,0 + ,26080 + ,16 + ,62 + ,1 + ,0 + ,36450 + ,16 + ,74 + ,0 + ,1 + ,37990 + ,15 + ,67 + ,1 + ,1 + ,40860 + ,18 + ,67 + ,1 + ,1 + ,28870 + ,15 + ,63 + ,0 + ,0 + ,40700 + ,16 + ,70 + ,1 + ,1 + ,39600 + ,17 + ,70 + ,0 + ,1 + ,42990 + ,16 + ,66 + ,0 + ,1 + ,29810 + ,16 + ,66 + ,0 + ,0 + ,22640 + ,15 + ,63 + ,1 + ,0 + ,44040 + ,18 + ,71 + ,1 + ,1 + ,22780 + ,15 + ,60 + ,1 + ,0 + ,45040 + ,16 + ,72 + ,0 + ,1 + ,56380 + ,17 + ,70 + ,0 + ,1 + ,48720 + ,16 + ,72 + ,1 + ,1 + ,42700 + ,16 + ,67 + ,1 + ,1 + ,37270 + ,15 + ,68 + ,1 + ,1 + ,28530 + ,18 + ,60 + ,0 + ,0 + ,27950 + ,16 + ,63 + ,1 + ,0 + ,32110 + ,15 + ,67 + ,0 + ,0) + ,dim=c(5 + ,651) + ,dimnames=list(c('FEV' + ,'Age' + ,'Ht' + ,'Smoker' + ,'gender') + ,1:651)) > y <- array(NA,dim=c(5,651),dimnames=list(c('FEV','Age','Ht','Smoker','gender'),1:651)) > 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 FEV Age Ht Smoker gender 1 17080 9 57 0 0 2 17240 8 68 0 0 3 17200 7 55 0 0 4 15580 9 53 0 1 5 18950 9 57 0 1 6 23360 8 61 0 0 7 19190 6 58 0 0 8 14150 6 56 0 0 9 19870 8 59 0 0 10 19420 9 60 0 0 11 16020 6 53 0 0 12 17350 8 54 0 1 13 21930 8 59 0 0 14 21180 8 61 0 1 15 22580 8 58 0 1 16 19320 7 53 0 1 17 14720 5 50 0 1 18 18780 6 53 0 0 19 23520 9 59 0 1 20 26040 9 62 0 1 21 14000 5 49 0 0 22 12560 5 53 0 0 23 25780 7 63 0 1 24 29880 9 65 0 0 25 14040 3 52 0 1 26 23480 9 60 0 1 27 17550 5 52 0 1 28 29800 8 60 0 0 29 21000 9 60 0 0 30 12820 5 49 0 0 31 30000 9 66 0 1 32 26730 8 60 0 0 33 20930 7 58 0 0 34 16120 5 52 0 0 35 21750 8 59 0 0 36 27250 9 59 0 1 37 20710 8 55 0 1 38 15470 8 57 0 1 39 20040 8 57 0 1 40 31350 9 60 0 0 41 24200 8 59 0 1 42 17760 5 51 0 1 43 19310 8 57 0 0 44 13430 5 50 0 0 45 20760 9 57 0 0 46 16240 7 54 0 1 47 13440 8 53 0 0 48 16500 6 55 0 1 49 27320 8 61 0 1 50 20170 5 55 0 1 51 27970 9 62 0 0 52 35560 9 62 0 1 53 17030 8 55 0 1 54 16340 6 54 0 1 55 25700 9 57 0 1 56 30160 9 63 0 0 57 24190 7 60 0 0 58 15690 4 50 0 0 59 16980 8 58 0 0 60 21230 8 60 0 1 61 24810 8 60 0 0 62 14810 6 51 0 0 63 15770 4 49 0 0 64 19400 8 59 0 1 65 17470 6 58 0 1 66 20690 9 58 0 1 67 16310 7 56 0 0 68 15360 5 52 0 0 69 25600 9 61 0 0 70 19620 8 57 0 1 71 25310 8 58 0 0 72 27150 9 60 0 1 73 24570 9 59 0 1 74 20900 9 60 0 1 75 17890 7 56 0 1 76 18580 5 53 0 1 77 14520 5 51 0 1 78 38420 9 69 0 1 79 17190 6 53 0 0 80 21110 7 57 0 0 81 16950 6 53 0 0 82 22110 8 63 0 1 83 17940 8 55 0 1 84 19170 7 58 0 0 85 21440 8 63 0 0 86 12530 7 52 0 1 87 26590 9 62 0 1 88 15800 5 53 0 1 89 21260 9 62 0 1 90 30290 9 62 0 0 91 29640 9 65 0 1 92 16110 7 58 0 1 93 22150 8 60 0 0 94 23880 8 60 0 0 95 21960 9 61 0 1 96 17510 9 58 0 1 97 21650 9 62 0 1 98 16820 7 55 0 1 99 15230 8 55 0 1 100 12920 8 52 0 0 101 16490 7 54 0 1 102 25880 9 63 0 1 103 25740 9 61 0 0 104 19790 6 56 0 1 105 23540 8 59 0 1 106 17180 6 55 0 1 107 17420 7 59 0 0 108 16030 7 51 0 0 109 26390 8 60 0 0 110 18290 7 54 0 0 111 20840 7 58 0 1 112 22200 7 58 0 1 113 14730 7 53 0 0 114 23410 8 61 0 0 115 16980 7 55 0 0 116 11960 5 47 0 0 117 18720 8 57 0 0 118 22190 7 55 0 1 119 24200 9 57 0 1 120 18270 7 55 0 0 121 14610 7 54 0 0 122 13380 6 53 0 1 123 20900 8 57 0 1 124 16970 8 59 0 0 125 15620 8 55 0 1 126 20400 9 56 0 0 127 16090 7 52 0 0 128 24580 8 61 0 0 129 26500 9 64 0 1 130 14290 8 58 0 1 131 16750 8 53 0 1 132 19470 9 57 0 0 133 20690 8 54 0 1 134 15720 6 52 0 1 135 13480 6 53 0 1 136 22880 8 62 0 0 137 17730 9 59 0 1 138 19050 7 58 0 1 139 24630 9 61 0 0 140 14310 6 51 0 1 141 26310 9 62 0 0 142 31140 9 65 0 1 143 21350 9 59 0 1 144 15270 6 53 0 1 145 22930 8 58 0 0 146 30420 9 66 0 0 147 29270 8 64 0 1 148 26650 8 64 0 0 149 23010 9 59 0 1 150 24600 9 64 0 1 151 25920 9 61 0 0 152 17500 7 55 0 0 153 17590 8 53 0 1 154 15360 6 48 0 1 155 22590 9 59 0 0 156 20480 9 65 0 0 157 25710 9 61 0 1 158 20460 7 56 0 1 159 17800 8 59 0 0 160 15520 5 54 0 0 161 19530 8 58 0 0 162 28930 9 65 0 1 163 17130 6 51 0 1 164 28510 9 60 0 0 165 16240 6 52 0 1 166 26310 8 59 0 1 167 18190 5 53 0 1 168 16580 7 53 0 1 169 21580 7 54 0 1 170 17890 4 52 0 1 171 30040 9 64 0 0 172 25030 8 63 0 1 173 19330 9 58 0 0 174 20910 9 59 0 0 175 23160 9 60 0 0 176 17040 5 51 0 0 177 16060 9 58 0 0 178 11650 7 47 0 1 179 21020 6 56 0 0 180 23200 9 57 0 0 181 22300 9 61 0 1 182 17160 9 56 0 1 183 17900 7 54 0 1 184 11460 5 50 0 0 185 21870 8 62 0 0 186 27170 9 62 0 1 187 17960 7 55 0 1 188 19530 9 58 1 1 189 13350 8 57 0 0 190 21190 9 57 0 1 191 16660 6 52 0 1 192 18260 6 53 0 1 193 27090 8 63 0 0 194 28710 9 65 0 1 195 10920 5 50 0 0 196 22620 6 58 0 1 197 21040 6 57 0 1 198 21660 9 58 0 0 199 16900 7 54 0 0 200 29730 9 60 0 1 201 21450 8 60 0 0 202 19710 5 58 0 1 203 20950 7 57 0 0 204 16970 6 55 0 0 205 24550 9 60 0 0 206 19200 7 57 0 1 207 21640 9 60 0 1 208 21300 9 59 0 0 209 29930 8 63 0 0 210 25290 9 59 0 0 211 17260 7 53 0 0 212 24420 9 62 0 0 213 11020 4 48 0 0 214 20560 9 63 0 0 215 18080 5 56 0 1 216 23050 8 65 0 0 217 19690 9 59 0 0 218 15560 8 59 0 0 219 10720 3 46 0 0 220 20420 9 62 0 1 221 15120 8 53 0 0 222 14230 6 50 0 1 223 36810 9 68 0 1 224 19910 8 60 0 1 225 18970 8 56 0 1 226 13700 7 55 0 0 227 13380 6 52 0 0 228 20160 8 56 0 1 229 26390 9 63 0 0 230 13890 4 48 0 0 231 16120 7 57 0 1 232 21350 8 59 0 0 233 26810 8 61 0 1 234 32230 9 65 0 0 235 17960 6 55 0 0 236 20100 8 55 0 1 237 15230 6 51 0 0 238 17440 8 53 0 1 239 24850 9 64 0 0 240 23350 8 59 0 0 241 14150 7 54 0 0 242 20760 9 61 0 1 243 24350 8 60 0 1 244 17280 7 57 0 0 245 28500 9 63 0 0 246 18440 8 57 0 0 247 17540 9 62 0 0 248 13430 6 52 0 0 249 23030 8 57 0 1 250 22460 9 64 0 1 251 24760 8 63 0 0 252 32390 9 65 0 1 253 24570 9 62 0 1 254 23820 8 62 0 0 255 16400 7 55 0 0 256 15890 5 51 0 0 257 20560 7 54 0 1 258 22260 8 57 0 1 259 18860 9 56 0 0 260 28330 9 62 0 1 261 17150 6 53 0 1 262 26310 8 59 0 1 263 25500 7 56 0 1 264 19120 9 59 0 0 265 18770 7 53 0 0 266 19350 7 53 0 0 267 15390 5 50 0 0 268 28030 9 60 0 1 269 29230 9 64 0 1 270 23580 8 61 0 0 271 20940 8 58 0 1 272 18550 9 60 0 1 273 15350 6 55 0 0 274 21350 7 56 0 1 275 19300 5 51 0 1 276 21820 9 60 0 0 277 13590 5 51 0 1 278 20020 7 58 0 0 279 16990 6 54 0 1 280 25000 8 57 0 1 281 23660 7 58 0 0 282 20690 8 60 0 0 283 14180 4 49 0 0 284 23330 8 57 0 0 285 15140 5 52 0 1 286 17580 8 52 0 0 287 25350 7 60 0 1 288 25640 7 58 0 0 289 24870 9 64 0 0 290 15910 9 57 0 0 291 16240 8 53 0 1 292 27980 9 62 0 1 293 16910 6 53 0 1 294 19990 8 57 0 0 295 18690 9 57 0 1 296 10040 4 48 0 1 297 14270 6 50 0 1 298 18260 7 51 0 1 299 26880 9 60 0 0 300 16570 8 56 0 1 301 16720 6 54 0 0 302 20150 8 58 0 0 303 23710 7 56 0 0 304 21150 5 50 0 1 305 23280 8 60 0 0 306 14950 7 57 0 0 307 28840 11 69 0 1 308 23280 10 64 0 1 309 33810 14 63 0 1 310 21700 11 58 0 0 311 34700 11 67 0 1 312 30580 12 61 0 0 313 18110 10 57 0 1 314 25240 11 64 0 1 315 26420 10 61 0 0 316 37410 14 69 0 1 317 43360 13 70 0 1 318 48420 14 72 0 1 319 45500 12 71 0 1 320 28410 12 63 0 0 321 31660 10 62 0 0 322 38160 13 64 0 0 323 25610 10 62 0 1 324 36540 11 65 0 0 325 24810 10 61 0 1 326 26650 11 63 0 0 327 32030 10 66 0 1 328 35490 13 68 0 1 329 22360 14 66 1 0 330 32220 11 72 0 1 331 31110 10 66 0 1 332 34900 11 67 0 0 333 31470 13 64 0 0 334 25200 10 61 0 0 335 22920 10 63 0 1 336 28890 12 64 0 0 337 22460 10 61 0 1 338 19370 10 62 0 1 339 26460 10 60 0 1 340 29570 11 65 0 1 341 40070 11 67 0 1 342 23860 11 62 0 0 343 32510 10 66 0 1 344 27620 11 60 0 0 345 30110 11 64 0 0 346 43050 13 69 0 1 347 39060 13 67 0 1 348 35830 11 67 0 1 349 32360 11 66 0 0 350 34360 14 63 0 1 351 30580 11 61 0 1 352 30070 10 62 0 1 353 34890 10 67 0 1 354 28640 10 60 0 0 355 34280 14 64 1 0 356 28190 13 62 0 0 357 22500 10 58 0 0 358 46830 14 69 0 1 359 23520 10 62 0 1 360 31080 11 65 0 1 361 39940 13 67 0 1 362 43930 12 69 0 1 363 32080 13 61 1 0 364 25920 10 65 0 1 365 31930 13 70 0 1 366 16940 11 60 1 1 367 39570 14 72 1 1 368 23460 11 59 0 0 369 47890 13 69 1 1 370 35150 11 68 0 1 371 27540 11 66 0 0 372 27200 10 66 0 1 373 24630 11 65 0 1 374 26330 11 62 0 0 375 30480 10 66 0 0 376 31110 11 68 0 1 377 37450 13 68 0 0 378 23840 12 64 1 0 379 20940 10 59 0 1 380 31830 10 66 0 0 381 30740 14 65 1 0 382 39770 11 71 0 1 383 33540 10 63 0 1 384 34110 11 64 0 0 385 23870 10 66 1 0 386 31710 11 63 0 0 387 38870 13 68 0 1 388 26460 13 62 0 0 389 25040 10 60 0 0 390 35870 11 65 0 1 391 38450 11 69 0 1 392 29710 12 65 0 1 393 28910 10 61 0 0 394 18230 10 57 0 0 395 24170 11 63 0 1 396 21750 10 58 0 0 397 27350 11 63 0 0 398 42730 14 73 0 1 399 29760 13 66 0 1 400 38350 12 70 1 0 401 40650 11 67 0 1 402 23180 11 59 0 0 403 35960 11 68 0 1 404 33950 14 67 0 0 405 27510 12 63 0 0 406 26730 10 65 0 0 407 25560 12 62 0 0 408 25420 11 62 0 0 409 26080 10 66 0 1 410 23540 11 62 0 0 411 25990 13 63 1 0 412 14580 10 57 0 0 413 37950 10 69 0 1 414 24910 11 59 0 0 415 30600 13 62 0 0 416 25450 10 65 0 1 417 29930 11 67 0 1 418 33050 10 65 0 0 419 47560 13 68 1 1 420 37740 11 67 0 0 421 28550 10 65 0 1 422 29880 11 70 0 1 423 24980 11 60 0 1 424 31690 14 64 0 0 425 28870 11 63 0 1 426 27040 13 61 0 0 427 35150 11 64 0 0 428 34250 11 66 0 1 429 22870 10 61 0 0 430 24340 13 65 0 0 431 23650 10 64 0 0 432 30860 13 68 1 0 433 26960 10 66 0 1 434 28680 12 62 0 0 435 28130 10 62 0 0 436 43090 14 69 1 1 437 32550 12 66 0 0 438 34130 10 66 1 0 439 45930 11 69 0 1 440 41110 14 71 0 1 441 19160 12 61 0 1 442 18580 10 58 0 1 443 29750 10 63 1 0 444 33500 10 69 0 1 445 29010 10 60 0 1 446 22410 12 64 0 1 447 42250 13 74 0 1 448 32230 11 65 0 0 449 52240 12 70 0 1 450 40730 11 67 0 1 451 40800 12 65 0 1 452 26060 11 65 0 0 453 31690 11 63 1 0 454 44110 12 68 0 1 455 37910 12 69 0 1 456 30890 13 68 0 1 457 24650 11 60 0 1 458 33430 12 68 1 1 459 32000 10 65 0 1 460 29130 12 64 0 1 461 48770 13 73 0 1 462 23580 10 59 0 0 463 32790 12 71 0 1 464 25810 10 66 0 1 465 23470 12 62 0 0 466 26910 10 67 0 0 467 28270 11 63 0 0 468 18730 10 53 0 1 469 37510 12 72 1 1 470 25380 14 71 0 0 471 27580 10 66 0 1 472 30500 10 60 0 0 473 30790 12 60 0 0 474 22010 10 61 0 1 475 18580 10 59 0 1 476 22160 13 68 1 0 477 34030 12 62 0 0 478 35010 12 65 0 0 479 25780 11 63 0 0 480 30780 13 66 1 0 481 31860 12 67 1 0 482 16650 10 57 0 1 483 20810 11 63 0 0 484 29740 11 62 0 0 485 32970 13 65 1 0 486 40730 12 69 0 1 487 44480 13 69 0 1 488 39840 13 71 0 1 489 22500 10 58 0 0 490 27520 12 64 0 0 491 23040 12 67 1 1 492 36800 14 67 0 1 493 31020 11 64 1 0 494 28620 10 61 0 0 495 26770 13 67 1 0 496 30230 11 68 0 0 497 36810 11 68 0 0 498 32550 13 67 0 0 499 36920 12 67 0 1 500 23560 10 61 0 0 501 45910 10 70 0 1 502 30820 12 64 0 0 503 32970 13 65 1 0 504 32580 11 63 0 0 505 22160 10 61 0 1 506 32470 11 66 0 1 507 43240 11 68 0 1 508 23620 11 61 0 0 509 25630 11 63 0 0 510 32060 11 64 0 1 511 35850 14 70 0 1 512 47200 12 72 0 1 513 33310 13 66 0 0 514 50830 13 74 0 1 515 34980 10 68 1 1 516 24170 12 61 0 0 517 23640 10 61 0 1 518 23410 10 61 0 1 519 27590 12 62 1 0 520 29530 11 67 1 0 521 32310 12 63 0 1 522 30780 11 68 0 1 523 33690 11 71 0 1 524 35290 12 71 0 1 525 28660 12 62 0 0 526 28910 14 62 0 0 527 30220 11 62 0 0 528 31270 10 62 0 1 529 28660 11 61 0 0 530 26050 12 63 0 0 531 30560 13 63 0 0 532 25690 12 63 0 0 533 25010 11 62 0 0 534 33200 11 66 0 1 535 21230 11 65 0 1 536 37800 14 70 0 1 537 38470 11 66 0 1 538 37850 13 63 1 0 539 39240 12 68 0 1 540 21320 10 59 0 1 541 27520 12 69 0 1 542 24490 13 63 0 0 543 34560 10 63 0 1 544 30730 10 66 0 0 545 26880 10 62 0 0 546 33290 10 68 0 1 547 42710 14 73 0 1 548 35300 12 64 0 1 549 29280 11 66 0 1 550 26890 11 62 0 0 551 23320 12 57 0 1 552 29340 14 64 0 0 553 22760 14 66 1 1 554 31100 10 65 0 1 555 28940 11 67 0 1 556 46370 11 72 1 1 557 24350 10 65 0 0 558 28380 10 63 0 0 559 30350 12 62 0 0 560 48310 12 71 0 1 561 28120 11 61 0 1 562 27140 12 66 0 0 563 30860 10 62 0 0 564 35190 12 66 0 0 565 42320 13 71 0 1 566 27700 10 62 0 1 567 33410 12 66 0 0 568 30900 10 65 0 1 569 25310 13 61 0 1 570 28220 12 70 0 1 571 30380 10 65 1 0 572 29350 12 66 0 1 573 25680 10 64 0 0 574 23870 11 61 0 1 575 24990 12 65 0 1 576 41300 11 67 0 1 577 30010 12 64 0 0 578 31320 10 60 0 0 579 35770 13 64 0 0 580 32220 12 61 0 0 581 32800 11 66 0 1 582 26590 11 64 0 1 583 28220 11 62 0 0 584 21400 11 61 0 0 585 42030 12 71 0 1 586 29970 14 65 0 0 587 31200 11 61 1 0 588 25620 11 63 0 0 589 30820 12 65 0 0 590 38060 14 68 0 1 591 33390 11 69 1 1 592 31520 13 62 1 0 593 24580 11 60 0 0 594 23910 10 60 0 1 595 31410 13 61 0 0 596 25790 12 63 0 0 597 31040 11 68 1 0 598 40450 13 69 1 1 599 47630 14 68 1 1 600 21000 10 58 0 1 601 30690 11 65 1 0 602 27850 11 69 0 1 603 42840 15 70 0 1 604 45060 15 71 1 1 605 29060 18 66 0 0 606 51020 19 72 0 1 607 35190 19 66 1 0 608 36880 16 68 1 1 609 44290 17 70 0 1 610 42790 15 68 0 1 611 45000 15 70 0 1 612 26350 15 64 0 0 613 26790 15 66 1 0 614 21980 15 62 1 0 615 33450 19 66 1 0 616 30820 18 65 0 0 617 33870 16 67 0 0 618 30820 17 67 1 1 619 29030 16 63 1 0 620 30040 15 64 1 0 621 57930 15 69 0 1 622 39850 15 71 0 1 623 42200 18 68 0 1 624 47240 17 71 0 1 625 37310 15 67 0 1 626 34060 17 69 1 1 627 35000 17 62 0 0 628 36740 16 68 0 0 629 56330 17 73 0 1 630 31220 15 64 1 0 631 33300 15 69 1 0 632 26080 16 62 1 0 633 36450 16 74 0 1 634 37990 15 67 1 1 635 40860 18 67 1 1 636 28870 15 63 0 0 637 40700 16 70 1 1 638 39600 17 70 0 1 639 42990 16 66 0 1 640 29810 16 66 0 0 641 22640 15 63 1 0 642 44040 18 71 1 1 643 22780 15 60 1 0 644 45040 16 72 0 1 645 56380 17 70 0 1 646 48720 16 72 1 1 647 42700 16 67 1 1 648 37270 15 68 1 1 649 28530 18 60 0 0 650 27950 16 63 1 0 651 32110 15 67 0 0 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Age Ht Smoker gender -43911.5 672.5 1025.4 -784.9 1578.0 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -13955.8 -2552.9 -37.2 2628.9 19423.3 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -43911.52 2244.90 -19.561 < 2e-16 *** Age 672.49 95.15 7.067 4.11e-12 *** Ht 1025.40 47.60 21.541 < 2e-16 *** Smoker -784.92 596.54 -1.316 0.189 gender 1577.99 334.85 4.713 3.00e-06 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 4147 on 646 degrees of freedom Multiple R-squared: 0.769, Adjusted R-squared: 0.7675 F-statistic: 537.5 on 4 and 646 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,] 5.563147e-01 8.873706e-01 0.4436853 [2,] 4.065641e-01 8.131281e-01 0.5934359 [3,] 2.682019e-01 5.364039e-01 0.7317981 [4,] 1.663742e-01 3.327484e-01 0.8336258 [5,] 9.947359e-02 1.989472e-01 0.9005264 [6,] 9.330633e-02 1.866127e-01 0.9066937 [7,] 7.422211e-02 1.484442e-01 0.9257779 [8,] 7.859224e-02 1.571845e-01 0.9214078 [9,] 5.259950e-02 1.051990e-01 0.9474005 [10,] 3.485864e-02 6.971728e-02 0.9651414 [11,] 2.675245e-02 5.350489e-02 0.9732476 [12,] 2.588096e-02 5.176192e-02 0.9741190 [13,] 3.194714e-02 6.389428e-02 0.9680529 [14,] 1.971675e-02 3.943350e-02 0.9802832 [15,] 1.688055e-02 3.376110e-02 0.9831195 [16,] 1.625636e-02 3.251272e-02 0.9837436 [17,] 5.658476e-02 1.131695e-01 0.9434152 [18,] 4.029448e-02 8.058896e-02 0.9597055 [19,] 2.761421e-02 5.522841e-02 0.9723858 [20,] 1.962470e-02 3.924939e-02 0.9803753 [21,] 1.151633e-01 2.303266e-01 0.8848367 [22,] 8.808537e-02 1.761707e-01 0.9119146 [23,] 6.580483e-02 1.316097e-01 0.9341952 [24,] 6.036546e-02 1.207309e-01 0.9396345 [25,] 8.190856e-02 1.638171e-01 0.9180914 [26,] 6.343755e-02 1.268751e-01 0.9365625 [27,] 4.913859e-02 9.827717e-02 0.9508614 [28,] 3.649001e-02 7.298003e-02 0.9635100 [29,] 4.103357e-02 8.206714e-02 0.9589664 [30,] 3.048772e-02 6.097545e-02 0.9695123 [31,] 4.582099e-02 9.164199e-02 0.9541790 [32,] 3.475852e-02 6.951704e-02 0.9652415 [33,] 1.059798e-01 2.119597e-01 0.8940202 [34,] 8.924417e-02 1.784883e-01 0.9107558 [35,] 7.982307e-02 1.596461e-01 0.9201769 [36,] 6.366212e-02 1.273242e-01 0.9363379 [37,] 4.981680e-02 9.963361e-02 0.9501832 [38,] 3.828965e-02 7.657929e-02 0.9617104 [39,] 3.211118e-02 6.422237e-02 0.9678888 [40,] 3.310712e-02 6.621424e-02 0.9668929 [41,] 2.702816e-02 5.405632e-02 0.9729718 [42,] 2.605016e-02 5.210033e-02 0.9739498 [43,] 2.162438e-02 4.324876e-02 0.9783756 [44,] 2.166229e-02 4.332459e-02 0.9783377 [45,] 8.469398e-02 1.693880e-01 0.9153060 [46,] 7.708426e-02 1.541685e-01 0.9229157 [47,] 6.453189e-02 1.290638e-01 0.9354681 [48,] 6.309809e-02 1.261962e-01 0.9369019 [49,] 6.812201e-02 1.362440e-01 0.9318780 [50,] 5.842251e-02 1.168450e-01 0.9415775 [51,] 5.597676e-02 1.119535e-01 0.9440232 [52,] 5.804217e-02 1.160843e-01 0.9419578 [53,] 5.203858e-02 1.040772e-01 0.9479614 [54,] 4.456611e-02 8.913221e-02 0.9554339 [55,] 3.590124e-02 7.180247e-02 0.9640988 [56,] 3.680726e-02 7.361452e-02 0.9631927 [57,] 3.567552e-02 7.135104e-02 0.9643245 [58,] 3.452515e-02 6.905030e-02 0.9654748 [59,] 2.896235e-02 5.792469e-02 0.9710377 [60,] 2.598190e-02 5.196380e-02 0.9740181 [61,] 2.074673e-02 4.149346e-02 0.9792533 [62,] 1.693913e-02 3.387826e-02 0.9830609 [63,] 1.373945e-02 2.747891e-02 0.9862605 [64,] 1.430749e-02 2.861498e-02 0.9856925 [65,] 1.314912e-02 2.629824e-02 0.9868509 [66,] 1.047554e-02 2.095108e-02 0.9895245 [67,] 9.832962e-03 1.966592e-02 0.9901670 [68,] 8.101900e-03 1.620380e-02 0.9918981 [69,] 6.877716e-03 1.375543e-02 0.9931223 [70,] 5.287399e-03 1.057480e-02 0.9947126 [71,] 9.353268e-03 1.870654e-02 0.9906467 [72,] 7.463815e-03 1.492763e-02 0.9925362 [73,] 5.889525e-03 1.177905e-02 0.9941105 [74,] 4.604182e-03 9.208365e-03 0.9953958 [75,] 5.038006e-03 1.007601e-02 0.9949620 [76,] 4.018780e-03 8.037560e-03 0.9959812 [77,] 3.245392e-03 6.490784e-03 0.9967546 [78,] 3.418395e-03 6.836790e-03 0.9965816 [79,] 3.483638e-03 6.967277e-03 0.9965164 [80,] 2.697221e-03 5.394442e-03 0.9973028 [81,] 2.055848e-03 4.111695e-03 0.9979442 [82,] 2.338228e-03 4.676456e-03 0.9976618 [83,] 2.983579e-03 5.967158e-03 0.9970164 [84,] 2.368587e-03 4.737173e-03 0.9976314 [85,] 3.064751e-03 6.129503e-03 0.9969352 [86,] 2.384455e-03 4.768909e-03 0.9976155 [87,] 1.840123e-03 3.680245e-03 0.9981599 [88,] 1.648273e-03 3.296545e-03 0.9983517 [89,] 1.909573e-03 3.819146e-03 0.9980904 [90,] 1.953472e-03 3.906945e-03 0.9980465 [91,] 1.576616e-03 3.153233e-03 0.9984234 [92,] 1.557389e-03 3.114777e-03 0.9984426 [93,] 1.409342e-03 2.818684e-03 0.9985907 [94,] 1.091237e-03 2.182475e-03 0.9989088 [95,] 8.352214e-04 1.670443e-03 0.9991648 [96,] 6.565790e-04 1.313158e-03 0.9993434 [97,] 4.974491e-04 9.948982e-04 0.9995026 [98,] 3.823632e-04 7.647263e-04 0.9996176 [99,] 2.882816e-04 5.765631e-04 0.9997117 [100,] 3.240324e-04 6.480648e-04 0.9996760 [101,] 2.535551e-04 5.071101e-04 0.9997464 [102,] 2.416032e-04 4.832064e-04 0.9997584 [103,] 1.828067e-04 3.656135e-04 0.9998172 [104,] 1.344422e-04 2.688843e-04 0.9998656 [105,] 1.009570e-04 2.019140e-04 0.9998990 [106,] 7.872806e-05 1.574561e-04 0.9999213 [107,] 5.729663e-05 1.145933e-04 0.9999427 [108,] 4.279945e-05 8.559889e-05 0.9999572 [109,] 3.392834e-05 6.785669e-05 0.9999661 [110,] 2.584456e-05 5.168913e-05 0.9999742 [111,] 2.624355e-05 5.248710e-05 0.9999738 [112,] 2.420098e-05 4.840197e-05 0.9999758 [113,] 1.732470e-05 3.464939e-05 0.9999827 [114,] 1.478774e-05 2.957548e-05 0.9999852 [115,] 1.319564e-05 2.639127e-05 0.9999868 [116,] 9.364888e-06 1.872978e-05 0.9999906 [117,] 1.288950e-05 2.577901e-05 0.9999871 [118,] 1.201779e-05 2.403559e-05 0.9999880 [119,] 8.648694e-06 1.729739e-05 0.9999914 [120,] 6.244297e-06 1.248859e-05 0.9999938 [121,] 4.461509e-06 8.923017e-06 0.9999955 [122,] 3.229688e-06 6.459376e-06 0.9999968 [123,] 9.022577e-06 1.804515e-05 0.9999910 [124,] 6.413724e-06 1.282745e-05 0.9999936 [125,] 4.669087e-06 9.338174e-06 0.9999953 [126,] 4.074964e-06 8.149927e-06 0.9999959 [127,] 2.873787e-06 5.747575e-06 0.9999971 [128,] 2.474479e-06 4.948957e-06 0.9999975 [129,] 1.883472e-06 3.766944e-06 0.9999981 [130,] 2.590997e-06 5.181994e-06 0.9999974 [131,] 2.025099e-06 4.050197e-06 0.9999980 [132,] 1.427548e-06 2.855095e-06 0.9999986 [133,] 9.938341e-07 1.987668e-06 0.9999990 [134,] 7.229568e-07 1.445914e-06 0.9999993 [135,] 6.712670e-07 1.342534e-06 0.9999993 [136,] 4.931730e-07 9.863460e-07 0.9999995 [137,] 3.469252e-07 6.938505e-07 0.9999997 [138,] 2.623063e-07 5.246126e-07 0.9999997 [139,] 1.949523e-07 3.899046e-07 0.9999998 [140,] 1.588985e-07 3.177969e-07 0.9999998 [141,] 1.080558e-07 2.161116e-07 0.9999999 [142,] 7.430094e-08 1.486019e-07 0.9999999 [143,] 6.428092e-08 1.285618e-07 0.9999999 [144,] 4.721425e-08 9.442850e-08 1.0000000 [145,] 3.215284e-08 6.430568e-08 1.0000000 [146,] 2.194812e-08 4.389623e-08 1.0000000 [147,] 2.294290e-08 4.588580e-08 1.0000000 [148,] 1.531980e-08 3.063961e-08 1.0000000 [149,] 5.341125e-08 1.068225e-07 0.9999999 [150,] 3.843537e-08 7.687075e-08 1.0000000 [151,] 2.701671e-08 5.403343e-08 1.0000000 [152,] 3.067238e-08 6.134477e-08 1.0000000 [153,] 2.131269e-08 4.262538e-08 1.0000000 [154,] 1.529919e-08 3.059838e-08 1.0000000 [155,] 1.082712e-08 2.165424e-08 1.0000000 [156,] 8.952322e-09 1.790464e-08 1.0000000 [157,] 1.449938e-08 2.899875e-08 1.0000000 [158,] 9.938992e-09 1.987798e-08 1.0000000 [159,] 1.066246e-08 2.132491e-08 1.0000000 [160,] 8.361246e-09 1.672249e-08 1.0000000 [161,] 5.592068e-09 1.118414e-08 1.0000000 [162,] 6.024474e-09 1.204895e-08 1.0000000 [163,] 5.468503e-09 1.093701e-08 1.0000000 [164,] 5.140619e-09 1.028124e-08 1.0000000 [165,] 3.620362e-09 7.240724e-09 1.0000000 [166,] 2.758449e-09 5.516898e-09 1.0000000 [167,] 1.940859e-09 3.881719e-09 1.0000000 [168,] 1.283890e-09 2.567781e-09 1.0000000 [169,] 1.245866e-09 2.491732e-09 1.0000000 [170,] 1.993109e-09 3.986217e-09 1.0000000 [171,] 1.337350e-09 2.674699e-09 1.0000000 [172,] 1.078453e-09 2.156907e-09 1.0000000 [173,] 9.211799e-10 1.842360e-09 1.0000000 [174,] 7.361467e-10 1.472293e-09 1.0000000 [175,] 6.399203e-10 1.279841e-09 1.0000000 [176,] 4.216487e-10 8.432975e-10 1.0000000 [177,] 3.084678e-10 6.169355e-10 1.0000000 [178,] 2.595806e-10 5.191611e-10 1.0000000 [179,] 1.903824e-10 3.807648e-10 1.0000000 [180,] 1.256293e-10 2.512586e-10 1.0000000 [181,] 8.477927e-11 1.695585e-10 1.0000000 [182,] 2.450428e-10 4.900856e-10 1.0000000 [183,] 1.622350e-10 3.244700e-10 1.0000000 [184,] 1.103764e-10 2.207527e-10 1.0000000 [185,] 7.976877e-11 1.595375e-10 1.0000000 [186,] 5.660727e-11 1.132145e-10 1.0000000 [187,] 3.823706e-11 7.647412e-11 1.0000000 [188,] 2.883351e-11 5.766702e-11 1.0000000 [189,] 2.006762e-11 4.013524e-11 1.0000000 [190,] 1.318713e-11 2.637426e-11 1.0000000 [191,] 8.518418e-12 1.703684e-11 1.0000000 [192,] 5.479069e-12 1.095814e-11 1.0000000 [193,] 1.227070e-11 2.454141e-11 1.0000000 [194,] 8.341323e-12 1.668265e-11 1.0000000 [195,] 5.583735e-12 1.116747e-11 1.0000000 [196,] 3.775715e-12 7.551431e-12 1.0000000 [197,] 2.458323e-12 4.916645e-12 1.0000000 [198,] 1.691552e-12 3.383105e-12 1.0000000 [199,] 1.134296e-12 2.268591e-12 1.0000000 [200,] 8.650435e-13 1.730087e-12 1.0000000 [201,] 5.575458e-13 1.115092e-12 1.0000000 [202,] 7.084073e-13 1.416815e-12 1.0000000 [203,] 6.367157e-13 1.273431e-12 1.0000000 [204,] 4.286692e-13 8.573383e-13 1.0000000 [205,] 2.723161e-13 5.446322e-13 1.0000000 [206,] 1.855079e-13 3.710157e-13 1.0000000 [207,] 3.298526e-13 6.597053e-13 1.0000000 [208,] 2.128036e-13 4.256071e-13 1.0000000 [209,] 2.630825e-13 5.261649e-13 1.0000000 [210,] 2.048781e-13 4.097561e-13 1.0000000 [211,] 5.220070e-13 1.044014e-12 1.0000000 [212,] 4.516800e-13 9.033599e-13 1.0000000 [213,] 7.631463e-13 1.526293e-12 1.0000000 [214,] 5.114072e-13 1.022814e-12 1.0000000 [215,] 3.303118e-13 6.606235e-13 1.0000000 [216,] 6.917327e-13 1.383465e-12 1.0000000 [217,] 7.044679e-13 1.408936e-12 1.0000000 [218,] 4.597834e-13 9.195667e-13 1.0000000 [219,] 5.209782e-13 1.041956e-12 1.0000000 [220,] 3.584734e-13 7.169469e-13 1.0000000 [221,] 2.286388e-13 4.572775e-13 1.0000000 [222,] 1.480654e-13 2.961307e-13 1.0000000 [223,] 1.570332e-13 3.140663e-13 1.0000000 [224,] 1.935326e-13 3.870652e-13 1.0000000 [225,] 1.228239e-13 2.456478e-13 1.0000000 [226,] 9.588019e-14 1.917604e-13 1.0000000 [227,] 1.295358e-13 2.590716e-13 1.0000000 [228,] 8.371244e-14 1.674249e-13 1.0000000 [229,] 5.586025e-14 1.117205e-13 1.0000000 [230,] 3.954830e-14 7.909659e-14 1.0000000 [231,] 2.495882e-14 4.991765e-14 1.0000000 [232,] 1.789673e-14 3.579347e-14 1.0000000 [233,] 1.223340e-14 2.446681e-14 1.0000000 [234,] 9.931446e-15 1.986289e-14 1.0000000 [235,] 1.134925e-14 2.269849e-14 1.0000000 [236,] 7.313957e-15 1.462791e-14 1.0000000 [237,] 5.553134e-15 1.110627e-14 1.0000000 [238,] 4.570841e-15 9.141682e-15 1.0000000 [239,] 3.044556e-15 6.089112e-15 1.0000000 [240,] 1.325265e-14 2.650530e-14 1.0000000 [241,] 8.925160e-15 1.785032e-14 1.0000000 [242,] 6.774596e-15 1.354919e-14 1.0000000 [243,] 1.143305e-14 2.286610e-14 1.0000000 [244,] 7.307420e-15 1.461484e-14 1.0000000 [245,] 8.236901e-15 1.647380e-14 1.0000000 [246,] 5.576099e-15 1.115220e-14 1.0000000 [247,] 3.546446e-15 7.092892e-15 1.0000000 [248,] 2.344739e-15 4.689479e-15 1.0000000 [249,] 1.925479e-15 3.850958e-15 1.0000000 [250,] 1.628500e-15 3.257000e-15 1.0000000 [251,] 1.098641e-15 2.197282e-15 1.0000000 [252,] 6.774026e-16 1.354805e-15 1.0000000 [253,] 5.521536e-16 1.104307e-15 1.0000000 [254,] 3.466018e-16 6.932036e-16 1.0000000 [255,] 3.518843e-16 7.037686e-16 1.0000000 [256,] 7.923744e-16 1.584749e-15 1.0000000 [257,] 6.505133e-16 1.301027e-15 1.0000000 [258,] 5.792761e-16 1.158552e-15 1.0000000 [259,] 5.960058e-16 1.192012e-15 1.0000000 [260,] 5.626877e-16 1.125375e-15 1.0000000 [261,] 6.633911e-16 1.326782e-15 1.0000000 [262,] 4.628907e-16 9.257814e-16 1.0000000 [263,] 2.849449e-16 5.698897e-16 1.0000000 [264,] 1.794121e-16 3.588243e-16 1.0000000 [265,] 3.239820e-16 6.479639e-16 1.0000000 [266,] 2.367409e-16 4.734818e-16 1.0000000 [267,] 1.631834e-16 3.263668e-16 1.0000000 [268,] 2.636739e-16 5.273477e-16 1.0000000 [269,] 1.684131e-16 3.368263e-16 1.0000000 [270,] 1.111289e-16 2.222577e-16 1.0000000 [271,] 6.868188e-17 1.373638e-16 1.0000000 [272,] 4.286268e-17 8.572536e-17 1.0000000 [273,] 5.145662e-17 1.029132e-16 1.0000000 [274,] 4.711355e-17 9.422710e-17 1.0000000 [275,] 3.272517e-17 6.545035e-17 1.0000000 [276,] 3.483986e-17 6.967972e-17 1.0000000 [277,] 3.552150e-17 7.104299e-17 1.0000000 [278,] 2.320878e-17 4.641755e-17 1.0000000 [279,] 1.957010e-17 3.914019e-17 1.0000000 [280,] 1.394545e-17 2.789090e-17 1.0000000 [281,] 2.497321e-17 4.994642e-17 1.0000000 [282,] 1.732185e-17 3.464370e-17 1.0000000 [283,] 1.910718e-17 3.821436e-17 1.0000000 [284,] 1.191629e-17 2.383259e-17 1.0000000 [285,] 8.881598e-18 1.776320e-17 1.0000000 [286,] 5.760248e-18 1.152050e-17 1.0000000 [287,] 3.532291e-18 7.064582e-18 1.0000000 [288,] 2.548875e-18 5.097750e-18 1.0000000 [289,] 1.958142e-18 3.916285e-18 1.0000000 [290,] 1.373428e-18 2.746856e-18 1.0000000 [291,] 1.606592e-18 3.213183e-18 1.0000000 [292,] 1.880355e-18 3.760711e-18 1.0000000 [293,] 1.601426e-18 3.202852e-18 1.0000000 [294,] 1.089852e-18 2.179703e-18 1.0000000 [295,] 6.635263e-19 1.327053e-18 1.0000000 [296,] 1.350932e-18 2.701865e-18 1.0000000 [297,] 1.713239e-17 3.426477e-17 1.0000000 [298,] 1.089004e-17 2.178008e-17 1.0000000 [299,] 1.340213e-17 2.680425e-17 1.0000000 [300,] 1.917310e-17 3.834620e-17 1.0000000 [301,] 2.607296e-17 5.214593e-17 1.0000000 [302,] 6.124716e-17 1.224943e-16 1.0000000 [303,] 3.786803e-17 7.573606e-17 1.0000000 [304,] 4.035995e-17 8.071990e-17 1.0000000 [305,] 8.093236e-17 1.618647e-16 1.0000000 [306,] 7.059522e-17 1.411904e-16 1.0000000 [307,] 6.587377e-17 1.317475e-16 1.0000000 [308,] 4.916140e-17 9.832281e-17 1.0000000 [309,] 5.298781e-17 1.059756e-16 1.0000000 [310,] 3.739065e-16 7.478129e-16 1.0000000 [311,] 9.027786e-15 1.805557e-14 1.0000000 [312,] 7.635527e-14 1.527105e-13 1.0000000 [313,] 5.097287e-14 1.019457e-13 1.0000000 [314,] 9.937771e-14 1.987554e-13 1.0000000 [315,] 5.781977e-13 1.156395e-12 1.0000000 [316,] 3.979301e-13 7.958602e-13 1.0000000 [317,] 1.119520e-12 2.239040e-12 1.0000000 [318,] 7.663259e-13 1.532652e-12 1.0000000 [319,] 5.171716e-13 1.034343e-12 1.0000000 [320,] 3.672713e-13 7.345427e-13 1.0000000 [321,] 2.638293e-13 5.276587e-13 1.0000000 [322,] 6.305860e-13 1.261172e-12 1.0000000 [323,] 1.027813e-12 2.055626e-12 1.0000000 [324,] 7.066316e-13 1.413263e-12 1.0000000 [325,] 6.759839e-13 1.351968e-12 1.0000000 [326,] 4.894694e-13 9.789388e-13 1.0000000 [327,] 3.268200e-13 6.536400e-13 1.0000000 [328,] 4.285183e-13 8.570366e-13 1.0000000 [329,] 2.849991e-13 5.699982e-13 1.0000000 [330,] 2.593729e-13 5.187459e-13 1.0000000 [331,] 8.097451e-13 1.619490e-12 1.0000000 [332,] 5.804880e-13 1.160976e-12 1.0000000 [333,] 4.049333e-13 8.098666e-13 1.0000000 [334,] 1.365977e-12 2.731954e-12 1.0000000 [335,] 1.093748e-12 2.187497e-12 1.0000000 [336,] 8.073483e-13 1.614697e-12 1.0000000 [337,] 7.134804e-13 1.426961e-12 1.0000000 [338,] 5.114100e-13 1.022820e-12 1.0000000 [339,] 1.429064e-12 2.858127e-12 1.0000000 [340,] 1.918961e-12 3.837922e-12 1.0000000 [341,] 1.706394e-12 3.412789e-12 1.0000000 [342,] 1.243557e-12 2.487114e-12 1.0000000 [343,] 1.139273e-12 2.278547e-12 1.0000000 [344,] 1.110462e-12 2.220924e-12 1.0000000 [345,] 9.394064e-13 1.878813e-12 1.0000000 [346,] 7.890124e-13 1.578025e-12 1.0000000 [347,] 9.571044e-13 1.914209e-12 1.0000000 [348,] 2.424500e-12 4.849001e-12 1.0000000 [349,] 1.637825e-12 3.275650e-12 1.0000000 [350,] 1.152166e-12 2.304333e-12 1.0000000 [351,] 8.056422e-12 1.611284e-11 1.0000000 [352,] 7.635557e-12 1.527111e-11 1.0000000 [353,] 5.214073e-12 1.042815e-11 1.0000000 [354,] 7.213833e-12 1.442767e-11 1.0000000 [355,] 2.404704e-11 4.809408e-11 1.0000000 [356,] 4.932800e-11 9.865600e-11 1.0000000 [357,] 5.368805e-11 1.073761e-10 1.0000000 [358,] 8.726377e-11 1.745275e-10 1.0000000 [359,] 1.837376e-10 3.674753e-10 1.0000000 [360,] 1.600493e-10 3.200986e-10 1.0000000 [361,] 1.155555e-10 2.311110e-10 1.0000000 [362,] 2.613393e-09 5.226786e-09 1.0000000 [363,] 1.940916e-09 3.881831e-09 1.0000000 [364,] 1.774473e-09 3.548947e-09 1.0000000 [365,] 1.860926e-09 3.721852e-09 1.0000000 [366,] 3.254024e-09 6.508049e-09 1.0000000 [367,] 2.345376e-09 4.690753e-09 1.0000000 [368,] 1.680136e-09 3.360271e-09 1.0000000 [369,] 1.523975e-09 3.047950e-09 1.0000000 [370,] 1.300372e-09 2.600744e-09 1.0000000 [371,] 1.497162e-09 2.994324e-09 1.0000000 [372,] 1.339484e-09 2.678967e-09 1.0000000 [373,] 1.007342e-09 2.014683e-09 1.0000000 [374,] 7.128220e-10 1.425644e-09 1.0000000 [375,] 6.055950e-10 1.211190e-09 1.0000000 [376,] 7.562637e-10 1.512527e-09 1.0000000 [377,] 9.469629e-10 1.893926e-09 1.0000000 [378,] 1.144836e-09 2.289672e-09 1.0000000 [379,] 1.102286e-09 2.204571e-09 1.0000000 [380,] 9.470344e-10 1.894069e-09 1.0000000 [381,] 7.355461e-10 1.471092e-09 1.0000000 [382,] 5.401194e-10 1.080239e-09 1.0000000 [383,] 5.984562e-10 1.196912e-09 1.0000000 [384,] 5.242469e-10 1.048494e-09 1.0000000 [385,] 4.165790e-10 8.331579e-10 1.0000000 [386,] 4.099375e-10 8.198751e-10 1.0000000 [387,] 3.485279e-10 6.970557e-10 1.0000000 [388,] 4.198740e-10 8.397479e-10 1.0000000 [389,] 3.054189e-10 6.108377e-10 1.0000000 [390,] 2.150413e-10 4.300825e-10 1.0000000 [391,] 1.633905e-10 3.267809e-10 1.0000000 [392,] 1.667992e-10 3.335984e-10 1.0000000 [393,] 1.613137e-10 3.226274e-10 1.0000000 [394,] 3.495948e-10 6.991896e-10 1.0000000 [395,] 2.520856e-10 5.041713e-10 1.0000000 [396,] 1.873528e-10 3.747055e-10 1.0000000 [397,] 1.307783e-10 2.615566e-10 1.0000000 [398,] 9.319728e-11 1.863946e-10 1.0000000 [399,] 7.376382e-11 1.475276e-10 1.0000000 [400,] 5.561261e-11 1.112252e-10 1.0000000 [401,] 3.996106e-11 7.992212e-11 1.0000000 [402,] 5.575713e-11 1.115143e-10 1.0000000 [403,] 4.835188e-11 9.670377e-11 1.0000000 [404,] 3.642998e-11 7.285996e-11 1.0000000 [405,] 5.852289e-11 1.170458e-10 1.0000000 [406,] 5.074022e-11 1.014804e-10 1.0000000 [407,] 3.718104e-11 7.436209e-11 1.0000000 [408,] 2.931538e-11 5.863076e-11 1.0000000 [409,] 3.698110e-11 7.396221e-11 1.0000000 [410,] 3.451863e-11 6.903725e-11 1.0000000 [411,] 3.271052e-11 6.542104e-11 1.0000000 [412,] 6.590425e-10 1.318085e-09 1.0000000 [413,] 8.887272e-10 1.777454e-09 1.0000000 [414,] 6.901545e-10 1.380309e-09 1.0000000 [415,] 1.378349e-09 2.756698e-09 1.0000000 [416,] 9.924098e-10 1.984820e-09 1.0000000 [417,] 7.042062e-10 1.408412e-09 1.0000000 [418,] 4.947935e-10 9.895870e-10 1.0000000 [419,] 3.483765e-10 6.967529e-10 1.0000000 [420,] 5.571752e-10 1.114350e-09 1.0000000 [421,] 4.178349e-10 8.356698e-10 1.0000000 [422,] 3.185934e-10 6.371867e-10 1.0000000 [423,] 6.167757e-10 1.233551e-09 1.0000000 [424,] 6.547020e-10 1.309404e-09 1.0000000 [425,] 5.496234e-10 1.099247e-09 1.0000000 [426,] 6.403088e-10 1.280618e-09 1.0000000 [427,] 4.646709e-10 9.293417e-10 1.0000000 [428,] 3.573410e-10 7.146820e-10 1.0000000 [429,] 5.192185e-10 1.038437e-09 1.0000000 [430,] 3.664809e-10 7.329618e-10 1.0000000 [431,] 3.829846e-10 7.659693e-10 1.0000000 [432,] 2.149568e-09 4.299136e-09 1.0000000 [433,] 1.608131e-09 3.216263e-09 1.0000000 [434,] 4.925912e-09 9.851823e-09 1.0000000 [435,] 5.255267e-09 1.051053e-08 1.0000000 [436,] 4.757202e-09 9.514404e-09 1.0000000 [437,] 3.674419e-09 7.348838e-09 1.0000000 [438,] 3.377604e-09 6.755208e-09 1.0000000 [439,] 1.053951e-08 2.107902e-08 1.0000000 [440,] 8.212794e-09 1.642559e-08 1.0000000 [441,] 6.481168e-09 1.296234e-08 1.0000000 [442,] 2.119654e-07 4.239308e-07 0.9999998 [443,] 3.551544e-07 7.103087e-07 0.9999996 [444,] 8.742091e-07 1.748418e-06 0.9999991 [445,] 8.294589e-07 1.658918e-06 0.9999992 [446,] 9.144949e-07 1.828990e-06 0.9999991 [447,] 2.178837e-06 4.357674e-06 0.9999978 [448,] 1.693142e-06 3.386285e-06 0.9999983 [449,] 2.076025e-06 4.152050e-06 0.9999979 [450,] 1.607054e-06 3.214108e-06 0.9999984 [451,] 1.245252e-06 2.490505e-06 0.9999988 [452,] 9.491952e-07 1.898390e-06 0.9999991 [453,] 7.600360e-07 1.520072e-06 0.9999992 [454,] 1.236869e-06 2.473738e-06 0.9999988 [455,] 9.499727e-07 1.899945e-06 0.9999991 [456,] 1.360089e-06 2.720178e-06 0.9999986 [457,] 1.915461e-06 3.830921e-06 0.9999981 [458,] 1.804548e-06 3.609097e-06 0.9999982 [459,] 1.883096e-06 3.766192e-06 0.9999981 [460,] 1.415133e-06 2.830267e-06 0.9999986 [461,] 1.102320e-06 2.204640e-06 0.9999989 [462,] 8.775423e-07 1.755085e-06 0.9999991 [463,] 1.253613e-05 2.507225e-05 0.9999875 [464,] 1.316928e-05 2.633857e-05 0.9999868 [465,] 2.078487e-05 4.156974e-05 0.9999792 [466,] 2.701555e-05 5.403110e-05 0.9999730 [467,] 2.688031e-05 5.376061e-05 0.9999731 [468,] 3.137283e-05 6.274567e-05 0.9999686 [469,] 1.676100e-04 3.352201e-04 0.9998324 [470,] 2.440308e-04 4.880617e-04 0.9997560 [471,] 2.443911e-04 4.887821e-04 0.9997556 [472,] 2.029115e-04 4.058230e-04 0.9997971 [473,] 1.619416e-04 3.238832e-04 0.9998381 [474,] 1.276333e-04 2.552666e-04 0.9998724 [475,] 1.385318e-04 2.770635e-04 0.9998615 [476,] 2.032425e-04 4.064850e-04 0.9997968 [477,] 1.838360e-04 3.676719e-04 0.9998162 [478,] 1.581112e-04 3.162224e-04 0.9998419 [479,] 1.488541e-04 2.977082e-04 0.9998511 [480,] 2.099193e-04 4.198386e-04 0.9997901 [481,] 1.706507e-04 3.413014e-04 0.9998293 [482,] 1.398394e-04 2.796787e-04 0.9998602 [483,] 1.155582e-04 2.311163e-04 0.9998844 [484,] 4.221337e-04 8.442673e-04 0.9995779 [485,] 3.401579e-04 6.803157e-04 0.9996598 [486,] 3.054923e-04 6.109845e-04 0.9996945 [487,] 3.004259e-04 6.008517e-04 0.9996996 [488,] 3.782888e-04 7.565776e-04 0.9996217 [489,] 3.518710e-04 7.037419e-04 0.9996481 [490,] 3.151631e-04 6.303263e-04 0.9996848 [491,] 2.544675e-04 5.089350e-04 0.9997455 [492,] 2.125704e-04 4.251408e-04 0.9997874 [493,] 1.695585e-04 3.391171e-04 0.9998304 [494,] 3.604819e-04 7.209638e-04 0.9996395 [495,] 2.913275e-04 5.826550e-04 0.9997087 [496,] 2.495079e-04 4.990157e-04 0.9997505 [497,] 2.733424e-04 5.466848e-04 0.9997267 [498,] 2.691221e-04 5.382442e-04 0.9997309 [499,] 2.129392e-04 4.258785e-04 0.9997871 [500,] 3.724269e-04 7.448537e-04 0.9996276 [501,] 3.031210e-04 6.062420e-04 0.9996969 [502,] 2.491055e-04 4.982110e-04 0.9997509 [503,] 2.002710e-04 4.005420e-04 0.9997997 [504,] 1.968464e-04 3.936929e-04 0.9998032 [505,] 2.664711e-04 5.329423e-04 0.9997335 [506,] 2.101771e-04 4.203543e-04 0.9997898 [507,] 3.221960e-04 6.443921e-04 0.9996778 [508,] 2.586590e-04 5.173180e-04 0.9997413 [509,] 2.096906e-04 4.193813e-04 0.9997903 [510,] 1.809666e-04 3.619331e-04 0.9998190 [511,] 1.599289e-04 3.198577e-04 0.9998401 [512,] 1.275151e-04 2.550303e-04 0.9998725 [513,] 1.014495e-04 2.028991e-04 0.9998986 [514,] 8.230838e-05 1.646168e-04 0.9999177 [515,] 8.473754e-05 1.694751e-04 0.9999153 [516,] 9.697498e-05 1.939500e-04 0.9999030 [517,] 1.005986e-04 2.011971e-04 0.9998994 [518,] 7.989613e-05 1.597923e-04 0.9999201 [519,] 6.108396e-05 1.221679e-04 0.9999389 [520,] 5.792604e-05 1.158521e-04 0.9999421 [521,] 5.265925e-05 1.053185e-04 0.9999473 [522,] 4.870027e-05 9.740054e-05 0.9999513 [523,] 3.893271e-05 7.786542e-05 0.9999611 [524,] 3.068253e-05 6.136507e-05 0.9999693 [525,] 2.493555e-05 4.987111e-05 0.9999751 [526,] 1.898059e-05 3.796118e-05 0.9999810 [527,] 1.407829e-05 2.815658e-05 0.9999859 [528,] 6.137992e-05 1.227598e-04 0.9999386 [529,] 5.127358e-05 1.025472e-04 0.9999487 [530,] 5.686904e-05 1.137381e-04 0.9999431 [531,] 1.905036e-04 3.810072e-04 0.9998095 [532,] 1.636448e-04 3.272895e-04 0.9998364 [533,] 1.394423e-04 2.788846e-04 0.9998606 [534,] 4.657059e-04 9.314118e-04 0.9995343 [535,] 4.586763e-04 9.173525e-04 0.9995413 [536,] 5.237319e-04 1.047464e-03 0.9994763 [537,] 4.044943e-04 8.089887e-04 0.9995955 [538,] 3.203417e-04 6.406833e-04 0.9996797 [539,] 2.619243e-04 5.238486e-04 0.9997381 [540,] 2.158409e-04 4.316818e-04 0.9997842 [541,] 1.979538e-04 3.959075e-04 0.9998020 [542,] 1.955620e-04 3.911241e-04 0.9998044 [543,] 1.491321e-04 2.982642e-04 0.9998509 [544,] 1.128077e-04 2.256155e-04 0.9998872 [545,] 8.622044e-05 1.724409e-04 0.9999138 [546,] 4.604482e-04 9.208964e-04 0.9995396 [547,] 3.558492e-04 7.116984e-04 0.9996442 [548,] 4.685296e-04 9.370591e-04 0.9995315 [549,] 6.146621e-04 1.229324e-03 0.9993853 [550,] 6.560380e-04 1.312076e-03 0.9993440 [551,] 5.177861e-04 1.035572e-03 0.9994822 [552,] 4.697030e-04 9.394060e-04 0.9995303 [553,] 8.249810e-04 1.649962e-03 0.9991750 [554,] 6.333466e-04 1.266693e-03 0.9993667 [555,] 6.416498e-04 1.283300e-03 0.9993584 [556,] 7.396190e-04 1.479238e-03 0.9992604 [557,] 6.775385e-04 1.355077e-03 0.9993225 [558,] 5.354293e-04 1.070859e-03 0.9994646 [559,] 4.040543e-04 8.081087e-04 0.9995959 [560,] 3.209746e-04 6.419493e-04 0.9996790 [561,] 2.417051e-04 4.834102e-04 0.9997583 [562,] 2.245464e-04 4.490928e-04 0.9997755 [563,] 1.053964e-03 2.107929e-03 0.9989460 [564,] 8.925915e-04 1.785183e-03 0.9991074 [565,] 1.061664e-03 2.123327e-03 0.9989383 [566,] 8.579606e-04 1.715921e-03 0.9991420 [567,] 8.969425e-04 1.793885e-03 0.9991031 [568,] 2.379487e-03 4.758975e-03 0.9976205 [569,] 2.645481e-03 5.290962e-03 0.9973545 [570,] 2.048975e-03 4.097950e-03 0.9979510 [571,] 3.922171e-03 7.844342e-03 0.9960778 [572,] 5.181280e-03 1.036256e-02 0.9948187 [573,] 8.014074e-03 1.602815e-02 0.9919859 [574,] 6.555522e-03 1.311104e-02 0.9934445 [575,] 7.936503e-03 1.587301e-02 0.9920635 [576,] 6.955297e-03 1.391059e-02 0.9930447 [577,] 6.020651e-03 1.204130e-02 0.9939793 [578,] 4.806177e-03 9.612353e-03 0.9951938 [579,] 3.736587e-03 7.473173e-03 0.9962634 [580,] 7.882089e-03 1.576418e-02 0.9921179 [581,] 6.171335e-03 1.234267e-02 0.9938287 [582,] 4.974935e-03 9.949870e-03 0.9950251 [583,] 3.961160e-03 7.922320e-03 0.9960388 [584,] 3.498982e-03 6.997964e-03 0.9965010 [585,] 5.119498e-03 1.023900e-02 0.9948805 [586,] 4.320891e-03 8.641783e-03 0.9956791 [587,] 3.726111e-03 7.452223e-03 0.9962739 [588,] 5.370675e-03 1.074135e-02 0.9946293 [589,] 4.113293e-03 8.226587e-03 0.9958867 [590,] 3.319439e-03 6.638878e-03 0.9966806 [591,] 2.713370e-03 5.426741e-03 0.9972866 [592,] 9.272002e-03 1.854400e-02 0.9907280 [593,] 8.105773e-03 1.621155e-02 0.9918942 [594,] 1.235501e-02 2.471001e-02 0.9876450 [595,] 3.366302e-02 6.732603e-02 0.9663370 [596,] 2.744601e-02 5.489203e-02 0.9725540 [597,] 2.536577e-02 5.073154e-02 0.9746342 [598,] 3.047384e-02 6.094768e-02 0.9695262 [599,] 2.638609e-02 5.277218e-02 0.9736139 [600,] 2.309726e-02 4.619451e-02 0.9769027 [601,] 1.903655e-02 3.807309e-02 0.9809635 [602,] 1.503883e-02 3.007765e-02 0.9849612 [603,] 1.180121e-02 2.360243e-02 0.9881988 [604,] 9.261364e-03 1.852273e-02 0.9907386 [605,] 9.125366e-03 1.825073e-02 0.9908746 [606,] 7.540212e-03 1.508042e-02 0.9924598 [607,] 7.083576e-03 1.416715e-02 0.9929164 [608,] 5.472754e-03 1.094551e-02 0.9945272 [609,] 4.413701e-03 8.827402e-03 0.9955863 [610,] 3.101669e-03 6.203339e-03 0.9968983 [611,] 7.163315e-03 1.432663e-02 0.9928367 [612,] 5.181480e-03 1.036296e-02 0.9948185 [613,] 3.977141e-03 7.954282e-03 0.9960229 [614,] 1.228905e-01 2.457809e-01 0.8771095 [615,] 1.093490e-01 2.186979e-01 0.8906510 [616,] 9.579312e-02 1.915862e-01 0.9042069 [617,] 7.560365e-02 1.512073e-01 0.9243964 [618,] 6.360374e-02 1.272075e-01 0.9363963 [619,] 1.223904e-01 2.447808e-01 0.8776096 [620,] 1.183970e-01 2.367939e-01 0.8816030 [621,] 9.945096e-02 1.989019e-01 0.9005490 [622,] 2.203682e-01 4.407364e-01 0.7796318 [623,] 2.150684e-01 4.301368e-01 0.7849316 [624,] 2.147465e-01 4.294929e-01 0.7852535 [625,] 1.660237e-01 3.320474e-01 0.8339763 [626,] 5.123349e-01 9.753302e-01 0.4876651 [627,] 4.356767e-01 8.713535e-01 0.5643233 [628,] 3.651906e-01 7.303812e-01 0.6348094 [629,] 2.912347e-01 5.824695e-01 0.7087653 [630,] 2.373092e-01 4.746183e-01 0.7626908 [631,] 4.867398e-01 9.734796e-01 0.5132602 [632,] 3.992097e-01 7.984194e-01 0.6007903 [633,] 3.229024e-01 6.458048e-01 0.6770976 [634,] 2.398189e-01 4.796378e-01 0.7601811 [635,] 3.073279e-01 6.146558e-01 0.6926721 [636,] 2.426154e-01 4.852308e-01 0.7573846 > postscript(file="/var/www/html/rcomp/tmp/1n73f1292085873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2n73f1292085873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3n73f1292085873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4xy2i1292085873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5xy2i1292085873.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 = 651 Frequency = 1 1 2 3 4 5 -3508.885708 -13955.834059 6.900395 -2485.257352 -3216.871153 6 7 8 9 10 -658.009907 -406.820354 -3396.013454 -2097.203007 -4245.096058 11 12 13 14 15 1550.196896 -1068.171201 -37.203007 -4415.995352 60.214998 16 17 18 19 20 2599.721851 2420.911403 4310.196896 -697.678053 -1253.888403 21 22 23 24 25 4304.300298 -1237.313502 -1194.312651 1087.886691 1035.083706 26 27 28 29 30 -1763.081503 3200.104503 6807.393543 -2665.096058 3124.300298 31 32 33 34 35 -1395.502204 3737.393543 660.690044 3348.089948 -217.203007 36 37 38 39 40 3032.321947 1266.425349 -6024.381551 -1454.381551 7684.903942 41 42 43 44 45 654.811548 4435.507953 -606.396107 2708.896848 171.114292 46 47 48 49 50 -1505.681600 -2374.782306 -1598.595449 1724.004648 2743.894153 51 52 53 54 55 2254.097041 8266.111597 -2413.574651 -733.191998 3533.128847 56 57 58 59 60 3418.693591 1869.883144 5641.386449 -3961.799557 -3340.591902 61 62 63 64 65 1817.393543 2391.003797 6746.789899 -4145.188452 -3704.805799 66 67 68 69 70 -2502.274603 -1908.503055 2588.089948 909.500491 -1874.381551 71 72 73 74 75 4368.200443 1906.918497 352.321947 -4343.081503 -1906.488500 76 77 78 79 80 3204.701053 1195.507953 3948.287445 2720.196896 1866.093494 81 82 83 84 85 2480.196896 -5536.802252 -1503.574651 -1099.309956 -4628.816808 86 87 88 89 90 -3164.874699 -703.888403 424.701053 -6033.888403 4574.097041 91 92 93 94 95 -730.098754 -5737.295400 -842.606457 887.393543 -4308.484953 96 97 98 99 100 -5682.274603 -5643.888403 -1951.085050 -4213.574651 -1869.378856 101 102 103 104 105 -1255.681600 -2439.291854 1049.500491 666.001101 -5.188452 106 107 108 109 110 -918.595449 -3874.713406 2938.514195 3397.393543 2122.303845 111 112 113 114 115 -1007.295400 352.704600 -412.292705 -608.009907 -213.099605 116 117 118 119 120 4315.107199 -1196.396107 3418.914950 2033.128847 1076.900395 121 122 123 124 125 -1557.696155 -2667.788548 -594.381551 -4997.203007 -3823.574651 126 127 128 129 130 836.517742 1973.110745 561.990093 -2844.695304 -8229.785002 131 132 133 134 135 -642.767751 -1118.885708 2271.828799 697.614902 -2567.788548 136 137 138 139 140 -2163.413358 -6487.678053 -2797.295400 -60.499509 313.018352 141 142 143 144 145 594.097041 769.901246 -2867.678053 -777.788548 1988.200443 146 147 148 149 150 602.483241 597.794297 -444.220258 -1207.678053 -4744.695304 151 152 153 154 155 1229.500491 306.900395 197.232249 4439.228703 -49.692608 156 157 158 159 160 -8312.113309 -558.484953 663.511500 -4167.203007 697.283047 161 162 163 164 165 -1411.799557 -1440.098754 3133.018352 4844.903942 1217.614902 166 167 168 169 170 2764.811548 2814.701053 -140.278149 3834.318400 4212.594104 171 172 173 174 175 2273.290141 -2616.802252 -2284.289158 -1729.692608 -505.096058 176 177 178 179 180 5293.493398 -5554.289158 1082.142552 3473.986546 2611.114292 181 182 183 184 185 -3968.484953 -3981.467702 154.318400 738.896848 -3173.413358 186 187 188 189 190 -123.888403 -811.085050 -2877.355103 -6566.396107 -976.871153 191 192 193 194 195 1637.614902 2212.211452 1021.183192 -1660.098754 198.896848 196 197 198 199 200 1445.194201 890.597651 45.710842 732.303845 4486.918497 201 202 203 204 205 -1542.606457 -792.316198 1706.093494 449.389996 884.903942 206 207 208 209 210 -1621.891950 -3603.081503 -1339.692608 3861.183192 2650.307392 211 212 213 214 215 2117.707295 -1295.902959 3022.193350 -6181.306409 -371.509298 216 217 218 219 220 -5069.623708 -2949.692608 -6407.203007 5445.489851 -6873.888403 221 222 223 224 225 -694.782306 1258.421802 3363.690896 -4660.591902 -1498.978101 226 227 228 229 230 -3493.099605 -64.399653 -308.978101 -351.306409 5892.193350 231 232 233 234 235 -4701.891950 -617.203007 1214.004648 3437.886691 1439.389996 236 237 238 239 240 656.425349 2811.003797 47.232249 -2916.709859 1382.796993 241 242 243 244 245 -2017.696155 -5508.484953 -220.591902 -1963.906506 1758.693591 246 247 248 249 250 -1476.396107 -8175.902959 -14.399653 1535.618449 -6884.695304 251 252 253 254 255 -1308.816808 2019.901246 -2723.888403 -1223.413358 -793.099605 256 257 258 259 260 4143.493398 2814.318400 765.618449 -703.482258 1036.111597 261 262 263 264 265 1102.211452 2764.811548 5703.511500 -3519.692608 3627.707295 266 267 268 269 270 4207.707295 4668.896848 2786.918497 -114.695304 -438.009907 271 272 273 274 275 -1579.785002 -6693.081503 -1170.610004 1553.511500 5975.507953 276 277 278 279 280 -1845.096058 265.507953 -249.309956 -83.191998 3505.618449 281 282 283 284 285 3390.690044 -2302.606457 5156.789899 3413.603893 790.104503 286 287 288 289 290 2790.621144 1451.897699 5370.690044 -2896.709859 -4678.885708 291 292 293 294 295 -1152.767751 686.111597 862.211452 73.603893 -3476.871153 296 297 298 299 300 464.207905 1298.421802 3590.528751 3214.903942 -3898.978101 301 302 303 304 305 1224.793446 -791.799557 5491.496945 8850.911403 287.393543 306 307 308 309 310 -4293.906506 -6976.691757 -6737.184905 2128.260140 -1259.268361 311 312 313 314 315 934.115143 3872.031688 -4729.360754 -5449.674506 1057.010890 316 317 318 319 320 -424.160561 5172.925590 7509.629089 6960.011741 -348.775213 321 322 323 324 325 5271.607440 7703.331736 -2356.378005 6402.907488 -2130.974555 326 327 328 329 330 -1436.285611 -37.991805 -646.267509 -10035.045266 -6672.902108 331 332 333 334 335 -957.991805 2712.100588 1013.331736 -162.989110 -6071.781455 336 337 338 339 340 -894.178663 -4480.974555 -8596.378005 544.428896 -2145.077956 341 342 343 344 345 6304.115143 -3200.882161 442.008195 2609.924739 998.310938 346 347 348 349 350 5888.329040 3949.135941 2064.115143 1197.504038 2678.260140 351 352 353 354 355 2966.535844 2103.621995 1796.604744 4302.414340 3935.761635 356 357 358 359 360 -215.861364 213.221241 8995.839439 -4446.378005 -635.077956 361 362 363 364 365 4829.135941 7440.818642 5484.461586 -5122.588355 -6257.074410 366 367 368 369 370 -8863.141206 -555.451411 -524.671811 11513.248540 358.711693 371 372 373 374 375 -3622.495962 -4867.991805 -7085.077956 -730.882161 -10.006361 376 377 378 379 380 -3681.288307 2891.717935 -5159.259163 -3950.167654 1339.993639 381 382 383 384 385 -629.641816 1902.501343 4548.218545 4998.310938 -5835.086861 386 387 388 389 390 3623.714389 2733.732491 -1945.861364 702.414340 4154.922044 391 392 393 394 395 2633.308243 -2677.567558 3547.010890 -3031.375309 -5494.271056 396 397 398 399 400 -536.778759 -736.285611 794.225639 -4325.460609 3198.320136 401 402 403 404 405 6884.115143 -804.671811 1168.711693 -255.368216 -1248.775213 406 407 408 409 410 -2734.602911 -2173.371762 -1640.882161 -5987.991805 -3520.882161 411 412 413 414 415 -2656.345314 -6681.375309 2805.797844 925.328189 2194.138636 416 417 418 419 420 -5592.588355 -3835.884857 3585.397089 12208.651990 5552.100588 421 422 423 424 425 -2492.588355 -6962.095207 -1608.060706 560.842135 -794.271056 426 427 428 429 430 -340.457914 6038.310938 1509.518593 -2492.989110 -7142.071714 431 432 433 434 435 -4789.199460 -2913.362565 -5107.991805 946.628238 1741.607440 436 437 438 439 440 6040.758939 715.014437 4424.913139 10113.308243 1225.032539 441 442 443 444 445 -9125.953757 -5284.764204 3121.123490 -1644.202156 3094.428896 446 447 448 449 450 -8952.164107 -38.688210 2092.907488 14725.415192 6964.115143 451 452 453 454 455 8412.432442 -4077.092512 4388.633888 8646.222092 1420.818642 456 457 458 459 460 -5246.267509 -1938.060706 -1248.858408 957.411645 -2232.164107 461 462 463 464 465 7506.715240 267.817790 -5749.988259 -6257.991805 -4263.371762 466 467 468 469 470 -4605.409811 183.714389 -7.746953 -1270.472209 -12926.982016 471 472 473 474 475 -4487.991805 6162.414340 5107.435138 -4930.974555 -6310.167654 476 477 478 479 480 -11613.362565 6296.628238 4200.417887 -2306.285611 -942.555665 481 482 483 484 485 -215.469514 -6189.360754 -7276.285611 2679.117839 2272.847786 486 487 488 489 490 4240.818642 7318.329040 627.522140 213.221241 -2264.178663 491 492 493 494 495 -10613.454958 1016.646340 2693.230438 3257.010890 -5977.959115 496 497 498 499 500 -2983.302862 3596.697138 -982.878615 2481.625542 -1802.989110 501 502 503 504 505 9740.394394 1035.821337 2272.847786 4493.714389 -4780.974555 506 507 508 509 510 -270.481407 8448.711693 -2415.478711 -2456.285611 1370.325494 511 512 513 514 515 -3009.564011 7634.608291 802.524836 8541.311790 1646.120794 516 517 518 519 520 -2537.968312 -3300.974555 -3530.974555 641.547737 -1872.979912 521 522 523 524 525 1973.239343 -4011.288307 -4177.498657 -3249.988259 926.628238 526 527 528 529 530 -168.350965 3159.117839 3303.621995 2624.521289 -2708.775213 531 532 533 534 535 1128.735186 -3068.775213 -2050.882161 459.518593 -10485.077956 536 537 538 539 540 -1059.564011 5729.518593 9203.654686 3776.222092 -3570.167654 541 542 543 544 545 -8969.181358 -4941.264814 5568.218545 239.993639 491.607440 546 547 548 549 550 -828.798706 774.225639 3937.835893 -3460.481407 -170.882161 551 552 553 554 555 -864.339956 -1789.157865 -11213.030710 57.411645 -4825.884857 556 557 558 559 560 8262.017392 -5114.602911 966.203990 2616.628238 9770.011741 561 562 563 564 565 506.535844 -4694.985563 4471.607440 3355.014437 3107.522140 566 567 568 569 570 -266.378005 1575.014437 -142.588355 -3648.443358 -9294.584808 571 572 573 574 575 1700.316589 -4062.971008 -2759.199460 -3743.464156 -7397.567558 576 577 578 579 580 7534.115143 225.821337 6982.414340 5313.331736 5512.031688 581 582 583 584 585 59.518593 -4099.674506 1159.117839 -4635.478711 3490.011741 586 587 588 589 590 -2184.561315 5949.440789 -2466.285611 10.417887 1251.242889 591 592 593 594 595 -1641.772257 3899.058136 -430.075261 -2005.571104 4029.542086 596 597 598 599 600 -2968.775213 -1388.383362 4073.248540 11606.162389 -2864.764204 601 602 603 604 605 1337.826988 -7966.691757 3307.946388 5287.462438 -6809.923171 606 607 608 609 610 6747.181083 -567.493272 -488.816813 3412.967185 5308.753288 611 612 613 614 615 5467.946388 -5451.647467 -6277.534867 -6985.921066 -2307.493272 616 617 618 619 620 -4024.519720 -1680.347418 -6195.902964 -1633.814118 -976.727967 621 622 623 624 625 19423.349838 -707.457062 2701.284485 5337.563735 854.156738 626 627 628 629 630 -5006.709865 3904.180231 164.249132 12376.756835 203.272033 631 632 633 634 635 -2843.745218 -3558.410668 -7856.157014 2319.076238 3171.607434 636 637 638 639 640 -1906.244016 1280.376286 -1277.032815 6887.070587 -4714.943968 641 642 643 644 645 -7351.324517 2249.993634 -4135.114166 2784.649886 15502.967185 646 647 648 649 650 7249.569386 6356.586637 573.672788 -1187.502470 -2713.814118 651 -2767.857817 > postscript(file="/var/www/html/rcomp/tmp/6xy2i1292085873.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 = 651 Frequency = 1 lag(myerror, k = 1) myerror 0 -3508.885708 NA 1 -13955.834059 -3508.885708 2 6.900395 -13955.834059 3 -2485.257352 6.900395 4 -3216.871153 -2485.257352 5 -658.009907 -3216.871153 6 -406.820354 -658.009907 7 -3396.013454 -406.820354 8 -2097.203007 -3396.013454 9 -4245.096058 -2097.203007 10 1550.196896 -4245.096058 11 -1068.171201 1550.196896 12 -37.203007 -1068.171201 13 -4415.995352 -37.203007 14 60.214998 -4415.995352 15 2599.721851 60.214998 16 2420.911403 2599.721851 17 4310.196896 2420.911403 18 -697.678053 4310.196896 19 -1253.888403 -697.678053 20 4304.300298 -1253.888403 21 -1237.313502 4304.300298 22 -1194.312651 -1237.313502 23 1087.886691 -1194.312651 24 1035.083706 1087.886691 25 -1763.081503 1035.083706 26 3200.104503 -1763.081503 27 6807.393543 3200.104503 28 -2665.096058 6807.393543 29 3124.300298 -2665.096058 30 -1395.502204 3124.300298 31 3737.393543 -1395.502204 32 660.690044 3737.393543 33 3348.089948 660.690044 34 -217.203007 3348.089948 35 3032.321947 -217.203007 36 1266.425349 3032.321947 37 -6024.381551 1266.425349 38 -1454.381551 -6024.381551 39 7684.903942 -1454.381551 40 654.811548 7684.903942 41 4435.507953 654.811548 42 -606.396107 4435.507953 43 2708.896848 -606.396107 44 171.114292 2708.896848 45 -1505.681600 171.114292 46 -2374.782306 -1505.681600 47 -1598.595449 -2374.782306 48 1724.004648 -1598.595449 49 2743.894153 1724.004648 50 2254.097041 2743.894153 51 8266.111597 2254.097041 52 -2413.574651 8266.111597 53 -733.191998 -2413.574651 54 3533.128847 -733.191998 55 3418.693591 3533.128847 56 1869.883144 3418.693591 57 5641.386449 1869.883144 58 -3961.799557 5641.386449 59 -3340.591902 -3961.799557 60 1817.393543 -3340.591902 61 2391.003797 1817.393543 62 6746.789899 2391.003797 63 -4145.188452 6746.789899 64 -3704.805799 -4145.188452 65 -2502.274603 -3704.805799 66 -1908.503055 -2502.274603 67 2588.089948 -1908.503055 68 909.500491 2588.089948 69 -1874.381551 909.500491 70 4368.200443 -1874.381551 71 1906.918497 4368.200443 72 352.321947 1906.918497 73 -4343.081503 352.321947 74 -1906.488500 -4343.081503 75 3204.701053 -1906.488500 76 1195.507953 3204.701053 77 3948.287445 1195.507953 78 2720.196896 3948.287445 79 1866.093494 2720.196896 80 2480.196896 1866.093494 81 -5536.802252 2480.196896 82 -1503.574651 -5536.802252 83 -1099.309956 -1503.574651 84 -4628.816808 -1099.309956 85 -3164.874699 -4628.816808 86 -703.888403 -3164.874699 87 424.701053 -703.888403 88 -6033.888403 424.701053 89 4574.097041 -6033.888403 90 -730.098754 4574.097041 91 -5737.295400 -730.098754 92 -842.606457 -5737.295400 93 887.393543 -842.606457 94 -4308.484953 887.393543 95 -5682.274603 -4308.484953 96 -5643.888403 -5682.274603 97 -1951.085050 -5643.888403 98 -4213.574651 -1951.085050 99 -1869.378856 -4213.574651 100 -1255.681600 -1869.378856 101 -2439.291854 -1255.681600 102 1049.500491 -2439.291854 103 666.001101 1049.500491 104 -5.188452 666.001101 105 -918.595449 -5.188452 106 -3874.713406 -918.595449 107 2938.514195 -3874.713406 108 3397.393543 2938.514195 109 2122.303845 3397.393543 110 -1007.295400 2122.303845 111 352.704600 -1007.295400 112 -412.292705 352.704600 113 -608.009907 -412.292705 114 -213.099605 -608.009907 115 4315.107199 -213.099605 116 -1196.396107 4315.107199 117 3418.914950 -1196.396107 118 2033.128847 3418.914950 119 1076.900395 2033.128847 120 -1557.696155 1076.900395 121 -2667.788548 -1557.696155 122 -594.381551 -2667.788548 123 -4997.203007 -594.381551 124 -3823.574651 -4997.203007 125 836.517742 -3823.574651 126 1973.110745 836.517742 127 561.990093 1973.110745 128 -2844.695304 561.990093 129 -8229.785002 -2844.695304 130 -642.767751 -8229.785002 131 -1118.885708 -642.767751 132 2271.828799 -1118.885708 133 697.614902 2271.828799 134 -2567.788548 697.614902 135 -2163.413358 -2567.788548 136 -6487.678053 -2163.413358 137 -2797.295400 -6487.678053 138 -60.499509 -2797.295400 139 313.018352 -60.499509 140 594.097041 313.018352 141 769.901246 594.097041 142 -2867.678053 769.901246 143 -777.788548 -2867.678053 144 1988.200443 -777.788548 145 602.483241 1988.200443 146 597.794297 602.483241 147 -444.220258 597.794297 148 -1207.678053 -444.220258 149 -4744.695304 -1207.678053 150 1229.500491 -4744.695304 151 306.900395 1229.500491 152 197.232249 306.900395 153 4439.228703 197.232249 154 -49.692608 4439.228703 155 -8312.113309 -49.692608 156 -558.484953 -8312.113309 157 663.511500 -558.484953 158 -4167.203007 663.511500 159 697.283047 -4167.203007 160 -1411.799557 697.283047 161 -1440.098754 -1411.799557 162 3133.018352 -1440.098754 163 4844.903942 3133.018352 164 1217.614902 4844.903942 165 2764.811548 1217.614902 166 2814.701053 2764.811548 167 -140.278149 2814.701053 168 3834.318400 -140.278149 169 4212.594104 3834.318400 170 2273.290141 4212.594104 171 -2616.802252 2273.290141 172 -2284.289158 -2616.802252 173 -1729.692608 -2284.289158 174 -505.096058 -1729.692608 175 5293.493398 -505.096058 176 -5554.289158 5293.493398 177 1082.142552 -5554.289158 178 3473.986546 1082.142552 179 2611.114292 3473.986546 180 -3968.484953 2611.114292 181 -3981.467702 -3968.484953 182 154.318400 -3981.467702 183 738.896848 154.318400 184 -3173.413358 738.896848 185 -123.888403 -3173.413358 186 -811.085050 -123.888403 187 -2877.355103 -811.085050 188 -6566.396107 -2877.355103 189 -976.871153 -6566.396107 190 1637.614902 -976.871153 191 2212.211452 1637.614902 192 1021.183192 2212.211452 193 -1660.098754 1021.183192 194 198.896848 -1660.098754 195 1445.194201 198.896848 196 890.597651 1445.194201 197 45.710842 890.597651 198 732.303845 45.710842 199 4486.918497 732.303845 200 -1542.606457 4486.918497 201 -792.316198 -1542.606457 202 1706.093494 -792.316198 203 449.389996 1706.093494 204 884.903942 449.389996 205 -1621.891950 884.903942 206 -3603.081503 -1621.891950 207 -1339.692608 -3603.081503 208 3861.183192 -1339.692608 209 2650.307392 3861.183192 210 2117.707295 2650.307392 211 -1295.902959 2117.707295 212 3022.193350 -1295.902959 213 -6181.306409 3022.193350 214 -371.509298 -6181.306409 215 -5069.623708 -371.509298 216 -2949.692608 -5069.623708 217 -6407.203007 -2949.692608 218 5445.489851 -6407.203007 219 -6873.888403 5445.489851 220 -694.782306 -6873.888403 221 1258.421802 -694.782306 222 3363.690896 1258.421802 223 -4660.591902 3363.690896 224 -1498.978101 -4660.591902 225 -3493.099605 -1498.978101 226 -64.399653 -3493.099605 227 -308.978101 -64.399653 228 -351.306409 -308.978101 229 5892.193350 -351.306409 230 -4701.891950 5892.193350 231 -617.203007 -4701.891950 232 1214.004648 -617.203007 233 3437.886691 1214.004648 234 1439.389996 3437.886691 235 656.425349 1439.389996 236 2811.003797 656.425349 237 47.232249 2811.003797 238 -2916.709859 47.232249 239 1382.796993 -2916.709859 240 -2017.696155 1382.796993 241 -5508.484953 -2017.696155 242 -220.591902 -5508.484953 243 -1963.906506 -220.591902 244 1758.693591 -1963.906506 245 -1476.396107 1758.693591 246 -8175.902959 -1476.396107 247 -14.399653 -8175.902959 248 1535.618449 -14.399653 249 -6884.695304 1535.618449 250 -1308.816808 -6884.695304 251 2019.901246 -1308.816808 252 -2723.888403 2019.901246 253 -1223.413358 -2723.888403 254 -793.099605 -1223.413358 255 4143.493398 -793.099605 256 2814.318400 4143.493398 257 765.618449 2814.318400 258 -703.482258 765.618449 259 1036.111597 -703.482258 260 1102.211452 1036.111597 261 2764.811548 1102.211452 262 5703.511500 2764.811548 263 -3519.692608 5703.511500 264 3627.707295 -3519.692608 265 4207.707295 3627.707295 266 4668.896848 4207.707295 267 2786.918497 4668.896848 268 -114.695304 2786.918497 269 -438.009907 -114.695304 270 -1579.785002 -438.009907 271 -6693.081503 -1579.785002 272 -1170.610004 -6693.081503 273 1553.511500 -1170.610004 274 5975.507953 1553.511500 275 -1845.096058 5975.507953 276 265.507953 -1845.096058 277 -249.309956 265.507953 278 -83.191998 -249.309956 279 3505.618449 -83.191998 280 3390.690044 3505.618449 281 -2302.606457 3390.690044 282 5156.789899 -2302.606457 283 3413.603893 5156.789899 284 790.104503 3413.603893 285 2790.621144 790.104503 286 1451.897699 2790.621144 287 5370.690044 1451.897699 288 -2896.709859 5370.690044 289 -4678.885708 -2896.709859 290 -1152.767751 -4678.885708 291 686.111597 -1152.767751 292 862.211452 686.111597 293 73.603893 862.211452 294 -3476.871153 73.603893 295 464.207905 -3476.871153 296 1298.421802 464.207905 297 3590.528751 1298.421802 298 3214.903942 3590.528751 299 -3898.978101 3214.903942 300 1224.793446 -3898.978101 301 -791.799557 1224.793446 302 5491.496945 -791.799557 303 8850.911403 5491.496945 304 287.393543 8850.911403 305 -4293.906506 287.393543 306 -6976.691757 -4293.906506 307 -6737.184905 -6976.691757 308 2128.260140 -6737.184905 309 -1259.268361 2128.260140 310 934.115143 -1259.268361 311 3872.031688 934.115143 312 -4729.360754 3872.031688 313 -5449.674506 -4729.360754 314 1057.010890 -5449.674506 315 -424.160561 1057.010890 316 5172.925590 -424.160561 317 7509.629089 5172.925590 318 6960.011741 7509.629089 319 -348.775213 6960.011741 320 5271.607440 -348.775213 321 7703.331736 5271.607440 322 -2356.378005 7703.331736 323 6402.907488 -2356.378005 324 -2130.974555 6402.907488 325 -1436.285611 -2130.974555 326 -37.991805 -1436.285611 327 -646.267509 -37.991805 328 -10035.045266 -646.267509 329 -6672.902108 -10035.045266 330 -957.991805 -6672.902108 331 2712.100588 -957.991805 332 1013.331736 2712.100588 333 -162.989110 1013.331736 334 -6071.781455 -162.989110 335 -894.178663 -6071.781455 336 -4480.974555 -894.178663 337 -8596.378005 -4480.974555 338 544.428896 -8596.378005 339 -2145.077956 544.428896 340 6304.115143 -2145.077956 341 -3200.882161 6304.115143 342 442.008195 -3200.882161 343 2609.924739 442.008195 344 998.310938 2609.924739 345 5888.329040 998.310938 346 3949.135941 5888.329040 347 2064.115143 3949.135941 348 1197.504038 2064.115143 349 2678.260140 1197.504038 350 2966.535844 2678.260140 351 2103.621995 2966.535844 352 1796.604744 2103.621995 353 4302.414340 1796.604744 354 3935.761635 4302.414340 355 -215.861364 3935.761635 356 213.221241 -215.861364 357 8995.839439 213.221241 358 -4446.378005 8995.839439 359 -635.077956 -4446.378005 360 4829.135941 -635.077956 361 7440.818642 4829.135941 362 5484.461586 7440.818642 363 -5122.588355 5484.461586 364 -6257.074410 -5122.588355 365 -8863.141206 -6257.074410 366 -555.451411 -8863.141206 367 -524.671811 -555.451411 368 11513.248540 -524.671811 369 358.711693 11513.248540 370 -3622.495962 358.711693 371 -4867.991805 -3622.495962 372 -7085.077956 -4867.991805 373 -730.882161 -7085.077956 374 -10.006361 -730.882161 375 -3681.288307 -10.006361 376 2891.717935 -3681.288307 377 -5159.259163 2891.717935 378 -3950.167654 -5159.259163 379 1339.993639 -3950.167654 380 -629.641816 1339.993639 381 1902.501343 -629.641816 382 4548.218545 1902.501343 383 4998.310938 4548.218545 384 -5835.086861 4998.310938 385 3623.714389 -5835.086861 386 2733.732491 3623.714389 387 -1945.861364 2733.732491 388 702.414340 -1945.861364 389 4154.922044 702.414340 390 2633.308243 4154.922044 391 -2677.567558 2633.308243 392 3547.010890 -2677.567558 393 -3031.375309 3547.010890 394 -5494.271056 -3031.375309 395 -536.778759 -5494.271056 396 -736.285611 -536.778759 397 794.225639 -736.285611 398 -4325.460609 794.225639 399 3198.320136 -4325.460609 400 6884.115143 3198.320136 401 -804.671811 6884.115143 402 1168.711693 -804.671811 403 -255.368216 1168.711693 404 -1248.775213 -255.368216 405 -2734.602911 -1248.775213 406 -2173.371762 -2734.602911 407 -1640.882161 -2173.371762 408 -5987.991805 -1640.882161 409 -3520.882161 -5987.991805 410 -2656.345314 -3520.882161 411 -6681.375309 -2656.345314 412 2805.797844 -6681.375309 413 925.328189 2805.797844 414 2194.138636 925.328189 415 -5592.588355 2194.138636 416 -3835.884857 -5592.588355 417 3585.397089 -3835.884857 418 12208.651990 3585.397089 419 5552.100588 12208.651990 420 -2492.588355 5552.100588 421 -6962.095207 -2492.588355 422 -1608.060706 -6962.095207 423 560.842135 -1608.060706 424 -794.271056 560.842135 425 -340.457914 -794.271056 426 6038.310938 -340.457914 427 1509.518593 6038.310938 428 -2492.989110 1509.518593 429 -7142.071714 -2492.989110 430 -4789.199460 -7142.071714 431 -2913.362565 -4789.199460 432 -5107.991805 -2913.362565 433 946.628238 -5107.991805 434 1741.607440 946.628238 435 6040.758939 1741.607440 436 715.014437 6040.758939 437 4424.913139 715.014437 438 10113.308243 4424.913139 439 1225.032539 10113.308243 440 -9125.953757 1225.032539 441 -5284.764204 -9125.953757 442 3121.123490 -5284.764204 443 -1644.202156 3121.123490 444 3094.428896 -1644.202156 445 -8952.164107 3094.428896 446 -38.688210 -8952.164107 447 2092.907488 -38.688210 448 14725.415192 2092.907488 449 6964.115143 14725.415192 450 8412.432442 6964.115143 451 -4077.092512 8412.432442 452 4388.633888 -4077.092512 453 8646.222092 4388.633888 454 1420.818642 8646.222092 455 -5246.267509 1420.818642 456 -1938.060706 -5246.267509 457 -1248.858408 -1938.060706 458 957.411645 -1248.858408 459 -2232.164107 957.411645 460 7506.715240 -2232.164107 461 267.817790 7506.715240 462 -5749.988259 267.817790 463 -6257.991805 -5749.988259 464 -4263.371762 -6257.991805 465 -4605.409811 -4263.371762 466 183.714389 -4605.409811 467 -7.746953 183.714389 468 -1270.472209 -7.746953 469 -12926.982016 -1270.472209 470 -4487.991805 -12926.982016 471 6162.414340 -4487.991805 472 5107.435138 6162.414340 473 -4930.974555 5107.435138 474 -6310.167654 -4930.974555 475 -11613.362565 -6310.167654 476 6296.628238 -11613.362565 477 4200.417887 6296.628238 478 -2306.285611 4200.417887 479 -942.555665 -2306.285611 480 -215.469514 -942.555665 481 -6189.360754 -215.469514 482 -7276.285611 -6189.360754 483 2679.117839 -7276.285611 484 2272.847786 2679.117839 485 4240.818642 2272.847786 486 7318.329040 4240.818642 487 627.522140 7318.329040 488 213.221241 627.522140 489 -2264.178663 213.221241 490 -10613.454958 -2264.178663 491 1016.646340 -10613.454958 492 2693.230438 1016.646340 493 3257.010890 2693.230438 494 -5977.959115 3257.010890 495 -2983.302862 -5977.959115 496 3596.697138 -2983.302862 497 -982.878615 3596.697138 498 2481.625542 -982.878615 499 -1802.989110 2481.625542 500 9740.394394 -1802.989110 501 1035.821337 9740.394394 502 2272.847786 1035.821337 503 4493.714389 2272.847786 504 -4780.974555 4493.714389 505 -270.481407 -4780.974555 506 8448.711693 -270.481407 507 -2415.478711 8448.711693 508 -2456.285611 -2415.478711 509 1370.325494 -2456.285611 510 -3009.564011 1370.325494 511 7634.608291 -3009.564011 512 802.524836 7634.608291 513 8541.311790 802.524836 514 1646.120794 8541.311790 515 -2537.968312 1646.120794 516 -3300.974555 -2537.968312 517 -3530.974555 -3300.974555 518 641.547737 -3530.974555 519 -1872.979912 641.547737 520 1973.239343 -1872.979912 521 -4011.288307 1973.239343 522 -4177.498657 -4011.288307 523 -3249.988259 -4177.498657 524 926.628238 -3249.988259 525 -168.350965 926.628238 526 3159.117839 -168.350965 527 3303.621995 3159.117839 528 2624.521289 3303.621995 529 -2708.775213 2624.521289 530 1128.735186 -2708.775213 531 -3068.775213 1128.735186 532 -2050.882161 -3068.775213 533 459.518593 -2050.882161 534 -10485.077956 459.518593 535 -1059.564011 -10485.077956 536 5729.518593 -1059.564011 537 9203.654686 5729.518593 538 3776.222092 9203.654686 539 -3570.167654 3776.222092 540 -8969.181358 -3570.167654 541 -4941.264814 -8969.181358 542 5568.218545 -4941.264814 543 239.993639 5568.218545 544 491.607440 239.993639 545 -828.798706 491.607440 546 774.225639 -828.798706 547 3937.835893 774.225639 548 -3460.481407 3937.835893 549 -170.882161 -3460.481407 550 -864.339956 -170.882161 551 -1789.157865 -864.339956 552 -11213.030710 -1789.157865 553 57.411645 -11213.030710 554 -4825.884857 57.411645 555 8262.017392 -4825.884857 556 -5114.602911 8262.017392 557 966.203990 -5114.602911 558 2616.628238 966.203990 559 9770.011741 2616.628238 560 506.535844 9770.011741 561 -4694.985563 506.535844 562 4471.607440 -4694.985563 563 3355.014437 4471.607440 564 3107.522140 3355.014437 565 -266.378005 3107.522140 566 1575.014437 -266.378005 567 -142.588355 1575.014437 568 -3648.443358 -142.588355 569 -9294.584808 -3648.443358 570 1700.316589 -9294.584808 571 -4062.971008 1700.316589 572 -2759.199460 -4062.971008 573 -3743.464156 -2759.199460 574 -7397.567558 -3743.464156 575 7534.115143 -7397.567558 576 225.821337 7534.115143 577 6982.414340 225.821337 578 5313.331736 6982.414340 579 5512.031688 5313.331736 580 59.518593 5512.031688 581 -4099.674506 59.518593 582 1159.117839 -4099.674506 583 -4635.478711 1159.117839 584 3490.011741 -4635.478711 585 -2184.561315 3490.011741 586 5949.440789 -2184.561315 587 -2466.285611 5949.440789 588 10.417887 -2466.285611 589 1251.242889 10.417887 590 -1641.772257 1251.242889 591 3899.058136 -1641.772257 592 -430.075261 3899.058136 593 -2005.571104 -430.075261 594 4029.542086 -2005.571104 595 -2968.775213 4029.542086 596 -1388.383362 -2968.775213 597 4073.248540 -1388.383362 598 11606.162389 4073.248540 599 -2864.764204 11606.162389 600 1337.826988 -2864.764204 601 -7966.691757 1337.826988 602 3307.946388 -7966.691757 603 5287.462438 3307.946388 604 -6809.923171 5287.462438 605 6747.181083 -6809.923171 606 -567.493272 6747.181083 607 -488.816813 -567.493272 608 3412.967185 -488.816813 609 5308.753288 3412.967185 610 5467.946388 5308.753288 611 -5451.647467 5467.946388 612 -6277.534867 -5451.647467 613 -6985.921066 -6277.534867 614 -2307.493272 -6985.921066 615 -4024.519720 -2307.493272 616 -1680.347418 -4024.519720 617 -6195.902964 -1680.347418 618 -1633.814118 -6195.902964 619 -976.727967 -1633.814118 620 19423.349838 -976.727967 621 -707.457062 19423.349838 622 2701.284485 -707.457062 623 5337.563735 2701.284485 624 854.156738 5337.563735 625 -5006.709865 854.156738 626 3904.180231 -5006.709865 627 164.249132 3904.180231 628 12376.756835 164.249132 629 203.272033 12376.756835 630 -2843.745218 203.272033 631 -3558.410668 -2843.745218 632 -7856.157014 -3558.410668 633 2319.076238 -7856.157014 634 3171.607434 2319.076238 635 -1906.244016 3171.607434 636 1280.376286 -1906.244016 637 -1277.032815 1280.376286 638 6887.070587 -1277.032815 639 -4714.943968 6887.070587 640 -7351.324517 -4714.943968 641 2249.993634 -7351.324517 642 -4135.114166 2249.993634 643 2784.649886 -4135.114166 644 15502.967185 2784.649886 645 7249.569386 15502.967185 646 6356.586637 7249.569386 647 573.672788 6356.586637 648 -1187.502470 573.672788 649 -2713.814118 -1187.502470 650 -2767.857817 -2713.814118 651 NA -2767.857817 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -13955.834059 -3508.885708 [2,] 6.900395 -13955.834059 [3,] -2485.257352 6.900395 [4,] -3216.871153 -2485.257352 [5,] -658.009907 -3216.871153 [6,] -406.820354 -658.009907 [7,] -3396.013454 -406.820354 [8,] -2097.203007 -3396.013454 [9,] -4245.096058 -2097.203007 [10,] 1550.196896 -4245.096058 [11,] -1068.171201 1550.196896 [12,] -37.203007 -1068.171201 [13,] -4415.995352 -37.203007 [14,] 60.214998 -4415.995352 [15,] 2599.721851 60.214998 [16,] 2420.911403 2599.721851 [17,] 4310.196896 2420.911403 [18,] -697.678053 4310.196896 [19,] -1253.888403 -697.678053 [20,] 4304.300298 -1253.888403 [21,] -1237.313502 4304.300298 [22,] -1194.312651 -1237.313502 [23,] 1087.886691 -1194.312651 [24,] 1035.083706 1087.886691 [25,] -1763.081503 1035.083706 [26,] 3200.104503 -1763.081503 [27,] 6807.393543 3200.104503 [28,] -2665.096058 6807.393543 [29,] 3124.300298 -2665.096058 [30,] -1395.502204 3124.300298 [31,] 3737.393543 -1395.502204 [32,] 660.690044 3737.393543 [33,] 3348.089948 660.690044 [34,] -217.203007 3348.089948 [35,] 3032.321947 -217.203007 [36,] 1266.425349 3032.321947 [37,] -6024.381551 1266.425349 [38,] -1454.381551 -6024.381551 [39,] 7684.903942 -1454.381551 [40,] 654.811548 7684.903942 [41,] 4435.507953 654.811548 [42,] -606.396107 4435.507953 [43,] 2708.896848 -606.396107 [44,] 171.114292 2708.896848 [45,] -1505.681600 171.114292 [46,] -2374.782306 -1505.681600 [47,] -1598.595449 -2374.782306 [48,] 1724.004648 -1598.595449 [49,] 2743.894153 1724.004648 [50,] 2254.097041 2743.894153 [51,] 8266.111597 2254.097041 [52,] -2413.574651 8266.111597 [53,] -733.191998 -2413.574651 [54,] 3533.128847 -733.191998 [55,] 3418.693591 3533.128847 [56,] 1869.883144 3418.693591 [57,] 5641.386449 1869.883144 [58,] -3961.799557 5641.386449 [59,] -3340.591902 -3961.799557 [60,] 1817.393543 -3340.591902 [61,] 2391.003797 1817.393543 [62,] 6746.789899 2391.003797 [63,] -4145.188452 6746.789899 [64,] -3704.805799 -4145.188452 [65,] -2502.274603 -3704.805799 [66,] -1908.503055 -2502.274603 [67,] 2588.089948 -1908.503055 [68,] 909.500491 2588.089948 [69,] -1874.381551 909.500491 [70,] 4368.200443 -1874.381551 [71,] 1906.918497 4368.200443 [72,] 352.321947 1906.918497 [73,] -4343.081503 352.321947 [74,] -1906.488500 -4343.081503 [75,] 3204.701053 -1906.488500 [76,] 1195.507953 3204.701053 [77,] 3948.287445 1195.507953 [78,] 2720.196896 3948.287445 [79,] 1866.093494 2720.196896 [80,] 2480.196896 1866.093494 [81,] -5536.802252 2480.196896 [82,] -1503.574651 -5536.802252 [83,] -1099.309956 -1503.574651 [84,] -4628.816808 -1099.309956 [85,] -3164.874699 -4628.816808 [86,] -703.888403 -3164.874699 [87,] 424.701053 -703.888403 [88,] -6033.888403 424.701053 [89,] 4574.097041 -6033.888403 [90,] -730.098754 4574.097041 [91,] -5737.295400 -730.098754 [92,] -842.606457 -5737.295400 [93,] 887.393543 -842.606457 [94,] -4308.484953 887.393543 [95,] -5682.274603 -4308.484953 [96,] -5643.888403 -5682.274603 [97,] -1951.085050 -5643.888403 [98,] -4213.574651 -1951.085050 [99,] -1869.378856 -4213.574651 [100,] -1255.681600 -1869.378856 [101,] -2439.291854 -1255.681600 [102,] 1049.500491 -2439.291854 [103,] 666.001101 1049.500491 [104,] -5.188452 666.001101 [105,] -918.595449 -5.188452 [106,] -3874.713406 -918.595449 [107,] 2938.514195 -3874.713406 [108,] 3397.393543 2938.514195 [109,] 2122.303845 3397.393543 [110,] -1007.295400 2122.303845 [111,] 352.704600 -1007.295400 [112,] -412.292705 352.704600 [113,] -608.009907 -412.292705 [114,] -213.099605 -608.009907 [115,] 4315.107199 -213.099605 [116,] -1196.396107 4315.107199 [117,] 3418.914950 -1196.396107 [118,] 2033.128847 3418.914950 [119,] 1076.900395 2033.128847 [120,] -1557.696155 1076.900395 [121,] -2667.788548 -1557.696155 [122,] -594.381551 -2667.788548 [123,] -4997.203007 -594.381551 [124,] -3823.574651 -4997.203007 [125,] 836.517742 -3823.574651 [126,] 1973.110745 836.517742 [127,] 561.990093 1973.110745 [128,] -2844.695304 561.990093 [129,] -8229.785002 -2844.695304 [130,] -642.767751 -8229.785002 [131,] -1118.885708 -642.767751 [132,] 2271.828799 -1118.885708 [133,] 697.614902 2271.828799 [134,] -2567.788548 697.614902 [135,] -2163.413358 -2567.788548 [136,] -6487.678053 -2163.413358 [137,] -2797.295400 -6487.678053 [138,] -60.499509 -2797.295400 [139,] 313.018352 -60.499509 [140,] 594.097041 313.018352 [141,] 769.901246 594.097041 [142,] -2867.678053 769.901246 [143,] -777.788548 -2867.678053 [144,] 1988.200443 -777.788548 [145,] 602.483241 1988.200443 [146,] 597.794297 602.483241 [147,] -444.220258 597.794297 [148,] -1207.678053 -444.220258 [149,] -4744.695304 -1207.678053 [150,] 1229.500491 -4744.695304 [151,] 306.900395 1229.500491 [152,] 197.232249 306.900395 [153,] 4439.228703 197.232249 [154,] -49.692608 4439.228703 [155,] -8312.113309 -49.692608 [156,] -558.484953 -8312.113309 [157,] 663.511500 -558.484953 [158,] -4167.203007 663.511500 [159,] 697.283047 -4167.203007 [160,] -1411.799557 697.283047 [161,] -1440.098754 -1411.799557 [162,] 3133.018352 -1440.098754 [163,] 4844.903942 3133.018352 [164,] 1217.614902 4844.903942 [165,] 2764.811548 1217.614902 [166,] 2814.701053 2764.811548 [167,] -140.278149 2814.701053 [168,] 3834.318400 -140.278149 [169,] 4212.594104 3834.318400 [170,] 2273.290141 4212.594104 [171,] -2616.802252 2273.290141 [172,] -2284.289158 -2616.802252 [173,] -1729.692608 -2284.289158 [174,] -505.096058 -1729.692608 [175,] 5293.493398 -505.096058 [176,] -5554.289158 5293.493398 [177,] 1082.142552 -5554.289158 [178,] 3473.986546 1082.142552 [179,] 2611.114292 3473.986546 [180,] -3968.484953 2611.114292 [181,] -3981.467702 -3968.484953 [182,] 154.318400 -3981.467702 [183,] 738.896848 154.318400 [184,] -3173.413358 738.896848 [185,] -123.888403 -3173.413358 [186,] -811.085050 -123.888403 [187,] -2877.355103 -811.085050 [188,] -6566.396107 -2877.355103 [189,] -976.871153 -6566.396107 [190,] 1637.614902 -976.871153 [191,] 2212.211452 1637.614902 [192,] 1021.183192 2212.211452 [193,] -1660.098754 1021.183192 [194,] 198.896848 -1660.098754 [195,] 1445.194201 198.896848 [196,] 890.597651 1445.194201 [197,] 45.710842 890.597651 [198,] 732.303845 45.710842 [199,] 4486.918497 732.303845 [200,] -1542.606457 4486.918497 [201,] -792.316198 -1542.606457 [202,] 1706.093494 -792.316198 [203,] 449.389996 1706.093494 [204,] 884.903942 449.389996 [205,] -1621.891950 884.903942 [206,] -3603.081503 -1621.891950 [207,] -1339.692608 -3603.081503 [208,] 3861.183192 -1339.692608 [209,] 2650.307392 3861.183192 [210,] 2117.707295 2650.307392 [211,] -1295.902959 2117.707295 [212,] 3022.193350 -1295.902959 [213,] -6181.306409 3022.193350 [214,] -371.509298 -6181.306409 [215,] -5069.623708 -371.509298 [216,] -2949.692608 -5069.623708 [217,] -6407.203007 -2949.692608 [218,] 5445.489851 -6407.203007 [219,] -6873.888403 5445.489851 [220,] -694.782306 -6873.888403 [221,] 1258.421802 -694.782306 [222,] 3363.690896 1258.421802 [223,] -4660.591902 3363.690896 [224,] -1498.978101 -4660.591902 [225,] -3493.099605 -1498.978101 [226,] -64.399653 -3493.099605 [227,] -308.978101 -64.399653 [228,] -351.306409 -308.978101 [229,] 5892.193350 -351.306409 [230,] -4701.891950 5892.193350 [231,] -617.203007 -4701.891950 [232,] 1214.004648 -617.203007 [233,] 3437.886691 1214.004648 [234,] 1439.389996 3437.886691 [235,] 656.425349 1439.389996 [236,] 2811.003797 656.425349 [237,] 47.232249 2811.003797 [238,] -2916.709859 47.232249 [239,] 1382.796993 -2916.709859 [240,] -2017.696155 1382.796993 [241,] -5508.484953 -2017.696155 [242,] -220.591902 -5508.484953 [243,] -1963.906506 -220.591902 [244,] 1758.693591 -1963.906506 [245,] -1476.396107 1758.693591 [246,] -8175.902959 -1476.396107 [247,] -14.399653 -8175.902959 [248,] 1535.618449 -14.399653 [249,] -6884.695304 1535.618449 [250,] -1308.816808 -6884.695304 [251,] 2019.901246 -1308.816808 [252,] -2723.888403 2019.901246 [253,] -1223.413358 -2723.888403 [254,] -793.099605 -1223.413358 [255,] 4143.493398 -793.099605 [256,] 2814.318400 4143.493398 [257,] 765.618449 2814.318400 [258,] -703.482258 765.618449 [259,] 1036.111597 -703.482258 [260,] 1102.211452 1036.111597 [261,] 2764.811548 1102.211452 [262,] 5703.511500 2764.811548 [263,] -3519.692608 5703.511500 [264,] 3627.707295 -3519.692608 [265,] 4207.707295 3627.707295 [266,] 4668.896848 4207.707295 [267,] 2786.918497 4668.896848 [268,] -114.695304 2786.918497 [269,] -438.009907 -114.695304 [270,] -1579.785002 -438.009907 [271,] -6693.081503 -1579.785002 [272,] -1170.610004 -6693.081503 [273,] 1553.511500 -1170.610004 [274,] 5975.507953 1553.511500 [275,] -1845.096058 5975.507953 [276,] 265.507953 -1845.096058 [277,] -249.309956 265.507953 [278,] -83.191998 -249.309956 [279,] 3505.618449 -83.191998 [280,] 3390.690044 3505.618449 [281,] -2302.606457 3390.690044 [282,] 5156.789899 -2302.606457 [283,] 3413.603893 5156.789899 [284,] 790.104503 3413.603893 [285,] 2790.621144 790.104503 [286,] 1451.897699 2790.621144 [287,] 5370.690044 1451.897699 [288,] -2896.709859 5370.690044 [289,] -4678.885708 -2896.709859 [290,] -1152.767751 -4678.885708 [291,] 686.111597 -1152.767751 [292,] 862.211452 686.111597 [293,] 73.603893 862.211452 [294,] -3476.871153 73.603893 [295,] 464.207905 -3476.871153 [296,] 1298.421802 464.207905 [297,] 3590.528751 1298.421802 [298,] 3214.903942 3590.528751 [299,] -3898.978101 3214.903942 [300,] 1224.793446 -3898.978101 [301,] -791.799557 1224.793446 [302,] 5491.496945 -791.799557 [303,] 8850.911403 5491.496945 [304,] 287.393543 8850.911403 [305,] -4293.906506 287.393543 [306,] -6976.691757 -4293.906506 [307,] -6737.184905 -6976.691757 [308,] 2128.260140 -6737.184905 [309,] -1259.268361 2128.260140 [310,] 934.115143 -1259.268361 [311,] 3872.031688 934.115143 [312,] -4729.360754 3872.031688 [313,] -5449.674506 -4729.360754 [314,] 1057.010890 -5449.674506 [315,] -424.160561 1057.010890 [316,] 5172.925590 -424.160561 [317,] 7509.629089 5172.925590 [318,] 6960.011741 7509.629089 [319,] -348.775213 6960.011741 [320,] 5271.607440 -348.775213 [321,] 7703.331736 5271.607440 [322,] -2356.378005 7703.331736 [323,] 6402.907488 -2356.378005 [324,] -2130.974555 6402.907488 [325,] -1436.285611 -2130.974555 [326,] -37.991805 -1436.285611 [327,] -646.267509 -37.991805 [328,] -10035.045266 -646.267509 [329,] -6672.902108 -10035.045266 [330,] -957.991805 -6672.902108 [331,] 2712.100588 -957.991805 [332,] 1013.331736 2712.100588 [333,] -162.989110 1013.331736 [334,] -6071.781455 -162.989110 [335,] -894.178663 -6071.781455 [336,] -4480.974555 -894.178663 [337,] -8596.378005 -4480.974555 [338,] 544.428896 -8596.378005 [339,] -2145.077956 544.428896 [340,] 6304.115143 -2145.077956 [341,] -3200.882161 6304.115143 [342,] 442.008195 -3200.882161 [343,] 2609.924739 442.008195 [344,] 998.310938 2609.924739 [345,] 5888.329040 998.310938 [346,] 3949.135941 5888.329040 [347,] 2064.115143 3949.135941 [348,] 1197.504038 2064.115143 [349,] 2678.260140 1197.504038 [350,] 2966.535844 2678.260140 [351,] 2103.621995 2966.535844 [352,] 1796.604744 2103.621995 [353,] 4302.414340 1796.604744 [354,] 3935.761635 4302.414340 [355,] -215.861364 3935.761635 [356,] 213.221241 -215.861364 [357,] 8995.839439 213.221241 [358,] -4446.378005 8995.839439 [359,] -635.077956 -4446.378005 [360,] 4829.135941 -635.077956 [361,] 7440.818642 4829.135941 [362,] 5484.461586 7440.818642 [363,] -5122.588355 5484.461586 [364,] -6257.074410 -5122.588355 [365,] -8863.141206 -6257.074410 [366,] -555.451411 -8863.141206 [367,] -524.671811 -555.451411 [368,] 11513.248540 -524.671811 [369,] 358.711693 11513.248540 [370,] -3622.495962 358.711693 [371,] -4867.991805 -3622.495962 [372,] -7085.077956 -4867.991805 [373,] -730.882161 -7085.077956 [374,] -10.006361 -730.882161 [375,] -3681.288307 -10.006361 [376,] 2891.717935 -3681.288307 [377,] -5159.259163 2891.717935 [378,] -3950.167654 -5159.259163 [379,] 1339.993639 -3950.167654 [380,] -629.641816 1339.993639 [381,] 1902.501343 -629.641816 [382,] 4548.218545 1902.501343 [383,] 4998.310938 4548.218545 [384,] -5835.086861 4998.310938 [385,] 3623.714389 -5835.086861 [386,] 2733.732491 3623.714389 [387,] -1945.861364 2733.732491 [388,] 702.414340 -1945.861364 [389,] 4154.922044 702.414340 [390,] 2633.308243 4154.922044 [391,] -2677.567558 2633.308243 [392,] 3547.010890 -2677.567558 [393,] -3031.375309 3547.010890 [394,] -5494.271056 -3031.375309 [395,] -536.778759 -5494.271056 [396,] -736.285611 -536.778759 [397,] 794.225639 -736.285611 [398,] -4325.460609 794.225639 [399,] 3198.320136 -4325.460609 [400,] 6884.115143 3198.320136 [401,] -804.671811 6884.115143 [402,] 1168.711693 -804.671811 [403,] -255.368216 1168.711693 [404,] -1248.775213 -255.368216 [405,] -2734.602911 -1248.775213 [406,] -2173.371762 -2734.602911 [407,] -1640.882161 -2173.371762 [408,] -5987.991805 -1640.882161 [409,] -3520.882161 -5987.991805 [410,] -2656.345314 -3520.882161 [411,] -6681.375309 -2656.345314 [412,] 2805.797844 -6681.375309 [413,] 925.328189 2805.797844 [414,] 2194.138636 925.328189 [415,] -5592.588355 2194.138636 [416,] -3835.884857 -5592.588355 [417,] 3585.397089 -3835.884857 [418,] 12208.651990 3585.397089 [419,] 5552.100588 12208.651990 [420,] -2492.588355 5552.100588 [421,] -6962.095207 -2492.588355 [422,] -1608.060706 -6962.095207 [423,] 560.842135 -1608.060706 [424,] -794.271056 560.842135 [425,] -340.457914 -794.271056 [426,] 6038.310938 -340.457914 [427,] 1509.518593 6038.310938 [428,] -2492.989110 1509.518593 [429,] -7142.071714 -2492.989110 [430,] -4789.199460 -7142.071714 [431,] -2913.362565 -4789.199460 [432,] -5107.991805 -2913.362565 [433,] 946.628238 -5107.991805 [434,] 1741.607440 946.628238 [435,] 6040.758939 1741.607440 [436,] 715.014437 6040.758939 [437,] 4424.913139 715.014437 [438,] 10113.308243 4424.913139 [439,] 1225.032539 10113.308243 [440,] -9125.953757 1225.032539 [441,] -5284.764204 -9125.953757 [442,] 3121.123490 -5284.764204 [443,] -1644.202156 3121.123490 [444,] 3094.428896 -1644.202156 [445,] -8952.164107 3094.428896 [446,] -38.688210 -8952.164107 [447,] 2092.907488 -38.688210 [448,] 14725.415192 2092.907488 [449,] 6964.115143 14725.415192 [450,] 8412.432442 6964.115143 [451,] -4077.092512 8412.432442 [452,] 4388.633888 -4077.092512 [453,] 8646.222092 4388.633888 [454,] 1420.818642 8646.222092 [455,] -5246.267509 1420.818642 [456,] -1938.060706 -5246.267509 [457,] -1248.858408 -1938.060706 [458,] 957.411645 -1248.858408 [459,] -2232.164107 957.411645 [460,] 7506.715240 -2232.164107 [461,] 267.817790 7506.715240 [462,] -5749.988259 267.817790 [463,] -6257.991805 -5749.988259 [464,] -4263.371762 -6257.991805 [465,] -4605.409811 -4263.371762 [466,] 183.714389 -4605.409811 [467,] -7.746953 183.714389 [468,] -1270.472209 -7.746953 [469,] -12926.982016 -1270.472209 [470,] -4487.991805 -12926.982016 [471,] 6162.414340 -4487.991805 [472,] 5107.435138 6162.414340 [473,] -4930.974555 5107.435138 [474,] -6310.167654 -4930.974555 [475,] -11613.362565 -6310.167654 [476,] 6296.628238 -11613.362565 [477,] 4200.417887 6296.628238 [478,] -2306.285611 4200.417887 [479,] -942.555665 -2306.285611 [480,] -215.469514 -942.555665 [481,] -6189.360754 -215.469514 [482,] -7276.285611 -6189.360754 [483,] 2679.117839 -7276.285611 [484,] 2272.847786 2679.117839 [485,] 4240.818642 2272.847786 [486,] 7318.329040 4240.818642 [487,] 627.522140 7318.329040 [488,] 213.221241 627.522140 [489,] -2264.178663 213.221241 [490,] -10613.454958 -2264.178663 [491,] 1016.646340 -10613.454958 [492,] 2693.230438 1016.646340 [493,] 3257.010890 2693.230438 [494,] -5977.959115 3257.010890 [495,] -2983.302862 -5977.959115 [496,] 3596.697138 -2983.302862 [497,] -982.878615 3596.697138 [498,] 2481.625542 -982.878615 [499,] -1802.989110 2481.625542 [500,] 9740.394394 -1802.989110 [501,] 1035.821337 9740.394394 [502,] 2272.847786 1035.821337 [503,] 4493.714389 2272.847786 [504,] -4780.974555 4493.714389 [505,] -270.481407 -4780.974555 [506,] 8448.711693 -270.481407 [507,] -2415.478711 8448.711693 [508,] -2456.285611 -2415.478711 [509,] 1370.325494 -2456.285611 [510,] -3009.564011 1370.325494 [511,] 7634.608291 -3009.564011 [512,] 802.524836 7634.608291 [513,] 8541.311790 802.524836 [514,] 1646.120794 8541.311790 [515,] -2537.968312 1646.120794 [516,] -3300.974555 -2537.968312 [517,] -3530.974555 -3300.974555 [518,] 641.547737 -3530.974555 [519,] -1872.979912 641.547737 [520,] 1973.239343 -1872.979912 [521,] -4011.288307 1973.239343 [522,] -4177.498657 -4011.288307 [523,] -3249.988259 -4177.498657 [524,] 926.628238 -3249.988259 [525,] -168.350965 926.628238 [526,] 3159.117839 -168.350965 [527,] 3303.621995 3159.117839 [528,] 2624.521289 3303.621995 [529,] -2708.775213 2624.521289 [530,] 1128.735186 -2708.775213 [531,] -3068.775213 1128.735186 [532,] -2050.882161 -3068.775213 [533,] 459.518593 -2050.882161 [534,] -10485.077956 459.518593 [535,] -1059.564011 -10485.077956 [536,] 5729.518593 -1059.564011 [537,] 9203.654686 5729.518593 [538,] 3776.222092 9203.654686 [539,] -3570.167654 3776.222092 [540,] -8969.181358 -3570.167654 [541,] -4941.264814 -8969.181358 [542,] 5568.218545 -4941.264814 [543,] 239.993639 5568.218545 [544,] 491.607440 239.993639 [545,] -828.798706 491.607440 [546,] 774.225639 -828.798706 [547,] 3937.835893 774.225639 [548,] -3460.481407 3937.835893 [549,] -170.882161 -3460.481407 [550,] -864.339956 -170.882161 [551,] -1789.157865 -864.339956 [552,] -11213.030710 -1789.157865 [553,] 57.411645 -11213.030710 [554,] -4825.884857 57.411645 [555,] 8262.017392 -4825.884857 [556,] -5114.602911 8262.017392 [557,] 966.203990 -5114.602911 [558,] 2616.628238 966.203990 [559,] 9770.011741 2616.628238 [560,] 506.535844 9770.011741 [561,] -4694.985563 506.535844 [562,] 4471.607440 -4694.985563 [563,] 3355.014437 4471.607440 [564,] 3107.522140 3355.014437 [565,] -266.378005 3107.522140 [566,] 1575.014437 -266.378005 [567,] -142.588355 1575.014437 [568,] -3648.443358 -142.588355 [569,] -9294.584808 -3648.443358 [570,] 1700.316589 -9294.584808 [571,] -4062.971008 1700.316589 [572,] -2759.199460 -4062.971008 [573,] -3743.464156 -2759.199460 [574,] -7397.567558 -3743.464156 [575,] 7534.115143 -7397.567558 [576,] 225.821337 7534.115143 [577,] 6982.414340 225.821337 [578,] 5313.331736 6982.414340 [579,] 5512.031688 5313.331736 [580,] 59.518593 5512.031688 [581,] -4099.674506 59.518593 [582,] 1159.117839 -4099.674506 [583,] -4635.478711 1159.117839 [584,] 3490.011741 -4635.478711 [585,] -2184.561315 3490.011741 [586,] 5949.440789 -2184.561315 [587,] -2466.285611 5949.440789 [588,] 10.417887 -2466.285611 [589,] 1251.242889 10.417887 [590,] -1641.772257 1251.242889 [591,] 3899.058136 -1641.772257 [592,] -430.075261 3899.058136 [593,] -2005.571104 -430.075261 [594,] 4029.542086 -2005.571104 [595,] -2968.775213 4029.542086 [596,] -1388.383362 -2968.775213 [597,] 4073.248540 -1388.383362 [598,] 11606.162389 4073.248540 [599,] -2864.764204 11606.162389 [600,] 1337.826988 -2864.764204 [601,] -7966.691757 1337.826988 [602,] 3307.946388 -7966.691757 [603,] 5287.462438 3307.946388 [604,] -6809.923171 5287.462438 [605,] 6747.181083 -6809.923171 [606,] -567.493272 6747.181083 [607,] -488.816813 -567.493272 [608,] 3412.967185 -488.816813 [609,] 5308.753288 3412.967185 [610,] 5467.946388 5308.753288 [611,] -5451.647467 5467.946388 [612,] -6277.534867 -5451.647467 [613,] -6985.921066 -6277.534867 [614,] -2307.493272 -6985.921066 [615,] -4024.519720 -2307.493272 [616,] -1680.347418 -4024.519720 [617,] -6195.902964 -1680.347418 [618,] -1633.814118 -6195.902964 [619,] -976.727967 -1633.814118 [620,] 19423.349838 -976.727967 [621,] -707.457062 19423.349838 [622,] 2701.284485 -707.457062 [623,] 5337.563735 2701.284485 [624,] 854.156738 5337.563735 [625,] -5006.709865 854.156738 [626,] 3904.180231 -5006.709865 [627,] 164.249132 3904.180231 [628,] 12376.756835 164.249132 [629,] 203.272033 12376.756835 [630,] -2843.745218 203.272033 [631,] -3558.410668 -2843.745218 [632,] -7856.157014 -3558.410668 [633,] 2319.076238 -7856.157014 [634,] 3171.607434 2319.076238 [635,] -1906.244016 3171.607434 [636,] 1280.376286 -1906.244016 [637,] -1277.032815 1280.376286 [638,] 6887.070587 -1277.032815 [639,] -4714.943968 6887.070587 [640,] -7351.324517 -4714.943968 [641,] 2249.993634 -7351.324517 [642,] -4135.114166 2249.993634 [643,] 2784.649886 -4135.114166 [644,] 15502.967185 2784.649886 [645,] 7249.569386 15502.967185 [646,] 6356.586637 7249.569386 [647,] 573.672788 6356.586637 [648,] -1187.502470 573.672788 [649,] -2713.814118 -1187.502470 [650,] -2767.857817 -2713.814118 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -13955.834059 -3508.885708 2 6.900395 -13955.834059 3 -2485.257352 6.900395 4 -3216.871153 -2485.257352 5 -658.009907 -3216.871153 6 -406.820354 -658.009907 7 -3396.013454 -406.820354 8 -2097.203007 -3396.013454 9 -4245.096058 -2097.203007 10 1550.196896 -4245.096058 11 -1068.171201 1550.196896 12 -37.203007 -1068.171201 13 -4415.995352 -37.203007 14 60.214998 -4415.995352 15 2599.721851 60.214998 16 2420.911403 2599.721851 17 4310.196896 2420.911403 18 -697.678053 4310.196896 19 -1253.888403 -697.678053 20 4304.300298 -1253.888403 21 -1237.313502 4304.300298 22 -1194.312651 -1237.313502 23 1087.886691 -1194.312651 24 1035.083706 1087.886691 25 -1763.081503 1035.083706 26 3200.104503 -1763.081503 27 6807.393543 3200.104503 28 -2665.096058 6807.393543 29 3124.300298 -2665.096058 30 -1395.502204 3124.300298 31 3737.393543 -1395.502204 32 660.690044 3737.393543 33 3348.089948 660.690044 34 -217.203007 3348.089948 35 3032.321947 -217.203007 36 1266.425349 3032.321947 37 -6024.381551 1266.425349 38 -1454.381551 -6024.381551 39 7684.903942 -1454.381551 40 654.811548 7684.903942 41 4435.507953 654.811548 42 -606.396107 4435.507953 43 2708.896848 -606.396107 44 171.114292 2708.896848 45 -1505.681600 171.114292 46 -2374.782306 -1505.681600 47 -1598.595449 -2374.782306 48 1724.004648 -1598.595449 49 2743.894153 1724.004648 50 2254.097041 2743.894153 51 8266.111597 2254.097041 52 -2413.574651 8266.111597 53 -733.191998 -2413.574651 54 3533.128847 -733.191998 55 3418.693591 3533.128847 56 1869.883144 3418.693591 57 5641.386449 1869.883144 58 -3961.799557 5641.386449 59 -3340.591902 -3961.799557 60 1817.393543 -3340.591902 61 2391.003797 1817.393543 62 6746.789899 2391.003797 63 -4145.188452 6746.789899 64 -3704.805799 -4145.188452 65 -2502.274603 -3704.805799 66 -1908.503055 -2502.274603 67 2588.089948 -1908.503055 68 909.500491 2588.089948 69 -1874.381551 909.500491 70 4368.200443 -1874.381551 71 1906.918497 4368.200443 72 352.321947 1906.918497 73 -4343.081503 352.321947 74 -1906.488500 -4343.081503 75 3204.701053 -1906.488500 76 1195.507953 3204.701053 77 3948.287445 1195.507953 78 2720.196896 3948.287445 79 1866.093494 2720.196896 80 2480.196896 1866.093494 81 -5536.802252 2480.196896 82 -1503.574651 -5536.802252 83 -1099.309956 -1503.574651 84 -4628.816808 -1099.309956 85 -3164.874699 -4628.816808 86 -703.888403 -3164.874699 87 424.701053 -703.888403 88 -6033.888403 424.701053 89 4574.097041 -6033.888403 90 -730.098754 4574.097041 91 -5737.295400 -730.098754 92 -842.606457 -5737.295400 93 887.393543 -842.606457 94 -4308.484953 887.393543 95 -5682.274603 -4308.484953 96 -5643.888403 -5682.274603 97 -1951.085050 -5643.888403 98 -4213.574651 -1951.085050 99 -1869.378856 -4213.574651 100 -1255.681600 -1869.378856 101 -2439.291854 -1255.681600 102 1049.500491 -2439.291854 103 666.001101 1049.500491 104 -5.188452 666.001101 105 -918.595449 -5.188452 106 -3874.713406 -918.595449 107 2938.514195 -3874.713406 108 3397.393543 2938.514195 109 2122.303845 3397.393543 110 -1007.295400 2122.303845 111 352.704600 -1007.295400 112 -412.292705 352.704600 113 -608.009907 -412.292705 114 -213.099605 -608.009907 115 4315.107199 -213.099605 116 -1196.396107 4315.107199 117 3418.914950 -1196.396107 118 2033.128847 3418.914950 119 1076.900395 2033.128847 120 -1557.696155 1076.900395 121 -2667.788548 -1557.696155 122 -594.381551 -2667.788548 123 -4997.203007 -594.381551 124 -3823.574651 -4997.203007 125 836.517742 -3823.574651 126 1973.110745 836.517742 127 561.990093 1973.110745 128 -2844.695304 561.990093 129 -8229.785002 -2844.695304 130 -642.767751 -8229.785002 131 -1118.885708 -642.767751 132 2271.828799 -1118.885708 133 697.614902 2271.828799 134 -2567.788548 697.614902 135 -2163.413358 -2567.788548 136 -6487.678053 -2163.413358 137 -2797.295400 -6487.678053 138 -60.499509 -2797.295400 139 313.018352 -60.499509 140 594.097041 313.018352 141 769.901246 594.097041 142 -2867.678053 769.901246 143 -777.788548 -2867.678053 144 1988.200443 -777.788548 145 602.483241 1988.200443 146 597.794297 602.483241 147 -444.220258 597.794297 148 -1207.678053 -444.220258 149 -4744.695304 -1207.678053 150 1229.500491 -4744.695304 151 306.900395 1229.500491 152 197.232249 306.900395 153 4439.228703 197.232249 154 -49.692608 4439.228703 155 -8312.113309 -49.692608 156 -558.484953 -8312.113309 157 663.511500 -558.484953 158 -4167.203007 663.511500 159 697.283047 -4167.203007 160 -1411.799557 697.283047 161 -1440.098754 -1411.799557 162 3133.018352 -1440.098754 163 4844.903942 3133.018352 164 1217.614902 4844.903942 165 2764.811548 1217.614902 166 2814.701053 2764.811548 167 -140.278149 2814.701053 168 3834.318400 -140.278149 169 4212.594104 3834.318400 170 2273.290141 4212.594104 171 -2616.802252 2273.290141 172 -2284.289158 -2616.802252 173 -1729.692608 -2284.289158 174 -505.096058 -1729.692608 175 5293.493398 -505.096058 176 -5554.289158 5293.493398 177 1082.142552 -5554.289158 178 3473.986546 1082.142552 179 2611.114292 3473.986546 180 -3968.484953 2611.114292 181 -3981.467702 -3968.484953 182 154.318400 -3981.467702 183 738.896848 154.318400 184 -3173.413358 738.896848 185 -123.888403 -3173.413358 186 -811.085050 -123.888403 187 -2877.355103 -811.085050 188 -6566.396107 -2877.355103 189 -976.871153 -6566.396107 190 1637.614902 -976.871153 191 2212.211452 1637.614902 192 1021.183192 2212.211452 193 -1660.098754 1021.183192 194 198.896848 -1660.098754 195 1445.194201 198.896848 196 890.597651 1445.194201 197 45.710842 890.597651 198 732.303845 45.710842 199 4486.918497 732.303845 200 -1542.606457 4486.918497 201 -792.316198 -1542.606457 202 1706.093494 -792.316198 203 449.389996 1706.093494 204 884.903942 449.389996 205 -1621.891950 884.903942 206 -3603.081503 -1621.891950 207 -1339.692608 -3603.081503 208 3861.183192 -1339.692608 209 2650.307392 3861.183192 210 2117.707295 2650.307392 211 -1295.902959 2117.707295 212 3022.193350 -1295.902959 213 -6181.306409 3022.193350 214 -371.509298 -6181.306409 215 -5069.623708 -371.509298 216 -2949.692608 -5069.623708 217 -6407.203007 -2949.692608 218 5445.489851 -6407.203007 219 -6873.888403 5445.489851 220 -694.782306 -6873.888403 221 1258.421802 -694.782306 222 3363.690896 1258.421802 223 -4660.591902 3363.690896 224 -1498.978101 -4660.591902 225 -3493.099605 -1498.978101 226 -64.399653 -3493.099605 227 -308.978101 -64.399653 228 -351.306409 -308.978101 229 5892.193350 -351.306409 230 -4701.891950 5892.193350 231 -617.203007 -4701.891950 232 1214.004648 -617.203007 233 3437.886691 1214.004648 234 1439.389996 3437.886691 235 656.425349 1439.389996 236 2811.003797 656.425349 237 47.232249 2811.003797 238 -2916.709859 47.232249 239 1382.796993 -2916.709859 240 -2017.696155 1382.796993 241 -5508.484953 -2017.696155 242 -220.591902 -5508.484953 243 -1963.906506 -220.591902 244 1758.693591 -1963.906506 245 -1476.396107 1758.693591 246 -8175.902959 -1476.396107 247 -14.399653 -8175.902959 248 1535.618449 -14.399653 249 -6884.695304 1535.618449 250 -1308.816808 -6884.695304 251 2019.901246 -1308.816808 252 -2723.888403 2019.901246 253 -1223.413358 -2723.888403 254 -793.099605 -1223.413358 255 4143.493398 -793.099605 256 2814.318400 4143.493398 257 765.618449 2814.318400 258 -703.482258 765.618449 259 1036.111597 -703.482258 260 1102.211452 1036.111597 261 2764.811548 1102.211452 262 5703.511500 2764.811548 263 -3519.692608 5703.511500 264 3627.707295 -3519.692608 265 4207.707295 3627.707295 266 4668.896848 4207.707295 267 2786.918497 4668.896848 268 -114.695304 2786.918497 269 -438.009907 -114.695304 270 -1579.785002 -438.009907 271 -6693.081503 -1579.785002 272 -1170.610004 -6693.081503 273 1553.511500 -1170.610004 274 5975.507953 1553.511500 275 -1845.096058 5975.507953 276 265.507953 -1845.096058 277 -249.309956 265.507953 278 -83.191998 -249.309956 279 3505.618449 -83.191998 280 3390.690044 3505.618449 281 -2302.606457 3390.690044 282 5156.789899 -2302.606457 283 3413.603893 5156.789899 284 790.104503 3413.603893 285 2790.621144 790.104503 286 1451.897699 2790.621144 287 5370.690044 1451.897699 288 -2896.709859 5370.690044 289 -4678.885708 -2896.709859 290 -1152.767751 -4678.885708 291 686.111597 -1152.767751 292 862.211452 686.111597 293 73.603893 862.211452 294 -3476.871153 73.603893 295 464.207905 -3476.871153 296 1298.421802 464.207905 297 3590.528751 1298.421802 298 3214.903942 3590.528751 299 -3898.978101 3214.903942 300 1224.793446 -3898.978101 301 -791.799557 1224.793446 302 5491.496945 -791.799557 303 8850.911403 5491.496945 304 287.393543 8850.911403 305 -4293.906506 287.393543 306 -6976.691757 -4293.906506 307 -6737.184905 -6976.691757 308 2128.260140 -6737.184905 309 -1259.268361 2128.260140 310 934.115143 -1259.268361 311 3872.031688 934.115143 312 -4729.360754 3872.031688 313 -5449.674506 -4729.360754 314 1057.010890 -5449.674506 315 -424.160561 1057.010890 316 5172.925590 -424.160561 317 7509.629089 5172.925590 318 6960.011741 7509.629089 319 -348.775213 6960.011741 320 5271.607440 -348.775213 321 7703.331736 5271.607440 322 -2356.378005 7703.331736 323 6402.907488 -2356.378005 324 -2130.974555 6402.907488 325 -1436.285611 -2130.974555 326 -37.991805 -1436.285611 327 -646.267509 -37.991805 328 -10035.045266 -646.267509 329 -6672.902108 -10035.045266 330 -957.991805 -6672.902108 331 2712.100588 -957.991805 332 1013.331736 2712.100588 333 -162.989110 1013.331736 334 -6071.781455 -162.989110 335 -894.178663 -6071.781455 336 -4480.974555 -894.178663 337 -8596.378005 -4480.974555 338 544.428896 -8596.378005 339 -2145.077956 544.428896 340 6304.115143 -2145.077956 341 -3200.882161 6304.115143 342 442.008195 -3200.882161 343 2609.924739 442.008195 344 998.310938 2609.924739 345 5888.329040 998.310938 346 3949.135941 5888.329040 347 2064.115143 3949.135941 348 1197.504038 2064.115143 349 2678.260140 1197.504038 350 2966.535844 2678.260140 351 2103.621995 2966.535844 352 1796.604744 2103.621995 353 4302.414340 1796.604744 354 3935.761635 4302.414340 355 -215.861364 3935.761635 356 213.221241 -215.861364 357 8995.839439 213.221241 358 -4446.378005 8995.839439 359 -635.077956 -4446.378005 360 4829.135941 -635.077956 361 7440.818642 4829.135941 362 5484.461586 7440.818642 363 -5122.588355 5484.461586 364 -6257.074410 -5122.588355 365 -8863.141206 -6257.074410 366 -555.451411 -8863.141206 367 -524.671811 -555.451411 368 11513.248540 -524.671811 369 358.711693 11513.248540 370 -3622.495962 358.711693 371 -4867.991805 -3622.495962 372 -7085.077956 -4867.991805 373 -730.882161 -7085.077956 374 -10.006361 -730.882161 375 -3681.288307 -10.006361 376 2891.717935 -3681.288307 377 -5159.259163 2891.717935 378 -3950.167654 -5159.259163 379 1339.993639 -3950.167654 380 -629.641816 1339.993639 381 1902.501343 -629.641816 382 4548.218545 1902.501343 383 4998.310938 4548.218545 384 -5835.086861 4998.310938 385 3623.714389 -5835.086861 386 2733.732491 3623.714389 387 -1945.861364 2733.732491 388 702.414340 -1945.861364 389 4154.922044 702.414340 390 2633.308243 4154.922044 391 -2677.567558 2633.308243 392 3547.010890 -2677.567558 393 -3031.375309 3547.010890 394 -5494.271056 -3031.375309 395 -536.778759 -5494.271056 396 -736.285611 -536.778759 397 794.225639 -736.285611 398 -4325.460609 794.225639 399 3198.320136 -4325.460609 400 6884.115143 3198.320136 401 -804.671811 6884.115143 402 1168.711693 -804.671811 403 -255.368216 1168.711693 404 -1248.775213 -255.368216 405 -2734.602911 -1248.775213 406 -2173.371762 -2734.602911 407 -1640.882161 -2173.371762 408 -5987.991805 -1640.882161 409 -3520.882161 -5987.991805 410 -2656.345314 -3520.882161 411 -6681.375309 -2656.345314 412 2805.797844 -6681.375309 413 925.328189 2805.797844 414 2194.138636 925.328189 415 -5592.588355 2194.138636 416 -3835.884857 -5592.588355 417 3585.397089 -3835.884857 418 12208.651990 3585.397089 419 5552.100588 12208.651990 420 -2492.588355 5552.100588 421 -6962.095207 -2492.588355 422 -1608.060706 -6962.095207 423 560.842135 -1608.060706 424 -794.271056 560.842135 425 -340.457914 -794.271056 426 6038.310938 -340.457914 427 1509.518593 6038.310938 428 -2492.989110 1509.518593 429 -7142.071714 -2492.989110 430 -4789.199460 -7142.071714 431 -2913.362565 -4789.199460 432 -5107.991805 -2913.362565 433 946.628238 -5107.991805 434 1741.607440 946.628238 435 6040.758939 1741.607440 436 715.014437 6040.758939 437 4424.913139 715.014437 438 10113.308243 4424.913139 439 1225.032539 10113.308243 440 -9125.953757 1225.032539 441 -5284.764204 -9125.953757 442 3121.123490 -5284.764204 443 -1644.202156 3121.123490 444 3094.428896 -1644.202156 445 -8952.164107 3094.428896 446 -38.688210 -8952.164107 447 2092.907488 -38.688210 448 14725.415192 2092.907488 449 6964.115143 14725.415192 450 8412.432442 6964.115143 451 -4077.092512 8412.432442 452 4388.633888 -4077.092512 453 8646.222092 4388.633888 454 1420.818642 8646.222092 455 -5246.267509 1420.818642 456 -1938.060706 -5246.267509 457 -1248.858408 -1938.060706 458 957.411645 -1248.858408 459 -2232.164107 957.411645 460 7506.715240 -2232.164107 461 267.817790 7506.715240 462 -5749.988259 267.817790 463 -6257.991805 -5749.988259 464 -4263.371762 -6257.991805 465 -4605.409811 -4263.371762 466 183.714389 -4605.409811 467 -7.746953 183.714389 468 -1270.472209 -7.746953 469 -12926.982016 -1270.472209 470 -4487.991805 -12926.982016 471 6162.414340 -4487.991805 472 5107.435138 6162.414340 473 -4930.974555 5107.435138 474 -6310.167654 -4930.974555 475 -11613.362565 -6310.167654 476 6296.628238 -11613.362565 477 4200.417887 6296.628238 478 -2306.285611 4200.417887 479 -942.555665 -2306.285611 480 -215.469514 -942.555665 481 -6189.360754 -215.469514 482 -7276.285611 -6189.360754 483 2679.117839 -7276.285611 484 2272.847786 2679.117839 485 4240.818642 2272.847786 486 7318.329040 4240.818642 487 627.522140 7318.329040 488 213.221241 627.522140 489 -2264.178663 213.221241 490 -10613.454958 -2264.178663 491 1016.646340 -10613.454958 492 2693.230438 1016.646340 493 3257.010890 2693.230438 494 -5977.959115 3257.010890 495 -2983.302862 -5977.959115 496 3596.697138 -2983.302862 497 -982.878615 3596.697138 498 2481.625542 -982.878615 499 -1802.989110 2481.625542 500 9740.394394 -1802.989110 501 1035.821337 9740.394394 502 2272.847786 1035.821337 503 4493.714389 2272.847786 504 -4780.974555 4493.714389 505 -270.481407 -4780.974555 506 8448.711693 -270.481407 507 -2415.478711 8448.711693 508 -2456.285611 -2415.478711 509 1370.325494 -2456.285611 510 -3009.564011 1370.325494 511 7634.608291 -3009.564011 512 802.524836 7634.608291 513 8541.311790 802.524836 514 1646.120794 8541.311790 515 -2537.968312 1646.120794 516 -3300.974555 -2537.968312 517 -3530.974555 -3300.974555 518 641.547737 -3530.974555 519 -1872.979912 641.547737 520 1973.239343 -1872.979912 521 -4011.288307 1973.239343 522 -4177.498657 -4011.288307 523 -3249.988259 -4177.498657 524 926.628238 -3249.988259 525 -168.350965 926.628238 526 3159.117839 -168.350965 527 3303.621995 3159.117839 528 2624.521289 3303.621995 529 -2708.775213 2624.521289 530 1128.735186 -2708.775213 531 -3068.775213 1128.735186 532 -2050.882161 -3068.775213 533 459.518593 -2050.882161 534 -10485.077956 459.518593 535 -1059.564011 -10485.077956 536 5729.518593 -1059.564011 537 9203.654686 5729.518593 538 3776.222092 9203.654686 539 -3570.167654 3776.222092 540 -8969.181358 -3570.167654 541 -4941.264814 -8969.181358 542 5568.218545 -4941.264814 543 239.993639 5568.218545 544 491.607440 239.993639 545 -828.798706 491.607440 546 774.225639 -828.798706 547 3937.835893 774.225639 548 -3460.481407 3937.835893 549 -170.882161 -3460.481407 550 -864.339956 -170.882161 551 -1789.157865 -864.339956 552 -11213.030710 -1789.157865 553 57.411645 -11213.030710 554 -4825.884857 57.411645 555 8262.017392 -4825.884857 556 -5114.602911 8262.017392 557 966.203990 -5114.602911 558 2616.628238 966.203990 559 9770.011741 2616.628238 560 506.535844 9770.011741 561 -4694.985563 506.535844 562 4471.607440 -4694.985563 563 3355.014437 4471.607440 564 3107.522140 3355.014437 565 -266.378005 3107.522140 566 1575.014437 -266.378005 567 -142.588355 1575.014437 568 -3648.443358 -142.588355 569 -9294.584808 -3648.443358 570 1700.316589 -9294.584808 571 -4062.971008 1700.316589 572 -2759.199460 -4062.971008 573 -3743.464156 -2759.199460 574 -7397.567558 -3743.464156 575 7534.115143 -7397.567558 576 225.821337 7534.115143 577 6982.414340 225.821337 578 5313.331736 6982.414340 579 5512.031688 5313.331736 580 59.518593 5512.031688 581 -4099.674506 59.518593 582 1159.117839 -4099.674506 583 -4635.478711 1159.117839 584 3490.011741 -4635.478711 585 -2184.561315 3490.011741 586 5949.440789 -2184.561315 587 -2466.285611 5949.440789 588 10.417887 -2466.285611 589 1251.242889 10.417887 590 -1641.772257 1251.242889 591 3899.058136 -1641.772257 592 -430.075261 3899.058136 593 -2005.571104 -430.075261 594 4029.542086 -2005.571104 595 -2968.775213 4029.542086 596 -1388.383362 -2968.775213 597 4073.248540 -1388.383362 598 11606.162389 4073.248540 599 -2864.764204 11606.162389 600 1337.826988 -2864.764204 601 -7966.691757 1337.826988 602 3307.946388 -7966.691757 603 5287.462438 3307.946388 604 -6809.923171 5287.462438 605 6747.181083 -6809.923171 606 -567.493272 6747.181083 607 -488.816813 -567.493272 608 3412.967185 -488.816813 609 5308.753288 3412.967185 610 5467.946388 5308.753288 611 -5451.647467 5467.946388 612 -6277.534867 -5451.647467 613 -6985.921066 -6277.534867 614 -2307.493272 -6985.921066 615 -4024.519720 -2307.493272 616 -1680.347418 -4024.519720 617 -6195.902964 -1680.347418 618 -1633.814118 -6195.902964 619 -976.727967 -1633.814118 620 19423.349838 -976.727967 621 -707.457062 19423.349838 622 2701.284485 -707.457062 623 5337.563735 2701.284485 624 854.156738 5337.563735 625 -5006.709865 854.156738 626 3904.180231 -5006.709865 627 164.249132 3904.180231 628 12376.756835 164.249132 629 203.272033 12376.756835 630 -2843.745218 203.272033 631 -3558.410668 -2843.745218 632 -7856.157014 -3558.410668 633 2319.076238 -7856.157014 634 3171.607434 2319.076238 635 -1906.244016 3171.607434 636 1280.376286 -1906.244016 637 -1277.032815 1280.376286 638 6887.070587 -1277.032815 639 -4714.943968 6887.070587 640 -7351.324517 -4714.943968 641 2249.993634 -7351.324517 642 -4135.114166 2249.993634 643 2784.649886 -4135.114166 644 15502.967185 2784.649886 645 7249.569386 15502.967185 646 6356.586637 7249.569386 647 573.672788 6356.586637 648 -1187.502470 573.672788 649 -2713.814118 -1187.502470 650 -2767.857817 -2713.814118 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7q8jl1292085873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/8jh151292085873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/9jh151292085873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/rcomp/tmp/10jh151292085873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/11frze1292085873.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12jrx21292085873.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13fjdb1292085873.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/14ijbh1292085873.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/153kan1292085873.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/16p38b1292085873.tab") + } > > try(system("convert tmp/1n73f1292085873.ps tmp/1n73f1292085873.png",intern=TRUE)) character(0) > try(system("convert tmp/2n73f1292085873.ps tmp/2n73f1292085873.png",intern=TRUE)) character(0) > try(system("convert tmp/3n73f1292085873.ps tmp/3n73f1292085873.png",intern=TRUE)) character(0) > try(system("convert tmp/4xy2i1292085873.ps tmp/4xy2i1292085873.png",intern=TRUE)) character(0) > try(system("convert tmp/5xy2i1292085873.ps tmp/5xy2i1292085873.png",intern=TRUE)) character(0) > try(system("convert tmp/6xy2i1292085873.ps tmp/6xy2i1292085873.png",intern=TRUE)) character(0) > try(system("convert tmp/7q8jl1292085873.ps tmp/7q8jl1292085873.png",intern=TRUE)) character(0) > try(system("convert tmp/8jh151292085873.ps tmp/8jh151292085873.png",intern=TRUE)) character(0) > try(system("convert tmp/9jh151292085873.ps tmp/9jh151292085873.png",intern=TRUE)) character(0) > try(system("convert tmp/10jh151292085873.ps tmp/10jh151292085873.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 24.189 3.016 55.977