Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 1 1 1418 1418 56 56 30 30 145 145 3 3 210907 1 2 2 869 869 56 56 28 28 101 101 4 4 120982 1 3 3 1530 1530 54 54 38 38 98 98 12 12 176508 1 4 4 2172 2172 89 89 30 30 132 132 2 2 179321 1 5 5 901 901 40 40 22 22 60 60 1 1 123185 1 6 6 463 463 25 25 26 26 38 38 3 3 52746 1 7 7 3201 3201 92 92 25 25 144 144 0 0 385534 1 8 8 371 371 18 18 18 18 5 5 0 0 33170 1 9 9 1583 1583 44 44 26 26 84 84 5 5 149061 1 10 10 1439 1439 33 33 25 25 79 79 0 0 165446 1 11 11 1764 1764 84 84 38 38 127 127 0 0 237213 1 12 12 1495 1495 88 88 44 44 78 78 7 7 173326 1 13 13 1373 1373 55 55 30 30 60 60 7 7 133131 1 14 14 2187 2187 60 60 40 40 131 131 3 3 258873 1 15 15 1491 1491 66 66 34 34 84 84 9 9 180083 1 16 16 4041 4041 154 154 47 47 133 133 0 0 324799 1 17 17 1706 1706 53 53 30 30 150 150 4 4 230964 1 18 18 2152 2152 119 119 31 31 91 91 3 3 236785 1 19 19 1036 1036 41 41 23 23 132 132 0 0 135473 1 20 20 1882 1882 61 61 36 36 136 136 7 7 202925 1 21 21 1929 1929 58 58 36 36 124 124 0 0 215147 1 22 22 2242 2242 75 75 30 30 118 118 1 1 344297 1 23 23 1220 1220 33 33 25 25 70 70 5 5 153935 1 24 24 1289 1289 40 40 39 39 107 107 7 7 132943 1 25 25 2515 2515 92 92 34 34 119 119 0 0 174724 1 26 26 2147 2147 100 100 31 31 89 89 0 0 174415 1 27 27 2352 2352 112 112 31 31 112 112 5 5 225548 1 28 28 1638 1638 73 73 33 33 108 108 0 0 223632 1 29 29 1222 1222 40 40 25 25 52 52 0 0 124817 1 30 30 1812 1812 45 45 33 33 112 112 0 0 221698 1 31 31 1677 1677 60 60 35 35 116 116 3 3 210767 1 32 32 1579 1579 62 62 42 42 123 123 4 4 170266 1 33 33 1731 1731 75 75 43 43 125 125 1 1 260561 1 34 34 807 807 31 31 30 30 27 27 4 4 84853 1 35 35 2452 2452 77 77 33 33 162 162 2 2 294424 1 36 36 1940 1940 46 46 32 32 64 64 0 0 215641 1 37 37 2662 2662 99 99 36 36 92 92 0 0 325107 1 38 38 1499 1499 66 66 28 28 83 83 2 2 167542 1 39 39 865 865 30 30 14 14 41 41 1 1 106408 1 40 40 2527 2527 146 146 32 32 120 120 2 2 265769 1 41 41 2747 2747 67 67 30 30 105 105 10 10 269651 1 42 42 1324 1324 56 56 35 35 79 79 6 6 149112 1 43 43 1383 1383 58 58 28 28 70 70 5 5 152871 1 44 44 1179 1179 34 34 28 28 55 55 4 4 111665 1 45 45 2099 2099 61 61 39 39 39 39 1 1 116408 1 46 46 4308 4308 119 119 34 34 67 67 2 2 362301 1 47 47 918 918 42 42 26 26 21 21 2 2 78800 1 48 48 1831 1831 66 66 39 39 127 127 0 0 183167 1 49 49 3373 3373 89 89 39 39 152 152 8 8 277965 1 50 50 1713 1713 44 44 33 33 113 113 3 3 150629 1 51 51 1438 1438 66 66 28 28 99 99 0 0 168809 1 52 52 496 496 24 24 4 4 7 7 0 0 24188 1 53 53 2253 2253 259 259 39 39 141 141 8 8 329267 1 54 54 744 744 17 17 18 18 21 21 5 5 65029 1 55 55 1161 1161 64 64 14 14 35 35 3 3 101097 1 56 56 2352 2352 41 41 29 29 109 109 1 1 218946 1 57 57 2144 2144 68 68 44 44 133 133 5 5 244052 1 58 58 2694 2694 132 132 28 28 230 230 5 5 233328 1 59 59 1973 1973 105 105 35 35 166 166 0 0 256462 1 60 60 1769 1769 71 71 28 28 68 68 12 12 206161 1 61 61 3148 3148 112 112 38 38 147 147 8 8 311473 1 62 62 2474 2474 94 94 23 23 179 179 8 8 235800 1 63 63 2084 2084 82 82 36 36 61 61 8 8 177939 1 64 64 1954 1954 70 70 32 32 101 101 8 8 207176 1 65 65 1226 1226 57 57 29 29 108 108 2 2 196553 1 66 66 1389 1389 53 53 25 25 90 90 0 0 174184 1 67 67 1496 1496 103 103 27 27 114 114 5 5 143246 1 68 68 2269 2269 121 121 36 36 103 103 8 8 187559 1 69 69 1833 1833 62 62 28 28 142 142 2 2 187681 1 70 70 1268 1268 52 52 23 23 79 79 5 5 119016 1 71 71 1943 1943 52 52 40 40 88 88 12 12 182192 1 72 72 893 893 32 32 23 23 25 25 6 6 73566 1 73 73 1762 1762 62 62 40 40 83 83 7 7 194979 1 74 74 1403 1403 45 45 28 28 113 113 2 2 167488 1 75 75 1425 1425 46 46 34 34 118 118 0 0 143756 1 76 76 1857 1857 63 63 33 33 110 110 4 4 275541 1 77 77 1840 1840 75 75 28 28 129 129 3 3 243199 1 78 78 1502 1502 88 88 34 34 51 51 6 6 182999 1 79 79 1441 1441 46 46 30 30 93 93 2 2 135649 1 80 80 1420 1420 53 53 33 33 76 76 0 0 152299 1 81 81 1416 1416 37 37 22 22 49 49 1 1 120221 1 82 82 2970 2970 90 90 38 38 118 118 0 0 346485 1 83 83 1317 1317 63 63 26 26 38 38 5 5 145790 1 84 84 1644 1644 78 78 35 35 141 141 2 2 193339 1 85 85 870 870 25 25 8 8 58 58 0 0 80953 1 86 86 1654 1654 45 45 24 24 27 27 0 0 122774 1 87 87 1054 1054 46 46 29 29 91 91 5 5 130585 1 88 88 3004 3004 144 144 29 29 63 63 1 1 286468 1 89 89 2008 2008 82 82 45 45 56 56 0 0 241066 1 90 90 2547 2547 91 91 37 37 144 144 1 1 148446 1 91 91 1885 1885 71 71 33 33 73 73 1 1 204713 1 92 92 1626 1626 63 63 33 33 168 168 2 2 182079 1 93 93 1468 1468 53 53 25 25 64 64 6 6 140344 1 94 94 2445 2445 62 62 32 32 97 97 1 1 220516 1 95 95 1964 1964 63 63 29 29 117 117 4 4 243060 1 96 96 1381 1381 32 32 28 28 100 100 2 2 162765 1 97 97 1369 1369 39 39 28 28 149 149 3 3 182613 1 98 98 1659 1659 62 62 31 31 187 187 0 0 232138 1 99 99 2888 2888 117 117 52 52 127 127 10 10 265318 1 100 100 2845 2845 92 92 24 24 245 245 9 9 310839 1 101 101 1982 1982 93 93 41 41 87 87 7 7 225060 1 102 102 1904 1904 54 54 33 33 177 177 0 0 232317 1 103 103 1391 1391 144 144 32 32 49 49 0 0 144966 1 104 104 602 602 14 14 19 19 49 49 4 4 43287 1 105 105 1743 1743 61 61 20 20 73 73 4 4 155754 1 106 106 1559 1559 109 109 31 31 177 177 0 0 164709 1 107 107 2014 2014 38 38 31 31 94 94 0 0 201940 1 108 108 2143 2143 73 73 32 32 117 117 0 0 235454 1 109 109 874 874 50 50 23 23 55 55 0 0 99466 1 110 110 1281 1281 72 72 30 30 58 58 0 0 100750 1 111 111 1401 1401 50 50 31 31 95 95 4 4 224549 1 112 112 1944 1944 71 71 42 42 129 129 0 0 243511 1 113 113 391 391 10 10 1 1 11 11 0 0 22938 1 114 114 1605 1605 65 65 32 32 101 101 0 0 152474 1 115 115 530 530 25 25 11 11 28 28 4 4 61857 1 116 116 1386 1386 41 41 36 36 89 89 0 0 132487 1 117 117 2395 2395 86 86 31 31 193 193 1 1 317394 1 118 118 387 387 16 16 0 0 4 4 0 0 21054 1 119 119 1742 1742 42 42 24 24 84 84 5 5 209641 1 120 120 449 449 19 19 8 8 39 39 0 0 31414 1 121 121 2699 2699 95 95 33 33 101 101 2 2 244749 1 122 122 1606 1606 49 49 40 40 82 82 7 7 184510 1 123 123 1204 1204 64 64 38 38 36 36 8 8 128423 1 124 124 1138 1138 38 38 24 24 75 75 2 2 97839 1 125 125 568 568 34 34 8 8 16 16 0 0 38214 1 126 126 1459 1459 32 32 35 35 55 55 2 2 151101 1 127 127 2158 2158 65 65 43 43 131 131 0 0 272458 1 128 128 1111 1111 52 52 43 43 131 131 0 0 172494 1 129 129 2833 2833 65 65 41 41 144 144 3 3 328107 1 130 130 1955 1955 83 83 38 38 139 139 0 0 250579 1 131 131 2922 2922 95 95 45 45 211 211 3 3 351067 1 132 132 1002 1002 29 29 31 31 78 78 0 0 158015 1 133 133 956 956 33 33 28 28 39 39 0 0 85439 1 134 134 2186 2186 247 247 31 31 90 90 4 4 229242 1 135 135 3604 3604 139 139 40 40 166 166 4 4 351619 1 136 136 1035 1035 29 29 30 30 12 12 11 11 84207 1 137 137 3261 3261 110 110 37 37 133 133 0 0 324598 1 138 138 1587 1587 67 67 30 30 69 69 4 4 131069 1 139 139 1424 1424 42 42 35 35 119 119 0 0 204271 1 140 140 1701 1701 65 65 32 32 119 119 1 1 165543 1 141 141 1249 1249 94 94 27 27 65 65 0 0 141722 1 142 142 3352 3352 95 95 31 31 101 101 9 9 299775 1 143 143 1641 1641 67 67 31 31 196 196 1 1 195838 1 144 144 2035 2035 63 63 21 21 15 15 3 3 173260 1 145 145 2312 2312 83 83 39 39 136 136 10 10 254488 1 146 146 1369 1369 45 45 41 41 89 89 5 5 104389 1 147 147 2201 2201 70 70 32 32 123 123 2 2 199476 1 148 148 1900 1900 83 83 39 39 163 163 1 1 224330 1 149 149 207 207 10 10 0 0 5 5 0 0 14688 1 150 150 1645 1645 70 70 30 30 96 96 2 2 181633 1 151 151 2429 2429 103 103 37 37 151 151 1 1 271856 1 152 152 151 151 5 5 0 0 6 6 0 0 7199 1 153 153 474 474 20 20 5 5 13 13 0 0 46660 1 154 154 141 141 5 5 1 1 3 3 0 0 17547 1 155 155 872 872 34 34 32 32 23 23 0 0 95227 1 156 156 1318 1318 48 48 24 24 57 57 2 2 152601 0 157 0 1018 0 40 0 17 0 14 0 0 0 98146 0 158 0 1383 0 43 0 11 0 43 0 3 0 79619 0 159 0 1314 0 31 0 24 0 20 0 6 0 59194 0 160 0 1335 0 42 0 22 0 72 0 0 0 139942 0 161 0 1403 0 46 0 12 0 87 0 2 0 118612 0 162 0 910 0 33 0 19 0 21 0 0 0 72880 0 163 0 616 0 18 0 13 0 56 0 2 0 65475 0 164 0 1407 0 55 0 17 0 59 0 1 0 99643 0 165 0 771 0 35 0 15 0 82 0 1 0 71965 0 166 0 766 0 59 0 16 0 43 0 2 0 77272 0 167 0 473 0 19 0 24 0 25 0 1 0 49289 0 168 0 1376 0 66 0 15 0 38 0 0 0 135131 0 169 0 1232 0 60 0 17 0 25 0 1 0 108446 0 170 0 1521 0 36 0 18 0 38 0 3 0 89746 0 171 0 572 0 25 0 20 0 12 0 0 0 44296 0 172 0 1059 0 47 0 16 0 29 0 0 0 77648 0 173 0 1544 0 54 0 16 0 47 0 0 0 181528 0 174 0 1230 0 53 0 18 0 45 0 0 0 134019 0 175 0 1206 0 40 0 22 0 40 0 1 0 124064 0 176 0 1205 0 40 0 8 0 30 0 4 0 92630 0 177 0 1255 0 39 0 17 0 41 0 0 0 121848 0 178 0 613 0 14 0 18 0 25 0 0 0 52915 0 179 0 721 0 45 0 16 0 23 0 0 0 81872 0 180 0 1109 0 36 0 23 0 14 0 7 0 58981 0 181 0 740 0 28 0 22 0 16 0 2 0 53515 0 182 0 1126 0 44 0 13 0 26 0 0 0 60812 0 183 0 728 0 30 0 13 0 21 0 7 0 56375 0 184 0 689 0 22 0 16 0 27 0 3 0 65490 0 185 0 592 0 17 0 16 0 9 0 0 0 80949 0 186 0 995 0 31 0 20 0 33 0 0 0 76302 0 187 0 1613 0 55 0 22 0 42 0 6 0 104011 0 188 0 2048 0 54 0 17 0 68 0 2 0 98104 0 189 0 705 0 21 0 18 0 32 0 0 0 67989 0 190 0 301 0 14 0 17 0 6 0 0 0 30989 0 191 0 1803 0 81 0 12 0 67 0 3 0 135458 0 192 0 799 0 35 0 7 0 33 0 0 0 73504 0 193 0 861 0 43 0 17 0 77 0 1 0 63123 0 194 0 1186 0 46 0 14 0 46 0 1 0 61254 0 195 0 1451 0 30 0 23 0 30 0 0 0 74914 0 196 0 628 0 23 0 17 0 0 0 1 0 31774 0 197 0 1161 0 38 0 14 0 36 0 0 0 81437 0 198 0 1463 0 54 0 15 0 46 0 0 0 87186 0 199 0 742 0 20 0 17 0 18 0 0 0 50090 0 200 0 979 0 53 0 21 0 48 0 0 0 65745 0 201 0 675 0 45 0 18 0 29 0 0 0 56653 0 202 0 1241 0 39 0 18 0 28 0 0 0 158399 0 203 0 676 0 20 0 17 0 34 0 0 0 46455 0 204 0 1049 0 24 0 17 0 33 0 0 0 73624 0 205 0 620 0 31 0 16 0 34 0 0 0 38395 0 206 0 1081 0 35 0 15 0 33 0 0 0 91899 0 207 0 1688 0 151 0 21 0 80 0 0 0 139526 0 208 0 736 0 52 0 16 0 32 0 0 0 52164 0 209 0 617 0 30 0 14 0 30 0 2 0 51567 0 210 0 812 0 31 0 15 0 41 0 0 0 70551 0 211 0 1051 0 29 0 17 0 41 0 1 0 84856 0 212 0 1656 0 57 0 15 0 51 0 1 0 102538 0 213 0 705 0 40 0 15 0 18 0 0 0 86678 0 214 0 945 0 44 0 10 0 34 0 0 0 85709 0 215 0 554 0 25 0 6 0 31 0 0 0 34662 0 216 0 1597 0 77 0 22 0 39 0 0 0 150580 0 217 0 982 0 35 0 21 0 54 0 0 0 99611 0 218 0 222 0 11 0 1 0 14 0 0 0 19349 0 219 0 1212 0 63 0 18 0 24 0 1 0 99373 0 220 0 1143 0 44 0 17 0 24 0 0 0 86230 0 221 0 435 0 19 0 4 0 8 0 0 0 30837 0 222 0 532 0 13 0 10 0 26 0 0 0 31706 0 223 0 882 0 42 0 16 0 19 0 0 0 89806 0 224 0 608 0 38 0 16 0 11 0 1 0 62088 0 225 0 459 0 29 0 9 0 14 0 0 0 40151 0 226 0 578 0 20 0 16 0 1 0 0 0 27634 0 227 0 826 0 27 0 17 0 39 0 0 0 76990 0 228 0 509 0 20 0 7 0 5 0 0 0 37460 0 229 0 717 0 19 0 15 0 37 0 0 0 54157 0 230 0 637 0 37 0 14 0 32 0 0 0 49862 0 231 0 857 0 26 0 14 0 38 0 0 0 84337 0 232 0 830 0 42 0 18 0 47 0 0 0 64175 0 233 0 652 0 49 0 12 0 47 0 0 0 59382 0 234 0 707 0 30 0 16 0 37 0 0 0 119308 0 235 0 954 0 49 0 21 0 51 0 0 0 76702 0 236 0 1461 0 67 0 19 0 45 0 1 0 103425 0 237 0 672 0 28 0 16 0 21 0 0 0 70344 0 238 0 778 0 19 0 1 0 1 0 0 0 43410 0 239 0 1141 0 49 0 16 0 42 0 1 0 104838 0 240 0 680 0 27 0 10 0 26 0 0 0 62215 0 241 0 1090 0 30 0 19 0 21 0 6 0 69304 0 242 0 616 0 22 0 12 0 4 0 3 0 53117 0 243 0 285 0 12 0 2 0 10 0 1 0 19764 0 244 0 1145 0 31 0 14 0 43 0 2 0 86680 0 245 0 733 0 20 0 17 0 34 0 0 0 84105 0 246 0 888 0 20 0 19 0 31 0 0 0 77945 0 247 0 849 0 39 0 14 0 19 0 0 0 89113 0 248 0 1182 0 29 0 11 0 34 0 3 0 91005 0 249 0 528 0 16 0 4 0 6 0 1 0 40248 0 250 0 642 0 27 0 16 0 11 0 0 0 64187 0 251 0 947 0 21 0 20 0 24 0 0 0 50857 0 252 0 819 0 19 0 12 0 16 0 1 0 56613 0 253 0 757 0 35 0 15 0 72 0 0 0 62792 0 254 0 894 0 14 0 16 0 21 0 0 0 72535 0 255 0 1192 0 63 0 11 0 28 0 0 0 101645 0 256 0 186 0 17 0 0 0 0 0 0 0 7176 0 257 0 1793 0 76 0 17 0 47 0 0 0 96560 0 258 0 829 0 34 0 13 0 32 0 0 0 101011 0 259 0 2702 0 107 0 20 0 65 0 0 0 175824 0 260 0 4691 0 168 0 21 0 123 0 1 0 341570 0 261 0 1112 0 43 0 16 0 26 0 1 0 103597 0 262 0 937 0 41 0 20 0 48 0 0 0 112611 0 263 0 1290 0 34 0 21 0 37 0 0 0 85574 0 264 0 2146 0 75 0 18 0 60 0 1 0 220801 0 265 0 1590 0 61 0 17 0 39 0 1 0 92661 0 266 0 1590 0 55 0 20 0 64 0 0 0 133328 0 267 0 1210 0 77 0 12 0 26 0 0 0 61361 0 268 0 2072 0 75 0 17 0 64 0 4 0 125930 0 269 0 834 0 32 0 10 0 25 0 4 0 82316 0 270 0 1105 0 53 0 13 0 26 0 3 0 102010 0 271 0 1272 0 42 0 22 0 76 0 0 0 101523 0 272 0 761 0 35 0 9 0 2 0 5 0 41566 0 273 0 1988 0 66 0 25 0 36 0 0 0 99923 0 274 0 620 0 19 0 13 0 23 0 0 0 22648 0 275 0 800 0 45 0 13 0 14 0 0 0 46698 0 276 0 1684 0 65 0 19 0 78 0 0 0 131698 0 277 0 1050 0 35 0 18 0 14 0 0 0 91735 0 278 0 1639 0 36 0 16 0 56 0 1 0 133368 0 279 0 1254 0 31 0 14 0 29 0 2 0 135781 0 280 0 1335 0 67 0 7 0 35 0 4 0 74408 0 281 0 1597 0 66 0 17 0 13 0 0 0 81240 0 282 0 961 0 32 0 18 0 21 0 0 0 92499 0 283 0 1577 0 30 0 13 0 40 0 0 0 136084 0 284 0 946 0 64 0 20 0 61 0 0 0 116048 0 285 0 1926 0 81 0 18 0 49 0 0 0 250047 0 286 0 1417 0 118 0 16 0 57 0 0 0 120445 0 287 0 1060 0 18 0 13 0 50 0 0 0 98866 0 288 0 1421 0 62 0 14 0 39 0 1 0 108043 0 289 0 1502 0 37 0 22 0 24 0 1 0 79863
Names of X columns:
Pop t pop_t Pageviews Pageviews_p Logins Logins_p Compendiums_Reviewed Compendiums_Reviewed_p Totblogs Totblogs_p Shared_Compendiums Shared_Compendiums_p Time_in_rfc
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