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(2550,2867,3458,2961,3163,2880,3331,3062,3534,3622,4464,5411,2564,2820,3508,3088,3299,2939,3320,3418,3604,3495,4163,4882,2211,3260,2992,2425,2707,3244,3965,3315,3333,3583,4021,4904,2252,2952,3573,3048,3059,2731,3563,3092,3478,3478,4308,5029,2075,3264,3308,3688,3136,2824,3644,4694,2914,3686,4358,5587,2265,3685,3754,3708,3210,3517,3905,3670,4221,4404,5086,5725,2367,3819,4067,4022,3937,4365,4290) > 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] 79 > (np <- floor(n / par1)) [1] 6 > 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] 7 7 7 7 7 7 7 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 2550 2564 2211 2252 2075 2265 2367 [2,] 2867 2820 3260 2952 3264 3685 3819 [3,] 3458 3508 2992 3573 3308 3754 4067 [4,] 2961 3088 2425 3048 3688 3708 4022 [5,] 3163 3299 2707 3059 3136 3210 3937 [6,] 2880 2939 3244 2731 2824 3517 4365 [7,] 3331 3320 3965 3563 3644 3905 4290 [8,] 3062 3418 3315 3092 4694 3670 NA [9,] 3534 3604 3333 3478 2914 4221 NA [10,] 3622 3495 3583 3478 3686 4404 NA [11,] 4464 4163 4021 4308 4358 5086 NA [12,] 5411 4882 4904 5029 5587 5725 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 317 256 1049 700 1189 1420 1452 [2,] 591 688 -268 621 44 69 248 [3,] -497 -420 -567 -525 380 -46 -45 [4,] 202 211 282 11 -552 -498 -85 [5,] -283 -360 537 -328 -312 307 428 [6,] 451 381 721 832 820 388 -75 [7,] -269 98 -650 -471 1050 -235 NA [8,] 472 186 18 386 -1780 551 NA [9,] 88 -109 250 0 772 183 NA [10,] 842 668 438 830 672 682 NA [11,] 947 719 883 721 1229 639 NA [12,] -2847 -2671 -2652 -2954 -3322 -3358 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/1zk3l1272299400.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/2zk3l1272299400.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/39b2o1272299400.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/49b2o1272299400.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] [,8] [,9] [,10] [,11] [1,] 2075.0 2820.0 2992.0 2425.0 3059.0 2731.0 3320 3062.0 3333 3478.0 4021 [2,] 2231.5 2909.5 3383.0 3004.5 3097.5 2852.0 3447 3092.0 3333 3495.0 4163 [3,] 2265.0 3260.0 3508.0 3088.0 3163.0 2939.0 3644 3366.5 3506 3602.5 4333 [4,] 2458.5 3474.5 3663.5 3698.0 3254.5 3380.5 3935 3670.0 3604 3686.0 4464 [5,] 2564.0 3819.0 4067.0 4022.0 3299.0 3517.0 4290 3670.0 3604 3686.0 4464 [,12] [1,] 4882 [2,] 4904 [3,] 5220 [4,] 5587 [5,] 5725 $n [1] 7 7 7 7 7 7 7 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2129.439 2922.591 3340.49 2673.853 3069.242 2623.388 3352.574 2993.671 [2,] 2400.561 3597.409 3675.51 3502.147 3256.758 3254.612 3935.426 3739.329 [,9] [,10] [,11] [,12] [1,] 3331.196 3479.299 4138.845 4779.443 [2,] 3680.804 3725.701 4527.155 5660.557 $out [1] 2707 3937 4365 4694 2914 4221 4404 5086 $group [1] 5 5 6 8 9 9 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(2075, 2231.5, 2265, 2458.5, 2564, 2820, 2909.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/59b2o1272299400.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] [,7] [,8] [,9] [,10] [,11] [1,] 256.0 -268.0 -567.0 -552.0 -360.0 -75.0 -650 18 -109.0 438 639 [2,] 508.5 56.5 -511.0 -291.5 -320.0 384.5 -471 18 0.0 668 719 [3,] 1049.0 248.0 -420.0 11.0 -283.0 451.0 -252 286 135.5 677 802 [4,] 1304.5 606.0 -45.5 206.5 367.5 770.5 98 472 250.0 830 947 [5,] 1452.0 688.0 380.0 282.0 537.0 832.0 98 551 250.0 842 1229 [,12] [1,] -3358.0 [2,] -3322.0 [3,] -2900.5 [4,] -2671.0 [5,] -2652.0 $n [1] 7 7 7 7 7 7 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 573.6416 -80.15254 -697.9891 -286.3976 -693.5639 220.4870 -619.0234 [2,] 1524.3584 576.15254 -142.0109 308.3976 127.5639 681.5130 115.0234 [,8] [,9] [,10] [,11] [,12] [1,] -6.844664 -25.75807 572.5048 654.9326 -3320.416 [2,] 578.844664 296.75807 781.4952 949.0674 -2480.584 $out [1] 1050 -1780 772 $group [1] 7 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(256, 508.5, 1049, 1304.5, 1452, -268, 56.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/62kj91272299400.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] [,6] [,7] [1,] 2550.0 2564.0 2211.0 2252 2075 3210.0 3819.0 [2,] 2920.5 3013.5 2849.5 3000 3025 3593.5 3878.0 [3,] 3247.0 3369.0 3287.5 3285 3476 3731.0 4022.0 [4,] 3578.0 3556.0 3774.0 3568 4023 4312.5 4178.5 [5,] 4464.0 4163.0 4904.0 4308 4694 5086.0 4365.0 $n [1] 12 12 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 2947.11 3121.562 2865.829 3025.931 3020.806 3403.059 3842.546 [2,] 3546.89 3616.438 3709.171 3544.069 3931.194 4058.941 4201.454 $out [1] 5411 4882 5029 5587 2265 5725 2367 $group [1] 1 2 4 5 6 6 7 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(2550, 2920.5, 3247, 3578, 4464, 2564, 3013.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/72kj91272299400.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,] 3214.286 2939.00 3116.250 [2,] 3227.000 3125.50 3184.000 [3,] 3518.429 3436.25 3427.625 [4,] 3714.095 3623.25 3642.250 [5,] 4400.000 4333.00 4318.375 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 3296.261 3209.223 3218.614 [2,] 3740.596 3663.277 3636.636 $out [1] 2326.286 5256.333 2265.000 5220.000 2345.000 5239.125 $group [1] 1 1 2 2 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(3214.28571428571, 3227, 3518.42857142857, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1zk3l1272299400.ps tmp/1zk3l1272299400.png",intern=TRUE)) character(0) > try(system("convert tmp/2zk3l1272299400.ps tmp/2zk3l1272299400.png",intern=TRUE)) character(0) > try(system("convert tmp/39b2o1272299400.ps tmp/39b2o1272299400.png",intern=TRUE)) character(0) > try(system("convert tmp/49b2o1272299400.ps tmp/49b2o1272299400.png",intern=TRUE)) character(0) > try(system("convert tmp/59b2o1272299400.ps tmp/59b2o1272299400.png",intern=TRUE)) character(0) > try(system("convert tmp/62kj91272299400.ps tmp/62kj91272299400.png",intern=TRUE)) character(0) > try(system("convert tmp/72kj91272299400.ps tmp/72kj91272299400.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.328 1.082 2.191