Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1418 210907 396 79 30 869 120982 297 58 28 1530 176508 559 60 38 2172 179321 967 108 30 901 123185 270 49 22 463 52746 143 0 26 3201 385534 1562 121 25 371 33170 109 1 18 1192 101645 371 20 11 1583 149061 656 43 26 1439 165446 511 69 25 1764 237213 655 78 38 1495 173326 465 86 44 1373 133131 525 44 30 2187 258873 885 104 40 1491 180083 497 63 34 4041 324799 1436 158 47 1706 230964 612 102 30 2152 236785 865 77 31 1036 135473 385 82 23 1882 202925 567 115 36 1929 215147 639 101 36 2242 344297 963 80 30 1220 153935 398 50 25 1289 132943 410 83 39 2515 174724 966 123 34 2147 174415 801 73 31 2352 225548 892 81 31 1638 223632 513 105 33 1222 124817 469 47 25 1812 221698 683 105 33 1677 210767 643 94 35 1579 170266 535 44 42 1731 260561 625 114 43 807 84853 264 38 30 2452 294424 992 107 33 829 101011 238 30 13 1940 215641 818 71 32 2662 325107 937 84 36 186 7176 70 0 0 1499 167542 507 59 28 865 106408 260 33 14 1793 96560 503 42 17 2527 265769 927 96 32 2747 269651 1269 106 30 1324 149112 537 56 35 2702 175824 910 57 20 1383 152871 532 59 28 1179 111665 345 39 28 2099 116408 918 34 39 4308 362301 1635 76 34 918 78800 330 20 26 1831 183167 557 91 39 3373 277965 1178 115 39 1713 150629 740 85 33 1438 168809 452 76 28 496 24188 218 8 4 2253 329267 764 79 39 744 65029 255 21 18 1161 101097 454 30 14 2352 218946 866 76 29 2144 244052 574 101 44 4691 341570 1276 94 21 1112 103597 379 27 16 2694 233328 825 92 28 1973 256462 798 123 35 1769 206161 663 75 28 3148 311473 1069 128 38 2474 235800 921 105 23 2084 177939 858 55 36 1954 207176 711 56 32 1226 196553 503 41 29 1389 174184 382 72 25 1496 143246 464 67 27 2269 187559 717 75 36 1833 187681 690 114 28 1268 119016 462 118 23 1943 182192 657 77 40 893 73566 385 22 23 1762 194979 577 66 40 1403 167488 619 69 28 1425 143756 479 105 34 1857 275541 817 116 33 1840 243199 752 88 28 1502 182999 430 73 34 1441 135649 451 99 30 1420 152299 537 62 33 1416 120221 519 53 22 2970 346485 1000 118 38 1317 145790 637 30 26 1644 193339 465 100 35 870 80953 437 49 8 1654 122774 711 24 24 1054 130585 299 67 29 937 112611 248 46 20 3004 286468 1162 57 29 2008 241066 714 75 45 2547 148446 905 135 37 1885 204713 649 68 33 1626 182079 512 124 33 1468 140344 472 33 25 2445 220516 905 98 32 1964 243060 786 58 29 1381 162765 489 68 28 1369 182613 479 81 28 1659 232138 617 131 31 2888 265318 925 110 52 1290 85574 351 37 21 2845 310839 1144 130 24 1982 225060 669 93 41 1904 232317 707 118 33 1391 144966 458 39 32 602 43287 214 13 19 1743 155754 599 74 20 1559 164709 572 81 31 2014 201940 897 109 31 2143 235454 819 151 32 2146 220801 720 51 18 874 99466 273 28 23 1590 92661 508 40 17 1590 133328 506 56 20 1210 61361 451 27 12 2072 125930 699 37 17 1281 100750 407 83 30 1401 224549 465 54 31 834 82316 245 27 10 1105 102010 370 28 13 1272 101523 316 59 22 1944 243511 603 133 42 391 22938 154 12 1 761 41566 229 0 9 1605 152474 577 106 32 530 61857 192 23 11 1988 99923 617 44 25 1386 132487 411 71 36 2395 317394 975 116 31 387 21054 146 4 0 1742 209641 705 62 24 620 22648 184 12 13 449 31414 200 18 8 800 46698 274 14 13 1684 131698 502 60 19 1050 91735 382 7 18 2699 244749 964 98 33 1606 184510 537 64 40 1502 79863 438 29 22 1204 128423 369 32 38 1138 97839 417 25 24 568 38214 276 16 8 1459 151101 514 48 35 2158 272458 822 100 43 1111 172494 389 46 43 1421 108043 466 45 14 2833 328107 1255 129 41 1955 250579 694 130 38 2922 351067 1024 136 45 1002 158015 400 59 31 1060 98866 397 25 13 956 85439 350 32 28 2186 229242 719 63 31 3604 351619 1277 95 40 1035 84207 356 14 30 1417 120445 457 36 16 3261 324598 1402 113 37 1587 131069 600 47 30 1424 204271 480 92 35 1701 165543 595 70 32 1249 141722 436 19 27 946 116048 230 50 20 1926 250047 651 41 18 3352 299775 1367 91 31 1641 195838 564 111 31 2035 173260 716 41 21 2312 254488 747 120 39 1369 104389 467 135 41 1577 136084 671 27 13 2201 199476 861 87 32 961 92499 319 25 18 1900 224330 612 131 39 1254 135781 433 45 14 1335 74408 434 29 7 1597 81240 503 58 17 207 14688 85 4 0 1645 181633 564 47 30 2429 271856 824 109 37 151 7199 74 7 0 474 46660 259 12 5 141 17547 69 0 1 1639 133368 535 37 16 872 95227 239 37 32 1318 152601 438 46 24 1018 98146 459 15 17 1383 79619 426 42 11 1314 59194 288 7 24 1335 139942 498 54 22 1403 118612 454 54 12 910 72880 376 14 19 616 65475 225 16 13 1407 99643 555 33 17 771 71965 252 32 15 766 77272 208 21 16 473 49289 130 15 24 1376 135131 481 38 15 1232 108446 389 22 17 1521 89746 565 28 18 572 44296 173 10 20 1059 77648 278 31 16 1544 181528 609 32 16 1230 134019 422 32 18 1206 124064 445 43 22 1205 92630 387 27 8 1255 121848 339 37 17 613 52915 181 20 18 721 81872 245 32 16 1109 58981 384 0 23 740 53515 212 5 22 1126 60812 399 26 13 728 56375 229 10 13 689 65490 224 27 16 592 80949 203 11 16 995 76302 333 29 20 1613 104011 384 25 22 2048 98104 636 55 17 705 67989 185 23 18 301 30989 93 5 17 1803 135458 581 43 12 799 73504 248 23 7 861 63123 304 34 17 1186 61254 344 36 14 1451 74914 407 35 23 628 31774 170 0 17 1161 81437 312 37 14 1463 87186 507 28 15 742 50090 224 16 17 979 65745 340 26 21 675 56653 168 38 18 1241 158399 443 23 18 676 46455 204 22 17 1049 73624 367 30 17 620 38395 210 16 16 1081 91899 335 18 15 1688 139526 364 28 21 736 52164 178 32 16 617 51567 206 21 14 812 70551 279 23 15 1051 84856 387 29 17 1656 102538 490 50 15 705 86678 238 12 15 945 85709 343 21 10 554 34662 232 18 6 1597 150580 530 27 22 982 99611 291 41 21 222 19349 67 13 1 1212 99373 397 12 18 1143 86230 467 21 17 435 30837 178 8 4 532 31706 175 26 10 882 89806 299 27 16 608 62088 154 13 16 459 40151 106 16 9 578 27634 189 2 16 826 76990 194 42 17 509 37460 135 5 7 717 54157 201 37 15 637 49862 207 17 14 857 84337 280 38 14 830 64175 260 37 18 652 59382 227 29 12 707 119308 239 32 16 954 76702 333 35 21 1461 103425 428 17 19 672 70344 230 20 16 778 43410 292 7 1 1141 104838 350 46 16 680 62215 186 24 10 1090 69304 326 40 19 616 53117 155 3 12 285 19764 75 10 2 1145 86680 361 37 14 733 84105 261 17 17 888 77945 299 28 19 849 89113 300 19 14 1182 91005 450 29 11 528 40248 183 8 4 642 64187 238 10 16 947 50857 165 15 20 819 56613 234 15 12 757 62792 176 28 15 894 72535 329 17 16
Names of X columns:
pageviews time_in_rfc compendium_views_info blogged_computations compendiums_reviewed
Sample Range:
(leave blank to include all observations)
From:
To:
Column Number of Endogenous Series
(?)
Fixed Seasonal Effects
0
Do not include Seasonal Dummies
Include Seasonal Dummies
Type of Equation
0
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation