Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
162687 0 48 21 82 20465 6200 23975 39 37 201906 1 58 20 80 33629 10265 85634 46 43 7215 0 0 0 0 1423 603 1929 0 0 146367 0 67 27 84 25629 8874 36294 54 54 257045 0 83 31 124 54002 20323 72255 93 86 524450 1 136 36 140 151036 26258 189748 198 181 188294 1 65 23 88 33287 10165 61834 42 42 195674 0 86 30 115 31172 8247 68167 59 59 177020 0 62 30 109 28113 8683 38462 49 46 325899 1 71 27 108 57803 16957 101219 83 77 121844 2 50 24 63 49830 8058 43270 49 49 203938 0 88 30 118 52143 20488 76183 83 79 113213 0 61 22 71 21055 7945 31476 39 37 220751 4 79 28 112 47007 13448 62157 93 92 172905 4 56 18 63 28735 5389 46261 31 31 156326 3 54 22 86 59147 6185 50063 29 28 145178 0 81 37 148 78950 24369 64483 104 103 89171 5 13 15 54 13497 70 2341 2 2 172624 0 74 34 134 46154 17327 48149 46 48 39790 0 18 18 57 53249 3878 12743 27 25 87927 0 31 15 59 10726 3149 18743 16 16 241285 0 99 30 113 83700 20517 97057 108 106 195820 1 38 25 96 40400 2570 17675 36 35 146946 1 59 34 96 33797 5162 33106 33 33 159763 1 54 21 78 36205 5299 53311 46 45 207078 0 63 21 80 30165 7233 42754 65 64 212394 0 66 25 93 58534 15657 59056 80 73 201536 0 90 31 109 44663 15329 101621 81 78 394662 0 72 31 115 92556 14881 118120 69 63 217892 0 61 20 79 40078 16318 79572 69 69 182286 0 61 28 103 34711 9556 42744 37 36 181740 2 61 22 71 31076 10462 65931 45 41 137978 4 53 17 66 74608 7192 38575 62 59 255929 0 118 25 100 58092 4362 28795 33 33 236489 1 73 25 100 42009 14349 94440 77 76 0 0 0 0 0 0 0 0 0 0 230761 0 54 31 121 36022 10881 38229 34 27 132807 3 54 14 51 23333 8022 31972 44 44 157118 9 46 35 119 53349 13073 40071 43 43 253254 0 83 34 136 92596 26641 132480 117 104 269329 2 106 22 84 49598 14426 62797 125 120 161273 0 44 34 136 44093 15604 40429 49 44 107181 2 27 23 84 84205 9184 45545 76 71 195891 1 64 24 92 63369 5989 57568 81 78 139667 2 71 26 103 60132 11270 39019 111 106 171101 2 44 23 85 37403 13958 53866 61 61 81407 1 23 35 106 24460 7162 38345 56 53 247563 0 78 24 96 46456 13275 50210 54 51 239807 1 60 31 124 66616 21224 80947 47 46 172743 8 73 30 106 41554 10615 43461 55 55 48188 0 12 22 82 22346 2102 14812 14 14 169355 0 104 23 87 30874 12396 37819 44 44 315622 0 83 27 97 68701 18717 102738 115 113 241518 0 57 30 107 35728 9724 54509 57 55 195583 1 67 33 126 29010 9863 62956 48 46 159913 8 44 12 43 23110 8374 55411 40 39 220241 0 53 26 96 38844 8030 50611 51 51 101694 1 26 26 100 27084 7509 26692 32 31 157258 0 67 23 91 35139 14146 60056 36 36 202536 10 36 38 136 57476 7768 25155 47 47 173505 6 56 32 128 33277 13823 42840 51 53 150518 0 52 21 83 31141 7230 39358 37 38 141491 11 54 22 74 61281 10170 47241 52 52 125612 3 57 26 96 25820 7573 49611 42 37 166049 0 27 28 102 23284 5753 41833 11 11 124197 0 58 33 122 35378 9791 48930 47 45 195043 8 76 36 144 74990 19365 110600 59 59 138708 2 93 25 90 29653 9422 52235 82 82 116552 0 59 25 97 64622 12310 53986 49 49 31970 0 5 21 78 4157 1283 4105 6 6 258158 3 57 19 72 29245 6372 59331 83 81 151184 1 42 12 45 50008 5413 47796 56 56 135926 2 88 30 120 52338 10837 38302 114 105 119629 1 53 21 59 13310 3394 14063 46 46 171518 0 81 39 150 92901 12964 54414 46 46 108949 2 35 32 117 10956 3495 9903 2 2 183471 1 102 28 123 34241 11580 53987 51 51 159966 0 71 29 114 75043 9970 88937 96 95 93786 0 28 21 75 21152 4911 21928 20 18 84971 0 34 31 114 42249 10138 29487 57 55 88882 0 54 26 94 42005 14697 35334 49 48 304603 0 49 29 116 41152 8464 57596 51 48 75101 1 30 23 86 14399 4204 29750 40 39 145043 0 57 25 90 28263 10226 41029 40 40 95827 0 54 22 87 17215 3456 12416 36 36 173924 0 38 26 99 48140 8895 51158 64 60 241957 0 63 33 132 62897 22557 79935 117 114 115367 0 58 24 96 22883 6900 26552 40 39 118408 7 46 24 91 41622 8620 25807 46 45 164078 0 46 21 77 40715 7820 50620 61 59 158931 5 51 28 104 65897 12112 61467 59 59 184139 1 87 28 100 76542 13178 65292 94 93 152856 0 39 25 94 37477 7028 55516 36 35 144014 0 28 15 60 53216 6616 42006 51 47 62535 0 26 13 46 40911 9570 26273 39 36 245196 0 52 36 135 57021 14612 90248 62 59 199841 0 96 27 99 73116 11219 61476 79 79 19349 0 13 1 2 3895 786 9604 14 14 247280 3 43 24 96 46609 11252 45108 45 42 159408 0 42 31 109 29351 9289 47232 43 41 72128 0 30 4 15 2325 593 3439 8 8 104253 0 59 21 68 31747 6562 30553 41 41 151090 0 73 27 102 32665 8208 24751 25 24 137382 1 39 23 84 19249 7488 34458 22 22 87448 1 36 12 46 15292 4574 24649 18 18 27676 0 2 16 59 5842 522 2342 3 1 165507 0 102 29 116 33994 12840 52739 54 53 132148 1 30 26 29 13018 1350 6245 6 6 0 0 0 0 0 0 0 0 0 0 95778 0 46 25 91 98177 10623 35381 50 49 109001 0 25 21 76 37941 5322 19595 33 33 158833 0 59 24 86 31032 7987 50848 54 50 147690 1 60 21 84 32683 10566 39443 63 64 89887 0 36 21 65 34545 1900 27023 56 53 3616 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 199005 0 45 23 84 27525 10698 61022 49 48 160930 0 79 33 114 66856 14884 63528 90 90 177948 2 30 32 132 28549 6852 34835 51 46 136061 0 43 23 92 38610 6873 37172 29 29 43410 0 7 1 3 2781 4 13 1 1 184277 1 80 29 109 41211 9188 62548 68 64 108858 0 32 20 81 22698 5141 31334 29 29 141744 8 81 33 121 41194 4260 20839 27 27 60493 3 3 12 48 32689 443 5084 4 4 19764 1 10 2 8 5752 2416 9927 10 10 177559 3 47 21 80 26757 9831 53229 47 47 140281 0 35 28 107 22527 5953 29877 44 44 164249 0 54 35 140 44810 9435 37310 53 51 11796 0 1 2 8 0 0 0 0 0 10674 0 0 0 0 0 0 0 0 0 151322 0 46 18 56 100674 7642 50067 40 38 6836 0 0 1 4 0 0 0 0 0 174712 6 51 21 70 57786 6837 47708 57 57 5118 0 5 0 0 0 0 0 0 0 40248 1 8 4 14 5444 775 6012 6 6 0 0 0 0 0 0 0 0 0 0 127628 0 38 29 104 28470 8191 27749 24 22 88837 0 21 26 89 61849 1661 47555 34 34 7131 1 0 0 0 0 0 0 0 0 9056 0 0 4 12 2179 548 1336 10 10 87957 1 18 19 60 8019 3080 11017 16 16 144470 0 53 22 84 39644 13400 55184 93 93 111408 1 17 22 88 23494 8181 43485 28 22
Names of X columns:
timeRFC compshared blogged reviewedcomp submfeedback characters revisions seconds inclhyperlinks inclblog
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