Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1587 252101 62 438 85 3 92 34 131 104 165119 951 134577 59 330 58 4 58 30 117 111 107269 1669 198520 62 609 51 14 62 38 146 93 93497 2283 189326 94 1015 131 2 108 34 132 119 100269 992 137449 43 294 44 1 55 25 80 57 91627 577 65295 27 164 42 3 8 31 117 80 47552 3916 439387 103 1912 94 0 134 29 112 107 233933 381 33186 19 111 46 0 1 18 67 22 6853 1790 178368 51 698 71 5 64 30 116 103 104380 1606 186657 38 556 65 0 77 29 107 72 98431 1886 261949 96 711 74 0 86 38 140 123 156949 1645 191051 95 495 55 7 93 49 186 164 81817 1433 138866 57 544 55 7 44 33 109 100 59238 2370 296878 66 959 98 3 106 46 159 143 101138 1645 192648 72 540 41 9 63 38 146 79 107158 4206 333462 162 1486 115 0 160 52 201 183 155499 1780 243571 58 635 45 4 104 32 124 123 156274 2352 263451 130 940 80 3 86 35 131 81 121777 1262 155679 48 452 37 3 93 25 96 74 105037 2026 227053 70 617 50 7 119 42 163 158 118661 2109 240028 63 695 93 0 107 40 151 133 131187 2451 388549 90 1046 93 1 86 35 128 128 145026 1253 156540 34 405 56 5 50 25 89 84 107016 1431 148421 43 477 57 9 92 46 184 184 87242 2631 177732 97 1012 138 0 123 36 136 127 91699 2293 191441 105 842 67 0 81 35 134 128 110087 2653 249893 122 994 88 5 93 38 146 118 145447 1709 236812 76 530 54 0 113 35 130 125 143307 1360 142329 45 515 47 0 52 28 105 89 61678 2051 259667 53 766 121 0 113 37 142 122 210080 1921 231625 65 734 44 3 112 40 155 151 165005 1623 176062 67 551 73 4 44 42 154 122 97806 1933 286683 79 718 49 1 123 44 169 162 184471 849 87485 33 280 36 4 38 33 125 121 27786 2641 322865 83 1055 61 2 111 35 135 132 184458 2235 247082 51 950 77 0 77 37 139 110 98765 2951 346011 106 1038 69 0 92 39 145 135 178441 1666 191653 74 552 63 2 74 32 124 80 100619 917 114673 31 275 36 1 33 17 55 46 58391 2745 284224 161 986 39 2 105 34 131 127 151672 2897 284195 72 1336 34 10 108 33 125 103 124437 1413 155363 59 565 65 6 66 35 128 95 79929 1535 177306 67 571 78 5 69 32 107 100 123064 1452 144571 49 404 67 5 62 35 130 102 50466 2369 140319 73 985 82 1 50 45 73 45 100991 4908 405267 135 1851 780 2 91 38 138 122 79367 918 78800 42 330 57 2 20 26 82 66 56968 2085 201970 69 611 72 0 101 45 173 159 106257 3655 302674 99 1249 112 9 129 44 169 153 178412 1923 164733 50 812 61 3 93 40 145 131 98520 1616 194221 68 501 39 0 89 33 134 113 153670 496 24188 24 218 20 0 8 4 12 7 15049 2328 342263 279 785 73 8 79 41 151 147 174478 744 65029 17 255 21 5 21 18 67 61 25109 1161 101097 64 454 70 3 30 14 52 41 45824 2598 246088 46 944 118 1 86 33 121 108 116772 2276 273108 75 600 72 5 116 49 186 184 189150 3185 282220 160 977 196 5 106 32 120 115 194404 2135 273495 119 863 58 0 127 37 135 132 185881 1863 214872 74 690 67 12 75 32 123 113 67508 3524 335121 124 1176 60 9 138 41 158 141 188597 2760 267171 106 1013 138 11 114 25 90 65 203618 2187 187938 88 890 69 9 55 40 157 87 87232 2161 229512 78 777 45 8 67 35 135 121 110875 1302 209798 61 521 54 2 45 33 125 112 144756 1540 201345 60 409 55 0 88 28 110 81 129825 1618 163833 114 493 46 6 67 31 121 116 92189 2403 204250 129 757 84 8 75 40 151 132 121158 1961 197813 67 736 71 2 114 32 123 104 96219 1402 132955 60 511 56 5 123 25 92 80 84128 2292 216092 59 789 55 13 86 42 162 145 97960 893 73566 32 385 39 6 22 23 88 67 23824 1935 213198 67 644 52 7 67 42 163 159 103515 1586 181713 49 664 94 2 77 38 133 90 91313 1494 148698 49 505 57 0 105 34 132 120 85407 1995 300103 70 878 82 4 119 38 144 126 95871 1904 251437 78 769 42 3 88 32 124 118 143846 1711 197295 101 499 45 6 78 37 140 112 155387 1746 158163 55 546 52 2 112 34 132 123 74429 1473 155529 57 551 63 0 66 33 122 98 74004 1538 132672 41 565 38 1 58 25 97 78 71987 3324 377205 100 1086 108 0 132 40 155 119 150629 1356 145905 66 649 45 5 30 26 99 99 68580 1931 223701 86 540 53 2 100 40 106 81 119855 870 80953 25 437 31 0 49 8 28 27 55792 1716 130805 47 732 169 0 26 27 101 77 25157 1091 135082 48 308 60 5 67 32 120 118 90895 3187 300805 156 1237 271 1 57 33 127 122 117510 2240 271806 95 783 84 0 95 50 178 103 144774 2677 150949 96 933 63 1 139 37 141 129 77529 2038 225805 79 710 54 1 73 33 122 69 103123 1783 197389 68 563 65 2 134 34 127 121 104669 1560 156583 56 508 80 6 37 28 102 81 82414 2540 222599 66 936 84 1 98 32 124 119 82390 2118 261601 70 838 115 4 58 32 124 116 128446 1521 178489 35 523 60 3 78 32 124 123 111542 1460 200657 43 500 62 3 88 31 111 111 136048 1866 259084 68 691 57 0 142 35 129 100 197257 3312 313075 130 1060 121 11 127 58 223 221 162079 3135 346933 100 1232 69 12 139 27 102 95 206286 2232 246440 104 735 60 8 108 45 174 153 109858 2061 252444 58 757 81 0 128 37 141 118 182125 1787 159965 159 574 100 0 62 32 122 50 74168 602 43287 14 214 43 4 13 19 71 64 19630 1973 172239 68 661 72 4 89 22 81 34 88634 1739 183738 120 631 61 0 83 35 131 76 128321 2282 227681 43 1015 101 0 116 36 139 112 118936 2339 260464 81 893 50 0 157 36 137 115 127044 923 106288 54 293 32 0 28 23 91 69 178377 1389 109632 77 446 74 0 83 36 142 108 69581 1673 268905 58 538 54 4 72 36 133 130 168019 2082 266805 78 627 65 0 134 42 155 110 113598 398 23623 11 156 9 0 12 1 0 0 5841 1605 152474 65 577 45 0 106 32 123 83 93116 530 61857 25 192 25 4 23 11 32 30 24610 1503 144889 43 437 102 0 83 40 149 106 60611 2636 346600 99 1054 59 1 126 34 128 91 226620 387 21054 16 146 2 0 4 0 0 0 6622 1849 224051 45 751 56 5 71 27 99 69 121996 449 31414 19 200 22 0 18 8 25 9 13155 2912 261043 105 1050 146 2 98 35 132 123 154158 1740 197819 57 590 63 7 66 41 155 143 78489 1401 154984 73 430 91 12 44 40 151 125 22007 1257 112933 45 467 46 2 29 28 103 81 72530 568 38214 34 276 52 0 16 8 27 21 13983 1512 158671 33 528 98 2 56 35 131 124 73397 2379 302148 70 898 105 0 112 47 178 168 143878 1180 177918 55 411 57 0 46 46 177 149 119956 3071 350552 70 1362 126 3 129 42 163 147 181558 2186 275578 91 743 120 0 139 48 187 145 208236 3065 368746 106 1069 104 3 136 49 182 172 237085 1127 172464 31 431 44 0 66 35 135 126 110297 1045 94381 35 380 48 0 42 32 118 89 61394 2419 243875 279 788 143 4 70 36 140 137 81420 3839 382487 153 1367 146 4 97 42 158 149 191154 1427 114525 40 449 91 14 49 35 132 121 11798 3406 335681 119 1461 129 0 113 37 136 133 135724 1730 147989 72 651 67 4 55 34 123 93 68614 1507 216638 45 494 74 0 100 36 134 119 139926 1907 192862 72 667 168 1 80 36 129 102 105203 1552 184818 107 510 69 0 29 32 125 45 80338 3611 336707 105 1472 99 9 95 33 128 104 121376 1925 215836 76 675 61 1 114 35 129 111 124922 2035 173260 63 716 37 3 41 21 79 78 10901 2503 271773 89 814 51 10 128 40 154 120 135471 1603 130908 52 556 121 5 142 49 188 176 66395 2272 204009 75 887 48 2 88 33 122 109 134041 2086 245514 92 663 52 1 147 39 144 132 153554 2 1 0 0 0 9 0 0 0 0 0 207 14688 10 85 0 0 4 0 0 0 7953 5 98 1 0 0 0 0 0 0 0 0 8 455 2 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1777 195765 75 607 51 2 56 33 120 78 98922 2781 326038 121 934 98 1 121 42 168 104 165395 0 0 0 0 0 0 0 0 0 0 0 4 203 4 0 0 0 0 0 0 0 0 151 7199 5 74 0 0 7 0 0 0 4245 474 46660 20 259 7 0 12 5 15 13 21509 141 17547 5 69 3 0 0 1 4 4 7670 976 107465 38 267 80 0 37 38 133 65 15167 29 969 2 0 0 0 0 0 0 0 0 1542 173102 58 517 43 2 47 28 101 55 63891
Names of X columns:
Pageviews TIMErfc #logins CCviews Cviews Shares BloggedComputations RVC Feedbackm Feedback+120 TimeCompendium
Sample Range:
(leave blank to include all observations)
From:
To:
Column Number of Endogenous Series
(?)
Fixed Seasonal Effects
-0.1
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
1
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