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