R version 2.8.0 (2008-10-20) 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. 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(23,15,25,18,21,19,15,22,19,20,26,26,21,18,19,19,18,19,24,28,20,29,27,18,19,24,21,22,25,19,15,34,23,19,26,15,15,17,30,19,28,23,23,21,18,19,24,15,20,24,9,20,20,10,44,20,20,20,11,21,21,19,21,17,16,14,19,21,16,19,19,16,24,29,21,20,19,23,18,19,23,19,21,26,13,23,16,17,30,19,22,14,14,21,21,33,23,30,19,21,25,18,29,25,21,16,17,23,26,18,19,28,20,29,19,18,25,15,24,12,11,19,25,12,15,25,14,19,23,19,24,20,16,13,20,30,18,22,21,25,18,25,44,12,28,17,26,18,21,24,20,24,28,20,33,19,19,25,35) > 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] 159 > (np <- floor(n / par1)) [1] 13 > 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] 14 14 14 13 13 13 13 13 13 13 13 13 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 23 21 19 15 20 21 24 13 23 26 11 16 28 [2,] 15 18 24 17 24 19 29 23 30 18 19 13 17 [3,] 25 19 21 30 9 21 21 16 19 19 25 20 26 [4,] 18 19 22 19 20 17 20 17 21 28 12 30 18 [5,] 21 18 25 28 20 16 19 30 25 20 15 18 21 [6,] 19 19 19 23 10 14 23 19 18 29 25 22 24 [7,] 15 24 15 23 44 19 18 22 29 19 14 21 20 [8,] 22 28 34 21 20 21 19 14 25 18 19 25 24 [9,] 19 20 23 18 20 16 23 14 21 25 23 18 28 [10,] 20 29 19 19 20 19 19 21 16 15 19 25 20 [11,] 26 27 26 24 11 19 21 21 17 24 24 44 33 [12,] 26 18 15 15 21 16 26 33 23 12 20 12 19 [,14] [1,] 19 [2,] 25 [3,] 35 [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] -8 -3 5 2 4 -2 5 10 7 -8 8 -3 -11 [2,] 10 1 -3 13 -15 2 -8 -7 -11 1 6 7 9 [3,] -7 0 1 -11 11 -4 -1 1 2 9 -13 10 -8 [4,] 3 -1 3 9 0 -1 -1 13 4 -8 3 -12 3 [5,] -2 1 -6 -5 -10 -2 4 -11 -7 9 10 4 3 [6,] -4 5 -4 0 34 5 -5 3 11 -10 -11 -1 -4 [7,] 7 4 19 -2 -24 2 1 -8 -4 -1 5 4 4 [8,] -3 -8 -11 -3 0 -5 4 0 -4 7 4 -7 4 [9,] 1 9 -4 1 0 3 -4 7 -5 -10 -4 7 -8 [10,] 6 -2 7 5 -9 0 2 0 1 9 5 19 13 [11,] 0 -9 -11 -9 10 -3 5 12 6 -12 -4 -32 -14 [12,] -5 1 0 5 0 8 -13 -10 3 -1 -4 16 0 [,14] [1,] 6 [2,] 10 [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] 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/freestat/rcomp/tmp/1d48m1289569512.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/freestat/rcomp/tmp/2d48m1289569512.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/freestat/rcomp/tmp/3owq71289569512.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/freestat/rcomp/tmp/4owq71289569512.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,] 11.0 13 16 17 15 14 14 14 14 19 17 12 [2,] 16.0 17 19 18 18 19 18 19 18 19 21 15 [3,] 20.5 19 21 19 20 19 20 21 20 19 24 19 [4,] 23.0 24 25 21 25 23 23 25 23 20 26 23 [5,] 28.0 30 30 22 30 29 29 34 28 21 33 33 $n [1] 14 14 14 13 13 13 13 13 13 13 13 13 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 17.54409 16.04409 18.46636 17.68536 16.93251 17.24715 17.80893 18.37072 [2,] 23.45591 21.95591 23.53364 20.31464 23.06749 20.75285 22.19107 23.62928 [,9] [,10] [,11] [,12] [1,] 17.80893 18.56179 21.80893 15.49429 [2,] 22.19107 19.43821 26.19107 22.50571 $out [1] 9 35 28 12 30 10 44 29 16 15 25 11 44 $group [1] 3 3 4 4 4 6 7 10 10 10 10 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(11, 16, 20.5, 23, 28, 13, 17, 19, 24, 30, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5owq71289569512.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,] -11 -15.0 -13 -1 -11 -11 -8 -11 -10 -9 -32 -13 [2,] -3 -7.0 -7 -1 -6 -4 -2 -5 -4 0 -11 -4 [3,] 3 1.5 0 3 -2 -1 2 -3 0 5 -4 0 [4,] 6 9.0 2 3 4 5 4 4 3 7 5 3 [5,] 10 13.0 11 9 10 11 7 7 9 13 12 8 $n [1] 14 14 13 13 13 13 13 13 13 13 13 13 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -0.8004549 -5.256364 -3.943918 1.247147 -6.382132 -4.943918 -0.6292789 [2,] 6.8004549 8.256364 3.943918 4.752853 2.382132 2.943918 4.6292789 [,8] [,9] [,10] [,11] [,12] [1,] -6.9439184 -3.067492 1.932508 -11.011410 -3.067492 [2,] 0.9439184 3.067492 8.067492 3.011410 3.067492 $out [1] 13 -8 -12 34 19 -24 19 16 $group [1] 4 4 4 6 7 7 10 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-11, -3, 3, 6, 10, -15, -7, 1.5, 9, 13, -13, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6y57a1289569512.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] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 15.0 18.0 15.0 15.0 9.0 14 18.0 13.0 16.0 12.0 11.0 12.0 17.0 19 [2,] 18.5 18.5 19.0 17.5 15.5 16 19.0 15.0 18.5 18.0 14.5 17.0 19.5 22 [3,] 20.5 19.5 21.5 20.0 20.0 19 21.0 20.0 22.0 19.5 19.0 20.5 22.5 25 [4,] 24.0 25.5 24.5 23.5 20.5 20 23.5 22.5 25.0 25.5 23.5 25.0 27.0 30 [5,] 26.0 29.0 26.0 30.0 24.0 21 29.0 33.0 30.0 29.0 25.0 30.0 33.0 35 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 17.99141 16.30725 18.99141 17.26336 17.71947 17.17557 18.94752 16.5792 [2,] 23.00859 22.69275 24.00859 22.73664 22.28053 20.82443 23.05248 23.4208 [,9] [,10] [,11] [,12] [,13] [,14] [1,] 19.03531 16.0792 14.89504 16.85115 19.0792 17.70229 [2,] 24.96469 22.9208 23.10496 24.14885 25.9208 32.29771 $out [1] 34 44 44 $group [1] 3 5 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" NA Warning message: In bxp(list(stats = c(15, 18.5, 20.5, 24, 26, 18, 18.5, 19.5, 25.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7y57a1289569512.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,] 19.69231 19.00 19.0000 [2,] 20.07692 19.00 19.6875 [3,] 20.70055 20.00 20.5625 [4,] 21.81319 20.75 21.7500 [5,] 24.38462 21.00 23.5000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 19.90863 19.20181 19.62178 [2,] 21.49247 20.79819 21.50322 $out [1] 24 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(19.6923076923077, 20.0769230769231, 20.7005494505495, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1d48m1289569512.ps tmp/1d48m1289569512.png",intern=TRUE)) character(0) > try(system("convert tmp/2d48m1289569512.ps tmp/2d48m1289569512.png",intern=TRUE)) character(0) > try(system("convert tmp/3owq71289569512.ps tmp/3owq71289569512.png",intern=TRUE)) character(0) > try(system("convert tmp/4owq71289569512.ps tmp/4owq71289569512.png",intern=TRUE)) character(0) > try(system("convert tmp/5owq71289569512.ps tmp/5owq71289569512.png",intern=TRUE)) character(0) > try(system("convert tmp/6y57a1289569512.ps tmp/6y57a1289569512.png",intern=TRUE)) character(0) > try(system("convert tmp/7y57a1289569512.ps tmp/7y57a1289569512.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.946 1.476 2.163