Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 162556 1081 807 213118 230380558 6282929 1 29790 309 444 81767 25266003 4324047 1 87550 458 412 153198 70164684 4108272 0 84738 588 428 -26007 -15292116 -1212617 1 54660 299 312 126942 37955658 1485329 1 42634 156 166 157214 24525384 1779876 0 40949 481 263 129352 62218312 1367203 1 42312 323 237 234817 75845891 2519076 1 37704 452 228 60448 27322496 912684 1 16275 109 129 47818 5212162 1443586 0 25830 115 104 245546 28237790 1220017 0 12679 110 122 48020 5282200 984885 1 18014 239 393 -1710 -408690 1457425 0 43556 247 190 32648 8064056 -572920 1 24524 497 275 95350 47388950 929144 0 6532 103 62 151352 15589256 1151176 0 7123 109 102 288170 31410530 790090 1 20813 502 255 114337 57397174 774497 1 37597 248 234 37884 9395232 990576 0 17821 373 277 122844 45820812 454195 1 12988 119 73 82340 9798460 876607 1 22330 84 67 79801 6703284 711969 0 13326 102 103 165548 16885896 702380 0 16189 295 290 116384 34333280 264449 0 7146 105 83 134028 14072940 450033 0 15824 64 56 63838 4085632 541063 1 26088 267 224 74996 20023932 588864 0 11326 129 64 31080 4009320 -37216 0 8568 37 34 32168 1190216 783310 0 14416 361 139 49857 17998377 467359 1 3369 28 26 87161 2440508 688779 1 11819 85 70 106113 9019605 608419 1 6620 44 40 80570 3545080 696348 1 4519 49 42 102129 5004321 597793 0 2220 22 12 301670 6636740 821730 0 18562 155 211 102313 15858515 377934 0 10327 91 74 88577 8060507 651939 1 5336 81 80 112477 9110637 697458 1 2365 79 83 191778 15150462 700368 0 4069 145 131 79804 11571580 225986 0 7710 816 187 128294 104687904 348695 0 13718 61 56 96448 5883328 373683 0 4525 226 89 93811 21201286 501709 0 6869 105 88 117520 12339600 413743 0 4628 62 39 69159 4287858 379825 1 3653 24 24 101792 2443008 336260 1 1265 26 19 210568 5474768 636765 1 7489 322 149 136996 44112712 481231 0 4901 84 58 121920 10241280 469107 0 2284 33 41 76403 2521299 211928 1 3160 108 90 108094 11674152 563925 1 4150 150 136 134759 20213850 511939 1 7285 115 97 188873 21720395 521016 1 1134 162 63 146216 23686992 543856 1 4658 158 114 156608 24744064 329304 0 2384 97 77 61348 5950756 423262 0 3748 9 6 50350 453150 509665 0 5371 66 47 87720 5789520 455881 0 1285 107 51 99489 10645323 367772 1 9327 101 85 87419 8829319 406339 1 5565 47 43 94355 4434685 493408 0 1528 38 32 60326 2292388 232942 1 3122 34 25 94670 3218780 416002 1 7317 84 74 82425 6923700 337430 0 2675 79 54 59017 4662343 361517 0 13253 947 251 90829 86015063 360962 0 880 74 15 80791 5978534 235561 1 2053 53 44 100423 5322419 408247 0 1424 94 73 131116 12324904 450296 1 4036 63 85 100269 6316947 418799 1 3045 58 49 27330 1585140 247405 0 5119 49 38 39039 1912911 378519 0 1431 34 35 106885 3634090 326638 0 554 11 9 79285 872135 328233 0 1975 35 34 118881 4160835 386225 1 1286 17 17 77623 1319591 283662 0 1012 47 29 114768 5394096 370225 0 810 43 11 74015 3182645 269236 0 1280 117 52 69465 8127405 365732 1 666 171 13 117869 20155599 420383 0 1380 26 29 60982 1585532 345811 1 4608 73 64 90131 6579563 431809 0 876 59 33 138971 8199289 418876 0 814 18 15 39625 713250 297476 0 514 15 15 102725 1540875 416776 1 5692 72 68 64239 4625208 357257 0 3642 86 100 90262 7762532 458343 0 540 14 13 103960 1455440 388386 0 2099 64 45 106611 6823104 358934 0 567 11 14 103345 1136795 407560 0 2001 52 36 95551 4968652 392558 1 2949 41 40 82903 3399023 373177 0 2253 99 68 63593 6295707 428370 1 6533 75 29 126910 9518250 369419 0 1889 45 43 37527 1688715 358649 1 3055 43 30 60247 2590621 376641 0 272 8 9 112995 903960 467427 1 1414 198 22 70184 13896432 364885 0 2564 22 19 130140 2863080 436230 1 1383 11 9 73221 805431 329118 1 1261 33 31 76114 2511762 317365 0 975 23 19 90534 2082282 286849 0 3366 80 55 108479 8678320 376685 0 576 18 8 113761 2047698 407198 0 1306 28 26 68696 1923488 377772 0 746 23 29 71561 1645903 271483 1 3192 60 48 59831 3589860 153661 1 2045 20 16 97890 1957800 513294 0 5477 59 45 101481 5987379 324881 1 1932 36 20 72954 2626344 264512 0 936 30 22 67939 2038170 420968 1 3437 47 33 48022 2257034 129302 0 5131 71 44 86111 6113881 191521 1 2397 14 13 74020 1036280 268673 1 1389 9 6 57530 517770 353179 0 1503 39 35 56364 2198196 354624 0 402 26 8 84990 2209740 363713 0 2239 21 17 88590 1860390 456657 1 2234 16 11 77200 1235200 211742 0 837 69 21 61262 4227078 338381 0 10579 92 92 110309 10148428 418530 0 875 14 12 67000 938000 351483 0 1395 103 108 93099 9589197 372928 1 1659 29 25 107577 3119733 485538 1 2647 37 17 62920 2328040 279268 1 3294 23 23 75832 1744136 219060 0 0 0 0 60720 0 325560 0 94 7 10 60793 425551 325314 0 422 28 23 57935 1622180 322046 0 0 0 0 60720 0 325560 0 34 8 7 60630 485040 325599 0 1558 63 25 55637 3505131 377028 0 0 0 0 60720 0 325560 0 43 3 20 60887 182661 323850 0 0 0 0 60720 0 325560 0 316 9 4 60505 544545 331514 0 115 13 10 60945 792285 325632 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 389 14 7 58990 825860 322265 0 0 0 0 60720 0 325560 0 1002 15 11 56750 851250 325906 0 36 3 4 60894 182682 325985 0 460 15 15 63346 950190 346145 0 309 11 9 56535 621885 325898 0 0 0 0 60720 0 325560 0 9 6 7 60835 365010 325356 0 0 0 0 60720 0 325560 0 14 1 0 61016 61016 325930 0 520 10 7 58650 586500 318020 0 1766 73 46 60438 4411974 326389 0 0 0 5 60720 0 325560 0 458 11 7 58625 644875 302925 0 20 3 2 60938 182814 325540 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 98 2 2 61490 122980 326736 0 405 7 5 60845 425915 340580 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 483 27 7 60830 1642410 331828 0 454 51 24 63261 3226311 323299 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 757 19 18 45689 868091 387722 0 0 0 4 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 36 4 3 61564 246256 324598 0 0 0 0 60720 0 325560 0 203 9 9 61938 557442 328726 0 0 0 0 60720 0 325560 0 90 8 6 60951 487608 325043 0 0 0 0 60720 0 325560 1 71 1 0 60745 60745 325806 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 972 34 19 71642 2435828 387732 1 531 10 11 71641 716410 349729 0 604 38 8 55792 2120096 332202 1 283 10 15 71873 718730 305442 1 23 5 5 62555 312775 329537 1 638 14 11 60370 845180 327055 1 699 16 10 64873 1037968 356245 0 149 5 6 62041 310205 328451 0 226 5 5 65745 328725 307062 0 0 0 0 60720 0 325560 0 275 4 7 59500 238000 331345 0 0 0 0 60720 0 325560 0 141 6 7 61630 369780 331824 0 0 0 0 60720 0 325560 0 28 2 3 60890 121780 325685 1 0 0 0 60720 0 325560 1 2566 80 74 113521 9081680 404480 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 1 472 20 19 80045 1600900 318314 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 1 203 27 12 50804 1371708 311807 1 496 17 12 87390 1485630 337724 1 10 2 5 61656 123312 326431 1 63 4 2 65688 262752 327556 1 0 0 0 60720 0 325560 1 1136 32 26 48522 1552704 356850 1 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 267 32 11 57640 1844480 322741 0 474 20 10 61977 1239540 310902 0 534 7 5 62620 438340 324295 0 0 0 2 60720 0 325560 0 15 8 6 60831 486648 326156 0 397 28 7 60646 1698088 326960 0 0 0 2 60720 0 325560 0 1061 20 22 56225 1124500 333411 0 288 4 3 60510 242040 297761 0 0 0 0 60720 0 325560 0 3 2 1 60698 121396 325536 0 0 0 6 60720 0 325560 0 20 2 1 60805 121610 325762 0 278 26 22 61404 1596504 327957 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 192 4 2 65276 261104 318521 0 0 0 0 60720 0 325560 0 317 9 7 63915 575235 319775 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 368 17 13 61686 1048662 332128 0 0 0 0 60720 0 325560 0 2 1 0 60743 60743 325486 0 0 0 0 60720 0 325560 0 53 6 6 60349 362094 325838 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 94 3 3 61360 184080 331767 0 0 0 0 60720 0 325560 0 24 8 7 59818 478544 324523 0 2332 4 2 72680 290720 339995 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 131 11 15 61808 679888 319582 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 206 9 9 53110 477990 307245 1 0 0 0 60720 0 325560 0 167 2 1 64245 128490 317967 0 622 73 38 73007 5329511 331488 0 885 85 49 82732 7032220 335452 0 0 0 0 60720 0 325560 0 365 8 6 54820 438560 334184 0 364 35 26 47705 1669675 313213 1 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 226 12 13 72835 874020 348678 0 307 15 10 58856 882840 328727 1 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 188 11 9 77655 854205 387978 1 0 0 0 60720 0 325560 0 138 6 26 69817 418902 336704 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 125 12 19 60798 729576 322076 0 0 0 0 60720 0 325560 0 282 30 12 62452 1873560 334272 0 335 33 23 64175 2117775 338197 0 0 0 0 60720 0 325560 1 813 82 24 67440 5530080 321024 0 176 28 8 68136 1907808 322145 1 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 249 72 26 56726 4084272 323351 0 0 0 0 60720 0 325560 0 333 13 9 70811 920543 327748 0 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 30 4 3 62045 248180 328157 1 0 0 0 60720 0 325560 0 249 62 13 54323 3368026 311594 1 0 0 0 60720 0 325560 0 165 24 12 62841 1508184 335962 0 453 21 19 81125 1703625 372426 0 0 0 0 60720 0 325560 0 53 14 10 59506 833084 319844 1 382 21 9 59365 1246665 355822 0 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 9 60720 0 325560 1 30 4 4 60798 243192 324047 0 290 2 1 58790 117580 311464 0 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 366 53 14 61808 3275824 353417 0 2 9 12 60735 546615 325590 1 0 0 0 60720 0 325560 1 209 13 19 64016 832208 328576 0 384 22 17 54683 1203026 326126 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 365 83 32 87192 7236936 369376 1 0 0 0 60720 0 325560 1 49 8 14 64107 512856 332013 0 3 4 8 60761 243044 325871 0 133 14 4 65990 923860 342165 0 32 1 0 59988 59988 324967 0 368 17 20 61167 1039839 314832 0 1 6 5 60719 364314 325557 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 22 2 1 60722 121444 322649 1 0 0 0 60720 0 325560 1 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 96 5 4 60379 301895 324598 0 1 2 1 60727 121454 325567 0 0 0 0 60720 0 325560 0 81 7 4 60925 426475 324005 0 0 0 0 60720 0 325560 0 26 1 1 60896 60896 325748 0 125 13 10 59734 776542 323385 0 304 15 12 62969 944535 315409 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 119 6 3 59118 354708 312275 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 312 14 3 58598 820372 320576 0 60 10 7 61124 611240 325246 0 587 12 10 59595 715140 332961 0 135 2 1 62065 124130 323010 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 514 52 15 78780 4096560 345253 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 1 4 4 60722 242888 325559 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 1 58 3 2 61600 184800 319634 0 180 11 9 59635 655985 319951 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 448 40 7 59781 2391240 318519 0 227 9 7 76644 689796 343222 0 174 1 3 64820 64820 317234 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 121 24 11 56178 1348272 314025 0 607 11 7 60436 664796 320249 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 0 0 0 60720 0 325560 0 530 60 18 73433 4405980 349365 0 571 80 14 41477 3318160 289197 0 0 0 0 60720 0 325560 0 78 16 12 62700 1003200 329245 0 2489 40 29 67804 2712160 240869 0 131 6 3 59661 357966 327182 0 923 8 6 58620 468960 322876 0 72 3 3 60398 181194 323117 0 572 16 8 58580 937280 306351 0 397 10 10 62710 627100 335137 0 450 8 6 59325 474600 308271 0 622 7 8 60950 426650 301731 0 694 8 6 68060 544480 382409 1 3425 12 9 83620 1003440 279230 0 562 13 8 58456 759928 298731 0 4917 42 26 52811 2218062 243650 1 1442 118 239 121173 14298414 532682 0 529 9 7 63870 574830 319771 1 2126 138 41 21001 2898138 171493 0 1061 5 3 70415 352075 347262 0 776 9 8 64230 578070 343945 0 611 8 6 59190 473520 311874 1 1526 25 21 69351 1733775 302211 0 592 7 7 64270 449890 316708 0 1182 13 11 70694 919022 333463 0 621 16 11 68005 1088080 344282 0 989 11 12 58930 648230 319635 0 438 11 9 58320 641520 301186 0 726 3 3 69980 209940 300381 0 1303 61 57 69863 4261643 318765 1 6341 24 16 63255 1518120 286146 1 1164 17 15 57320 974440 306844 1 3310 33 32 75230 2482590 307705 0 1366 7 5 79420 555940 312448 0 965 3 2 73490 220470 299715 0 3256 66 23 35250 2326500 373399 1 1135 17 20 62285 1058845 299446 0 1270 26 24 69206 1799356 325586 0 661 3 1 65920 197760 291221 0 1013 2 1 69770 139540 261173 0 2844 67 74 72683 4869761 255027 1 11528 70 68 -14545 -1018150 -78375 0 6526 26 20 55830 1451580 -58143 0 2264 24 20 55174 1324176 227033 1 4461 94 79 67038 6301572 235098 0 3999 30 21 51252 1537560 21267 0 35624 223 244 157278 35072994 238675 0 9252 48 32 79510 3816480 197687 0 15236 90 86 77440 6969600 418341 0 18073 180 69 27284 4911120 -297706
Names of X columns:
Group Costs Trades Orders Dividends TrDiv 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