R version 2.12.0 (2010-10-15) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-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. 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(16,29,22,30,20,39,18,9.6,10.2,20.2,50,120,19.8,18,3,11,15,27,28,14,5.6,6.5,8.5,87.9,5.8,25.2,7.5,13.7,34,17,9,9.2,5,24,40,86.5,0.54,14,4.8,28,16,5.8,16,9.1,6,17,26,99.6,41,72,23,42,40,18,45,18,2,10,13.6,160) > 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] 60 > (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] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 16.0 19.8 5.8 0.54 41.0 NA [2,] 29.0 18.0 25.2 14.00 72.0 NA [3,] 22.0 3.0 7.5 4.80 23.0 NA [4,] 30.0 11.0 13.7 28.00 42.0 NA [5,] 20.0 15.0 34.0 16.00 40.0 NA [6,] 39.0 27.0 17.0 5.80 18.0 NA [7,] 18.0 28.0 9.0 16.00 45.0 NA [8,] 9.6 14.0 9.2 9.10 18.0 NA [9,] 10.2 5.6 5.0 6.00 2.0 NA [10,] 20.2 6.5 24.0 17.00 10.0 NA [11,] 50.0 8.5 40.0 26.00 13.6 NA [12,] 120.0 87.9 86.5 99.60 160.0 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 13.0 -1.8 19.40 13.46 31.0 NA [2,] -7.0 -15.0 -17.70 -9.20 -49.0 NA [3,] 8.0 8.0 6.20 23.20 19.0 NA [4,] -10.0 4.0 20.30 -12.00 -2.0 NA [5,] 19.0 12.0 -17.00 -10.20 -22.0 NA [6,] -21.0 1.0 -8.00 10.20 27.0 NA [7,] -8.4 -14.0 0.20 -6.90 -27.0 NA [8,] 0.6 -8.4 -4.20 -3.10 -16.0 NA [9,] 10.0 0.9 19.00 11.00 8.0 NA [10,] 29.8 2.0 16.00 9.00 3.6 NA [11,] 70.0 79.4 46.50 73.60 146.4 NA [12,] -100.2 -82.1 -85.96 -58.60 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/rcomp/tmp/1aqz21292078370.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/rcomp/tmp/23ig51292078370.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/rcomp/tmp/33ig51292078370.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/rcomp/tmp/43ig51292078370.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] [,9] [,10] [,11] [,12] [1,] 0.54 14.0 3.0 11.0 15 5.8 9 9.1 5.0 6.5 8.5 86.5 [2,] 5.80 18.0 4.8 13.7 16 17.0 16 9.2 5.0 10.0 13.6 87.9 [3,] 16.00 25.2 7.5 28.0 20 18.0 18 9.6 5.6 17.0 26.0 99.6 [4,] 19.80 29.0 22.0 30.0 34 27.0 28 14.0 6.0 20.2 40.0 120.0 [5,] 19.80 29.0 23.0 42.0 40 39.0 45 18.0 6.0 24.0 50.0 160.0 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 6.107635 17.42743 -4.653477 16.48246 7.281245 10.93403 9.52083 [2,] 25.892365 32.97257 19.653477 39.51754 32.718755 25.06597 26.47917 [,8] [,9] [,10] [,11] [,12] [1,] 6.208332 4.893403 9.792706 7.345827 76.91822 [2,] 12.991668 6.306597 24.207294 44.654173 122.28178 $out [1] 41.0 72.0 10.2 2.0 $group [1] 1 2 9 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(0.54, 5.8, 16, 19.8, 19.8, 14, 18, 25.2, 29, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/53ig51292078370.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] [,8] [,9] [,10] [,11] [,12] [1,] 13.00 -17.7 6.2 -12.0 -22.0 -21.0 -14.0 -16.0 8 2.0 70.0 -100.20 [2,] 13.00 -17.7 8.0 -10.0 -17.0 -8.0 -14.0 -8.4 8 3.6 70.0 -93.08 [3,] 13.46 -15.0 8.0 -2.0 -10.2 1.0 -8.4 -4.2 10 9.0 73.6 -84.03 [4,] 19.40 -9.2 19.0 4.0 12.0 10.2 -6.9 -3.1 11 16.0 79.4 -70.35 [5,] 19.40 -7.0 23.2 20.3 19.0 27.0 0.2 0.6 11 29.8 79.4 -58.60 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 8.937776 -21.006079 0.2274277 -11.892365 -30.69133 -11.86007 -13.416842 [2,] 17.982224 -8.993921 15.7725723 7.892365 10.29133 13.86007 -3.383158 [,8] [,9] [,10] [,11] [,12] [1,] -7.9449666 7.880208 0.2381912 66.95798 -101.9867 [2,] -0.4550334 12.119792 17.7618088 80.24202 -66.0733 $out [1] -1.8 31.0 -49.0 -27.0 0.9 19.0 46.5 146.4 $group [1] 1 1 2 7 9 9 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(13, 13, 13.46, 19.4, 19.4, -17.7, -17.7, -15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6wrf81292078370.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] [1,] 9.6 3.0 5.00 0.54 2.0 NA [2,] 17.0 7.5 8.25 5.90 15.8 NA [3,] 21.1 14.5 15.35 15.00 31.5 NA [4,] 34.5 23.4 29.60 21.50 43.5 NA [5,] 50.0 28.0 40.00 28.00 72.0 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 13.11813 7.247903 5.612122 7.884735 18.86584 NA [2,] 29.08187 21.752097 25.087878 22.115265 44.13416 NA $out [1] 120.0 87.9 86.5 99.6 160.0 $group [1] 1 2 3 4 5 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(9.6, 17, 21.1, 34.5, 50, 3, 7.5, 14.5, 23.4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7wrf81292078370.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,] 5.76 5.6 5.500 [2,] 13.80 12.8 13.100 [3,] 22.28 18.0 21.925 [4,] 26.31 25.6 24.250 [5,] 31.64 28.0 26.800 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 16.57411 12.16183 16.83941 [2,] 27.98589 23.83817 27.01059 $out [1] 110.80 99.60 103.95 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(5.76, 13.8, 22.28, 26.31, 31.64, 5.6, 12.8, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > try(system("convert tmp/1aqz21292078370.ps tmp/1aqz21292078370.png",intern=TRUE)) character(0) > try(system("convert tmp/23ig51292078370.ps tmp/23ig51292078370.png",intern=TRUE)) character(0) > try(system("convert tmp/33ig51292078370.ps tmp/33ig51292078370.png",intern=TRUE)) character(0) > try(system("convert tmp/43ig51292078370.ps tmp/43ig51292078370.png",intern=TRUE)) character(0) > try(system("convert tmp/53ig51292078370.ps tmp/53ig51292078370.png",intern=TRUE)) character(0) > try(system("convert tmp/6wrf81292078370.ps tmp/6wrf81292078370.png",intern=TRUE)) character(0) > try(system("convert tmp/7wrf81292078370.ps tmp/7wrf81292078370.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.520 0.890 2.379