Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
210907 56 396 79 30 115 112285 145 120982 56 297 58 28 109 84786 101 176508 54 559 60 38 146 83123 98 179321 89 967 108 30 116 101193 132 123185 40 270 49 22 68 38361 60 52746 25 143 0 26 101 68504 38 385534 92 1562 121 25 96 119182 144 33170 18 109 1 18 67 22807 5 101645 63 371 20 11 44 17140 28 149061 44 656 43 26 100 116174 84 165446 33 511 69 25 93 57635 79 237213 84 655 78 38 140 66198 127 173326 88 465 86 44 166 71701 78 133131 55 525 44 30 99 57793 60 258873 60 885 104 40 139 80444 131 180083 66 497 63 34 130 53855 84 324799 154 1436 158 47 181 97668 133 230964 53 612 102 30 116 133824 150 236785 119 865 77 31 116 101481 91 135473 41 385 82 23 88 99645 132 202925 61 567 115 36 139 114789 136 215147 58 639 101 36 135 99052 124 344297 75 963 80 30 108 67654 118 153935 33 398 50 25 89 65553 70 132943 40 410 83 39 156 97500 107 174724 92 966 123 34 129 69112 119 174415 100 801 73 31 118 82753 89 225548 112 892 81 31 118 85323 112 223632 73 513 105 33 125 72654 108 124817 40 469 47 25 95 30727 52 221698 45 683 105 33 126 77873 112 210767 60 643 94 35 135 117478 116 170266 62 535 44 42 154 74007 123 260561 75 625 114 43 165 90183 125 84853 31 264 38 30 113 61542 27 294424 77 992 107 33 127 101494 162 101011 34 238 30 13 52 27570 32 215641 46 818 71 32 121 55813 64 325107 99 937 84 36 136 79215 92 7176 17 70 0 0 0 1423 0 167542 66 507 59 28 108 55461 83 106408 30 260 33 14 46 31081 41 96560 76 503 42 17 54 22996 47 265769 146 927 96 32 124 83122 120 269651 67 1269 106 30 115 70106 105 149112 56 537 56 35 128 60578 79 175824 107 910 57 20 80 39992 65 152871 58 532 59 28 97 79892 70 111665 34 345 39 28 104 49810 55 116408 61 918 34 39 59 71570 39 362301 119 1635 76 34 125 100708 67 78800 42 330 20 26 82 33032 21 183167 66 557 91 39 149 82875 127 277965 89 1178 115 39 149 139077 152 150629 44 740 85 33 122 71595 113 168809 66 452 76 28 118 72260 99 24188 24 218 8 4 12 5950 7 329267 259 764 79 39 144 115762 141 65029 17 255 21 18 67 32551 21 101097 64 454 30 14 52 31701 35 218946 41 866 76 29 108 80670 109 244052 68 574 101 44 166 143558 133 341570 168 1276 94 21 80 117105 123 103597 43 379 27 16 60 23789 26 233328 132 825 92 28 107 120733 230 256462 105 798 123 35 127 105195 166 206161 71 663 75 28 107 73107 68 311473 112 1069 128 38 146 132068 147 235800 94 921 105 23 84 149193 179 177939 82 858 55 36 141 46821 61 207176 70 711 56 32 123 87011 101 196553 57 503 41 29 111 95260 108 174184 53 382 72 25 98 55183 90 143246 103 464 67 27 105 106671 114 187559 121 717 75 36 135 73511 103 187681 62 690 114 28 107 92945 142 119016 52 462 118 23 85 78664 79 182192 52 657 77 40 155 70054 88 73566 32 385 22 23 88 22618 25 194979 62 577 66 40 155 74011 83 167488 45 619 69 28 104 83737 113 143756 46 479 105 34 132 69094 118 275541 63 817 116 33 127 93133 110 243199 75 752 88 28 108 95536 129 182999 88 430 73 34 129 225920 51 135649 46 451 99 30 116 62133 93 152299 53 537 62 33 122 61370 76 120221 37 519 53 22 85 43836 49 346485 90 1000 118 38 147 106117 118 145790 63 637 30 26 99 38692 38 193339 78 465 100 35 87 84651 141 80953 25 437 49 8 28 56622 58 122774 45 711 24 24 90 15986 27 130585 46 299 67 29 109 95364 91 112611 41 248 46 20 78 26706 48 286468 144 1162 57 29 111 89691 63 241066 82 714 75 45 158 67267 56 148446 91 905 135 37 141 126846 144 204713 71 649 68 33 122 41140 73 182079 63 512 124 33 124 102860 168 140344 53 472 33 25 93 51715 64 220516 62 905 98 32 124 55801 97 243060 63 786 58 29 112 111813 117 162765 32 489 68 28 108 120293 100 182613 39 479 81 28 99 138599 149 232138 62 617 131 31 117 161647 187 265318 117 925 110 52 199 115929 127 85574 34 351 37 21 78 24266 37 310839 92 1144 130 24 91 162901 245 225060 93 669 93 41 158 109825 87 232317 54 707 118 33 126 129838 177 144966 144 458 39 32 122 37510 49 43287 14 214 13 19 71 43750 49 155754 61 599 74 20 75 40652 73 164709 109 572 81 31 115 87771 177 201940 38 897 109 31 119 85872 94 235454 73 819 151 32 124 89275 117 220801 75 720 51 18 72 44418 60 99466 50 273 28 23 91 192565 55 92661 61 508 40 17 45 35232 39 133328 55 506 56 20 78 40909 64 61361 77 451 27 12 39 13294 26 125930 75 699 37 17 68 32387 64 100750 72 407 83 30 119 140867 58 224549 50 465 54 31 117 120662 95 82316 32 245 27 10 39 21233 25 102010 53 370 28 13 50 44332 26 101523 42 316 59 22 88 61056 76 243511 71 603 133 42 155 101338 129 22938 10 154 12 1 0 1168 11 41566 35 229 0 9 36 13497 2 152474 65 577 106 32 123 65567 101 61857 25 192 23 11 32 25162 28 99923 66 617 44 25 99 32334 36 132487 41 411 71 36 136 40735 89 317394 86 975 116 31 117 91413 193 21054 16 146 4 0 0 855 4 209641 42 705 62 24 88 97068 84 22648 19 184 12 13 39 44339 23 31414 19 200 18 8 25 14116 39 46698 45 274 14 13 52 10288 14 131698 65 502 60 19 75 65622 78 91735 35 382 7 18 71 16563 14 244749 95 964 98 33 124 76643 101 184510 49 537 64 40 151 110681 82 79863 37 438 29 22 71 29011 24 128423 64 369 32 38 145 92696 36 97839 38 417 25 24 87 94785 75 38214 34 276 16 8 27 8773 16 151101 32 514 48 35 131 83209 55 272458 65 822 100 43 162 93815 131 172494 52 389 46 43 165 86687 131 108043 62 466 45 14 54 34553 39 328107 65 1255 129 41 159 105547 144 250579 83 694 130 38 147 103487 139 351067 95 1024 136 45 170 213688 211 158015 29 400 59 31 119 71220 78 98866 18 397 25 13 49 23517 50 85439 33 350 32 28 104 56926 39 229242 247 719 63 31 120 91721 90 351619 139 1277 95 40 150 115168 166 84207 29 356 14 30 112 111194 12 120445 118 457 36 16 59 51009 57 324598 110 1402 113 37 136 135777 133 131069 67 600 47 30 107 51513 69 204271 42 480 92 35 130 74163 119 165543 65 595 70 32 115 51633 119 141722 94 436 19 27 107 75345 65 116048 64 230 50 20 75 33416 61 250047 81 651 41 18 71 83305 49 299775 95 1367 91 31 120 98952 101 195838 67 564 111 31 116 102372 196 173260 63 716 41 21 79 37238 15 254488 83 747 120 39 150 103772 136 104389 45 467 135 41 156 123969 89 136084 30 671 27 13 51 27142 40 199476 70 861 87 32 118 135400 123 92499 32 319 25 18 71 21399 21 224330 83 612 131 39 144 130115 163 135781 31 433 45 14 47 24874 29 74408 67 434 29 7 28 34988 35 81240 66 503 58 17 68 45549 13 14688 10 85 4 0 0 6023 5 181633 70 564 47 30 110 64466 96 271856 103 824 109 37 147 54990 151 7199 5 74 7 0 0 1644 6 46660 20 259 12 5 15 6179 13 17547 5 69 0 1 4 3926 3 133368 36 535 37 16 64 32755 56 95227 34 239 37 32 111 34777 23 152601 48 438 46 24 85 73224 57 98146 40 459 15 17 68 27114 14 79619 43 426 42 11 40 20760 43 59194 31 288 7 24 80 37636 20 139942 42 498 54 22 88 65461 72 118612 46 454 54 12 48 30080 87 72880 33 376 14 19 76 24094 21 65475 18 225 16 13 51 69008 56 99643 55 555 33 17 67 54968 59 71965 35 252 32 15 59 46090 82 77272 59 208 21 16 61 27507 43 49289 19 130 15 24 76 10672 25 135131 66 481 38 15 60 34029 38 108446 60 389 22 17 68 46300 25 89746 36 565 28 18 71 24760 38 44296 25 173 10 20 76 18779 12 77648 47 278 31 16 62 21280 29 181528 54 609 32 16 61 40662 47 134019 53 422 32 18 67 28987 45 124064 40 445 43 22 88 22827 40 92630 40 387 27 8 30 18513 30 121848 39 339 37 17 64 30594 41 52915 14 181 20 18 68 24006 25 81872 45 245 32 16 64 27913 23 58981 36 384 0 23 91 42744 14 53515 28 212 5 22 88 12934 16 60812 44 399 26 13 52 22574 26 56375 30 229 10 13 49 41385 21 65490 22 224 27 16 62 18653 27 80949 17 203 11 16 61 18472 9 76302 31 333 29 20 76 30976 33 104011 55 384 25 22 88 63339 42 98104 54 636 55 17 66 25568 68 67989 21 185 23 18 71 33747 32 30989 14 93 5 17 68 4154 6 135458 81 581 43 12 48 19474 67 73504 35 248 23 7 25 35130 33 63123 43 304 34 17 68 39067 77 61254 46 344 36 14 41 13310 46 74914 30 407 35 23 90 65892 30 31774 23 170 0 17 66 4143 0 81437 38 312 37 14 54 28579 36 87186 54 507 28 15 59 51776 46 50090 20 224 16 17 60 21152 18 65745 53 340 26 21 77 38084 48 56653 45 168 38 18 68 27717 29 158399 39 443 23 18 72 32928 28 46455 20 204 22 17 67 11342 34 73624 24 367 30 17 64 19499 33 38395 31 210 16 16 63 16380 34 91899 35 335 18 15 59 36874 33 139526 151 364 28 21 84 48259 80 52164 52 178 32 16 64 16734 32 51567 30 206 21 14 56 28207 30 70551 31 279 23 15 54 30143 41 84856 29 387 29 17 67 41369 41 102538 57 490 50 15 58 45833 51 86678 40 238 12 15 59 29156 18 85709 44 343 21 10 40 35944 34 34662 25 232 18 6 22 36278 31 150580 77 530 27 22 83 45588 39 99611 35 291 41 21 81 45097 54 19349 11 67 13 1 2 3895 14 99373 63 397 12 18 72 28394 24 86230 44 467 21 17 61 18632 24 30837 19 178 8 4 15 2325 8 31706 13 175 26 10 32 25139 26 89806 42 299 27 16 62 27975 19 62088 38 154 13 16 58 14483 11 40151 29 106 16 9 36 13127 14 27634 20 189 2 16 59 5839 1 76990 27 194 42 17 68 24069 39 37460 20 135 5 7 21 3738 5 54157 19 201 37 15 55 18625 37 49862 37 207 17 14 54 36341 32 84337 26 280 38 14 55 24548 38 64175 42 260 37 18 72 21792 47 59382 49 227 29 12 41 26263 47 119308 30 239 32 16 61 23686 37 76702 49 333 35 21 67 49303 51 103425 67 428 17 19 76 25659 45 70344 28 230 20 16 64 28904 21 43410 19 292 7 1 3 2781 1 104838 49 350 46 16 63 29236 42 62215 27 186 24 10 40 19546 26 69304 30 326 40 19 69 22818 21 53117 22 155 3 12 48 32689 4 19764 12 75 10 2 8 5752 10 86680 31 361 37 14 52 22197 43 84105 20 261 17 17 66 20055 34 77945 20 299 28 19 76 25272 31 89113 39 300 19 14 43 82206 19 91005 29 450 29 11 39 32073 34 40248 16 183 8 4 14 5444 6 64187 27 238 10 16 61 20154 11 50857 21 165 15 20 71 36944 24 56613 19 234 15 12 44 8019 16 62792 35 176 28 15 60 30884 72 72535 14 329 17 16 64 19540 21
Names of X columns:
Timeinrfc Logins Compendiumsviewed Bloggedcomputations Reviews Feedbackmessages Compendiumsize Totalblogs
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