Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 162556 162556 213118 213118 6282929 1 29790 29790 81767 81767 4324047 1 87550 87550 153198 153198 4108272 0 84738 0 -26007 0 -1212617 1 54660 54660 126942 126942 1485329 1 42634 42634 157214 157214 1779876 0 40949 0 129352 0 1367203 1 42312 42312 234817 234817 2519076 1 37704 37704 60448 60448 912684 1 16275 16275 47818 47818 1443586 0 25830 0 245546 0 1220017 0 12679 0 48020 0 984885 1 18014 18014 -1710 -1710 1457425 0 43556 0 32648 0 -572920 1 24524 24524 95350 95350 929144 0 6532 0 151352 0 1151176 0 7123 0 288170 0 790090 1 20813 20813 114337 114337 774497 1 37597 37597 37884 37884 990576 0 17821 0 122844 0 454195 1 12988 12988 82340 82340 876607 1 22330 22330 79801 79801 711969 0 13326 0 165548 0 702380 0 16189 0 116384 0 264449 0 7146 0 134028 0 450033 0 15824 0 63838 0 541063 1 26088 26088 74996 74996 588864 0 11326 0 31080 0 -37216 0 8568 0 32168 0 783310 0 14416 0 49857 0 467359 1 3369 3369 87161 87161 688779 1 11819 11819 106113 106113 608419 1 6620 6620 80570 80570 696348 1 4519 4519 102129 102129 597793 0 2220 0 301670 0 821730 0 18562 0 102313 0 377934 0 10327 0 88577 0 651939 1 5336 5336 112477 112477 697458 1 2365 2365 191778 191778 700368 0 4069 0 79804 0 225986 0 7710 0 128294 0 348695 0 13718 0 96448 0 373683 0 4525 0 93811 0 501709 0 6869 0 117520 0 413743 0 4628 0 69159 0 379825 1 3653 3653 101792 101792 336260 1 1265 1265 210568 210568 636765 1 7489 7489 136996 136996 481231 0 4901 0 121920 0 469107 0 2284 0 76403 0 211928 1 3160 3160 108094 108094 563925 1 4150 4150 134759 134759 511939 1 7285 7285 188873 188873 521016 1 1134 1134 146216 146216 543856 1 4658 4658 156608 156608 329304 0 2384 0 61348 0 423262 0 3748 0 50350 0 509665 0 5371 0 87720 0 455881 0 1285 0 99489 0 367772 1 9327 9327 87419 87419 406339 1 5565 5565 94355 94355 493408 0 1528 0 60326 0 232942 1 3122 3122 94670 94670 416002 1 7317 7317 82425 82425 337430 0 2675 0 59017 0 361517 0 13253 0 90829 0 360962 0 880 0 80791 0 235561 1 2053 2053 100423 100423 408247 0 1424 0 131116 0 450296 1 4036 4036 100269 100269 418799 1 3045 3045 27330 27330 247405 0 5119 0 39039 0 378519 0 1431 0 106885 0 326638 0 554 0 79285 0 328233 0 1975 0 118881 0 386225 1 1286 1286 77623 77623 283662 0 1012 0 114768 0 370225 0 810 0 74015 0 269236 0 1280 0 69465 0 365732 1 666 666 117869 117869 420383 0 1380 0 60982 0 345811 1 4608 4608 90131 90131 431809 0 876 0 138971 0 418876 0 814 0 39625 0 297476 0 514 0 102725 0 416776 1 5692 5692 64239 64239 357257 0 3642 0 90262 0 458343 0 540 0 103960 0 388386 0 2099 0 106611 0 358934 0 567 0 103345 0 407560 0 2001 0 95551 0 392558 1 2949 2949 82903 82903 373177 0 2253 0 63593 0 428370 1 6533 6533 126910 126910 369419 0 1889 0 37527 0 358649 1 3055 3055 60247 60247 376641 0 272 0 112995 0 467427 1 1414 1414 70184 70184 364885 0 2564 0 130140 0 436230 1 1383 1383 73221 73221 329118 1 1261 1261 76114 76114 317365 0 975 0 90534 0 286849 0 3366 0 108479 0 376685 0 576 0 113761 0 407198 0 1306 0 68696 0 377772 0 746 0 71561 0 271483 1 3192 3192 59831 59831 153661 1 2045 2045 97890 97890 513294 0 5477 0 101481 0 324881 1 1932 1932 72954 72954 264512 0 936 0 67939 0 420968 1 3437 3437 48022 48022 129302 0 5131 0 86111 0 191521 1 2397 2397 74020 74020 268673 1 1389 1389 57530 57530 353179 0 1503 0 56364 0 354624 0 402 0 84990 0 363713 0 2239 0 88590 0 456657 1 2234 2234 77200 77200 211742 0 837 0 61262 0 338381 0 10579 0 110309 0 418530 0 875 0 67000 0 351483 0 1395 0 93099 0 372928 1 1659 1659 107577 107577 485538 1 2647 2647 62920 62920 279268 1 3294 3294 75832 75832 219060 0 0 0 60720 0 325560 0 94 0 60793 0 325314 0 422 0 57935 0 322046 0 0 0 60720 0 325560 0 34 0 60630 0 325599 0 1558 0 55637 0 377028 0 0 0 60720 0 325560 0 43 0 60887 0 323850 0 0 0 60720 0 325560 0 316 0 60505 0 331514 0 115 0 60945 0 325632 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 389 0 58990 0 322265 0 0 0 60720 0 325560 0 1002 0 56750 0 325906 0 36 0 60894 0 325985 0 460 0 63346 0 346145 0 309 0 56535 0 325898 0 0 0 60720 0 325560 0 9 0 60835 0 325356 0 0 0 60720 0 325560 0 14 0 61016 0 325930 0 520 0 58650 0 318020 0 1766 0 60438 0 326389 0 0 0 60720 0 325560 0 458 0 58625 0 302925 0 20 0 60938 0 325540 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 98 0 61490 0 326736 0 405 0 60845 0 340580 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 483 0 60830 0 331828 0 454 0 63261 0 323299 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 757 0 45689 0 387722 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 36 0 61564 0 324598 0 0 0 60720 0 325560 0 203 0 61938 0 328726 0 0 0 60720 0 325560 0 90 0 60951 0 325043 0 0 0 60720 0 325560 1 71 71 60745 60745 325806 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 972 0 71642 0 387732 1 531 531 71641 71641 349729 0 604 0 55792 0 332202 1 283 283 71873 71873 305442 1 23 23 62555 62555 329537 1 638 638 60370 60370 327055 1 699 699 64873 64873 356245 0 149 0 62041 0 328451 0 226 0 65745 0 307062 0 0 0 60720 0 325560 0 275 0 59500 0 331345 0 0 0 60720 0 325560 0 141 0 61630 0 331824 0 0 0 60720 0 325560 0 28 0 60890 0 325685 1 0 0 60720 60720 325560 1 2566 2566 113521 113521 404480 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 1 472 472 80045 80045 318314 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 1 203 203 50804 50804 311807 1 496 496 87390 87390 337724 1 10 10 61656 61656 326431 1 63 63 65688 65688 327556 1 0 0 60720 60720 325560 1 1136 1136 48522 48522 356850 1 0 0 60720 60720 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 267 0 57640 0 322741 0 474 0 61977 0 310902 0 534 0 62620 0 324295 0 0 0 60720 0 325560 0 15 0 60831 0 326156 0 397 0 60646 0 326960 0 0 0 60720 0 325560 0 1061 0 56225 0 333411 0 288 0 60510 0 297761 0 0 0 60720 0 325560 0 3 0 60698 0 325536 0 0 0 60720 0 325560 0 20 0 60805 0 325762 0 278 0 61404 0 327957 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 192 0 65276 0 318521 0 0 0 60720 0 325560 0 317 0 63915 0 319775 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 368 0 61686 0 332128 0 0 0 60720 0 325560 0 2 0 60743 0 325486 0 0 0 60720 0 325560 0 53 0 60349 0 325838 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 94 0 61360 0 331767 0 0 0 60720 0 325560 0 24 0 59818 0 324523 0 2332 0 72680 0 339995 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 0 131 0 61808 0 319582 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 0 206 0 53110 0 307245 1 0 0 60720 60720 325560 0 167 0 64245 0 317967 0 622 0 73007 0 331488 0 885 0 82732 0 335452 0 0 0 60720 0 325560 0 365 0 54820 0 334184 0 364 0 47705 0 313213 1 0 0 60720 60720 325560 0 0 0 60720 0 325560 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 0 226 0 72835 0 348678 0 307 0 58856 0 328727 1 0 0 60720 60720 325560 0 0 0 60720 0 325560 1 0 0 60720 60720 325560 0 188 0 77655 0 387978 1 0 0 60720 60720 325560 0 138 0 69817 0 336704 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 0 125 0 60798 0 322076 0 0 0 60720 0 325560 0 282 0 62452 0 334272 0 335 0 64175 0 338197 0 0 0 60720 0 325560 1 813 813 67440 67440 321024 0 176 0 68136 0 322145 1 0 0 60720 60720 325560 0 0 0 60720 0 325560 0 249 0 56726 0 323351 0 0 0 60720 0 325560 0 333 0 70811 0 327748 0 0 0 60720 0 325560 1 0 0 60720 60720 325560 0 30 0 62045 0 328157 1 0 0 60720 60720 325560 0 249 0 54323 0 311594 1 0 0 60720 60720 325560 0 165 0 62841 0 335962 0 453 0 81125 0 372426 0 0 0 60720 0 325560 0 53 0 59506 0 319844 1 382 382 59365 59365 355822 0 0 0 60720 0 325560 1 0 0 60720 60720 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 1 30 30 60798 60798 324047 0 290 0 58790 0 311464 0 0 0 60720 0 325560 1 0 0 60720 60720 325560 0 366 0 61808 0 353417 0 2 0 60735 0 325590 1 0 0 60720 60720 325560 1 209 209 64016 64016 328576 0 384 0 54683 0 326126 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 0 365 0 87192 0 369376 1 0 0 60720 60720 325560 1 49 49 64107 64107 332013 0 3 0 60761 0 325871 0 133 0 65990 0 342165 0 32 0 59988 0 324967 0 368 0 61167 0 314832 0 1 0 60719 0 325557 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 1 0 0 60720 60720 325560 0 0 0 60720 0 325560 0 22 0 60722 0 322649 1 0 0 60720 60720 325560 1 0 0 60720 60720 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 96 0 60379 0 324598 0 1 0 60727 0 325567 0 0 0 60720 0 325560 0 81 0 60925 0 324005 0 0 0 60720 0 325560 0 26 0 60896 0 325748 0 125 0 59734 0 323385 0 304 0 62969 0 315409 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 119 0 59118 0 312275 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 312 0 58598 0 320576 0 60 0 61124 0 325246 0 587 0 59595 0 332961 0 135 0 62065 0 323010 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 514 0 78780 0 345253 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 1 0 60722 0 325559 0 0 0 60720 0 325560 0 0 0 60720 0 325560 1 58 58 61600 61600 319634 0 180 0 59635 0 319951 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 448 0 59781 0 318519 0 227 0 76644 0 343222 0 174 0 64820 0 317234 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 121 0 56178 0 314025 0 607 0 60436 0 320249 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 0 0 60720 0 325560 0 530 0 73433 0 349365 0 571 0 41477 0 289197 0 0 0 60720 0 325560 0 78 0 62700 0 329245 0 2489 0 67804 0 240869 0 131 0 59661 0 327182 0 923 0 58620 0 322876 0 72 0 60398 0 323117 0 572 0 58580 0 306351 0 397 0 62710 0 335137 0 450 0 59325 0 308271 0 622 0 60950 0 301731 0 694 0 68060 0 382409 1 3425 3425 83620 83620 279230 0 562 0 58456 0 298731 0 4917 0 52811 0 243650 1 1442 1442 121173 121173 532682 0 529 0 63870 0 319771 1 2126 2126 21001 21001 171493 0 1061 0 70415 0 347262 0 776 0 64230 0 343945 0 611 0 59190 0 311874 1 1526 1526 69351 69351 302211 0 592 0 64270 0 316708 0 1182 0 70694 0 333463 0 621 0 68005 0 344282 0 989 0 58930 0 319635 0 438 0 58320 0 301186 0 726 0 69980 0 300381 0 1303 0 69863 0 318765 1 6341 6341 63255 63255 286146 1 1164 1164 57320 57320 306844 1 3310 3310 75230 75230 307705 0 1366 0 79420 0 312448 0 965 0 73490 0 299715 0 3256 0 35250 0 373399 1 1135 1135 62285 62285 299446 0 1270 0 69206 0 325586 0 661 0 65920 0 291221 0 1013 0 69770 0 261173 0 2844 0 72683 0 255027 1 11528 11528 -14545 -14545 -78375 0 6526 0 55830 0 -58143 0 2264 0 55174 0 227033 1 4461 4461 67038 67038 235098 0 3999 0 51252 0 21267 0 35624 0 157278 0 238675 0 9252 0 79510 0 197687 0 15236 0 77440 0 418341 0 18073 0 27284 0 -297706
Names of X columns:
Group Costs GrCosts Dividends GrDiv Wealth
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