R version 2.7.0 (2008-04-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(10511,10812,10738,10171,9721,9897,9828,9924,10371,10846,10413,10709,10662,10570,10297,10635,10872,10296,10383,10431,10574,10653,10805,10872,10625,10407,10463,10556,10646,10702,11353,11346,11451,11964,12574,13031,13812,14544,14931,14886,16005,17064,15168,16050,15839,15137,14954,15648,15305,15579,16348,15928,16171,15937,15713,15594,15683,16438,17032,17696,17745) > 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] 61 > (np <- floor(n / par1)) [1] 5 > 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 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 10511 10662 10625 13812 15305 17745 [2,] 10812 10570 10407 14544 15579 NA [3,] 10738 10297 10463 14931 16348 NA [4,] 10171 10635 10556 14886 15928 NA [5,] 9721 10872 10646 16005 16171 NA [6,] 9897 10296 10702 17064 15937 NA [7,] 9828 10383 11353 15168 15713 NA [8,] 9924 10431 11346 16050 15594 NA [9,] 10371 10574 11451 15839 15683 NA [10,] 10846 10653 11964 15137 16438 NA [11,] 10413 10805 12574 14954 17032 NA [12,] 10709 10872 13031 15648 17696 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 301 -92 -218 732 274 NA [2,] -74 -273 56 387 769 NA [3,] -567 338 93 -45 -420 NA [4,] -450 237 90 1119 243 NA [5,] 176 -576 56 1059 -234 NA [6,] -69 87 651 -1896 -224 NA [7,] 96 48 -7 882 -119 NA [8,] 447 143 105 -211 89 NA [9,] 475 79 513 -702 755 NA [10,] -433 152 610 -183 594 NA [11,] 296 67 457 694 664 NA [12,] -47 -247 781 -343 49 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/1p3tv1226236026.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/24una1226236026.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/392xk1226236026.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/41wx01226236026.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] [,12] [1,] 10511 10407 10297 10171 9721 9897 9828 9924 10371 10653 10413 10709 [2,] 10625 10570 10463 10556 10646 10296 10383 10431 10574 10846 10805 10872 [3,] 12237 10812 10738 10635 10872 10702 11353 11346 11451 11964 12574 13031 [4,] 15305 14544 14931 14886 16005 15937 15168 15594 15683 15137 14954 15648 [5,] 17745 15579 16348 15928 16171 17064 15713 16050 15839 16438 17032 17696 $n [1] 6 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 9218.249 8003.982 7580.922 7575.433 7085.344 6716.084 7971.931 [2,] 15255.751 13620.018 13895.078 13694.567 14658.656 14687.916 14734.069 [,8] [,9] [,10] [,11] [,12] [1,] 7697.837 7840.993 8931.99 9642.327 9656.29 [2,] 14994.163 15061.007 14996.01 15505.673 16405.71 $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(10511, 10625, 12237, 15305, 17745, 10407, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/52qh41226236026.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] [,12] [1,] -218 -273 -567 90 -576 -224 -119 89 79 -433 67 -343 [2,] -92 -74 -420 90 -234 -224 -7 89 79 -183 296 -247 [3,] 274 56 -45 237 56 -69 48 105 475 152 457 -47 [4,] 301 387 93 243 176 87 96 143 513 594 664 49 [5,] 732 769 338 243 176 87 96 143 755 610 694 49 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3.69281 -269.7414 -407.4845 128.8906 -233.7050 -288.7518 -24.77954 [2,] 551.69281 381.7414 317.4845 345.1094 345.7050 150.7518 120.77954 [,8] [,9] [,10] [,11] [,12] [1,] 66.84374 168.3367 -397.0262 196.9721 -256.1529 [2,] 143.15626 781.6633 701.0262 717.0279 162.1529 $out [1] -450 1119 1059 651 -1896 882 447 -211 -702 781 $group [1] 4 4 5 6 6 7 8 8 9 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-218, -92, 274, 301, 732, -273, -74, 56, 387, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6bmwj1226236026.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] [1,] 9721.0 10296.0 10407.0 13812.0 15305.0 17745 [2,] 9910.5 10407.0 10590.5 14908.5 15638.5 17745 [3,] 10392.0 10604.5 11024.0 15152.5 15932.5 17745 [4,] 10723.5 10733.5 11707.5 15922.0 16393.0 17745 [5,] 10846.0 10872.0 13031.0 17064.0 17032.0 17745 $n [1] 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 10021.19 10455.58 10514.53 14690.24 15588.37 17745 [2,] 10762.81 10753.42 11533.47 15614.76 16276.63 17745 $out [1] 17696 $group [1] 5 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(9721, 9910.5, 10392, 10723.5, 10846, 10296, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7l9411226236026.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,] 12382.4 10635.0 12557.00 [2,] 12522.2 10775.0 12748.25 [3,] 12731.1 11349.5 12935.50 [4,] 13058.8 12100.5 13122.50 [5,] 13591.2 13031.0 13325.50 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 12486.35 10744.93 12764.80 [2,] 12975.85 11954.07 13106.20 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(12382.4, 12522.2, 12731.1, 13058.8, 13591.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1p3tv1226236026.ps tmp/1p3tv1226236026.png") > system("convert tmp/24una1226236026.ps tmp/24una1226236026.png") > system("convert tmp/392xk1226236026.ps tmp/392xk1226236026.png") > system("convert tmp/41wx01226236026.ps tmp/41wx01226236026.png") > system("convert tmp/52qh41226236026.ps tmp/52qh41226236026.png") > system("convert tmp/6bmwj1226236026.ps tmp/6bmwj1226236026.png") > system("convert tmp/7l9411226236026.ps tmp/7l9411226236026.png") > > > proc.time() user system elapsed 2.556 1.570 2.816