R version 2.8.1 (2008-12-22) Copyright (C) 2008 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(11025.1,10853.8,12252.6,11839.4,11669.1,11601.4,11178.4,9516.4,12102.8,12989,11610.2,10205.5,11356.2,11307.1,12648.6,11947.2,11714.1,12192.5,11268.8,9097.4,12639.8,13040.1,11687.3,11191.7,11391.9,11793.1,13933.2,12778.1,11810.3,13698.4,11956.6,10723.8,13938.9,13979.8,13807.4,12973.9,12509.8,12934.1,14908.3,13772.1,13012.6,14049.9,11816.5,11593.2,14466.2,13615.9,14733.9,13880.7,13527.5,13584,16170.2,13260.6,14741.9,15486.5,13154.5,12621.2,15031.6,15452.4,15428,13105.9,14716.8,14180,16202.2,14392.4,15140.6,15960.1,14351.3,13230.2,15202.1,17056,16077.7,13348.2,16707.5,16792.6,16831.3,17804.5,16370.2,17602.5,17065.6,14427.9,17818.5,18027.6) > 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] 82 > (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 7 7 7 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 11025.1 11356.2 11391.9 12509.8 13527.5 14716.8 16707.5 [2,] 10853.8 11307.1 11793.1 12934.1 13584.0 14180.0 16792.6 [3,] 12252.6 12648.6 13933.2 14908.3 16170.2 16202.2 16831.3 [4,] 11839.4 11947.2 12778.1 13772.1 13260.6 14392.4 17804.5 [5,] 11669.1 11714.1 11810.3 13012.6 14741.9 15140.6 16370.2 [6,] 11601.4 12192.5 13698.4 14049.9 15486.5 15960.1 17602.5 [7,] 11178.4 11268.8 11956.6 11816.5 13154.5 14351.3 17065.6 [8,] 9516.4 9097.4 10723.8 11593.2 12621.2 13230.2 14427.9 [9,] 12102.8 12639.8 13938.9 14466.2 15031.6 15202.1 17818.5 [10,] 12989.0 13040.1 13979.8 13615.9 15452.4 17056.0 18027.6 [11,] 11610.2 11687.3 13807.4 14733.9 15428.0 16077.7 NA [12,] 10205.5 11191.7 12973.9 13880.7 13105.9 13348.2 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -171.3 -49.1 401.2 424.3 56.5 -536.8 85.1 [2,] 1398.8 1341.5 2140.1 1974.2 2586.2 2022.2 38.7 [3,] -413.2 -701.4 -1155.1 -1136.2 -2909.6 -1809.8 973.2 [4,] -170.3 -233.1 -967.8 -759.5 1481.3 748.2 -1434.3 [5,] -67.7 478.4 1888.1 1037.3 744.6 819.5 1232.3 [6,] -423.0 -923.7 -1741.8 -2233.4 -2332.0 -1608.8 -536.9 [7,] -1662.0 -2171.4 -1232.8 -223.3 -533.3 -1121.1 -2637.7 [8,] 2586.4 3542.4 3215.1 2873.0 2410.4 1971.9 3390.6 [9,] 886.2 400.3 40.9 -850.3 420.8 1853.9 209.1 [10,] -1378.8 -1352.8 -172.4 1118.0 -24.4 -978.3 NA [11,] -1404.7 -495.6 -833.5 -853.2 -2322.1 -2729.5 NA [12,] 1150.7 200.2 -464.1 -353.2 1610.9 3359.3 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/16u611240832030.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/2ymfo1240832030.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/3og6v1240832030.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/452k21240832030.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,] 11025.10 10853.8 12252.6 11839.40 11669.10 11601.40 11178.40 9097.4 [2,] 11374.05 11550.1 13290.9 12362.65 11762.20 12945.45 11542.65 10120.1 [3,] 12509.80 12934.1 14908.3 13260.60 13012.60 14049.90 11956.60 11593.2 [4,] 14122.15 13882.0 16186.2 14082.25 14941.25 15723.30 13752.90 12925.7 [5,] 16707.50 16792.6 16831.3 14392.40 16370.20 17602.50 17065.60 14427.9 [,9] [,10] [,11] [,12] [1,] 12102.80 12989.0 11610.20 10205.5 [2,] 13289.35 13328.0 11687.30 11191.7 [3,] 14466.20 13979.8 14270.65 13039.9 [4,] 15116.85 16254.2 15428.00 13348.2 [5,] 17818.50 18027.6 16077.70 13880.7 $n [1] 7 7 7 7 7 7 7 7 7 7 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 10868.68 11541.53 13179.27 12233.68 11114.12 12391.01 10636.67 9917.74 [2,] 14150.92 14326.67 16637.33 14287.52 14911.08 15708.79 13276.53 13268.66 [,9] [,10] [,11] [,12] [1,] 13374.85 12232.32 11857.78 11648.89 [2,] 15557.55 15727.28 16683.52 14430.91 $out [1] 17804.5 $group [1] 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(11025.1, 11374.05, 12509.8, 14122.15, 16707.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/536kb1240832030.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] [1,] -536.80 1341.50 -1809.80 -1434.30 -67.7 -2332.0 -2637.7 1971.90 40.9 [2,] -110.20 1370.15 -1482.45 -863.65 611.5 -1987.6 -1916.7 2498.40 125.0 [3,] 56.50 1974.20 -1136.20 -233.10 819.5 -1608.8 -1232.8 2873.00 400.3 [4,] 243.15 2081.15 -557.30 288.95 1134.8 -730.3 -827.2 3302.85 653.5 [5,] 424.30 2586.20 -413.20 1481.30 1888.1 -423.0 -223.3 3542.40 886.2 [,10] [,11] [,12] [1,] -1378.80 -2729.50 -464.10 [2,] -1352.80 -2322.10 -353.20 [3,] -575.35 -1128.95 675.45 [4,] -24.40 -833.50 1610.90 [5,] 1118.00 -495.60 3359.30 $n [1] 7 7 7 7 7 7 7 7 7 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -154.5149 1549.602 -1688.6847 -921.4141 506.9937 -2359.6393 -1883.4318 [2,] 267.5149 2398.798 -583.7153 455.2141 1132.0063 -857.9607 -582.1682 [,8] [,9] [,10] [,11] [,12] [1,] 2392.595 84.68833 -1432.2109 -2089.1451 -591.4579 [2,] 3353.405 715.91167 281.5109 -168.7549 1942.3579 $out [1] 38.7 -2909.6 973.2 -850.3 1853.9 $group [1] 2 3 3 9 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-536.799999999999, -110.200000000001, 56.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6hee81240832030.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,] 9516.40 11191.70 10723.8 11593.20 12621.20 13230.20 16370.20 [2,] 10939.45 11287.95 11801.7 12721.95 13207.55 14265.65 16707.50 [3,] 11605.80 11700.70 12876.0 13694.00 14162.95 14928.70 16948.45 [4,] 11971.10 12416.15 13870.3 14258.05 15440.20 16018.90 17804.50 [5,] 12989.00 13040.10 13979.8 14908.30 16170.20 17056.00 18027.60 $n [1] 12 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 11135.26 11186.12 11932.50 12993.37 13144.62 14129.03 16400.35 [2,] 12076.34 12215.28 13819.50 14394.63 15181.28 15728.37 17496.55 $out [1] 9097.4 14427.9 $group [1] 2 7 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(9516.4, 10939.45, 11605.8, 11971.1, 12989, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/79x0l1240832030.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,] 11601.44 11593.20 11522.90 [2,] 13001.89 12721.95 12681.91 [3,] 13589.51 13150.25 13287.09 [4,] 14413.66 14160.27 14268.74 [5,] 14880.11 14908.30 14791.10 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 12945.59 12494.22 12563.33 [2,] 14233.42 13806.28 14010.85 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(11601.4428571429, 13001.8928571429, 13589.5071428571, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/16u611240832030.ps tmp/16u611240832030.png") > system("convert tmp/2ymfo1240832030.ps tmp/2ymfo1240832030.png") > system("convert tmp/3og6v1240832030.ps tmp/3og6v1240832030.png") > system("convert tmp/452k21240832030.ps tmp/452k21240832030.png") > system("convert tmp/536kb1240832030.ps tmp/536kb1240832030.png") > system("convert tmp/6hee81240832030.ps tmp/6hee81240832030.png") > system("convert tmp/79x0l1240832030.ps tmp/79x0l1240832030.png") > > > proc.time() user system elapsed 1.180 0.942 2.724