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(10788500,10708433,10625700,10547958,10478617,10421137,10376133,10332785,10286570,10251250,10226419,10203008,10181245,10156637,10136811,10115603,10084475,10045158,10004486,9967379,9937697,9901664,9870234,9861823,9858308,9855372,9855520,9856303,9858982,9859242,9848382,9839534,9830358,9818227,9800700,9772419,9741720,9711115,9673162,9655549,9646032,9618756,9580991,9527807,9463667,9378113,9289770,9220578,9183948,9153489) > 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] 50 > (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 4 4 4 4 4 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] 10788500 10181245 9858308 9741720 9183948 [2,] 10708433 10156637 9855372 9711115 9153489 [3,] 10625700 10136811 9855520 9673162 NA [4,] 10547958 10115603 9856303 9655549 NA [5,] 10478617 10084475 9858982 9646032 NA [6,] 10421137 10045158 9859242 9618756 NA [7,] 10376133 10004486 9848382 9580991 NA [8,] 10332785 9967379 9839534 9527807 NA [9,] 10286570 9937697 9830358 9463667 NA [10,] 10251250 9901664 9818227 9378113 NA [11,] 10226419 9870234 9800700 9289770 NA [12,] 10203008 9861823 9772419 9220578 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] -80067 -24608 -2936 -30605 -30459 [2,] -82733 -19826 148 -37953 NA [3,] -77742 -21208 783 -17613 NA [4,] -69341 -31128 2679 -9517 NA [5,] -57480 -39317 260 -27276 NA [6,] -45004 -40672 -10860 -37765 NA [7,] -43348 -37107 -8848 -53184 NA [8,] -46215 -29682 -9176 -64140 NA [9,] -35320 -36033 -12131 -85554 NA [10,] -24831 -31430 -17527 -88343 NA [11,] -23411 -8411 -28281 -69192 NA [12,] -21763 -3515 -30699 -36630 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/17puc1290710361.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/2n1141290710361.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/3n1141290710361.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/4n1141290710361.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] [1,] 9183948 9153489 9673162 9655549 9646032 9618756 9580991 9527807 [2,] 9741720 9711115 9764341 9755926 9752507 9738999 9714686 9683670 [3,] 9858308 9855372 9996166 9985953 9971728 9952200 9926434 9903456 [4,] 10181245 10156637 10381256 10331780 10281546 10233148 10190310 10150082 [5,] 10788500 10708433 10625700 10547958 10478617 10421137 10376133 10332785 [,9] [,10] [,11] [,12] [1,] 9463667 9378113 9289770 9220578 [2,] 9647012 9598170 9545235 9496498 [3,] 9884028 9859946 9835467 9817121 [4,] 10112134 10076457 10048326 10032416 [5,] 10286570 10251250 10226419 10203008 $n [1] 5 5 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 9547741 9540567 9508803 9531028 9553788 9561823 9550692 9534991 [2,] 10168875 10170177 10483528 10440878 10389669 10342577 10302176 10271922 [,9] [,10] [,11] [,12] [1,] 9516582 9482099 9438025 9393747 [2,] 10251473 10237792 10232909 10240495 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(9183948, 9741720, 9858308, 10181245, 10788500, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5s7bi1290710361.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] [1,] -30605 -82733.0 -77742.0 -69341.0 -57480.0 -45004.0 -53184.0 -64140.0 [2,] -30605 -60343.0 -49475.0 -50234.5 -48398.5 -42838.0 -48266.0 -55177.5 [3,] -30459 -28889.5 -19410.5 -20322.5 -33296.5 -39218.5 -40227.5 -37948.5 [4,] -24608 -9839.0 -8415.0 -3419.0 -13508.0 -24312.5 -22977.5 -19429.0 [5,] -24608 148.0 783.0 2679.0 260.0 -10860.0 -8848.0 -9176.0 [,9] [,10] [,11] [,12] [1,] -85554.0 -88343.0 -69192.0 -36630.0 [2,] -60793.5 -59886.5 -48736.5 -33664.5 [3,] -35676.5 -28130.5 -25846.0 -26231.0 [4,] -23725.5 -21179.0 -15911.0 -12639.0 [5,] -12131.0 -17527.0 -8411.0 -3515.0 $n [1] 5 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -34696.47 -68787.66 -51847.9 -57306.75 -60859.995 -53853.65 -60205.42 [2,] -26221.53 11008.66 13026.9 16661.75 -5733.005 -24583.35 -20249.58 [,8] [,9] [,10] [,11] [,12] [1,] -66189.815 -64960.22 -58709.425 -51778.145 -42841.145 [2,] -9707.185 -6392.78 2448.425 86.145 -9620.855 $out [1] -80067 -2936 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-30605, -30605, -30459, -24608, -24608, -82733, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6s7bi1290710361.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,] 10203008 9861823 9800700 9220578 9153489 [2,] 10268910 9919680 9824292 9420890 9153489 [3,] 10398635 10024822 9851877 9599874 9168718 [4,] 10586829 10126207 9857306 9664356 9183948 [5,] 10788500 10181245 9859242 9741720 9183948 $n [1] 12 12 12 12 2 $conf [,1] [,2] [,3] [,4] [,5] [1,] 10253630 9930624 9836820 9488827 9134689 [2,] 10543640 10119020 9866934 9710920 9202748 $out [1] 9772419 $group [1] 3 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(10203008, 10268910, 10398635, 10586829, 10788500, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/73zal1290710361.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,] 9764457 9817121 9790789 [2,] 9858443 9856840 9865215 [3,] 9933877 9893742 9936671 [4,] 10001550 9961964 9981757 [5,] 10072798 9996166 10034482 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 9868605 9845794 9883515 [2,] 9999149 9941690 9989827 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(9764457, 9858443.25, 9933876.7, 10001549.875, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/17puc1290710361.ps tmp/17puc1290710361.png",intern=TRUE)) character(0) > try(system("convert tmp/2n1141290710361.ps tmp/2n1141290710361.png",intern=TRUE)) character(0) > try(system("convert tmp/3n1141290710361.ps tmp/3n1141290710361.png",intern=TRUE)) character(0) > try(system("convert tmp/4n1141290710361.ps tmp/4n1141290710361.png",intern=TRUE)) character(0) > try(system("convert tmp/5s7bi1290710361.ps tmp/5s7bi1290710361.png",intern=TRUE)) character(0) > try(system("convert tmp/6s7bi1290710361.ps tmp/6s7bi1290710361.png",intern=TRUE)) character(0) > try(system("convert tmp/73zal1290710361.ps tmp/73zal1290710361.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.239 0.985 3.226