R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- c(-719.405774320179,-1892.13972512203,948.584751735815,694.702053597555,-1016.98885981099,-985.243198088241,-439.006139106049,-8289.64839707302,478.341587615535,-4619.19009697141,-438.889269676146,1922.27169831854,-9316.27276546704,4934.24926529499,294.506079291393,6062.59748359645,1545.14497838096,-5388.19323530356,-8728.74410350888,-10701.1030591330,10938.7367915057,-12697.2927303432,-1368.16155419152,545.462343565687,8646.56018969786,-5152.59361283694,-5230.8539005193,-1548.47014455579,-337.152015819184,1309.74791688059,2666.25355012676,-5370.63361781142,-7476.69691278903,8086.32165260496,6571.98868771161,3149.48850681552,104.637480526355,-1952.46924181325,3591.38504257275,962.019554951216,-5699.66826957231,5791.29992251809,-3685.82799741182,-4065.74959426738,6283.84760394879,6963.92536654231,4239.8594212716,5856.07066535003,7938.63020854731,4533.55760400507,1779.15891762194,-3469.0240935556,1497.02166817019,4510.32580191698,-4112.36752504752,254.655960321710,2559.86399726922) > par1 = '12' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > par1 <- as.numeric(par1) > (n <- length(x)) [1] 57 > (np <- floor(n / par1)) [1] 4 > arr <- array(NA,dim=c(par1,np+1)) > darr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > dx <- diff(x) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + darr[j,ari[j]] <- dx[i] + if (j == par1) j = 0 + } > ari [1] 5 5 5 5 5 5 5 5 5 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] -719.4058 -9316.2728 8646.560 104.6375 7938.6302 [2,] -1892.1397 4934.2493 -5152.594 -1952.4692 4533.5576 [3,] 948.5848 294.5061 -5230.854 3591.3850 1779.1589 [4,] 694.7021 6062.5975 -1548.470 962.0196 -3469.0241 [5,] -1016.9889 1545.1450 -337.152 -5699.6683 1497.0217 [6,] -985.2432 -5388.1932 1309.748 5791.2999 4510.3258 [7,] -439.0061 -8728.7441 2666.254 -3685.8280 -4112.3675 [8,] -8289.6484 -10701.1031 -5370.634 -4065.7496 254.6560 [9,] 478.3416 10938.7368 -7476.697 6283.8476 2559.8640 [10,] -4619.1901 -12697.2927 8086.322 6963.9254 NA [11,] -438.8893 -1368.1616 6571.989 4239.8594 NA [12,] 1922.2717 545.4623 3149.489 5856.0707 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] -1172.73395 14250.522 -13799.15380 -2057.1067 -3405.073 [2,] 2840.72448 -4639.743 -78.26029 5543.8543 -2754.399 [3,] -253.88270 5768.091 3682.38376 -2629.3655 -5248.183 [4,] -1711.69091 -4517.453 1211.31813 -6661.6878 4966.046 [5,] 31.74566 -6933.338 1646.89993 11490.9682 3013.304 [6,] 546.23706 -3340.551 1356.50563 -9477.1279 -8622.693 [7,] -7850.64226 -1972.359 -8036.88717 -379.9216 4367.023 [8,] 8767.98998 21639.840 -2106.06329 10349.5972 2305.208 [9,] -5097.53168 -23636.030 15563.01857 680.0778 NA [10,] 4180.30083 11329.131 -1514.33296 -2724.0659 NA [11,] 2361.16097 1913.624 -3422.50018 1616.2112 NA [12,] -11238.54446 8101.098 -3044.85103 2082.5595 NA > arr.mean <- array(NA,dim=par1) > arr.median <- array(NA,dim=par1) > arr.midrange <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.mean[j] <- mean(arr[j,],na.rm=TRUE) + arr.median[j] <- median(arr[j,],na.rm=TRUE) + arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2 + } > overall.mean <- mean(x) > overall.median <- median(x) > overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2 > postscript(file="/var/www/html/rcomp/tmp/12cjk1259789365.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.mean,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2hmlv1259789365.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.median,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3zgqc1259789365.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.midrange,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4x03x1259789365.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(arr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -9316.2728 -5152.594 294.5061 -3469.0241 -1016.989 -5388.1932 -8728.7441 [2,] -719.4058 -1952.469 294.5061 -1548.4701 -1016.989 -985.2432 -4112.3675 [3,] 104.6375 -1892.140 948.5848 694.7021 -337.152 1309.7479 -3685.8280 [4,] 7938.6302 4533.558 1779.1589 962.0196 1497.022 4510.3258 -439.0061 [5,] 8646.5602 4934.249 3591.3850 962.0196 1545.145 5791.2999 2666.2536 [,8] [,9] [,10] [,11] [,12] [1,] -10701.1031 -7476.6969 -12697.293 -1368.1616 545.4623 [2,] -8289.6484 478.3416 -8658.241 -903.5254 1233.8670 [3,] -5370.6336 2559.8640 1172.368 1900.4851 2535.8801 [4,] -4065.7496 6283.8476 7525.124 5405.9241 4502.7796 [5,] 254.6560 10938.7368 8086.322 6571.9887 5856.0707 $n [1] 5 5 5 5 5 5 5 5 5 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -6013.109 -6475.150 -100.4672 -1079.204 -2113.546 -2573.407 -6281.416 [2,] 6222.384 2690.871 1997.6367 2468.608 1439.241 5192.903 -1090.240 [,8] [,9] [,10] [,11] [,12] [1,] -8355.230 -1542.292 -11612.49 -3083.98 -46.56082 [2,] -2386.037 6662.020 13957.23 6884.95 5118.32103 $out [1] -5230.854 6062.597 -5699.668 $group [1] 3 4 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-9316.27276546704, -719.405774320179, 104.637480526355, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5zz861259789365.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(darr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [1,] -3405.073 -4639.74319 -5248.1830 -6661.688 31.74566 -9477.1279 [2,] -3405.073 -2754.39869 -2629.3655 -4517.453 31.74566 -8622.6933 [3,] -2057.107 -78.26029 -253.8827 -1711.691 1646.89993 -3340.5509 [4,] -1172.734 2840.72448 3682.3838 1211.318 3013.30413 546.2371 [5,] -1172.734 5543.85428 5768.0914 4966.046 3013.30413 1356.5056 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -8036.8872 -2106.063 -23636.030 -2724.066 -3422.5002 -11238.5445 [2,] -7850.6423 2305.208 -14366.781 -2119.199 -903.1445 -7141.6977 [3,] -1972.3590 8767.990 -2208.727 1332.984 1764.9176 -481.1457 [4,] -379.9216 10349.597 8121.548 7754.716 2137.3924 5091.8287 [5,] 4367.0235 21639.840 15563.019 11329.131 2361.1610 8101.0978 $n [1] 5 5 5 5 5 5 5 5 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3634.4716 -4031.76 -4713.749 -5759.626 -459.8618 -9819.294 -7251.151 [2,] -479.7419 3875.24 4205.983 2336.244 3753.6616 3138.192 3306.433 [,8] [,9] [,10] [,11] [,12] [1,] 3083.845 -19974.51 -6467.409 -637.1066 -10145.63 [2,] 14452.135 15557.05 9133.377 4166.9417 9183.34 $out [1] 14250.522 -13799.154 -6933.338 11490.968 $group [1] 1 1 5 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-3405.07260454224, -3405.07260454224, -2057.10672233960, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6md5l1259789365.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(arr) > names(z) <- c(1:np) > (boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks')) $stats [,1] [,2] [,3] [,4] [,5] [1,] -1892.1397 -12697.2927 -7476.6969 -5699.668 -4112.3675 [2,] -1454.5643 -9022.5084 -5191.7238 -2819.149 254.6560 [3,] -579.2060 -536.8277 486.2980 2276.702 1779.1589 [4,] 586.5218 3239.6971 4860.7386 5823.685 4510.3258 [5,] 1922.2717 10938.7368 8646.5602 6963.925 7938.6302 $n [1] 12 12 12 12 9 $conf [,1] [,2] [,3] [,4] [,5] [1,] -1510.1590 -6129.702 -4098.698 -1665.352 -462.1605 [2,] 351.7471 5056.047 5071.294 6218.757 4020.4784 $out [1] -8289.648 -4619.190 $group [1] 1 1 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-1892.13972512203, -1454.56429246651, -579.205956713114, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/77tdk1259789365.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.mean,arr.median,arr.midrange)) > names(z) <- list('mean','median','midrange') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] -2859.9384 -3685.8280 -2275.68683 [2,] -684.4437 -1114.6459 -26.60445 [3,] 408.4606 821.6434 1163.68834 [4,] 1791.0146 1605.1165 2388.97195 [5,] 2868.3233 2559.8640 3609.61222 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -720.6126 -418.8585 61.92773 [2,] 1537.5337 2062.1453 2265.44895 $out [1] -5634.496 -5370.634 -6177.699 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-2859.9384429895, -684.443725886051, 408.460574473643, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/12cjk1259789365.ps tmp/12cjk1259789365.png") > system("convert tmp/2hmlv1259789365.ps tmp/2hmlv1259789365.png") > system("convert tmp/3zgqc1259789365.ps tmp/3zgqc1259789365.png") > system("convert tmp/4x03x1259789365.ps tmp/4x03x1259789365.png") > system("convert tmp/5zz861259789365.ps tmp/5zz861259789365.png") > system("convert tmp/6md5l1259789365.ps tmp/6md5l1259789365.png") > system("convert tmp/77tdk1259789365.ps tmp/77tdk1259789365.png") > > > proc.time() user system elapsed 1.203 0.952 1.527