Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
102 122 88 1 9 8976 102 918 88 792 9 99 114 106 1 3 10494 99 297 106 318 3 97 140 70 1 2 6790 97 194 70 140 2 82 143 70 1 21 5740 82 1722 70 1470 21 77 122 56 1 10 4312 77 770 56 560 10 65 127 50 1 2 3250 65 130 50 100 2 64 113 48 1 4 3072 64 256 48 192 4 62 118 71 1 4 4402 62 248 71 284 4 62 161 61 1 4 3782 62 248 61 244 4 62 134 66 1 5 4092 62 310 66 330 5 61 96 80 1 2 4880 61 122 80 160 2 59 104 37 1 3 2183 59 177 37 111 3 57 135 53 1 3 3021 57 171 53 159 3 56 110 39 1 6 2184 56 336 39 234 6 54 128 40 1 3 2160 54 162 40 120 3 54 142 59 1 4 3186 54 216 59 236 4 53 117 42 1 3 2226 53 159 42 126 3 52 94 33 1 13 1716 52 676 33 429 13 51 135 36 1 3 1836 51 153 36 108 3 51 121 57 1 4 2907 51 204 57 228 4 51 103 38 1 8 1938 51 408 38 304 8 50 118 98 1 8 4900 50 400 98 784 8 50 127 43 1 3 2150 50 150 43 129 3 50 116 73 1 4 3650 50 200 73 292 4 49 129 52 1 2 2548 49 98 52 104 2 49 115 53 1 5 2597 49 245 53 265 5 49 135 51 1 4 2499 49 196 51 204 4 48 133 32 1 3 1536 48 144 32 96 3 48 113 43 1 3 2064 48 144 43 129 3 47 111 53 1 2 2491 47 94 53 106 2 47 92 50 1 3 2350 47 141 50 150 3 46 118 50 1 3 2300 46 138 50 150 3 46 134 56 1 3 2576 46 138 56 168 3 45 106 53 1 5 2385 45 225 53 265 5 45 137 47 1 3 2115 45 135 47 141 3 45 100 42 1 4 1890 45 180 42 168 4 44 102 29 1 3 1276 44 132 29 87 3 43 134 54 1 4 2322 43 172 54 216 4 42 130 40 1 8 1680 42 336 40 320 8 42 144 41 1 3 1722 42 126 41 123 3 42 120 37 1 4 1554 42 168 37 148 4 42 91 25 1 2 1050 42 84 25 50 2 42 100 27 1 5 1134 42 210 27 135 5 42 134 61 1 4 2562 42 168 61 244 4 41 161 54 1 7 2214 41 287 54 378 7 41 128 35 1 3 1435 41 123 35 105 3 41 124 55 1 4 2255 41 164 55 220 4 41 115 47 1 6 1927 41 246 47 282 6 41 123 49 1 7 2009 41 287 49 343 7 41 117 38 1 20 1558 41 820 38 760 20 41 111 52 1 49 2132 41 2009 52 2548 49 40 146 35 1 3 1400 40 120 35 105 3 40 101 52 1 3 2080 40 120 52 156 3 40 131 54 1 6 2160 40 240 54 324 6 40 122 40 1 6 1600 40 240 40 240 6 40 78 52 1 4 2080 40 160 52 208 4 39 120 34 1 5 1326 39 195 34 170 5 39 115 51 1 4 1989 39 156 51 204 4 38 142 43 1 4 1634 38 152 43 172 4 38 94 40 1 31 1520 38 1178 40 1240 31 36 114 38 1 3 1368 36 108 38 114 3 36 108 33 1 3 1188 36 108 33 99 3 35 119 27 1 4 945 35 140 27 108 4 35 117 34 1 6 1190 35 210 34 204 6 35 86 44 1 5 1540 35 175 44 220 5 35 138 46 1 3 1610 35 105 46 138 3 34 119 50 1 3 1700 34 102 50 150 3 34 117 31 1 2 1054 34 68 31 62 2 34 117 33 1 3 1122 34 102 33 99 3 33 76 37 1 3 1221 33 99 37 111 3 33 119 48 1 16 1584 33 528 48 768 16 33 119 33 1 3 1089 33 99 33 99 3 32 124 40 1 3 1280 32 96 40 120 3 32 116 21 1 3 672 32 96 21 63 3 32 118 33 1 2 1056 32 64 33 66 2 31 102 41 1 5 1271 31 155 41 205 5 31 116 35 1 3 1085 31 93 35 105 3 30 103 60 1 4 1800 30 120 60 240 4 30 117 30 1 5 900 30 150 30 150 5 30 108 45 1 2 1350 30 60 45 90 2 30 122 26 1 3 780 30 90 26 78 3 29 90 41 1 3 1189 29 87 41 123 3 28 133 48 1 14 1344 28 392 48 672 14 28 116 10 1 8 280 28 224 10 80 8 27 110 35 1 4 945 27 108 35 140 4 27 90 23 1 4 621 27 108 23 92 4 27 74 29 1 3 783 27 81 29 87 3 26 75 17 1 4 442 26 104 17 68 4 25 107 35 1 3 875 25 75 35 105 3 25 90 50 1 5 1250 25 125 50 250 5 25 96 33 1 3 825 25 75 33 99 3 24 115 23 1 3 552 24 72 23 69 3 24 91 22 1 2 528 24 48 22 44 2 23 77 52 1 4 1196 23 92 52 208 4 23 108 38 1 31 874 23 713 38 1178 31 23 83 32 1 2 736 23 46 32 64 2 23 77 28 1 5 644 23 115 28 140 5 23 99 43 1 5 989 23 115 43 215 5 22 115 32 1 2 704 22 44 32 64 2 22 99 35 1 2 770 22 44 35 70 2 22 106 25 1 3 550 22 66 25 75 3 22 77 14 1 8 308 22 176 14 112 8 20 115 17 1 6 340 20 120 17 102 6 19 67 18 1 3 342 19 57 18 54 3 19 8 12 1 2 228 19 38 12 24 2 17 69 27 1 5 459 17 85 27 135 5 17 88 28 1 3 476 17 51 28 84 3 16 107 12 1 5 192 16 80 12 60 5 16 120 21 1 5 336 16 80 21 105 5 5 3 9 1 4 45 5 20 9 36 4 4 1 11 1 2 44 4 8 11 22 2 3 0 3 1 4 9 3 12 3 12 4 156 111 111 0 4 17316 0 624 0 444 0 109 69 137 0 8 14933 0 872 0 1096 0 104 116 112 0 3 11648 0 312 0 336 0 98 103 73 0 4 7154 0 392 0 292 0 78 139 99 0 3 7722 0 234 0 297 0 77 135 115 0 3 8855 0 231 0 345 0 73 113 95 0 3 6935 0 219 0 285 0 71 99 60 0 4 4260 0 284 0 240 0 67 76 94 0 4 6298 0 268 0 376 0 64 110 70 0 5 4480 0 320 0 350 0 62 121 87 0 2 5394 0 124 0 174 0 61 95 102 0 3 6222 0 183 0 306 0 58 66 69 0 4 4002 0 232 0 276 0 58 111 111 0 7 6438 0 406 0 777 0 56 77 55 0 3 3080 0 168 0 165 0 56 101 118 0 4 6608 0 224 0 472 0 52 108 90 0 3 4680 0 156 0 270 0 51 135 81 0 4 4131 0 204 0 324 0 51 70 88 0 4 4488 0 204 0 352 0 50 124 63 0 3 3150 0 150 0 189 0 49 92 84 0 6 4116 0 294 0 504 0 49 104 87 0 4 4263 0 196 0 348 0 48 113 78 0 4 3744 0 192 0 312 0 47 95 93 0 4 4371 0 188 0 372 0 47 89 69 0 4 3243 0 188 0 276 0 46 83 67 0 3 3082 0 138 0 201 0 45 96 61 0 6 2745 0 270 0 366 0 45 95 123 0 4 5535 0 180 0 492 0 45 110 91 0 6 4095 0 270 0 546 0 45 106 98 0 6 4410 0 270 0 588 0 44 78 38 0 2 1672 0 88 0 76 0 44 115 72 0 3 3168 0 132 0 216 0 44 74 59 0 3 2596 0 132 0 177 0 43 93 78 0 2 3354 0 86 0 156 0 43 88 58 0 4 2494 0 172 0 232 0 42 104 97 0 5 4074 0 210 0 485 0 41 86 69 0 3 2829 0 123 0 207 0 41 104 50 0 7 2050 0 287 0 350 0 40 99 66 0 4 2640 0 160 0 264 0 39 101 70 0 3 2730 0 117 0 210 0 39 53 65 0 4 2535 0 156 0 260 0 39 96 69 0 4 2691 0 156 0 276 0 39 58 49 0 3 1911 0 117 0 147 0 39 117 72 0 3 2808 0 117 0 216 0 39 82 74 0 4 2886 0 156 0 296 0 39 57 82 0 5 3198 0 195 0 410 0 38 71 61 0 3 2318 0 114 0 183 0 38 105 72 0 4 2736 0 152 0 288 0 38 60 77 0 6 2926 0 228 0 462 0 38 77 64 0 5 2432 0 190 0 320 0 37 73 23 0 7 851 0 259 0 161 0 37 78 39 0 5 1443 0 185 0 195 0 37 81 87 0 5 3219 0 185 0 435 0 36 101 46 0 3 1656 0 108 0 138 0 36 118 66 0 5 2376 0 180 0 330 0 36 59 57 0 4 2052 0 144 0 228 0 36 101 48 0 9 1728 0 324 0 432 0 36 22 75 0 3 2700 0 108 0 225 0 36 77 35 0 3 1260 0 108 0 105 0 35 100 53 0 3 1855 0 105 0 159 0 35 39 60 0 3 2100 0 105 0 180 0 34 42 20 0 15 680 0 510 0 300 0 34 80 66 0 4 2244 0 136 0 264 0 34 48 34 0 6 1156 0 204 0 204 0 34 131 80 0 3 2720 0 102 0 240 0 34 46 63 0 3 2142 0 102 0 189 0 33 89 46 0 3 1518 0 99 0 138 0 33 51 20 0 5 660 0 165 0 100 0 33 108 73 0 3 2409 0 99 0 219 0 33 86 57 0 4 1881 0 132 0 228 0 33 105 65 0 7 2145 0 231 0 455 0 33 85 70 0 4 2310 0 132 0 280 0 32 103 53 0 5 1696 0 160 0 265 0 32 83 60 0 7 1920 0 224 0 420 0 32 77 34 0 14 1088 0 448 0 476 0 32 26 18 0 25 576 0 800 0 450 0 32 73 49 0 6 1568 0 192 0 294 0 31 42 27 0 4 837 0 124 0 108 0 31 71 45 0 3 1395 0 93 0 135 0 31 105 9 0 62 279 0 1922 0 558 0 30 73 23 0 5 690 0 150 0 115 0 30 98 61 0 10 1830 0 300 0 610 0 29 108 67 0 4 1943 0 116 0 268 0 29 57 72 0 5 2088 0 145 0 360 0 29 37 58 0 5 1682 0 145 0 290 0 28 70 55 0 4 1540 0 112 0 220 0 28 73 33 0 10 924 0 280 0 330 0 28 47 40 0 5 1120 0 140 0 200 0 28 73 57 0 3 1596 0 84 0 171 0 28 91 61 0 3 1708 0 84 0 183 0 28 110 87 0 17 2436 0 476 0 1479 0 27 78 65 0 4 1755 0 108 0 260 0 27 92 85 0 6 2295 0 162 0 510 0 27 52 85 0 3 2295 0 81 0 255 0 26 88 54 0 4 1404 0 104 0 216 0 26 100 24 0 8 624 0 208 0 192 0 26 33 31 0 3 806 0 78 0 93 0 26 42 64 0 4 1664 0 104 0 256 0 25 81 70 0 4 1750 0 100 0 280 0 25 67 2 0 47 50 0 1175 0 94 0 24 8 27 0 8 648 0 192 0 216 0 24 46 29 0 3 696 0 72 0 87 0 24 83 68 0 5 1632 0 120 0 340 0 24 87 42 0 3 1008 0 72 0 126 0 24 82 78 0 4 1872 0 96 0 312 0 24 63 13 0 30 312 0 720 0 390 0 24 27 52 0 4 1248 0 96 0 208 0 23 14 25 0 12 575 0 276 0 300 0 23 83 38 0 8 874 0 184 0 304 0 23 168 40 0 3 920 0 69 0 120 0 23 67 42 0 8 966 0 184 0 336 0 23 21 40 0 4 920 0 92 0 160 0 23 55 74 0 3 1702 0 69 0 222 0 23 54 73 0 4 1679 0 92 0 292 0 22 118 56 0 4 1232 0 88 0 224 0 22 69 3 0 21 66 0 462 0 63 0 21 77 9 0 12 189 0 252 0 108 0 21 72 68 0 3 1428 0 63 0 204 0 21 53 28 0 3 588 0 63 0 84 0 21 40 36 0 4 756 0 84 0 144 0 20 102 38 0 6 760 0 120 0 228 0 20 25 55 0 17 1100 0 340 0 935 0 20 31 36 0 3 720 0 60 0 108 0 20 77 17 0 18 340 0 360 0 306 0 20 38 54 0 5 1080 0 100 0 270 0 19 23 57 0 5 1083 0 95 0 285 0 19 91 30 0 16 570 0 304 0 480 0 19 58 40 0 10 760 0 190 0 400 0 18 42 37 0 5 666 0 90 0 185 0 18 44 46 0 4 828 0 72 0 184 0 18 58 32 0 3 576 0 54 0 96 0 18 35 34 0 5 612 0 90 0 170 0 18 88 22 0 4 396 0 72 0 88 0 17 25 59 0 5 1003 0 85 0 295 0 17 39 32 0 10 544 0 170 0 320 0 16 48 18 0 12 288 0 192 0 216 0 16 64 28 0 4 448 0 64 0 112 0 15 65 34 0 9 510 0 135 0 306 0 15 95 29 0 12 435 0 180 0 348 0 15 29 24 0 10 360 0 150 0 240 0 15 2 24 0 9 360 0 135 0 216 0 14 83 23 0 17 322 0 238 0 391 0 13 11 43 0 6 559 0 78 0 258 0 13 16 28 0 3 364 0 39 0 84 0 12 9 19 0 4 228 0 48 0 76 0 11 46 16 0 19 176 0 209 0 304 0 11 41 40 0 3 440 0 33 0 120 0 10 14 14 0 9 140 0 90 0 126 0 10 63 19 0 7 190 0 70 0 133 0 10 9 22 0 4 220 0 40 0 88 0 10 0 8 0 3 80 0 30 0 24 0 9 58 31 0 46 279 0 414 0 1426 0 8 18 9 0 31 72 0 248 0 279 0 8 42 18 0 21 144 0 168 0 378 0 7 26 9 0 7 63 0 49 0 63 0 7 38 5 0 29 35 0 203 0 145 0 4 1 11 0 5 44 0 20 0 55 0
Names of X columns:
hours lfm blogs uk spr hours_blogs hours_uk hours_spr blogs_uk blogs_spr uk_spr
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