Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 1 1 1418 1418 210907 56 56 396 396 3 3 115 115 112285 112285 1 2 2 869 869 120982 56 56 297 297 4 4 109 109 84786 84786 1 3 3 1530 1530 176508 54 54 559 559 12 12 146 146 83123 83123 1 4 4 2172 2172 179321 89 89 967 967 2 2 116 116 101193 101193 1 5 5 901 901 123185 40 40 270 270 1 1 68 68 38361 38361 1 6 6 463 463 52746 25 25 143 143 3 3 101 101 68504 68504 1 7 7 3201 3201 385534 92 92 1562 1562 0 0 96 96 119182 119182 1 8 8 371 371 33170 18 18 109 109 0 0 67 67 22807 22807 1 9 9 1583 1583 149061 44 44 656 656 5 5 100 100 116174 116174 1 10 10 1439 1439 165446 33 33 511 511 0 0 93 93 57635 57635 1 11 11 1764 1764 237213 84 84 655 655 0 0 140 140 66198 66198 1 12 12 1495 1495 173326 88 88 465 465 7 7 166 166 71701 71701 1 13 13 1373 1373 133131 55 55 525 525 7 7 99 99 57793 57793 1 14 14 2187 2187 258873 60 60 885 885 3 3 139 139 80444 80444 1 15 15 1491 1491 180083 66 66 497 497 9 9 130 130 53855 53855 1 16 16 4041 4041 324799 154 154 1436 1436 0 0 181 181 97668 97668 1 17 17 1706 1706 230964 53 53 612 612 4 4 116 116 133824 133824 1 18 18 2152 2152 236785 119 119 865 865 3 3 116 116 101481 101481 1 19 19 1036 1036 135473 41 41 385 385 0 0 88 88 99645 99645 1 20 20 1882 1882 202925 61 61 567 567 7 7 139 139 114789 114789 1 21 21 1929 1929 215147 58 58 639 639 0 0 135 135 99052 99052 1 22 22 2242 2242 344297 75 75 963 963 1 1 108 108 67654 67654 1 23 23 1220 1220 153935 33 33 398 398 5 5 89 89 65553 65553 1 24 24 1289 1289 132943 40 40 410 410 7 7 156 156 97500 97500 1 25 25 2515 2515 174724 92 92 966 966 0 0 129 129 69112 69112 1 26 26 2147 2147 174415 100 100 801 801 0 0 118 118 82753 82753 1 27 27 2352 2352 225548 112 112 892 892 5 5 118 118 85323 85323 1 28 28 1638 1638 223632 73 73 513 513 0 0 125 125 72654 72654 1 29 29 1222 1222 124817 40 40 469 469 0 0 95 95 30727 30727 1 30 30 1812 1812 221698 45 45 683 683 0 0 126 126 77873 77873 1 31 31 1677 1677 210767 60 60 643 643 3 3 135 135 117478 117478 1 32 32 1579 1579 170266 62 62 535 535 4 4 154 154 74007 74007 1 33 33 1731 1731 260561 75 75 625 625 1 1 165 165 90183 90183 1 34 34 807 807 84853 31 31 264 264 4 4 113 113 61542 61542 1 35 35 2452 2452 294424 77 77 992 992 2 2 127 127 101494 101494 1 36 36 1940 1940 215641 46 46 818 818 0 0 121 121 55813 55813 1 37 37 2662 2662 325107 99 99 937 937 0 0 136 136 79215 79215 1 38 38 1499 1499 167542 66 66 507 507 2 2 108 108 55461 55461 1 39 39 865 865 106408 30 30 260 260 1 1 46 46 31081 31081 1 40 40 2527 2527 265769 146 146 927 927 2 2 124 124 83122 83122 1 41 41 2747 2747 269651 67 67 1269 1269 10 10 115 115 70106 70106 1 42 42 1324 1324 149112 56 56 537 537 6 6 128 128 60578 60578 1 43 43 1383 1383 152871 58 58 532 532 5 5 97 97 79892 79892 1 44 44 1179 1179 111665 34 34 345 345 4 4 104 104 49810 49810 1 45 45 2099 2099 116408 61 61 918 918 1 1 59 59 71570 71570 1 46 46 4308 4308 362301 119 119 1635 1635 2 2 125 125 100708 100708 1 47 47 918 918 78800 42 42 330 330 2 2 82 82 33032 33032 1 48 48 1831 1831 183167 66 66 557 557 0 0 149 149 82875 82875 1 49 49 3373 3373 277965 89 89 1178 1178 8 8 149 149 139077 139077 1 50 50 1713 1713 150629 44 44 740 740 3 3 122 122 71595 71595 1 51 51 1438 1438 168809 66 66 452 452 0 0 118 118 72260 72260 1 52 52 496 496 24188 24 24 218 218 0 0 12 12 5950 5950 1 53 53 2253 2253 329267 259 259 764 764 8 8 144 144 115762 115762 1 54 54 744 744 65029 17 17 255 255 5 5 67 67 32551 32551 1 55 55 1161 1161 101097 64 64 454 454 3 3 52 52 31701 31701 1 56 56 2352 2352 218946 41 41 866 866 1 1 108 108 80670 80670 1 57 57 2144 2144 244052 68 68 574 574 5 5 166 166 143558 143558 1 58 58 2694 2694 233328 132 132 825 825 5 5 107 107 120733 120733 1 59 59 1973 1973 256462 105 105 798 798 0 0 127 127 105195 105195 1 60 60 1769 1769 206161 71 71 663 663 12 12 107 107 73107 73107 1 61 61 3148 3148 311473 112 112 1069 1069 8 8 146 146 132068 132068 1 62 62 2474 2474 235800 94 94 921 921 8 8 84 84 149193 149193 1 63 63 2084 2084 177939 82 82 858 858 8 8 141 141 46821 46821 1 64 64 1954 1954 207176 70 70 711 711 8 8 123 123 87011 87011 1 65 65 1226 1226 196553 57 57 503 503 2 2 111 111 95260 95260 1 66 66 1389 1389 174184 53 53 382 382 0 0 98 98 55183 55183 1 67 67 1496 1496 143246 103 103 464 464 5 5 105 105 106671 106671 1 68 68 2269 2269 187559 121 121 717 717 8 8 135 135 73511 73511 1 69 69 1833 1833 187681 62 62 690 690 2 2 107 107 92945 92945 1 70 70 1268 1268 119016 52 52 462 462 5 5 85 85 78664 78664 1 71 71 1943 1943 182192 52 52 657 657 12 12 155 155 70054 70054 1 72 72 893 893 73566 32 32 385 385 6 6 88 88 22618 22618 1 73 73 1762 1762 194979 62 62 577 577 7 7 155 155 74011 74011 1 74 74 1403 1403 167488 45 45 619 619 2 2 104 104 83737 83737 1 75 75 1425 1425 143756 46 46 479 479 0 0 132 132 69094 69094 1 76 76 1857 1857 275541 63 63 817 817 4 4 127 127 93133 93133 1 77 77 1840 1840 243199 75 75 752 752 3 3 108 108 95536 95536 1 78 78 1502 1502 182999 88 88 430 430 6 6 129 129 225920 225920 1 79 79 1441 1441 135649 46 46 451 451 2 2 116 116 62133 62133 1 80 80 1420 1420 152299 53 53 537 537 0 0 122 122 61370 61370 1 81 81 1416 1416 120221 37 37 519 519 1 1 85 85 43836 43836 1 82 82 2970 2970 346485 90 90 1000 1000 0 0 147 147 106117 106117 1 83 83 1317 1317 145790 63 63 637 637 5 5 99 99 38692 38692 1 84 84 1644 1644 193339 78 78 465 465 2 2 87 87 84651 84651 1 85 85 870 870 80953 25 25 437 437 0 0 28 28 56622 56622 1 86 86 1654 1654 122774 45 45 711 711 0 0 90 90 15986 15986 1 87 87 1054 1054 130585 46 46 299 299 5 5 109 109 95364 95364 1 88 88 3004 3004 286468 144 144 1162 1162 1 1 111 111 89691 89691 1 89 89 2008 2008 241066 82 82 714 714 0 0 158 158 67267 67267 1 90 90 2547 2547 148446 91 91 905 905 1 1 141 141 126846 126846 1 91 91 1885 1885 204713 71 71 649 649 1 1 122 122 41140 41140 1 92 92 1626 1626 182079 63 63 512 512 2 2 124 124 102860 102860 1 93 93 1468 1468 140344 53 53 472 472 6 6 93 93 51715 51715 1 94 94 2445 2445 220516 62 62 905 905 1 1 124 124 55801 55801 1 95 95 1964 1964 243060 63 63 786 786 4 4 112 112 111813 111813 1 96 96 1381 1381 162765 32 32 489 489 2 2 108 108 120293 120293 1 97 97 1369 1369 182613 39 39 479 479 3 3 99 99 138599 138599 1 98 98 1659 1659 232138 62 62 617 617 0 0 117 117 161647 161647 1 99 99 2888 2888 265318 117 117 925 925 10 10 199 199 115929 115929 1 100 100 2845 2845 310839 92 92 1144 1144 9 9 91 91 162901 162901 1 101 101 1982 1982 225060 93 93 669 669 7 7 158 158 109825 109825 1 102 102 1904 1904 232317 54 54 707 707 0 0 126 126 129838 129838 1 103 103 1391 1391 144966 144 144 458 458 0 0 122 122 37510 37510 1 104 104 602 602 43287 14 14 214 214 4 4 71 71 43750 43750 1 105 105 1743 1743 155754 61 61 599 599 4 4 75 75 40652 40652 1 106 106 1559 1559 164709 109 109 572 572 0 0 115 115 87771 87771 1 107 107 2014 2014 201940 38 38 897 897 0 0 119 119 85872 85872 1 108 108 2143 2143 235454 73 73 819 819 0 0 124 124 89275 89275 1 109 109 874 874 99466 50 50 273 273 0 0 91 91 192565 192565 1 110 110 1281 1281 100750 72 72 407 407 0 0 119 119 140867 140867 1 111 111 1401 1401 224549 50 50 465 465 4 4 117 117 120662 120662 1 112 112 1944 1944 243511 71 71 603 603 0 0 155 155 101338 101338 1 113 113 391 391 22938 10 10 154 154 0 0 0 0 1168 1168 1 114 114 1605 1605 152474 65 65 577 577 0 0 123 123 65567 65567 1 115 115 530 530 61857 25 25 192 192 4 4 32 32 25162 25162 1 116 116 1386 1386 132487 41 41 411 411 0 0 136 136 40735 40735 1 117 117 2395 2395 317394 86 86 975 975 1 1 117 117 91413 91413 1 118 118 387 387 21054 16 16 146 146 0 0 0 0 855 855 1 119 119 1742 1742 209641 42 42 705 705 5 5 88 88 97068 97068 1 120 120 449 449 31414 19 19 200 200 0 0 25 25 14116 14116 1 121 121 2699 2699 244749 95 95 964 964 2 2 124 124 76643 76643 1 122 122 1606 1606 184510 49 49 537 537 7 7 151 151 110681 110681 1 123 123 1204 1204 128423 64 64 369 369 8 8 145 145 92696 92696 1 124 124 1138 1138 97839 38 38 417 417 2 2 87 87 94785 94785 1 125 125 568 568 38214 34 34 276 276 0 0 27 27 8773 8773 1 126 126 1459 1459 151101 32 32 514 514 2 2 131 131 83209 83209 1 127 127 2158 2158 272458 65 65 822 822 0 0 162 162 93815 93815 1 128 128 1111 1111 172494 52 52 389 389 0 0 165 165 86687 86687 1 129 129 2833 2833 328107 65 65 1255 1255 3 3 159 159 105547 105547 1 130 130 1955 1955 250579 83 83 694 694 0 0 147 147 103487 103487 1 131 131 2922 2922 351067 95 95 1024 1024 3 3 170 170 213688 213688 1 132 132 1002 1002 158015 29 29 400 400 0 0 119 119 71220 71220 1 133 133 956 956 85439 33 33 350 350 0 0 104 104 56926 56926 1 134 134 2186 2186 229242 247 247 719 719 4 4 120 120 91721 91721 1 135 135 3604 3604 351619 139 139 1277 1277 4 4 150 150 115168 115168 1 136 136 1035 1035 84207 29 29 356 356 11 11 112 112 111194 111194 1 137 137 3261 3261 324598 110 110 1402 1402 0 0 136 136 135777 135777 1 138 138 1587 1587 131069 67 67 600 600 4 4 107 107 51513 51513 1 139 139 1424 1424 204271 42 42 480 480 0 0 130 130 74163 74163 1 140 140 1701 1701 165543 65 65 595 595 1 1 115 115 51633 51633 1 141 141 1249 1249 141722 94 94 436 436 0 0 107 107 75345 75345 1 142 142 3352 3352 299775 95 95 1367 1367 9 9 120 120 98952 98952 1 143 143 1641 1641 195838 67 67 564 564 1 1 116 116 102372 102372 1 144 144 2035 2035 173260 63 63 716 716 3 3 79 79 37238 37238 1 145 145 2312 2312 254488 83 83 747 747 10 10 150 150 103772 103772 1 146 146 1369 1369 104389 45 45 467 467 5 5 156 156 123969 123969 1 147 147 2201 2201 199476 70 70 861 861 2 2 118 118 135400 135400 1 148 148 1900 1900 224330 83 83 612 612 1 1 144 144 130115 130115 1 149 149 207 207 14688 10 10 85 85 0 0 0 0 6023 6023 1 150 150 1645 1645 181633 70 70 564 564 2 2 110 110 64466 64466 1 151 151 2429 2429 271856 103 103 824 824 1 1 147 147 54990 54990 1 152 152 151 151 7199 5 5 74 74 0 0 0 0 1644 1644 1 153 153 474 474 46660 20 20 259 259 0 0 15 15 6179 6179 1 154 154 141 141 17547 5 5 69 69 0 0 4 4 3926 3926 1 155 155 872 872 95227 34 34 239 239 0 0 111 111 34777 34777 1 156 156 1318 1318 152601 48 48 438 438 2 2 85 85 73224 73224 0 157 0 1192 0 101645 63 0 371 0 0 0 44 0 17140 0 0 158 0 829 0 101011 34 0 238 0 0 0 52 0 27570 0 0 159 0 186 0 7176 17 0 70 0 0 0 0 0 1423 0 0 160 0 1793 0 96560 76 0 503 0 0 0 54 0 22996 0 0 161 0 2702 0 175824 107 0 910 0 0 0 80 0 39992 0 0 162 0 4691 0 341570 168 0 1276 0 1 0 80 0 117105 0 0 163 0 1112 0 103597 43 0 379 0 1 0 60 0 23789 0 0 164 0 937 0 112611 41 0 248 0 0 0 78 0 26706 0 0 165 0 1290 0 85574 34 0 351 0 0 0 78 0 24266 0 0 166 0 2146 0 220801 75 0 720 0 1 0 72 0 44418 0 0 167 0 1590 0 92661 61 0 508 0 1 0 45 0 35232 0 0 168 0 1590 0 133328 55 0 506 0 0 0 78 0 40909 0 0 169 0 1210 0 61361 77 0 451 0 0 0 39 0 13294 0 0 170 0 2072 0 125930 75 0 699 0 4 0 68 0 32387 0 0 171 0 834 0 82316 32 0 245 0 4 0 39 0 21233 0 0 172 0 1105 0 102010 53 0 370 0 3 0 50 0 44332 0 0 173 0 1272 0 101523 42 0 316 0 0 0 88 0 61056 0 0 174 0 761 0 41566 35 0 229 0 5 0 36 0 13497 0 0 175 0 1988 0 99923 66 0 617 0 0 0 99 0 32334 0 0 176 0 620 0 22648 19 0 184 0 0 0 39 0 44339 0 0 177 0 800 0 46698 45 0 274 0 0 0 52 0 10288 0 0 178 0 1684 0 131698 65 0 502 0 0 0 75 0 65622 0 0 179 0 1050 0 91735 35 0 382 0 0 0 71 0 16563 0 0 180 0 1502 0 79863 37 0 438 0 1 0 71 0 29011 0 0 181 0 1421 0 108043 62 0 466 0 1 0 54 0 34553 0 0 182 0 1060 0 98866 18 0 397 0 0 0 49 0 23517 0 0 183 0 1417 0 120445 118 0 457 0 0 0 59 0 51009 0 0 184 0 946 0 116048 64 0 230 0 0 0 75 0 33416 0 0 185 0 1926 0 250047 81 0 651 0 0 0 71 0 83305 0 0 186 0 1577 0 136084 30 0 671 0 0 0 51 0 27142 0 0 187 0 961 0 92499 32 0 319 0 0 0 71 0 21399 0 0 188 0 1254 0 135781 31 0 433 0 2 0 47 0 24874 0 0 189 0 1335 0 74408 67 0 434 0 4 0 28 0 34988 0 0 190 0 1597 0 81240 66 0 503 0 0 0 68 0 45549 0 0 191 0 1639 0 133368 36 0 535 0 1 0 64 0 32755 0 0 192 0 1018 0 98146 40 0 459 0 0 0 68 0 27114 0 0 193 0 1383 0 79619 43 0 426 0 3 0 40 0 20760 0 0 194 0 1314 0 59194 31 0 288 0 6 0 80 0 37636 0 0 195 0 1335 0 139942 42 0 498 0 0 0 88 0 65461 0 0 196 0 1403 0 118612 46 0 454 0 2 0 48 0 30080 0 0 197 0 910 0 72880 33 0 376 0 0 0 76 0 24094 0
Names of X columns:
Pop t pop_t pageviews pageviews_p time_in_rfc logins logins_p compendium_views_info compendium_views_info_p shared_compendiums shared_compendiums_p feedback_messages_p1 feedback_messages_p1_p totsize totsize_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