Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 1 1 210907 56 56 79 79 30 30 94 94 112285 112285 1 2 2 120982 56 56 58 58 28 28 103 103 84786 84786 1 3 3 176508 54 54 60 60 38 38 93 93 83123 83123 1 4 4 179321 89 89 108 108 30 30 103 103 101193 101193 1 5 5 123185 40 40 49 49 22 22 51 51 38361 38361 1 6 6 52746 25 25 0 0 26 26 70 70 68504 68504 1 7 7 385534 92 92 121 121 25 25 91 91 119182 119182 1 8 8 33170 18 18 1 1 18 18 22 22 22807 22807 1 9 9 149061 44 44 43 43 26 26 93 93 116174 116174 1 10 10 165446 33 33 69 69 25 25 60 60 57635 57635 1 11 11 237213 84 84 78 78 38 38 123 123 66198 66198 1 12 12 173326 88 88 86 86 44 44 148 148 71701 71701 1 13 13 133131 55 55 44 44 30 30 90 90 57793 57793 1 14 14 258873 60 60 104 104 40 40 124 124 80444 80444 1 15 15 180083 66 66 63 63 34 34 70 70 53855 53855 1 16 16 324799 154 154 158 158 47 47 168 168 97668 97668 1 17 17 230964 53 53 102 102 30 30 115 115 133824 133824 1 18 18 236785 119 119 77 77 31 31 71 71 101481 101481 1 19 19 135473 41 41 82 82 23 23 66 66 99645 99645 1 20 20 202925 61 61 115 115 36 36 134 134 114789 114789 1 21 21 215147 58 58 101 101 36 36 117 117 99052 99052 1 22 22 344297 75 75 80 80 30 30 108 108 67654 67654 1 23 23 153935 33 33 50 50 25 25 84 84 65553 65553 1 24 24 132943 40 40 83 83 39 39 156 156 97500 97500 1 25 25 174724 92 92 123 123 34 34 120 120 69112 69112 1 26 26 174415 100 100 73 73 31 31 114 114 82753 82753 1 27 27 225548 112 112 81 81 31 31 94 94 85323 85323 1 28 28 223632 73 73 105 105 33 33 120 120 72654 72654 1 29 29 124817 40 40 47 47 25 25 81 81 30727 30727 1 30 30 221698 45 45 105 105 33 33 110 110 77873 77873 1 31 31 210767 60 60 94 94 35 35 133 133 117478 117478 1 32 32 170266 62 62 44 44 42 42 122 122 74007 74007 1 33 33 260561 75 75 114 114 43 43 158 158 90183 90183 1 34 34 84853 31 31 38 38 30 30 109 109 61542 61542 1 35 35 294424 77 77 107 107 33 33 124 124 101494 101494 1 36 36 215641 46 46 71 71 32 32 92 92 55813 55813 1 37 37 325107 99 99 84 84 36 36 126 126 79215 79215 1 38 38 167542 66 66 59 59 28 28 70 70 55461 55461 1 39 39 106408 30 30 33 33 14 14 37 37 31081 31081 1 40 40 265769 146 146 96 96 32 32 120 120 83122 83122 1 41 41 269651 67 67 106 106 30 30 93 93 70106 70106 1 42 42 149112 56 56 56 56 35 35 95 95 60578 60578 1 43 43 152871 58 58 59 59 28 28 90 90 79892 79892 1 44 44 111665 34 34 39 39 28 28 80 80 49810 49810 1 45 45 116408 61 61 34 34 39 39 31 31 71570 71570 1 46 46 362301 119 119 76 76 34 34 110 110 100708 100708 1 47 47 78800 42 42 20 20 26 26 66 66 33032 33032 1 48 48 183167 66 66 91 91 39 39 138 138 82875 82875 1 49 49 277965 89 89 115 115 39 39 133 133 139077 139077 1 50 50 150629 44 44 85 85 33 33 113 113 71595 71595 1 51 51 168809 66 66 76 76 28 28 100 100 72260 72260 1 52 52 24188 24 24 8 8 4 4 7 7 5950 5950 1 53 53 329267 259 259 79 79 39 39 140 140 115762 115762 1 54 54 65029 17 17 21 21 18 18 61 61 32551 32551 1 55 55 101097 64 64 30 30 14 14 41 41 31701 31701 1 56 56 218946 41 41 76 76 29 29 96 96 80670 80670 1 57 57 244052 68 68 101 101 44 44 164 164 143558 143558 1 58 58 233328 132 132 92 92 28 28 102 102 120733 120733 1 59 59 256462 105 105 123 123 35 35 124 124 105195 105195 1 60 60 206161 71 71 75 75 28 28 99 99 73107 73107 1 61 61 311473 112 112 128 128 38 38 129 129 132068 132068 1 62 62 235800 94 94 105 105 23 23 62 62 149193 149193 1 63 63 177939 82 82 55 55 36 36 73 73 46821 46821 1 64 64 207176 70 70 56 56 32 32 114 114 87011 87011 1 65 65 196553 57 57 41 41 29 29 99 99 95260 95260 1 66 66 174184 53 53 72 72 25 25 70 70 55183 55183 1 67 67 143246 103 103 67 67 27 27 104 104 106671 106671 1 68 68 187559 121 121 75 75 36 36 116 116 73511 73511 1 69 69 187681 62 62 114 114 28 28 91 91 92945 92945 1 70 70 119016 52 52 118 118 23 23 74 74 78664 78664 1 71 71 182192 52 52 77 77 40 40 138 138 70054 70054 1 72 72 73566 32 32 22 22 23 23 67 67 22618 22618 1 73 73 194979 62 62 66 66 40 40 151 151 74011 74011 1 74 74 167488 45 45 69 69 28 28 72 72 83737 83737 1 75 75 143756 46 46 105 105 34 34 120 120 69094 69094 1 76 76 275541 63 63 116 116 33 33 115 115 93133 93133 1 77 77 243199 75 75 88 88 28 28 105 105 95536 95536 1 78 78 182999 88 88 73 73 34 34 104 104 225920 225920 1 79 79 135649 46 46 99 99 30 30 108 108 62133 62133 1 80 80 152299 53 53 62 62 33 33 98 98 61370 61370 1 81 81 120221 37 37 53 53 22 22 69 69 43836 43836 1 82 82 346485 90 90 118 118 38 38 111 111 106117 106117 1 83 83 145790 63 63 30 30 26 26 99 99 38692 38692 1 84 84 193339 78 78 100 100 35 35 71 71 84651 84651 1 85 85 80953 25 25 49 49 8 8 27 27 56622 56622 1 86 86 122774 45 45 24 24 24 24 69 69 15986 15986 1 87 87 130585 46 46 67 67 29 29 107 107 95364 95364 1 88 88 286468 144 144 57 57 29 29 107 107 89691 89691 1 89 89 241066 82 82 75 75 45 45 93 93 67267 67267 1 90 90 148446 91 91 135 135 37 37 129 129 126846 126846 1 91 91 204713 71 71 68 68 33 33 69 69 41140 41140 1 92 92 182079 63 63 124 124 33 33 118 118 102860 102860 1 93 93 140344 53 53 33 33 25 25 73 73 51715 51715 1 94 94 220516 62 62 98 98 32 32 119 119 55801 55801 1 95 95 243060 63 63 58 58 29 29 104 104 111813 111813 1 96 96 162765 32 32 68 68 28 28 107 107 120293 120293 1 97 97 182613 39 39 81 81 28 28 99 99 138599 138599 1 98 98 232138 62 62 131 131 31 31 90 90 161647 161647 1 99 99 265318 117 117 110 110 52 52 197 197 115929 115929 1 100 100 310839 92 92 130 130 24 24 85 85 162901 162901 1 101 101 225060 93 93 93 93 41 41 139 139 109825 109825 1 102 102 232317 54 54 118 118 33 33 106 106 129838 129838 1 103 103 144966 144 144 39 39 32 32 50 50 37510 37510 1 104 104 43287 14 14 13 13 19 19 64 64 43750 43750 1 105 105 155754 61 61 74 74 20 20 31 31 40652 40652 1 106 106 164709 109 109 81 81 31 31 63 63 87771 87771 1 107 107 201940 38 38 109 109 31 31 92 92 85872 85872 1 108 108 235454 73 73 151 151 32 32 106 106 89275 89275 1 109 109 99466 50 50 28 28 23 23 69 69 192565 192565 1 110 110 100750 72 72 83 83 30 30 93 93 140867 140867 1 111 111 224549 50 50 54 54 31 31 114 114 120662 120662 1 112 112 243511 71 71 133 133 42 42 110 110 101338 101338 1 113 113 22938 10 10 12 12 1 1 0 0 1168 1168 1 114 114 152474 65 65 106 106 32 32 83 83 65567 65567 1 115 115 61857 25 25 23 23 11 11 30 30 25162 25162 1 116 116 132487 41 41 71 71 36 36 98 98 40735 40735 1 117 117 317394 86 86 116 116 31 31 82 82 91413 91413 1 118 118 21054 16 16 4 4 0 0 0 0 855 855 1 119 119 209641 42 42 62 62 24 24 60 60 97068 97068 1 120 120 31414 19 19 18 18 8 8 9 9 14116 14116 1 121 121 244749 95 95 98 98 33 33 115 115 76643 76643 1 122 122 184510 49 49 64 64 40 40 140 140 110681 110681 1 123 123 128423 64 64 32 32 38 38 120 120 92696 92696 1 124 124 97839 38 38 25 25 24 24 66 66 94785 94785 1 125 125 38214 34 34 16 16 8 8 21 21 8773 8773 1 126 126 151101 32 32 48 48 35 35 124 124 83209 83209 1 127 127 272458 65 65 100 100 43 43 152 152 93815 93815 1 128 128 172494 52 52 46 46 43 43 139 139 86687 86687 1 129 129 328107 65 65 129 129 41 41 144 144 105547 105547 1 130 130 250579 83 83 130 130 38 38 120 120 103487 103487 1 131 131 351067 95 95 136 136 45 45 160 160 213688 213688 1 132 132 158015 29 29 59 59 31 31 114 114 71220 71220 1 133 133 85439 33 33 32 32 28 28 78 78 56926 56926 1 134 134 229242 247 247 63 63 31 31 119 119 91721 91721 1 135 135 351619 139 139 95 95 40 40 141 141 115168 115168 1 136 136 84207 29 29 14 14 30 30 101 101 111194 111194 1 137 137 324598 110 110 113 113 37 37 133 133 135777 135777 1 138 138 131069 67 67 47 47 30 30 83 83 51513 51513 1 139 139 204271 42 42 92 92 35 35 116 116 74163 74163 1 140 140 165543 65 65 70 70 32 32 90 90 51633 51633 1 141 141 141722 94 94 19 19 27 27 36 36 75345 75345 1 142 142 299775 95 95 91 91 31 31 97 97 98952 98952 1 143 143 195838 67 67 111 111 31 31 98 98 102372 102372 1 144 144 173260 63 63 41 41 21 21 78 78 37238 37238 1 145 145 254488 83 83 120 120 39 39 117 117 103772 103772 1 146 146 104389 45 45 135 135 41 41 148 148 123969 123969 1 147 147 199476 70 70 87 87 32 32 105 105 135400 135400 1 148 148 224330 83 83 131 131 39 39 132 132 130115 130115 1 149 149 14688 10 10 4 4 0 0 0 0 6023 6023 1 150 150 181633 70 70 47 47 30 30 73 73 64466 64466 1 151 151 271856 103 103 109 109 37 37 86 86 54990 54990 1 152 152 7199 5 5 7 7 0 0 0 0 1644 1644 1 153 153 46660 20 20 12 12 5 5 13 13 6179 6179 1 154 154 17547 5 5 0 0 1 1 4 4 3926 3926 1 155 155 95227 34 34 37 37 32 32 48 48 34777 34777 1 156 156 152601 48 48 46 46 24 24 46 46 73224 73224 0 157 0 101645 63 0 20 0 11 0 38 0 17140 0 0 158 0 101011 34 0 30 0 13 0 39 0 27570 0 0 159 0 7176 17 0 0 0 0 0 0 0 1423 0 0 160 0 96560 76 0 42 0 17 0 38 0 22996 0 0 161 0 175824 107 0 57 0 20 0 77 0 39992 0 0 162 0 341570 168 0 94 0 21 0 78 0 117105 0 0 163 0 103597 43 0 27 0 16 0 49 0 23789 0 0 164 0 112611 41 0 46 0 20 0 73 0 26706 0 0 165 0 85574 34 0 37 0 21 0 36 0 24266 0 0 166 0 220801 75 0 51 0 18 0 63 0 44418 0 0 167 0 92661 61 0 40 0 17 0 41 0 35232 0 0 168 0 133328 55 0 56 0 20 0 56 0 40909 0 0 169 0 61361 77 0 27 0 12 0 25 0 13294 0 0 170 0 125930 75 0 37 0 17 0 65 0 32387 0 0 171 0 82316 32 0 27 0 10 0 38 0 21233 0 0 172 0 102010 53 0 28 0 13 0 44 0 44332 0 0 173 0 101523 42 0 59 0 22 0 87 0 61056 0 0 174 0 41566 35 0 0 0 9 0 27 0 13497 0 0 175 0 99923 66 0 44 0 25 0 80 0 32334 0 0 176 0 22648 19 0 12 0 13 0 28 0 44339 0 0 177 0 46698 45 0 14 0 13 0 33 0 10288 0 0 178 0 131698 65 0 60 0 19 0 59 0 65622 0 0 179 0 91735 35 0 7 0 18 0 49 0 16563 0 0 180 0 79863 37 0 29 0 22 0 49 0 29011 0 0 181 0 108043 62 0 45 0 14 0 38 0 34553 0 0 182 0 98866 18 0 25 0 13 0 39 0 23517 0 0 183 0 120445 118 0 36 0 16 0 56 0 51009 0 0 184 0 116048 64 0 50 0 20 0 50 0 33416 0 0 185 0 250047 81 0 41 0 18 0 61 0 83305 0 0 186 0 136084 30 0 27 0 13 0 41 0 27142 0 0 187 0 92499 32 0 25 0 18 0 55 0 21399 0 0 188 0 135781 31 0 45 0 14 0 44 0 24874 0 0 189 0 74408 67 0 29 0 7 0 21 0 34988 0 0 190 0 81240 66 0 58 0 17 0 50 0 45549 0 0 191 0 133368 36 0 37 0 16 0 57 0 32755 0 0 192 0 98146 40 0 15 0 17 0 48 0 27114 0 0 193 0 79619 43 0 42 0 11 0 32 0 20760 0 0 194 0 59194 31 0 7 0 24 0 68 0 37636 0 0 195 0 139942 42 0 54 0 22 0 87 0 65461 0 0 196 0 118612 46 0 54 0 12 0 43 0 30080 0 0 197 0 72880 33 0 14 0 19 0 67 0 24094 0 0 198 0 65475 18 0 16 0 13 0 46 0 69008 0 0 199 0 99643 55 0 33 0 17 0 46 0 54968 0 0 200 0 71965 35 0 32 0 15 0 56 0 46090 0 0 201 0 77272 59 0 21 0 16 0 48 0 27507 0 0 202 0 49289 19 0 15 0 24 0 44 0 10672 0 0 203 0 135131 66 0 38 0 15 0 60 0 34029 0 0 204 0 108446 60 0 22 0 17 0 65 0 46300 0 0 205 0 89746 36 0 28 0 18 0 55 0 24760 0 0 206 0 44296 25 0 10 0 20 0 38 0 18779 0 0 207 0 77648 47 0 31 0 16 0 52 0 21280 0 0 208 0 181528 54 0 32 0 16 0 60 0 40662 0 0 209 0 134019 53 0 32 0 18 0 54 0 28987 0 0 210 0 124064 40 0 43 0 22 0 86 0 22827 0 0 211 0 92630 40 0 27 0 8 0 24 0 18513 0 0 212 0 121848 39 0 37 0 17 0 52 0 30594 0 0 213 0 52915 14 0 20 0 18 0 49 0 24006 0 0 214 0 81872 45 0 32 0 16 0 61 0 27913 0 0 215 0 58981 36 0 0 0 23 0 61 0 42744 0 0 216 0 53515 28 0 5 0 22 0 81 0 12934 0 0 217 0 60812 44 0 26 0 13 0 43 0 22574 0 0 218 0 56375 30 0 10 0 13 0 40 0 41385 0 0 219 0 65490 22 0 27 0 16 0 40 0 18653 0 0 220 0 80949 17 0 11 0 16 0 56 0 18472 0 0 221 0 76302 31 0 29 0 20 0 68 0 30976 0 0 222 0 104011 55 0 25 0 22 0 79 0 63339 0 0 223 0 98104 54 0 55 0 17 0 47 0 25568 0 0 224 0 67989 21 0 23 0 18 0 57 0 33747 0 0 225 0 30989 14 0 5 0 17 0 41 0 4154 0 0 226 0 135458 81 0 43 0 12 0 29 0 19474 0 0 227 0 73504 35 0 23 0 7 0 3 0 35130 0 0 228 0 63123 43 0 34 0 17 0 60 0 39067 0 0 229 0 61254 46 0 36 0 14 0 30 0 13310 0 0 230 0 74914 30 0 35 0 23 0 79 0 65892 0 0 231 0 31774 23 0 0 0 17 0 47 0 4143 0 0 232 0 81437 38 0 37 0 14 0 40 0 28579 0 0 233 0 87186 54 0 28 0 15 0 48 0 51776 0 0 234 0 50090 20 0 16 0 17 0 36 0 21152 0 0 235 0 65745 53 0 26 0 21 0 42 0 38084 0 0 236 0 56653 45 0 38 0 18 0 49 0 27717 0 0 237 0 158399 39 0 23 0 18 0 57 0 32928 0 0 238 0 46455 20 0 22 0 17 0 12 0 11342 0 0 239 0 73624 24 0 30 0 17 0 40 0 19499 0 0 240 0 38395 31 0 16 0 16 0 43 0 16380 0 0 241 0 91899 35 0 18 0 15 0 33 0 36874 0 0 242 0 139526 151 0 28 0 21 0 77 0 48259 0 0 243 0 52164 52 0 32 0 16 0 43 0 16734 0 0 244 0 51567 30 0 21 0 14 0 45 0 28207 0 0 245 0 70551 31 0 23 0 15 0 47 0 30143 0 0 246 0 84856 29 0 29 0 17 0 43 0 41369 0 0 247 0 102538 57 0 50 0 15 0 45 0 45833 0 0 248 0 86678 40 0 12 0 15 0 50 0 29156 0 0 249 0 85709 44 0 21 0 10 0 35 0 35944 0 0 250 0 34662 25 0 18 0 6 0 7 0 36278 0 0 251 0 150580 77 0 27 0 22 0 71 0 45588 0 0 252 0 99611 35 0 41 0 21 0 67 0 45097 0 0 253 0 19349 11 0 13 0 1 0 0 0 3895 0 0 254 0 99373 63 0 12 0 18 0 62 0 28394 0 0 255 0 86230 44 0 21 0 17 0 54 0 18632 0 0 256 0 30837 19 0 8 0 4 0 4 0 2325 0 0 257 0 31706 13 0 26 0 10 0 25 0 25139 0 0 258 0 89806 42 0 27 0 16 0 40 0 27975 0 0 259 0 62088 38 0 13 0 16 0 38 0 14483 0 0 260 0 40151 29 0 16 0 9 0 19 0 13127 0 0 261 0 27634 20 0 2 0 16 0 17 0 5839 0 0 262 0 76990 27 0 42 0 17 0 67 0 24069 0 0 263 0 37460 20 0 5 0 7 0 14 0 3738 0 0 264 0 54157 19 0 37 0 15 0 30 0 18625 0 0 265 0 49862 37 0 17 0 14 0 54 0 36341 0 0 266 0 84337 26 0 38 0 14 0 35 0 24548 0 0 267 0 64175 42 0 37 0 18 0 59 0 21792 0 0 268 0 59382 49 0 29 0 12 0 24 0 26263 0 0 269 0 119308 30 0 32 0 16 0 58 0 23686 0 0 270 0 76702 49 0 35 0 21 0 42 0 49303 0 0 271 0 103425 67 0 17 0 19 0 46 0 25659 0 0 272 0 70344 28 0 20 0 16 0 61 0 28904 0 0 273 0 43410 19 0 7 0 1 0 3 0 2781 0 0 274 0 104838 49 0 46 0 16 0 52 0 29236 0 0 275 0 62215 27 0 24 0 10 0 25 0 19546 0 0 276 0 69304 30 0 40 0 19 0 40 0 22818 0 0 277 0 53117 22 0 3 0 12 0 32 0 32689 0 0 278 0 19764 12 0 10 0 2 0 4 0 5752 0 0 279 0 86680 31 0 37 0 14 0 49 0 22197 0 0 280 0 84105 20 0 17 0 17 0 63 0 20055 0 0 281 0 77945 20 0 28 0 19 0 67 0 25272 0 0 282 0 89113 39 0 19 0 14 0 32 0 82206 0 0 283 0 91005 29 0 29 0 11 0 23 0 32073 0 0 284 0 40248 16 0 8 0 4 0 7 0 5444 0 0 285 0 64187 27 0 10 0 16 0 54 0 20154 0 0 286 0 50857 21 0 15 0 20 0 37 0 36944 0 0 287 0 56613 19 0 15 0 12 0 35 0 8019 0 0 288 0 62792 35 0 28 0 15 0 51 0 30884 0 0 289 0 72535 14 0 17 0 16 0 39 0 19540 0
Names of X columns:
Pop T Pop_t Time_in_rfc Logins Logins_P Blogged_computations Blogged_computations_P Compendiums_reviewed Compendiums_reviewed_P Feedback_messages_p120 Feedback_messages_p120_P Total_number_of_characters Total_number_of_characters_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