*The author of this computation has been verified*

R Software Module: /rwasp_multipleregression.wasp (opens new window with default values)

Title produced by software: Multiple Regression

Date of computation: Sat, 11 Dec 2010 13:42:31 +0000

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8.htm/, Retrieved Sat, 11 Dec 2010 14:40:29 +0100

BibTeX entries for LaTeX users:

@Manual{KEY, author = {{YOUR NAME}}, publisher = {Office for Research Development and Education}, title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8.htm/}, year = {2010}, } @Manual{R, title = {R: A Language and Environment for Statistical Computing}, author = {{R Development Core Team}}, organization = {R Foundation for Statistical Computing}, address = {Vienna, Austria}, year = {2010}, note = {{ISBN} 3-900051-07-0}, url = {http://www.R-project.org}, }

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

Dataseries X:

» Textbox « » Textfile « » CSV «

9700 0 9081 0 9084 0 9743 0 8587 0 9731 0 9563 0 9998 0 9437 0 10038 0 9918 0 9252 0 9737 0 9035 0 9133 0 9487 0 8700 0 9627 0 8947 0 9283 0 8829 0 9947 0 9628 0 9318 0 9605 0 8640 0 9214 0 9567 0 8547 0 9185 0 9470 0 9123 0 9278 0 10170 0 9434 0 9655 0 9429 0 8739 0 9552 0 9687 1 9019 1 9672 1 9206 1 9069 1 9788 1 10312 1 10105 1 9863 1 9656 1 9295 1 9946 1 9701 1 9049 1 10190 1 9706 1 9765 1 9893 1 9994 1 10433 1 10073 1 10112 1 9266 1 9820 1 10097 1 9115 1 10411 1 9678 1 10408 1 10153 1 10368 1 10581 1 10597 1 10680 1 9738 1 9556 1

Output produced by software:

Enter (or paste) a matrix (table) containing all data (time) series. Every column represents a different variable and must be delimited by a space or Tab. Every row represents a period in time (or category) and must be delimited by hard returns. The easiest way to enter data is to copy and paste a block of spreadsheet cells. Please, do not use commas or spaces to seperate groups of digits! Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 6 seconds R Server 'RServer@AstonUniversity' @ vre.aston.ac.uk Multiple Linear Regression - Estimated Regression Equation geboortes[t] = + 9430.3111111111 + 286.706944444442x[t] + 99.1618551587258M1[t] -638.203273809524M2[t] -284.711259920635M3[t] -37.555158730158M4[t] -920.277430555555M5[t] + 41.0002976190481M6[t] -338.555307539682M7[t] -164.444246031746M8[t] -214.333184523809M9[t] + 355.611210317461M10[t] + 228.722271825397M11[t] + 5.2222718253969t + e[t] Multiple Linear Regression - Ordinary Least Squares Variable Parameter S.D. T-STAT H0: parameter = 0 2-tail p-value 1-tail p-value (Intercept) 9430.3111111111 140.713171 67.018 0 0 x 286.706944444442 134.583356 2.1303 0.037185 0.018593 M1 99.1618551587258 158.573874 0.6253 0.534083 0.267042 M2 -638.203273809524 158.462983 -4.0275 0.000159 7.9e-05 M3 -284.711259920635 158.413321 -1.7973 0.077243 0.038622 M4 -37.555158730158 166.197782 -0.226 0.821983 0.410991 M5 -920.277430555555 165.759024 -5.5519 1e-06 0 M6 41.0002976190481 165.377825 0.2479 0.80503 0.402515 M7 -338.555307539682 165.054585 -2.0512 0.044552 0.022276 M8 -164.444246031746 164.789644 -0.9979 0.322268 0.161134 M9 -214.333184523809 164.583284 -1.3023 0.197717 0.098859 M10 355.611210317461 164.435725 2.1626 0.034501 0.01725 M11 228.722271825397 164.347126 1.3917 0.169067 0.084533 t 5.2222718253969 3.116074 1.6759 0.098874 0.049437 Multiple Linear Regression - Regression Statistics Multiple R 0.857999753809147 R-squared 0.736163577536557 Adjusted R-squared 0.679936143241069 F-TEST (value) 13.0926048246813 F-TEST (DF numerator) 13 F-TEST (DF denominator) 61 p-value 3.89466237038505e-13 Multiple Linear Regression - Residual Statistics Residual Standard Deviation 284.606401731326 Sum Squared Residuals 4941049.03829364 Multiple Linear Regression - Actuals, Interpolation, and Residuals Time or Index Actuals Interpolation Forecast Residuals Prediction Error 1 9700 9534.69523809526 165.304761904736 2 9081 8802.55238095238 278.44761904762 3 9084 9161.26666666666 -77.2666666666647 4 9743 9413.64503968254 329.354960317462 5 8587 8536.14503968254 50.8549603174615 6 9731 9502.64503968254 228.354960317461 7 9563 9128.3117063492 434.688293650795 8 9998 9307.64503968254 690.354960317462 9 9437 9262.97837301587 174.021626984128 10 10038 9838.14503968254 199.854960317462 11 9918 9716.47837301587 201.521626984128 12 9252 9492.97837301587 -240.978373015872 13 9737 9597.3625 139.637500000006 14 9035 8865.21964285714 169.780357142858 15 9133 9223.93392857143 -90.9339285714277 16 9487 9476.3123015873 10.6876984126986 17 8700 8598.8123015873 101.187698412699 18 9627 9565.3123015873 61.6876984126987 19 8947 9190.97896825397 -243.978968253968 20 9283 9370.3123015873 -87.3123015873014 21 8829 9325.64563492063 -496.645634920634 22 9947 9900.8123015873 46.1876984126987 23 9628 9779.14563492063 -151.145634920635 24 9318 9555.64563492063 -237.645634920635 25 9605 9660.02976190476 -55.0297619047573 26 8640 8927.8869047619 -287.886904761905 27 9214 9286.60119047619 -72.6011904761906 28 9567 9538.97956349206 28.0204365079358 29 8547 8661.47956349206 -114.479563492064 30 9185 9627.97956349206 -442.979563492064 31 9470 9253.64623015873 216.353769841269 32 9123 9432.97956349206 -309.979563492064 33 9278 9388.3128968254 -110.312896825398 34 10170 9963.47956349206 206.520436507936 35 9434 9841.8128968254 -407.812896825397 36 9655 9618.3128968254 36.6871031746026 37 9429 9722.69702380952 -293.69702380952 38 8739 8990.55416666667 -251.554166666667 39 9552 9349.26845238095 202.731547619047 40 9687 9888.35376984127 -201.353769841269 41 9019 9010.85376984127 8.14623015873069 42 9672 9977.35376984127 -305.353769841269 43 9206 9603.02043650794 -397.020436507936 44 9069 9782.35376984127 -713.35376984127 45 9788 9737.6871031746 50.3128968253973 46 10312 10312.8537698413 -0.853769841269287 47 10105 10191.1871031746 -86.1871031746026 48 9863 9967.6871031746 -104.687103174603 49 9656 10072.0712301587 -416.071230158725 50 9295 9339.92837301587 -44.9283730158728 51 9946 9698.64265873016 247.357341269842 52 9701 9951.02103174603 -250.021031746032 53 9049 9073.52103174603 -24.5210317460321 54 10190 10040.021031746 149.978968253968 55 9706 9665.6876984127 40.3123015873013 56 9765 9845.02103174603 -80.021031746032 57 9893 9800.35436507937 92.6456349206346 58 9994 10375.521031746 -381.521031746032 59 10433 10253.8543650794 179.145634920635 60 10073 10030.3543650794 42.6456349206347 61 10112 10134.7384920635 -22.7384920634879 62 9266 9402.59563492064 -136.595634920636 63 9820 9761.30992063492 58.6900793650788 64 10097 10013.6882936508 83.3117063492052 65 9115 9136.1882936508 -21.1882936507949 66 10411 10102.6882936508 308.311706349205 67 9678 9728.35496031746 -50.3549603174615 68 10408 9907.6882936508 500.311706349205 69 10153 9863.02162698413 289.978373015872 70 10368 10438.1882936508 -70.1882936507949 71 10581 10316.5216269841 264.478373015872 72 10597 10093.0216269841 503.978373015872 73 10680 10197.4057539683 482.594246031749 74 9738 9465.2628968254 272.737103174602 75 9556 9823.97718253968 -267.977182539684 Goldfeld-Quandt test for Heteroskedasticity p-values Alternative Hypothesis breakpoint index greater 2-sided less 17 0.128704094160692 0.257408188321383 0.871295905839308 18 0.0585261049506291 0.117052209901258 0.941473895049371 19 0.332367614469383 0.664735228938767 0.667632385530617 20 0.550884107448718 0.898231785102564 0.449115892551282 21 0.586756671353637 0.826486657292726 0.413243328646363 22 0.513449018504325 0.97310196299135 0.486550981495675 23 0.407871981551796 0.815743963103592 0.592128018448204 24 0.370593804769147 0.741187609538295 0.629406195230853 25 0.319036561008134 0.638073122016268 0.680963438991866 26 0.249508717541665 0.49901743508333 0.750491282458335 27 0.280486322768488 0.560972645536976 0.719513677231512 28 0.242372334677351 0.484744669354702 0.757627665322649 29 0.182497335099524 0.364994670199048 0.817502664900476 30 0.197747143157017 0.395494286314035 0.802252856842983 31 0.299959030274241 0.599918060548483 0.700040969725759 32 0.28770785260175 0.575415705203501 0.71229214739825 33 0.288768577780508 0.577537155561016 0.711231422219492 34 0.382424533003739 0.764849066007479 0.617575466996261 35 0.374207099953584 0.748414199907168 0.625792900046416 36 0.442915661496232 0.885831322992463 0.557084338503768 37 0.384864629026904 0.769729258053809 0.615135370973096 38 0.359014037214078 0.718028074428156 0.640985962785922 39 0.44630739378285 0.892614787565699 0.55369260621715 40 0.375181631546515 0.75036326309303 0.624818368453485 41 0.364361540384443 0.728723080768885 0.635638459615557 42 0.316459530508155 0.63291906101631 0.683540469491845 43 0.282369215667823 0.564738431335647 0.717630784332177 44 0.545537159609201 0.908925680781599 0.454462840390799 45 0.562852318887975 0.87429536222405 0.437147681112025 46 0.627185925495234 0.745628149009532 0.372814074504766 47 0.564436560349906 0.871126879300189 0.435563439650094 48 0.503167450735575 0.99366509852885 0.496832549264425 49 0.559425109492357 0.881149781015287 0.440574890507643 50 0.483790186403116 0.967580372806233 0.516209813596883 51 0.783017429210443 0.433965141579113 0.216982570789557 52 0.705302666891455 0.58939466621709 0.294697333108545 53 0.648188520995924 0.703622958008151 0.351811479004076 54 0.584657412834872 0.830685174330256 0.415342587165128 55 0.58036611538795 0.8392677692241 0.41963388461205 56 0.558480063448771 0.883039873102458 0.441519936551229 57 0.429348480632675 0.858696961265349 0.570651519367325 58 0.293459673420112 0.586919346840224 0.706540326579888 Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity Description # significant tests % significant tests OK/NOK 1% type I error level 0 0 OK 5% type I error level 0 0 OK 10% type I error level 0 0 OK

Charts produced by software:

http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/10effz1292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/10effz1292074941.ps (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/1pw051292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/1pw051292074941.ps (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/205iq1292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/205iq1292074941.ps (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/305iq1292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/305iq1292074941.ps (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/405iq1292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/405iq1292074941.ps (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/5sfzt1292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/5sfzt1292074941.ps (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/6sfzt1292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/6sfzt1292074941.ps (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/736ye1292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/736ye1292074941.ps (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/8effz1292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/8effz1292074941.ps (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/9effz1292074941.png (open in new window) http://www.freestatistics.org/blog/date/2010/Dec/11/t1292074826qptza6k4afxx5f8/9effz1292074941.ps (open in new window)

Parameters (Session):

par1 = 1 ; par2 = Include Monthly Dummies ; par3 = Linear Trend ;

Parameters (R input):

par1 = 1 ; par2 = Include Monthly Dummies ; par3 = Linear Trend ;

R code (references can be found in the software module):

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') }

Copyright

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa

Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.

Privacy Policy

We may request personal information to be submitted to our servers in order to be able to:

• personalize online software applications according to your needs
• enforce strict security rules with respect to the data that you upload (e.g. statistical data)
• manage user sessions of online applications
• alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.

FreeStatistics.org is powered by