Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 1 1 130631 130631 58 58 58 58 30 30 117 117 98956 98956 1 2 2 189326 189326 94 94 108 108 34 34 132 132 106816 106816 1 3 3 65295 65295 27 27 8 8 31 31 117 117 76173 76173 1 4 4 33186 33186 19 19 1 1 18 18 67 67 22807 22807 1 5 5 261949 261949 95 95 86 86 38 38 140 140 72117 72117 1 6 6 190794 190794 94 94 93 93 49 49 186 186 79738 79738 1 7 7 150733 150733 47 47 92 92 25 25 96 96 114029 114029 1 8 8 223226 223226 67 67 119 119 42 42 163 163 124550 124550 1 9 9 384138 384138 86 86 86 86 35 35 128 128 72875 72875 1 10 10 156540 156540 34 34 50 50 25 25 89 89 81964 81964 1 11 11 191441 191441 105 105 81 81 35 35 134 134 96740 96740 1 12 12 259667 259667 53 53 113 113 37 37 142 142 90694 90694 1 13 13 228871 228871 64 64 109 109 40 40 155 155 125369 125369 1 14 14 87485 87485 33 33 38 38 33 33 125 125 65702 65702 1 15 15 322865 322865 82 82 111 111 35 35 135 135 108179 108179 1 16 16 340093 340093 103 103 92 92 37 37 139 139 95066 95066 1 17 17 155363 155363 59 59 66 66 35 35 128 128 68966 68966 1 18 18 174198 174198 65 65 69 69 32 32 107 107 88766 88766 1 19 19 392666 392666 132 132 91 91 34 34 125 125 109249 109249 1 20 20 302674 302674 99 99 129 129 44 44 169 169 146648 146648 1 21 21 164733 164733 50 50 93 93 40 40 145 145 80613 80613 1 22 22 24188 24188 24 24 8 8 4 4 12 12 5950 5950 1 23 23 340411 340411 274 274 79 79 41 41 151 151 131106 131106 1 24 24 65029 65029 17 17 21 21 18 18 67 67 32551 32551 1 25 25 101097 101097 64 64 30 30 14 14 52 52 31701 31701 1 26 26 282220 282220 160 160 106 106 32 32 120 120 143950 143950 1 27 27 273495 273495 118 118 127 127 37 37 135 135 112368 112368 1 28 28 214872 214872 74 74 75 75 32 32 123 123 82124 82124 1 29 29 184474 184474 87 87 55 55 38 38 149 149 55062 55062 1 30 30 205675 205675 60 60 43 43 33 33 125 125 105612 105612 1 31 31 197760 197760 66 66 114 114 32 32 123 123 108461 108461 1 32 32 73566 73566 32 32 22 22 23 23 88 88 22618 22618 1 33 33 177949 177949 48 48 77 77 34 34 120 120 92059 92059 1 34 34 148698 148698 49 49 105 105 34 34 132 132 77993 77993 1 35 35 300103 300103 69 69 119 119 38 38 144 144 104155 104155 1 36 36 251437 251437 78 78 88 88 32 32 124 124 109840 109840 1 37 37 132672 132672 41 41 58 58 25 25 97 97 48194 48194 1 38 38 376465 376465 99 99 132 132 40 40 155 155 134796 134796 1 39 39 135042 135042 47 47 67 67 32 32 120 120 113402 113402 1 40 40 300074 300074 152 152 57 57 33 33 127 127 97967 97967 1 41 41 271757 271757 93 93 95 95 50 50 178 178 74844 74844 1 42 42 150949 150949 95 95 139 139 37 37 141 141 136051 136051 1 43 43 216802 216802 77 77 70 70 33 33 122 122 50548 50548 1 44 44 222599 222599 65 65 98 98 32 32 124 124 59938 59938 1 45 45 261601 261601 70 70 58 58 32 32 124 124 137639 137639 1 46 46 200657 200657 43 43 88 88 31 31 111 111 138599 138599 1 47 47 259084 259084 67 67 142 142 35 35 129 129 174110 174110 1 48 48 159965 159965 156 156 62 62 32 32 122 122 44244 44244 1 49 49 43287 43287 14 14 13 13 19 19 71 71 43750 43750 1 50 50 172212 172212 67 67 89 89 22 22 81 81 48029 48029 1 51 51 227681 227681 43 43 116 116 36 36 139 139 92288 92288 1 52 52 106288 106288 54 54 28 28 23 23 91 91 197426 197426 1 53 53 268905 268905 58 58 72 72 36 36 133 133 139206 139206 1 54 54 266568 266568 77 77 134 134 42 42 155 155 106271 106271 1 55 55 152474 152474 65 65 106 106 32 32 123 123 71764 71764 1 56 56 330910 330910 95 95 120 120 34 34 128 128 101817 101817 1 57 57 259747 259747 103 103 98 98 35 35 132 132 85008 85008 1 58 58 190495 190495 55 55 66 66 40 40 151 151 124254 124254 1 59 59 154984 154984 73 73 44 44 40 40 151 151 105793 105793 1 60 60 38214 38214 34 34 16 16 8 8 27 27 8773 8773 1 61 61 158671 158671 33 33 56 56 35 35 131 131 94747 94747 1 62 62 299775 299775 68 68 112 112 45 45 170 170 107549 107549 1 63 63 172464 172464 31 31 66 66 35 35 135 135 74783 74783 1 64 64 94381 94381 35 35 42 42 32 32 118 118 66089 66089 1 65 65 243875 243875 274 274 70 70 36 36 140 140 95684 95684 1 66 66 334926 334926 117 117 113 113 37 37 136 136 153824 153824 1 67 67 147979 147979 72 72 55 55 34 34 123 123 63995 63995 1 68 68 216638 216638 44 44 100 100 36 36 134 134 84891 84891 1 69 69 192853 192853 71 71 80 80 36 36 129 129 61263 61263 1 70 70 336678 336678 103 103 95 95 33 33 128 128 113587 113587 1 71 71 271773 271773 89 89 128 128 40 40 154 154 119906 119906 1 72 72 203606 203606 73 73 88 88 33 33 122 122 151611 151611 1 73 73 230177 230177 87 87 132 132 39 39 144 144 144645 144645 1 74 74 1 1 0 0 0 0 0 0 0 0 0 0 1 75 75 14688 14688 10 10 4 4 0 0 0 0 6023 6023 1 76 76 455 455 2 2 0 0 0 0 0 0 0 0 1 77 77 0 0 0 0 0 0 0 0 0 0 0 0 1 78 78 195765 195765 75 75 56 56 33 33 120 120 77457 77457 1 79 79 306514 306514 117 117 111 111 42 42 168 168 62464 62464 1 80 80 203 203 4 4 0 0 0 0 0 0 0 0 1 81 81 7199 7199 5 5 7 7 0 0 0 0 1644 1644 1 82 82 17547 17547 5 5 0 0 1 1 4 4 3926 3926 1 83 83 105044 105044 37 37 37 37 38 38 133 133 42087 42087 1 84 84 969 969 2 2 0 0 0 0 0 0 0 0 2 85 0 0 236496 0 61 0 85 0 34 0 131 0 124252 2 86 0 0 198514 0 62 0 62 0 38 0 146 0 98073 2 87 0 0 137449 0 43 0 55 0 25 0 80 0 41449 2 88 0 0 439387 0 103 0 134 0 29 0 112 0 177551 2 89 0 0 174859 0 50 0 63 0 30 0 116 0 126938 2 90 0 0 186657 0 38 0 77 0 29 0 107 0 61680 2 91 0 0 138866 0 57 0 44 0 33 0 109 0 57793 2 92 0 0 296878 0 65 0 106 0 46 0 159 0 91677 2 93 0 0 192648 0 71 0 63 0 38 0 146 0 64631 2 94 0 0 333348 0 161 0 160 0 52 0 201 0 106385 2 95 0 0 242212 0 57 0 104 0 32 0 124 0 161961 2 96 0 0 263451 0 130 0 86 0 35 0 131 0 112669 2 97 0 0 240028 0 63 0 107 0 40 0 151 0 105416 2 98 0 0 148421 0 43 0 92 0 46 0 184 0 104880 2 99 0 0 176502 0 96 0 123 0 36 0 136 0 76302 2 100 0 0 249735 0 120 0 93 0 38 0 146 0 93071 2 101 0 0 236812 0 76 0 113 0 35 0 130 0 78912 2 102 0 0 142329 0 45 0 52 0 28 0 105 0 35224 2 103 0 0 176054 0 66 0 44 0 42 0 154 0 80849 2 104 0 0 286683 0 79 0 123 0 44 0 169 0 104434 2 105 0 0 247013 0 50 0 77 0 37 0 139 0 63583 2 106 0 0 191653 0 72 0 74 0 32 0 124 0 62486 2 107 0 0 114673 0 31 0 33 0 17 0 55 0 31081 2 108 0 0 284210 0 160 0 105 0 34 0 131 0 94584 2 109 0 0 284195 0 72 0 108 0 33 0 125 0 87408 2 110 0 0 142986 0 48 0 62 0 35 0 130 0 57139 2 111 0 0 140319 0 73 0 50 0 45 0 73 0 90586 2 112 0 0 78800 0 42 0 20 0 26 0 82 0 33032 2 113 0 0 201970 0 69 0 101 0 45 0 173 0 96056 2 114 0 0 194221 0 68 0 89 0 33 0 134 0 87026 2 115 0 0 243889 0 45 0 86 0 33 0 121 0 91072 2 116 0 0 273003 0 74 0 116 0 49 0 186 0 159803 2 117 0 0 333165 0 122 0 138 0 41 0 158 0 144068 2 118 0 0 260981 0 105 0 114 0 25 0 90 0 162627 2 119 0 0 222366 0 76 0 67 0 34 0 131 0 95329 2 120 0 0 201345 0 60 0 88 0 28 0 110 0 62853 2 121 0 0 163043 0 110 0 67 0 31 0 121 0 125976 2 122 0 0 204250 0 128 0 75 0 40 0 151 0 79146 2 123 0 0 127260 0 57 0 119 0 25 0 92 0 99971 2 124 0 0 216092 0 59 0 86 0 42 0 162 0 77826 2 125 0 0 213198 0 67 0 67 0 42 0 163 0 84892 2 126 0 0 191971 0 99 0 75 0 37 0 140 0 238712 2 127 0 0 154651 0 53 0 112 0 34 0 132 0 67486 2 128 0 0 155473 0 56 0 66 0 33 0 122 0 68007 2 129 0 0 145869 0 65 0 30 0 26 0 99 0 38692 2 130 0 0 223666 0 85 0 100 0 40 0 106 0 93587 2 131 0 0 80953 0 25 0 49 0 8 0 28 0 56622 2 132 0 0 130789 0 46 0 26 0 27 0 101 0 15986 2 133 0 0 197389 0 67 0 134 0 34 0 127 0 112215 2 134 0 0 156583 0 56 0 37 0 28 0 102 0 59591 2 135 0 0 178489 0 35 0 78 0 32 0 124 0 143372 2 136 0 0 302789 0 129 0 127 0 52 0 199 0 135062 2 137 0 0 342025 0 98 0 139 0 27 0 102 0 175681 2 138 0 0 246440 0 104 0 108 0 45 0 174 0 130307 2 139 0 0 251306 0 56 0 128 0 37 0 141 0 139141 2 140 0 0 181781 0 117 0 83 0 35 0 131 0 95216 2 141 0 0 260464 0 80 0 157 0 36 0 137 0 94588 2 142 0 0 109632 0 76 0 83 0 36 0 142 0 151244 2 143 0 0 23623 0 11 0 12 0 1 0 0 0 1168 2 144 0 0 61857 0 25 0 23 0 11 0 32 0 25162 2 145 0 0 144889 0 43 0 83 0 40 0 149 0 45635 2 146 0 0 21054 0 16 0 4 0 0 0 0 0 855 2 147 0 0 223718 0 44 0 71 0 27 0 99 0 100174 2 148 0 0 31414 0 19 0 18 0 8 0 25 0 14116 2 149 0 0 112933 0 45 0 29 0 28 0 103 0 117129 2 150 0 0 172783 0 54 0 46 0 43 0 165 0 97392 2 151 0 0 348678 0 69 0 129 0 41 0 159 0 126893 2 152 0 0 266701 0 89 0 139 0 43 0 167 0 118850 2 153 0 0 358933 0 99 0 136 0 47 0 178 0 234853 2 154 0 0 382487 0 153 0 97 0 42 0 158 0 139537 2 155 0 0 111853 0 39 0 49 0 35 0 132 0 144253 2 156 0 0 173710 0 103 0 29 0 27 0 107 0 106221 2 157 0 0 212961 0 75 0 114 0 35 0 129 0 113864 2 158 0 0 173260 0 63 0 41 0 21 0 79 0 37238 2 159 0 0 127096 0 51 0 142 0 47 0 180 0 135096 2 160 0 0 98 0 1 0 0 0 0 0 0 0 0 2 161 0 0 0 0 0 0 0 0 0 0 0 0 0 2 162 0 0 0 0 0 0 0 0 0 0 0 0 0 2 163 0 0 46660 0 20 0 12 0 5 0 15 0 6179 2 164 0 0 165838 0 56 0 46 0 28 0 101 0 87656
Names of X columns:
Pop t pop_t Total_Time_spent_in_RFC Total_Time_spent_in_RFC_p Number_of_Logins Number_of_Logins_p Total_Number_of_Blogged_Computations Total_Number_of_Blogged_Computations_p Total_Number_of_Reviewed_Compendiums Total_Number_of_Reviewed_Compendiums_p Total_Number_of_submitted_Feedback_Messages_in_Peer_Reviews Total_Number_of_submitted_Feedback_Messages_in_Peer_Reviews_p Compendium_Writing:_total_number_of_characters Compendium_Writing:_total_number_of_characters_p
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