Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
0 65 146455 0 1 0 95556 0 114468 0 127 0 0 54 84944 0 4 0 54565 0 88594 0 90 0 0 58 113337 0 9 0 63016 0 74151 0 68 0 0 75 128655 0 2 0 79774 0 77921 0 111 0 0 41 74398 0 1 0 31258 0 53212 0 51 0 0 0 35523 0 2 0 52491 0 34956 0 33 0 1 111 293403 293403 0 0 91256 91256 149703 149703 123 123 0 1 32750 0 0 0 22807 0 6853 0 5 0 0 36 106539 0 5 0 77411 0 58907 0 63 0 0 60 130539 0 0 0 48821 0 67067 0 66 0 0 63 154991 0 0 0 52295 0 110563 0 99 0 1 71 126683 126683 7 7 63262 63262 58126 58126 72 72 1 38 100672 100672 6 6 50466 50466 57113 57113 55 55 1 76 179562 179562 3 3 62932 62932 77993 77993 116 116 1 61 125971 125971 4 4 38439 38439 68091 68091 71 71 0 125 234509 0 0 0 70817 0 124676 0 125 0 1 84 158980 158980 4 4 105965 105965 109522 109522 123 123 0 69 184217 0 3 0 73795 0 75865 0 74 0 1 77 107342 107342 0 0 82043 82043 79746 79746 116 116 1 95 141371 141371 5 5 74349 74349 77844 77844 117 117 1 78 154730 154730 0 0 82204 82204 98681 98681 98 98 1 76 264020 264020 1 1 55709 55709 105531 105531 101 101 1 40 90938 90938 3 3 37137 37137 51428 51428 43 43 0 81 101324 0 5 0 70780 0 65703 0 103 0 0 102 130232 0 0 0 55027 0 72562 0 107 0 0 70 137793 0 0 0 56699 0 81728 0 77 0 0 75 161678 0 4 0 65911 0 95580 0 87 0 1 93 151503 151503 0 0 56316 56316 98278 98278 99 99 1 42 105324 105324 0 0 26982 26982 46629 46629 46 46 1 95 175914 175914 0 0 54628 54628 115189 115189 96 96 0 87 181853 0 3 0 96750 0 124865 0 92 0 0 44 114928 0 4 0 53009 0 59392 0 96 0 1 84 190410 190410 1 1 64664 64664 127818 127818 96 96 1 28 61499 61499 4 4 36990 36990 17821 17821 15 15 1 87 223004 223004 1 1 85224 85224 154076 154076 147 147 1 71 167131 167131 0 0 37048 37048 64881 64881 56 56 1 68 233482 233482 0 0 59635 59635 136506 136506 81 81 1 50 121185 121185 2 2 42051 42051 66524 66524 69 69 1 30 78776 78776 1 1 26998 26998 45988 45988 34 34 1 86 188967 188967 2 2 63717 63717 107445 107445 98 98 1 75 199512 199512 8 8 55071 55071 102772 102772 82 82 1 46 102531 102531 5 5 40001 40001 46657 46657 64 64 1 52 118958 118958 3 3 54506 54506 97563 97563 61 61 0 31 68948 0 4 0 35838 0 36663 0 45 0 0 30 93125 0 1 0 50838 0 55369 0 37 0 0 70 277108 0 2 0 86997 0 77921 0 64 0 0 20 78800 0 2 0 33032 0 56968 0 21 0 0 84 157250 0 0 0 61704 0 77519 0 104 0 1 81 210554 210554 6 6 117986 117986 129805 129805 126 126 0 79 127324 0 3 0 56733 0 72761 0 104 0 1 70 114397 114397 0 0 55064 55064 81278 81278 87 87 1 8 24188 24188 0 0 5950 5950 15049 15049 7 7 1 67 246209 246209 6 6 84607 84607 113935 113935 130 130 0 21 65029 0 5 0 32551 0 25109 0 21 0 0 30 98030 0 3 0 31701 0 45824 0 35 0 1 70 173587 173587 1 1 71170 71170 89644 89644 97 97 0 87 172684 0 5 0 101773 0 109011 0 103 0 0 87 191381 0 5 0 101653 0 134245 0 210 0 1 112 191276 191276 0 0 81493 81493 136692 136692 151 151 1 54 134043 134043 9 9 55901 55901 50741 50741 57 57 1 96 233406 233406 6 6 109104 109104 149510 149510 117 117 1 93 195304 195304 6 6 114425 114425 147888 147888 152 152 1 49 127619 127619 5 5 36311 36311 54987 54987 52 52 0 49 162810 0 6 0 70027 0 74467 0 83 0 0 38 129100 0 2 0 73713 0 100033 0 87 0 1 64 108715 108715 0 0 40671 40671 85505 85505 80 80 0 62 106469 0 3 0 89041 0 62426 0 88 0 1 66 142069 142069 8 8 57231 57231 82932 82932 83 83 0 98 143937 0 2 0 78792 0 79169 0 140 0 1 97 84256 84256 5 5 59155 59155 65469 65469 76 76 0 56 118807 0 11 0 55827 0 63572 0 70 0 1 22 69471 69471 6 6 22618 22618 23824 23824 26 26 0 51 122433 0 5 0 58425 0 73831 0 66 0 1 56 131122 131122 1 1 65724 65724 63551 63551 89 89 0 94 94763 0 0 0 56979 0 56756 0 100 0 1 98 188780 188780 3 3 72369 72369 81399 81399 98 98 0 76 191467 0 3 0 79194 0 117881 0 109 0 0 57 105615 0 6 0 202316 0 70711 0 51 0 0 75 89318 0 1 0 44970 0 50495 0 82 0 0 48 107335 0 0 0 49319 0 53845 0 65 0 0 48 98599 0 1 0 36252 0 51390 0 46 0 0 109 260646 0 0 0 75741 0 104953 0 104 0 1 27 131876 131876 5 5 38417 38417 65983 65983 36 36 1 83 119291 119291 2 2 64102 64102 76839 76839 123 123 1 49 80953 80953 0 0 56622 56622 55792 55792 59 59 1 24 99768 99768 0 0 15430 15430 25155 25155 27 27 1 43 84572 84572 5 5 72571 72571 55291 55291 84 84 1 44 202373 202373 1 1 67271 67271 84279 84279 61 61 1 49 166790 166790 0 0 43460 43460 99692 99692 46 46 0 106 99946 0 1 0 99501 0 59633 0 125 0 1 42 116900 116900 1 1 28340 28340 63249 63249 58 58 0 108 142146 0 2 0 76013 0 82928 0 152 0 1 27 99246 99246 4 4 37361 37361 50000 50000 52 52 0 79 156833 0 1 0 48204 0 69455 0 85 0 1 49 175078 175078 4 4 76168 76168 84068 84068 95 95 0 64 130533 0 0 0 85168 0 76195 0 78 0 1 75 142339 142339 2 2 125410 125410 114634 114634 144 144 0 115 176789 0 0 0 123328 0 139357 0 149 0 1 92 181379 181379 7 7 83038 83038 110044 110044 101 101 0 106 228548 0 7 0 120087 0 155118 0 205 0 1 73 142141 142141 6 6 91939 91939 83061 83061 61 61 1 105 167845 167845 0 0 103646 103646 127122 127122 145 145 1 30 103012 103012 0 0 29467 29467 45653 45653 28 28 1 13 43287 43287 4 4 43750 43750 19630 19630 49 49 1 69 125366 125366 4 4 34497 34497 67229 67229 68 68 1 72 118372 118372 0 0 66477 66477 86060 86060 142 142 0 80 135171 0 0 0 71181 0 88003 0 82 0 0 106 175568 0 0 0 74482 0 95815 0 105 0 0 28 74112 0 0 0 174949 0 85499 0 52 0 0 70 88817 0 0 0 46765 0 27220 0 56 0 1 51 164767 164767 4 4 90257 90257 109882 109882 81 81 0 90 141933 0 0 0 51370 0 72579 0 100 0 0 12 22938 0 0 0 1168 0 5841 0 11 0 0 84 115199 0 0 0 51360 0 68369 0 87 0 0 23 61857 0 4 0 25162 0 24610 0 31 0 1 57 91185 91185 0 0 21067 21067 30995 30995 67 67 0 84 213765 0 1 0 58233 0 150662 0 150 0 1 4 21054 21054 0 0 855 855 6622 6622 4 4 0 56 167105 0 5 0 85903 0 93694 0 75 0 0 18 31414 0 0 0 14116 0 13155 0 39 0 1 86 178863 178863 1 1 57637 57637 111908 111908 88 88 0 39 126681 0 7 0 94137 0 57550 0 67 0 1 16 64320 64320 5 5 62147 62147 16356 16356 24 24 1 18 67746 67746 2 2 62832 62832 40174 40174 58 58 1 16 38214 38214 0 0 8773 8773 13983 13983 16 16 1 42 90961 90961 1 1 63785 63785 52316 52316 49 49 1 75 181510 181510 0 0 65196 65196 99585 99585 109 109 0 30 116775 0 0 0 73087 0 86271 0 124 0 0 104 223914 0 2 0 72631 0 131012 0 115 0 0 121 185139 0 0 0 86281 0 130274 0 128 0 0 106 242879 0 2 0 162365 0 159051 0 159 0 1 57 139144 139144 0 0 56530 56530 76506 76506 75 75 1 28 75812 75812 0 0 35606 35606 49145 49145 30 30 1 56 178218 178218 4 4 70111 70111 66398 66398 83 83 1 81 246834 246834 4 4 92046 92046 127546 127546 135 135 0 2 50999 0 8 0 63989 0 6802 0 8 0 0 88 223842 0 0 0 104911 0 99509 0 115 0 0 41 93577 0 4 0 43448 0 43106 0 60 0 1 83 155383 155383 0 0 60029 60029 108303 108303 99 99 1 55 111664 111664 1 1 38650 38650 64167 64167 98 98 1 3 75426 75426 0 0 47261 47261 8579 8579 36 36 1 54 243551 243551 9 9 73586 73586 97811 97811 93 93 1 89 136548 136548 0 0 83042 83042 84365 84365 158 158 1 41 173260 173260 3 3 37238 37238 10901 10901 16 16 0 94 185039 0 7 0 63958 0 91346 0 100 0 0 101 67507 0 5 0 78956 0 33660 0 49 0 0 70 139350 0 2 0 99518 0 93634 0 89 0 0 111 172964 0 1 0 111436 0 109348 0 153 0 1 0 0 0 9 9 0 0 0 0 0 0 1 4 14688 14688 0 0 6023 6023 7953 7953 5 5 1 0 98 98 0 0 0 0 0 0 0 0 1 0 455 455 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 0 0 0 1 42 128066 128066 2 2 42564 42564 63538 63538 80 80 0 97 176460 0 1 0 38885 0 108281 0 122 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 203 203 0 0 0 0 0 0 0 0 1 7 7199 7199 0 0 1644 1644 4245 4245 6 6 1 12 46660 46660 0 0 6179 6179 21509 21509 13 13 0 0 17547 0 0 0 3926 0 7670 0 3 0 1 37 73567 73567 0 0 23238 23238 10641 10641 18 18 0 0 969 0 0 0 0 0 0 0 0 0 0 39 101060 0 2 0 49288 0 41243 0 49 0
Names of X columns:
Pop BloggedComputations TotalTime TotalTimep Shared Sharedp Characters Charactersp Writing Writingp Hyperlinks Hyperlinksp
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
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation