Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
146283 56 396 81 3 79 30 144 145 98364 56 297 55 4 58 28 103 101 86146 54 559 50 12 60 38 98 98 96933 89 967 125 2 108 30 135 132 79234 40 270 40 1 49 22 61 60 42551 25 143 37 3 0 26 39 38 195663 92 1562 63 0 121 25 150 144 6853 18 109 44 0 1 18 5 5 21529 63 371 88 0 20 11 28 28 95757 44 656 66 5 43 26 84 84 85584 33 511 57 0 69 25 80 79 143983 84 655 74 0 78 38 130 127 75851 88 465 49 7 86 44 82 78 59238 55 525 52 7 44 30 60 60 93163 60 885 88 3 104 40 131 131 96037 66 497 36 9 63 34 84 84 151511 154 1436 108 0 158 47 140 133 136368 53 612 43 4 102 30 151 150 112642 119 865 75 3 77 31 91 91 94728 41 385 32 0 82 23 138 132 105499 61 567 44 7 115 36 150 136 121527 58 639 85 0 101 36 124 124 127766 75 963 86 1 80 30 119 118 98958 33 398 56 5 50 25 73 70 77900 40 410 50 7 83 39 110 107 85646 92 966 135 0 123 34 123 119 98579 100 801 63 0 73 31 90 89 130767 112 892 81 5 81 31 116 112 131741 73 513 52 0 105 33 113 108 53907 40 469 44 0 47 25 56 52 178812 45 683 113 0 105 33 115 112 146761 60 643 39 3 94 35 119 116 82036 62 535 73 4 44 42 129 123 163253 75 625 48 1 114 43 127 125 27032 31 264 33 4 38 30 27 27 171975 77 992 59 2 107 33 175 162 65990 34 238 41 0 30 13 35 32 86572 46 818 69 0 71 32 64 64 159676 99 937 64 0 84 36 96 92 1929 17 70 1 0 0 0 0 0 85371 66 507 59 2 59 28 84 83 58391 30 260 32 1 33 14 41 41 31580 76 503 129 0 42 17 47 47 136815 146 927 37 2 96 32 126 120 120642 67 1269 31 10 106 30 105 105 69107 56 537 65 6 56 35 80 79 50495 107 910 107 0 57 20 70 65 108016 58 532 74 5 59 28 73 70 46341 34 345 54 4 39 28 57 55 78348 61 918 76 1 34 39 40 39 79336 119 1635 715 2 76 34 68 67 56968 42 330 57 2 20 26 21 21 93176 66 557 66 0 91 39 127 127 161632 89 1178 106 8 115 39 154 152 87850 44 740 54 3 85 33 116 113 127969 66 452 32 0 76 28 102 99 15049 24 218 20 0 8 4 7 7 155135 259 764 71 8 79 39 148 141 25109 17 255 21 5 21 18 21 21 45824 64 454 70 3 30 14 35 35 102996 41 866 112 1 76 29 112 109 160604 68 574 66 5 101 44 137 133 158051 168 1276 190 1 94 21 135 123 44547 43 379 66 1 27 16 26 26 162647 132 825 165 5 92 28 230 230 174141 105 798 56 0 123 35 181 166 60622 71 663 61 12 75 28 71 68 179566 112 1069 53 8 128 38 147 147 184301 94 921 127 8 105 23 190 179 75661 82 858 63 8 55 36 64 61 96144 70 711 38 8 56 32 105 101 129847 57 503 50 2 41 29 107 108 117286 53 382 52 0 72 25 94 90 71180 103 464 42 5 67 27 116 114 109377 121 717 76 8 75 36 106 103 85298 62 690 67 2 114 28 143 142 73631 52 462 50 5 118 23 81 79 86767 52 657 53 12 77 40 89 88 23824 32 385 39 6 22 23 26 25 93487 62 577 50 7 66 40 84 83 82981 45 619 77 2 69 28 113 113 73815 46 479 57 0 105 34 120 118 94552 63 817 73 4 116 33 110 110 132190 75 752 34 3 88 28 134 129 128754 88 430 39 6 73 34 54 51 66363 46 451 46 2 99 30 96 93 67808 53 537 63 0 62 33 78 76 61724 37 519 35 1 53 22 51 49 131722 90 1000 106 0 118 38 121 118 68580 63 637 43 5 30 26 38 38 106175 78 465 47 2 100 35 145 141 55792 25 437 31 0 49 8 59 58 25157 45 711 162 0 24 24 27 27 76669 46 299 57 5 67 29 91 91 57283 41 248 36 0 46 20 48 48 105805 144 1162 263 1 57 29 68 63 129484 82 714 78 0 75 45 58 56 72413 91 905 63 1 135 37 150 144 87831 71 649 54 1 68 33 74 73 96971 63 512 63 2 124 33 181 168 71299 53 472 77 6 33 25 65 64 77494 62 905 79 1 98 32 97 97 120336 63 786 110 4 58 29 121 117 93913 32 489 56 2 68 28 99 100 136048 39 479 56 3 81 28 152 149 181248 62 617 43 0 131 31 188 187 146123 117 925 111 10 110 52 138 127 32036 34 351 71 0 37 21 40 37 186646 92 1144 62 9 130 24 254 245 102255 93 669 56 7 93 41 87 87 168237 54 707 74 0 118 33 178 177 64219 144 458 60 0 39 32 51 49 19630 14 214 43 4 13 19 49 49 76825 61 599 68 4 74 20 73 73 115338 109 572 53 0 81 31 176 177 109427 38 897 87 0 109 31 94 94 118168 73 819 46 0 151 32 120 117 84845 75 720 105 1 51 18 66 60 153197 50 273 32 0 28 23 56 55 29877 61 508 133 1 40 17 39 39 63506 55 506 79 0 56 20 66 64 22445 77 451 51 0 27 12 27 26 47695 75 699 207 4 37 17 65 64 68370 72 407 67 0 83 30 58 58 146304 50 465 47 4 54 31 98 95 38233 32 245 34 4 27 10 25 25 42071 53 370 66 3 28 13 26 26 50517 42 316 76 0 59 22 77 76 103950 71 603 65 0 133 42 130 129 5841 10 154 9 0 12 1 11 11 2341 35 229 42 5 0 9 2 2 84396 65 577 45 0 106 32 101 101 24610 25 192 25 4 23 11 31 28 35753 66 617 115 0 44 25 36 36 55515 41 411 97 0 71 36 120 89 209056 86 975 53 1 116 31 195 193 6622 16 146 2 0 4 0 4 4 115814 42 705 52 5 62 24 89 84 11609 19 184 44 0 12 13 24 23 13155 19 200 22 0 18 8 39 39 18274 45 274 35 0 14 13 14 14 72875 65 502 74 0 60 19 78 78 10112 35 382 103 0 7 18 15 14 142775 95 964 144 2 98 33 106 101 68847 49 537 60 7 64 40 83 82 17659 37 438 134 1 29 22 24 24 20112 64 369 89 8 32 38 37 36 61023 38 417 42 2 25 24 77 75 13983 34 276 52 0 16 8 16 16 65176 32 514 98 2 48 35 56 55 132432 65 822 99 0 100 43 132 131 112494 52 389 52 0 46 43 144 131 45109 62 466 29 1 45 14 40 39 170875 65 1255 125 3 129 41 153 144 180759 83 694 106 0 130 38 143 139 214921 95 1024 95 3 136 45 220 211 100226 29 400 40 0 59 31 79 78 32043 18 397 140 0 25 13 50 50 54454 33 350 43 0 32 28 39 39 78876 247 719 128 4 63 31 95 90 170745 139 1277 142 4 95 40 169 166 6940 29 356 73 11 14 30 12 12 49025 118 457 72 0 36 16 63 57 122037 110 1402 128 0 113 37 134 133 53782 67 600 61 4 47 30 69 69 127748 42 480 73 0 92 35 119 119 86839 65 595 148 1 70 32 119 119 44830 94 436 64 0 19 27 75 65 77395 64 230 45 0 50 20 63 61 89324 81 651 58 0 41 18 55 49 103300 95 1367 97 9 91 31 103 101 112283 67 564 50 1 111 31 197 196 10901 63 716 37 3 41 21 16 15 120691 83 747 50 10 120 39 140 136 58106 45 467 105 5 135 41 89 89 57140 30 671 69 0 27 13 40 40 122422 70 861 46 2 87 32 125 123 25899 32 319 57 0 25 18 21 21 139296 83 612 52 1 131 39 167 163 52678 31 433 98 2 45 14 32 29 23853 67 434 61 4 29 7 36 35 17306 66 503 89 0 58 17 13 13 7953 10 85 0 0 4 0 5 5 89455 70 564 48 2 47 30 96 96 147866 103 824 91 1 109 37 151 151 4245 5 74 0 0 7 0 6 6 21509 20 259 7 0 12 5 13 13 7670 5 69 3 0 0 1 3 3 66675 36 535 54 1 37 16 57 56 14336 34 239 70 0 37 32 23 23 53608 48 438 36 2 46 24 61 57 30059 40 459 37 0 15 17 21 14 29668 43 426 123 3 42 11 43 43 22097 31 288 247 6 7 24 20 20 96841 42 498 46 0 54 22 82 72 41907 46 454 72 2 54 12 90 87 27080 33 376 41 0 14 19 25 21 35885 18 225 24 2 16 13 60 56 41247 55 555 45 1 33 17 61 59 28313 35 252 33 1 32 15 85 82 36845 59 208 27 2 21 16 43 43 16548 19 130 36 1 15 24 25 25 36134 66 481 87 0 38 15 41 38 55764 60 389 90 1 22 17 26 25 28910 36 565 114 3 28 18 38 38 13339 25 173 31 0 10 20 12 12 25319 47 278 45 0 31 16 29 29 66956 54 609 69 0 32 16 49 47 47487 53 422 51 0 32 18 46 45 52785 40 445 34 1 43 22 41 40 44683 40 387 60 4 27 8 31 30 35619 39 339 45 0 37 17 41 41 21920 14 181 54 0 20 18 26 25 45608 45 245 25 0 32 16 23 23 7721 36 384 38 7 0 23 14 14 20634 28 212 52 2 5 22 16 16 29788 44 399 67 0 26 13 25 26 31931 30 229 74 7 10 13 21 21 37754 22 224 38 3 27 16 32 27 32505 17 203 30 0 11 16 9 9 40557 31 333 26 0 29 20 35 33 94238 55 384 67 6 25 22 42 42 44197 54 636 132 2 55 17 68 68 43228 21 185 42 0 23 18 32 32 4103 14 93 35 0 5 17 6 6 44144 81 581 118 3 43 12 68 67 32868 35 248 68 0 23 7 33 33 27640 43 304 43 1 34 17 84 77 14063 46 344 76 1 36 14 46 46 28990 30 407 64 0 35 23 30 30 4694 23 170 48 1 0 17 0 0 42648 38 312 64 0 37 14 36 36 64329 54 507 56 0 28 15 47 46 21928 20 224 71 0 16 17 20 18 25836 53 340 75 0 26 21 50 48 22779 45 168 39 0 38 18 30 29 40820 39 443 42 0 23 18 30 28 27530 20 204 39 0 22 17 34 34 32378 24 367 93 0 30 17 33 33 10824 31 210 38 0 16 16 34 34 39613 35 335 60 0 18 15 37 33 60865 151 364 71 0 28 21 83 80 19787 52 178 52 0 32 16 32 32 20107 30 206 27 2 21 14 30 30 36605 31 279 59 0 23 15 43 41 40961 29 387 40 1 29 17 41 41 48231 57 490 79 1 50 15 51 51 39725 40 238 44 0 12 15 19 18 21455 44 343 65 0 21 10 37 34 23430 25 232 10 0 18 6 33 31 62991 77 530 124 0 27 22 41 39 49363 35 291 81 0 41 21 54 54 9604 11 67 15 0 13 1 14 14 24552 63 397 92 1 12 18 25 24 31493 44 467 42 0 21 17 25 24 3439 19 178 10 0 8 4 8 8 19555 13 175 24 0 26 10 26 26 21228 42 299 64 0 27 16 20 19 23177 38 154 45 1 13 16 11 11 22094 29 106 22 0 16 9 14 14 2342 20 189 56 0 2 16 3 1 38798 27 194 94 0 42 17 40 39 3255 20 135 19 0 5 7 5 5 24261 19 201 35 0 37 15 38 37 18511 37 207 32 0 17 14 32 32 40798 26 280 35 0 38 14 41 38 28893 42 260 48 0 37 18 46 47 21425 49 227 49 0 29 12 47 47 50276 30 239 48 0 32 16 37 37 37643 49 333 62 0 35 21 51 51 30377 67 428 96 1 17 19 49 45 27126 28 230 45 0 20 16 21 21 13 19 292 63 0 7 1 1 1 42097 49 350 71 1 46 16 44 42 24451 27 186 26 0 24 10 26 26 14335 30 326 48 6 40 19 21 21 5084 22 155 29 3 3 12 4 4 9927 12 75 19 1 10 2 10 10 43527 31 361 45 2 37 14 43 43 27184 20 261 45 0 17 17 34 34 21610 20 299 67 0 28 19 32 31 20484 39 300 30 0 19 14 20 19 20156 29 450 36 3 29 11 34 34 6012 16 183 34 1 8 4 6 6 18475 27 238 36 0 10 16 12 11 12645 21 165 34 0 15 20 24 24 11017 19 234 37 1 15 12 16 16 37623 35 176 46 0 28 15 72 72 35873 14 329 44 0 17 16 27 21
Names of X columns:
totseconds logins compendium_views_info compendium_views_pr shared_compendiums blogged_computations compendiums_reviewed tothyperlinks totblogs
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