R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-redhat-linux-gnu (64-bit) 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. Natural language support but running in an English locale 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(5939520.00,89948768.00,80953652.00,85942882.00,8944937.00,82975432.00,24940816.00,21973899.00,37950221.00,45949881.00,85950373.00,48960313.00,81954506.00,24960419.00,65973338.00,22950513.00,54963528.00,90995659.00,91967517.00,28999053.00,96990529.00,38979852.00,81496957.00,74982424.00,70976192.00,90990000.00,12998850.00,92986156.00,67994976.00,91022206.00,87992489.00,421022698.00,11018942.00,79100042.00,65996442.00,51000620.00,12996871.00,44994249.00,99996135.00,91977037.00,63974211.00,15998036.00,65974265.00,33984410.00,45939098.00,67935827.00,66921032.00,89911836.00,71890975.00,72880342.00,28871286.00,41844334.00,82847667.00,24871401.00,3867451.00,99896846.00,41890361.00,45884264.00,69884586.00,95896400.00,39904491.00,81900399.00,27909863.00,88900470.00,89917101.00,2945005.00,4934411.00,61957264.00,31946515.00,3938309.00,52933321.00,21947613.00) > 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] 72 > (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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 5939520 81954506 70976192 12996871 71890975 39904491 NA [2,] 89948768 24960419 90990000 44994249 72880342 81900399 NA [3,] 80953652 65973338 12998850 99996135 28871286 27909863 NA [4,] 85942882 22950513 92986156 91977037 41844334 88900470 NA [5,] 8944937 54963528 67994976 63974211 82847667 89917101 NA [6,] 82975432 90995659 91022206 15998036 24871401 2945005 NA [7,] 24940816 91967517 87992489 65974265 3867451 4934411 NA [8,] 21973899 28999053 421022698 33984410 99896846 61957264 NA [9,] 37950221 96990529 11018942 45939098 41890361 31946515 NA [10,] 45949881 38979852 79100042 67935827 45884264 3938309 NA [11,] 85950373 81496957 65996442 66921032 69884586 52933321 NA [12,] 48960313 74982424 51000620 89911836 95896400 21947613 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 84009248 -56994087 20013808 31997378 989367 41995908 NA [2,] -8995116 41012919 -77991150 55001886 -44009056 -53990536 NA [3,] 4989230 -43022825 79987306 -8019098 12973048 60990607 NA [4,] -76997945 32013015 -24991180 -28002826 41003333 1016631 NA [5,] 74030495 36032131 23027230 -47976175 -57976266 -86972096 NA [6,] -58034616 971858 -3029717 49976229 -21003950 1989406 NA [7,] -2966917 -62968464 333030209 -31989855 96029395 57022853 NA [8,] 15976322 67991476 -410003756 11954688 -58006485 -30010749 NA [9,] 7999660 -58010677 68081100 21996729 3993903 -28008206 NA [10,] 40000492 42517105 -13103600 -1014795 24000322 48995012 NA [11,] -36990060 -6514533 -14995822 22990804 26011814 -30985708 NA [12,] 32994193 -4006232 -38003749 -18020861 -55991909 NA 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/wessaorg/rcomp/tmp/106791305713529.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/266ws1305713529.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/34eox1305713529.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4sxk71305713529.ps",horizontal=F,onefile=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] [1,] 5939520 24960419 12998850 22950513 54963528 2945005 3867451 21973899 [2,] 12996871 44994249 27909863 41844334 54963528 15998036 4934411 28999053 [3,] 55440342 77390370 47422312 87421676 65984594 53923416 45457540 47970837 [4,] 71890975 89948768 80953652 91977037 82847667 90995659 87992489 99896846 [5,] 81954506 90990000 99996135 92986156 89917101 91022206 91967517 99896846 [,9] [,10] [,11] [,12] [1,] 11018942 3938309 52933321 21947613 [2,] 31946515 38979852 65996442 48960313 [3,] 39920291 45917072 68402809 62991522 [4,] 45939098 67935827 81496957 89911836 [5,] 45939098 79100042 85950373 95896400 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 17451742 48393254 13207355 55084463 47998423 5547527 -8117602 2239471 [2,] 93428941 106387487 81637269 119758889 83970764 102299306 99032683 93702203 [,9] [,10] [,11] [,12] [1,] 30894623 27239533 58404476 36576467 [2,] 48945959 64594612 78401142 89406577 $out [1] 8944937 421022698 96990529 $group [1] 5 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(5939520, 12996871, 55440341.5, 71890975, 81954506, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/59xla1305713529.ps",horizontal=F,onefile=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] [1,] -56994087 -77991150 -43022825 -76997945 -86972096 -21003950 -62968464 [2,] 989367 -53990536 -8019098 -28002826 -57976266 -21003950 -31989855 [3,] 26005593 -26502086 8981139 -11987274 -12474472 -1028930 27027968 [4,] 41995908 41012919 60990607 32013015 36032131 1989406 96029395 [5,] 84009248 55001886 79987306 41003333 74030495 1989406 96029395 [,8] [,9] [,10] [,11] [,12] [1,] -58006485 -58010677 -13103600 -36990060 -55991909 [2,] -58006485 -28008206 -1014795 -30985708 -38003749 [3,] -9028030 5996782 32000407 -10755178 -18020861 [4,] 15976322 21996729 42517105 22990804 -4006232 [5,] 67991476 68081100 48995012 26011814 32994193 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -444950.4 -87782383 -35532350 -50699430 -73112925 -15860387 -55548583 [2,] 52456136.4 34778211 53494628 26724881 48163980 13802528 109604519 [,8] [,9] [,10] [,11] [,12] [1,] -56749331 -26258017 3920925 -45571771 -42043421 [2,] 38693270 38251580 60079889 24061416 6001699 $out [1] -58034616 49976229 333030209 -410003756 $group [1] 6 6 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-56994087, 989367, 26005593, 41995908, 84009248, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/69zo51305713529.ps",horizontal=F,onefile=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,] 5939520 22950513 11018942 12996871 3867451 2945005 NA [2,] 23457358 33989452 58498531 39489330 35357810 13441012 NA [3,] 47455097 70477881 75038117 64974238 57884425 35925503 NA [4,] 84459157 86475082 91006103 78923832 77864004 71928832 NA [5,] 89948768 96990529 92986156 99996135 99896846 89917101 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 19631767 46538833 60211195 46987897 38497064 9248816 NA [2,] 75278427 94416929 89865039 82960579 77271786 62602190 NA $out [1] 421022698 $group [1] 3 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(5939520, 23457357.5, 47455097, 84459157, 89948768, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/7i8hd1305713529.ps",horizontal=F,onefile=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,] 44289278 39920291 39187178 [2,] 47120894 46669692 48892310 [3,] 57112129 54681879 56966081 [4,] 69071407 67193701 69063285 [5,] 70766899 87421676 72410727 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 47100352 45320741 47765963 [2,] 67123905 64043017 66166198 $out [1] 111305695 $group [1] 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(44289277.6666667, 47120894.1666667, 57112128.6666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/106791305713529.ps tmp/106791305713529.png",intern=TRUE)) character(0) > try(system("convert tmp/266ws1305713529.ps tmp/266ws1305713529.png",intern=TRUE)) character(0) > try(system("convert tmp/34eox1305713529.ps tmp/34eox1305713529.png",intern=TRUE)) character(0) > try(system("convert tmp/4sxk71305713529.ps tmp/4sxk71305713529.png",intern=TRUE)) character(0) > try(system("convert tmp/59xla1305713529.ps tmp/59xla1305713529.png",intern=TRUE)) character(0) > try(system("convert tmp/69zo51305713529.ps tmp/69zo51305713529.png",intern=TRUE)) character(0) > try(system("convert tmp/7i8hd1305713529.ps tmp/7i8hd1305713529.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.58 0.47 5.46