Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 1418 1418 210907 56 56 3 3 79 79 30 30 115 115 1 869 869 120982 56 56 4 4 58 58 28 28 109 109 1 1530 1530 176508 54 54 12 12 60 60 38 38 146 146 1 2172 2172 179321 89 89 2 2 108 108 30 30 116 116 1 901 901 123185 40 40 1 1 49 49 22 22 68 68 1 463 463 52746 25 25 3 3 0 0 26 26 101 101 1 3201 3201 385534 92 92 0 0 121 121 25 25 96 96 1 371 371 33170 18 18 0 0 1 1 18 18 67 67 0 1192 0 101645 63 0 0 0 20 0 11 0 44 0 1 1583 1583 149061 44 44 5 5 43 43 26 26 100 100 1 1439 1439 165446 33 33 0 0 69 69 25 25 93 93 1 1764 1764 237213 84 84 0 0 78 78 38 38 140 140 1 1495 1495 173326 88 88 7 7 86 86 44 44 166 166 1 1373 1373 133131 55 55 7 7 44 44 30 30 99 99 1 2187 2187 258873 60 60 3 3 104 104 40 40 139 139 1 1491 1491 180083 66 66 9 9 63 63 34 34 130 130 1 4041 4041 324799 154 154 0 0 158 158 47 47 181 181 1 1706 1706 230964 53 53 4 4 102 102 30 30 116 116 1 2152 2152 236785 119 119 3 3 77 77 31 31 116 116 1 1036 1036 135473 41 41 0 0 82 82 23 23 88 88 1 1882 1882 202925 61 61 7 7 115 115 36 36 139 139 1 1929 1929 215147 58 58 0 0 101 101 36 36 135 135 1 2242 2242 344297 75 75 1 1 80 80 30 30 108 108 1 1220 1220 153935 33 33 5 5 50 50 25 25 89 89 1 1289 1289 132943 40 40 7 7 83 83 39 39 156 156 1 2515 2515 174724 92 92 0 0 123 123 34 34 129 129 1 2147 2147 174415 100 100 0 0 73 73 31 31 118 118 1 2352 2352 225548 112 112 5 5 81 81 31 31 118 118 1 1638 1638 223632 73 73 0 0 105 105 33 33 125 125 1 1222 1222 124817 40 40 0 0 47 47 25 25 95 95 1 1812 1812 221698 45 45 0 0 105 105 33 33 126 126 1 1677 1677 210767 60 60 3 3 94 94 35 35 135 135 1 1579 1579 170266 62 62 4 4 44 44 42 42 154 154 1 1731 1731 260561 75 75 1 1 114 114 43 43 165 165 1 807 807 84853 31 31 4 4 38 38 30 30 113 113 1 2452 2452 294424 77 77 2 2 107 107 33 33 127 127 0 829 0 101011 34 0 0 0 30 0 13 0 52 0 1 1940 1940 215641 46 46 0 0 71 71 32 32 121 121 1 2662 2662 325107 99 99 0 0 84 84 36 36 136 136 0 186 0 7176 17 0 0 0 0 0 0 0 0 0 1 1499 1499 167542 66 66 2 2 59 59 28 28 108 108 1 865 865 106408 30 30 1 1 33 33 14 14 46 46 0 1793 0 96560 76 0 0 0 42 0 17 0 54 0 1 2527 2527 265769 146 146 2 2 96 96 32 32 124 124 1 2747 2747 269651 67 67 10 10 106 106 30 30 115 115 1 1324 1324 149112 56 56 6 6 56 56 35 35 128 128 0 2702 0 175824 107 0 0 0 57 0 20 0 80 0 1 1383 1383 152871 58 58 5 5 59 59 28 28 97 97 1 1179 1179 111665 34 34 4 4 39 39 28 28 104 104 1 2099 2099 116408 61 61 1 1 34 34 39 39 59 59 1 4308 4308 362301 119 119 2 2 76 76 34 34 125 125 1 918 918 78800 42 42 2 2 20 20 26 26 82 82 1 1831 1831 183167 66 66 0 0 91 91 39 39 149 149 1 3373 3373 277965 89 89 8 8 115 115 39 39 149 149 1 1713 1713 150629 44 44 3 3 85 85 33 33 122 122 1 1438 1438 168809 66 66 0 0 76 76 28 28 118 118 1 496 496 24188 24 24 0 0 8 8 4 4 12 12 1 2253 2253 329267 259 259 8 8 79 79 39 39 144 144 1 744 744 65029 17 17 5 5 21 21 18 18 67 67 1 1161 1161 101097 64 64 3 3 30 30 14 14 52 52 1 2352 2352 218946 41 41 1 1 76 76 29 29 108 108 1 2144 2144 244052 68 68 5 5 101 101 44 44 166 166 0 4691 0 341570 168 0 1 0 94 0 21 0 80 0 0 1112 0 103597 43 0 1 0 27 0 16 0 60 0 1 2694 2694 233328 132 132 5 5 92 92 28 28 107 107 1 1973 1973 256462 105 105 0 0 123 123 35 35 127 127 1 1769 1769 206161 71 71 12 12 75 75 28 28 107 107 1 3148 3148 311473 112 112 8 8 128 128 38 38 146 146 1 2474 2474 235800 94 94 8 8 105 105 23 23 84 84 1 2084 2084 177939 82 82 8 8 55 55 36 36 141 141 1 1954 1954 207176 70 70 8 8 56 56 32 32 123 123 1 1226 1226 196553 57 57 2 2 41 41 29 29 111 111 1 1389 1389 174184 53 53 0 0 72 72 25 25 98 98 1 1496 1496 143246 103 103 5 5 67 67 27 27 105 105 1 2269 2269 187559 121 121 8 8 75 75 36 36 135 135 1 1833 1833 187681 62 62 2 2 114 114 28 28 107 107 1 1268 1268 119016 52 52 5 5 118 118 23 23 85 85 1 1943 1943 182192 52 52 12 12 77 77 40 40 155 155 1 893 893 73566 32 32 6 6 22 22 23 23 88 88 1 1762 1762 194979 62 62 7 7 66 66 40 40 155 155 1 1403 1403 167488 45 45 2 2 69 69 28 28 104 104 1 1425 1425 143756 46 46 0 0 105 105 34 34 132 132 1 1857 1857 275541 63 63 4 4 116 116 33 33 127 127 1 1840 1840 243199 75 75 3 3 88 88 28 28 108 108 1 1502 1502 182999 88 88 6 6 73 73 34 34 129 129 1 1441 1441 135649 46 46 2 2 99 99 30 30 116 116 1 1420 1420 152299 53 53 0 0 62 62 33 33 122 122 1 1416 1416 120221 37 37 1 1 53 53 22 22 85 85 1 2970 2970 346485 90 90 0 0 118 118 38 38 147 147 1 1317 1317 145790 63 63 5 5 30 30 26 26 99 99 1 1644 1644 193339 78 78 2 2 100 100 35 35 87 87 1 870 870 80953 25 25 0 0 49 49 8 8 28 28 1 1654 1654 122774 45 45 0 0 24 24 24 24 90 90 1 1054 1054 130585 46 46 5 5 67 67 29 29 109 109 0 937 0 112611 41 0 0 0 46 0 20 0 78 0 1 3004 3004 286468 144 144 1 1 57 57 29 29 111 111 1 2008 2008 241066 82 82 0 0 75 75 45 45 158 158 1 2547 2547 148446 91 91 1 1 135 135 37 37 141 141 1 1885 1885 204713 71 71 1 1 68 68 33 33 122 122 1 1626 1626 182079 63 63 2 2 124 124 33 33 124 124 1 1468 1468 140344 53 53 6 6 33 33 25 25 93 93 1 2445 2445 220516 62 62 1 1 98 98 32 32 124 124 1 1964 1964 243060 63 63 4 4 58 58 29 29 112 112 1 1381 1381 162765 32 32 2 2 68 68 28 28 108 108 1 1369 1369 182613 39 39 3 3 81 81 28 28 99 99 1 1659 1659 232138 62 62 0 0 131 131 31 31 117 117 1 2888 2888 265318 117 117 10 10 110 110 52 52 199 199 0 1290 0 85574 34 0 0 0 37 0 21 0 78 0 1 2845 2845 310839 92 92 9 9 130 130 24 24 91 91 1 1982 1982 225060 93 93 7 7 93 93 41 41 158 158 1 1904 1904 232317 54 54 0 0 118 118 33 33 126 126 1 1391 1391 144966 144 144 0 0 39 39 32 32 122 122 1 602 602 43287 14 14 4 4 13 13 19 19 71 71 1 1743 1743 155754 61 61 4 4 74 74 20 20 75 75 1 1559 1559 164709 109 109 0 0 81 81 31 31 115 115 1 2014 2014 201940 38 38 0 0 109 109 31 31 119 119 1 2143 2143 235454 73 73 0 0 151 151 32 32 124 124 0 2146 0 220801 75 0 1 0 51 0 18 0 72 0 1 874 874 99466 50 50 0 0 28 28 23 23 91 91 0 1590 0 92661 61 0 1 0 40 0 17 0 45 0 0 1590 0 133328 55 0 0 0 56 0 20 0 78 0 0 1210 0 61361 77 0 0 0 27 0 12 0 39 0 0 2072 0 125930 75 0 4 0 37 0 17 0 68 0 1 1281 1281 100750 72 72 0 0 83 83 30 30 119 119 1 1401 1401 224549 50 50 4 4 54 54 31 31 117 117 0 834 0 82316 32 0 4 0 27 0 10 0 39 0 0 1105 0 102010 53 0 3 0 28 0 13 0 50 0 0 1272 0 101523 42 0 0 0 59 0 22 0 88 0 1 1944 1944 243511 71 71 0 0 133 133 42 42 155 155 1 391 391 22938 10 10 0 0 12 12 1 1 0 0 0 761 0 41566 35 0 5 0 0 0 9 0 36 0 1 1605 1605 152474 65 65 0 0 106 106 32 32 123 123 1 530 530 61857 25 25 4 4 23 23 11 11 32 32 0 1988 0 99923 66 0 0 0 44 0 25 0 99 0 1 1386 1386 132487 41 41 0 0 71 71 36 36 136 136 1 2395 2395 317394 86 86 1 1 116 116 31 31 117 117 1 387 387 21054 16 16 0 0 4 4 0 0 0 0 1 1742 1742 209641 42 42 5 5 62 62 24 24 88 88 0 620 0 22648 19 0 0 0 12 0 13 0 39 0 1 449 449 31414 19 19 0 0 18 18 8 8 25 25 0 800 0 46698 45 0 0 0 14 0 13 0 52 0 0 1684 0 131698 65 0 0 0 60 0 19 0 75 0 0 1050 0 91735 35 0 0 0 7 0 18 0 71 0 1 2699 2699 244749 95 95 2 2 98 98 33 33 124 124 1 1606 1606 184510 49 49 7 7 64 64 40 40 151 151 0 1502 0 79863 37 0 1 0 29 0 22 0 71 0 1 1204 1204 128423 64 64 8 8 32 32 38 38 145 145 1 1138 1138 97839 38 38 2 2 25 25 24 24 87 87 1 568 568 38214 34 34 0 0 16 16 8 8 27 27 1 1459 1459 151101 32 32 2 2 48 48 35 35 131 131 1 2158 2158 272458 65 65 0 0 100 100 43 43 162 162 1 1111 1111 172494 52 52 0 0 46 46 43 43 165 165 0 1421 0 108043 62 0 1 0 45 0 14 0 54 0 1 2833 2833 328107 65 65 3 3 129 129 41 41 159 159 1 1955 1955 250579 83 83 0 0 130 130 38 38 147 147 1 2922 2922 351067 95 95 3 3 136 136 45 45 170 170 1 1002 1002 158015 29 29 0 0 59 59 31 31 119 119 0 1060 0 98866 18 0 0 0 25 0 13 0 49 0 1 956 956 85439 33 33 0 0 32 32 28 28 104 104 1 2186 2186 229242 247 247 4 4 63 63 31 31 120 120 1 3604 3604 351619 139 139 4 4 95 95 40 40 150 150 1 1035 1035 84207 29 29 11 11 14 14 30 30 112 112 0 1417 0 120445 118 0 0 0 36 0 16 0 59 0 1 3261 3261 324598 110 110 0 0 113 113 37 37 136 136 1 1587 1587 131069 67 67 4 4 47 47 30 30 107 107 1 1424 1424 204271 42 42 0 0 92 92 35 35 130 130 1 1701 1701 165543 65 65 1 1 70 70 32 32 115 115 1 1249 1249 141722 94 94 0 0 19 19 27 27 107 107 0 946 0 116048 64 0 0 0 50 0 20 0 75 0 0 1926 0 250047 81 0 0 0 41 0 18 0 71 0 1 3352 3352 299775 95 95 9 9 91 91 31 31 120 120 1 1641 1641 195838 67 67 1 1 111 111 31 31 116 116 1 2035 2035 173260 63 63 3 3 41 41 21 21 79 79 1 2312 2312 254488 83 83 10 10 120 120 39 39 150 150 1 1369 1369 104389 45 45 5 5 135 135 41 41 156 156 0 1577 0 136084 30 0 0 0 27 0 13 0 51 0 1 2201 2201 199476 70 70 2 2 87 87 32 32 118 118 0 961 0 92499 32 0 0 0 25 0 18 0 71 0 1 1900 1900 224330 83 83 1 1 131 131 39 39 144 144 0 1254 0 135781 31 0 2 0 45 0 14 0 47 0 0 1335 0 74408 67 0 4 0 29 0 7 0 28 0 0 1597 0 81240 66 0 0 0 58 0 17 0 68 0 1 207 207 14688 10 10 0 0 4 4 0 0 0 0 1 1645 1645 181633 70 70 2 2 47 47 30 30 110 110 1 2429 2429 271856 103 103 1 1 109 109 37 37 147 147 1 151 151 7199 5 5 0 0 7 7 0 0 0 0 1 474 474 46660 20 20 0 0 12 12 5 5 15 15 1 141 141 17547 5 5 0 0 0 0 1 1 4 4 0 1639 0 133368 36 0 1 0 37 0 16 0 64 0 1 872 872 95227 34 34 0 0 37 37 32 32 111 111 1 1318 1318 152601 48 48 2 2 46 46 24 24 85 85 0 1018 0 98146 40 0 0 0 15 0 17 0 68 0 0 1383 0 79619 43 0 3 0 42 0 11 0 40 0 0 1314 0 59194 31 0 6 0 7 0 24 0 80 0 0 1335 0 139942 42 0 0 0 54 0 22 0 88 0 0 1403 0 118612 46 0 2 0 54 0 12 0 48 0 0 910 0 72880 33 0 0 0 14 0 19 0 76 0 0 616 0 65475 18 0 2 0 16 0 13 0 51 0 0 1407 0 99643 55 0 1 0 33 0 17 0 67 0 0 771 0 71965 35 0 1 0 32 0 15 0 59 0 0 766 0 77272 59 0 2 0 21 0 16 0 61 0 0 473 0 49289 19 0 1 0 15 0 24 0 76 0 0 1376 0 135131 66 0 0 0 38 0 15 0 60 0 0 1232 0 108446 60 0 1 0 22 0 17 0 68 0 0 1521 0 89746 36 0 3 0 28 0 18 0 71 0 0 572 0 44296 25 0 0 0 10 0 20 0 76 0 0 1059 0 77648 47 0 0 0 31 0 16 0 62 0 0 1544 0 181528 54 0 0 0 32 0 16 0 61 0 0 1230 0 134019 53 0 0 0 32 0 18 0 67 0 0 1206 0 124064 40 0 1 0 43 0 22 0 88 0 0 1205 0 92630 40 0 4 0 27 0 8 0 30 0 0 1255 0 121848 39 0 0 0 37 0 17 0 64 0 0 613 0 52915 14 0 0 0 20 0 18 0 68 0 0 721 0 81872 45 0 0 0 32 0 16 0 64 0 0 1109 0 58981 36 0 7 0 0 0 23 0 91 0 0 740 0 53515 28 0 2 0 5 0 22 0 88 0 0 1126 0 60812 44 0 0 0 26 0 13 0 52 0 0 728 0 56375 30 0 7 0 10 0 13 0 49 0 0 689 0 65490 22 0 3 0 27 0 16 0 62 0 0 592 0 80949 17 0 0 0 11 0 16 0 61 0 0 995 0 76302 31 0 0 0 29 0 20 0 76 0 0 1613 0 104011 55 0 6 0 25 0 22 0 88 0 0 2048 0 98104 54 0 2 0 55 0 17 0 66 0 0 705 0 67989 21 0 0 0 23 0 18 0 71 0 0 301 0 30989 14 0 0 0 5 0 17 0 68 0 0 1803 0 135458 81 0 3 0 43 0 12 0 48 0 0 799 0 73504 35 0 0 0 23 0 7 0 25 0 0 861 0 63123 43 0 1 0 34 0 17 0 68 0 0 1186 0 61254 46 0 1 0 36 0 14 0 41 0 0 1451 0 74914 30 0 0 0 35 0 23 0 90 0 0 628 0 31774 23 0 1 0 0 0 17 0 66 0 0 1161 0 81437 38 0 0 0 37 0 14 0 54 0 0 1463 0 87186 54 0 0 0 28 0 15 0 59 0 0 742 0 50090 20 0 0 0 16 0 17 0 60 0 0 979 0 65745 53 0 0 0 26 0 21 0 77 0 0 675 0 56653 45 0 0 0 38 0 18 0 68 0 0 1241 0 158399 39 0 0 0 23 0 18 0 72 0 0 676 0 46455 20 0 0 0 22 0 17 0 67 0 0 1049 0 73624 24 0 0 0 30 0 17 0 64 0 0 620 0 38395 31 0 0 0 16 0 16 0 63 0 0 1081 0 91899 35 0 0 0 18 0 15 0 59 0 0 1688 0 139526 151 0 0 0 28 0 21 0 84 0 0 736 0 52164 52 0 0 0 32 0 16 0 64 0 0 617 0 51567 30 0 2 0 21 0 14 0 56 0 0 812 0 70551 31 0 0 0 23 0 15 0 54 0 0 1051 0 84856 29 0 1 0 29 0 17 0 67 0 0 1656 0 102538 57 0 1 0 50 0 15 0 58 0 0 705 0 86678 40 0 0 0 12 0 15 0 59 0 0 945 0 85709 44 0 0 0 21 0 10 0 40 0 0 554 0 34662 25 0 0 0 18 0 6 0 22 0 0 1597 0 150580 77 0 0 0 27 0 22 0 83 0 0 982 0 99611 35 0 0 0 41 0 21 0 81 0 0 222 0 19349 11 0 0 0 13 0 1 0 2 0 0 1212 0 99373 63 0 1 0 12 0 18 0 72 0 0 1143 0 86230 44 0 0 0 21 0 17 0 61 0 0 435 0 30837 19 0 0 0 8 0 4 0 15 0 0 532 0 31706 13 0 0 0 26 0 10 0 32 0 0 882 0 89806 42 0 0 0 27 0 16 0 62 0 0 608 0 62088 38 0 1 0 13 0 16 0 58 0 0 459 0 40151 29 0 0 0 16 0 9 0 36 0 0 578 0 27634 20 0 0 0 2 0 16 0 59 0 0 826 0 76990 27 0 0 0 42 0 17 0 68 0 0 509 0 37460 20 0 0 0 5 0 7 0 21 0 0 717 0 54157 19 0 0 0 37 0 15 0 55 0 0 637 0 49862 37 0 0 0 17 0 14 0 54 0 0 857 0 84337 26 0 0 0 38 0 14 0 55 0 0 830 0 64175 42 0 0 0 37 0 18 0 72 0 0 652 0 59382 49 0 0 0 29 0 12 0 41 0 0 707 0 119308 30 0 0 0 32 0 16 0 61 0 0 954 0 76702 49 0 0 0 35 0 21 0 67 0 0 1461 0 103425 67 0 1 0 17 0 19 0 76 0 0 672 0 70344 28 0 0 0 20 0 16 0 64 0 0 778 0 43410 19 0 0 0 7 0 1 0 3 0 0 1141 0 104838 49 0 1 0 46 0 16 0 63 0 0 680 0 62215 27 0 0 0 24 0 10 0 40 0 0 1090 0 69304 30 0 6 0 40 0 19 0 69 0 0 616 0 53117 22 0 3 0 3 0 12 0 48 0 0 285 0 19764 12 0 1 0 10 0 2 0 8 0 0 1145 0 86680 31 0 2 0 37 0 14 0 52 0 0 733 0 84105 20 0 0 0 17 0 17 0 66 0 0 888 0 77945 20 0 0 0 28 0 19 0 76 0 0 849 0 89113 39 0 0 0 19 0 14 0 43 0 0 1182 0 91005 29 0 3 0 29 0 11 0 39 0 0 528 0 40248 16 0 1 0 8 0 4 0 14 0 0 642 0 64187 27 0 0 0 10 0 16 0 61 0 0 947 0 50857 21 0 0 0 15 0 20 0 71 0 0 819 0 56613 19 0 1 0 15 0 12 0 44 0 0 757 0 62792 35 0 0 0 28 0 15 0 60 0 0 894 0 72535 14 0 0 0 17 0 16 0 64 0
Names of X columns:
pop pageviews pageviews_p time_in_rfc logins logins_p shared_compendiums shared_compendiums_p blogged_computations blogged_computations_p compendiums_reviewed compendiums_reviewed_p feedback_messages feedback_messages_p
Sample Range:
(leave blank to include all observations)
From:
To:
Column Number of Endogenous Series
(?)
Fixed Seasonal Effects
Do not include Seasonal Dummies
Do not include Seasonal Dummies
Include Seasonal Dummies
Type of Equation
No Linear Trend
No Linear Trend
Linear Trend
First Differences
Seasonal Differences (s)
First and Seasonal Differences (s)
Degree of Predetermination (lagged endogenous variables)
Degree of Seasonal Predetermination
Seasonality
12
1
2
3
4
5
6
7
8
9
10
11
12
Chart options
R Code
library(lattice) library(lmtest) n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test 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 k <- length(x[1,]) df <- as.data.frame(x) (mylm <- lm(df)) (mysum <- summary(mylm)) if (n > n25) { kp3 <- k + 3 nmkm3 <- n - k - 3 gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) numgqtests <- 0 numsignificant1 <- 0 numsignificant5 <- 0 numsignificant10 <- 0 for (mypoint in kp3:nmkm3) { j <- 0 numgqtests <- numgqtests + 1 for (myalt in c('greater', 'two.sided', 'less')) { j <- j + 1 gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value } if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 } gqarr } bitmap(file='test0.png') plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') points(x[,1]-mysum$resid) grid() dev.off() bitmap(file='test1.png') plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') grid() dev.off() bitmap(file='test2.png') hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') grid() dev.off() bitmap(file='test3.png') densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') dev.off() bitmap(file='test4.png') qqnorm(mysum$resid, main='Residual Normal Q-Q Plot') qqline(mysum$resid) grid() dev.off() (myerror <- as.ts(mysum$resid)) bitmap(file='test5.png') dum <- cbind(lag(myerror,k=1),myerror) dum dum1 <- dum[2:length(myerror),] dum1 z <- as.data.frame(dum1) z plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') lines(lowess(z)) abline(lm(z)) grid() dev.off() bitmap(file='test6.png') acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') grid() dev.off() bitmap(file='test7.png') pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') grid() dev.off() bitmap(file='test8.png') opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) plot(mylm, las = 1, sub='Residual Diagnostics') par(opar) dev.off() if (n > n25) { bitmap(file='test9.png') plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') grid() dev.off() } load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) a<-table.row.end(a) myeq <- colnames(x)[1] myeq <- paste(myeq, '[t] = ', sep='') for (i in 1:k){ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') myeq <- paste(myeq, 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='mytable1.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Variable',header=TRUE) a<-table.element(a,'Parameter',header=TRUE) a<-table.element(a,'S.D.',header=TRUE) a<-table.element(a,'T-STAT<br />H0: parameter = 0',header=TRUE) a<-table.element(a,'2-tail p-value',header=TRUE) a<-table.element(a,'1-tail p-value',header=TRUE) a<-table.row.end(a) for (i in 1:k){ a<-table.row.start(a) a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) a<-table.element(a,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='mytable2.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Multiple R',1,TRUE) a<-table.element(a, 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='mytable3.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<br />Forecast', 1, TRUE) a<-table.element(a, 'Residuals<br />Prediction Error', 1, TRUE) a<-table.row.end(a) for (i in 1:n) { a<-table.row.start(a) a<-table.element(a,i, 1, TRUE) a<-table.element(a,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='mytable4.tab') if (n > n25) { a<-table.start() a<-table.row.start(a) a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-values',header=TRUE) a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'breakpoint index',header=TRUE) a<-table.element(a,'greater',header=TRUE) a<-table.element(a,'2-sided',header=TRUE) a<-table.element(a,'less',header=TRUE) a<-table.row.end(a) for (mypoint in kp3:nmkm3) { a<-table.row.start(a) a<-table.element(a,mypoint,header=TRUE) a<-table.element(a,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='mytable5.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Description',header=TRUE) a<-table.element(a,'# significant tests',header=TRUE) a<-table.element(a,'% significant tests',header=TRUE) a<-table.element(a,'OK/NOK',header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'1% type I error level',header=TRUE) a<-table.element(a,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='mytable6.tab') }
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation