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(167.16,179.84,174.44,180.35,193.17,195.16,202.43,189.91,195.98,212.09,205.81,204.31,196.07,199.98,199.10,198.31,195.72,223.04,238.41,259.73,326.54,335.15,321.81,368.62,369.59,425.00,439.72,362.23,328.76,348.55,328.18,329.34,295.55,237.38,226.85,220.14,239.36,224.69,230.98,233.47,256.70,253.41,224.95,210.37,191.09,198.85,211.04,206.25,201.51,194.54,191.07,192.82,181.88,157.67,195.82,246.25,271.69,270.29)
> 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] 58
> (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 5 5 5 5 5 5 5 5 4 4
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]
 [1,] 167.16 196.07 369.59 239.36 201.51
 [2,] 179.84 199.98 425.00 224.69 194.54
 [3,] 174.44 199.10 439.72 230.98 191.07
 [4,] 180.35 198.31 362.23 233.47 192.82
 [5,] 193.17 195.72 328.76 256.70 181.88
 [6,] 195.16 223.04 348.55 253.41 157.67
 [7,] 202.43 238.41 328.18 224.95 195.82
 [8,] 189.91 259.73 329.34 210.37 246.25
 [9,] 195.98 326.54 295.55 191.09 271.69
[10,] 212.09 335.15 237.38 198.85 270.29
[11,] 205.81 321.81 226.85 211.04     NA
[12,] 204.31 368.62 220.14 206.25     NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]
 [1,]  12.68   3.91  55.41 -14.67  -6.97
 [2,]  -5.40  -0.88  14.72   6.29  -3.47
 [3,]   5.91  -0.79 -77.49   2.49   1.75
 [4,]  12.82  -2.59 -33.47  23.23 -10.94
 [5,]   1.99  27.32  19.79  -3.29 -24.21
 [6,]   7.27  15.37 -20.37 -28.46  38.15
 [7,] -12.52  21.32   1.16 -14.58  50.43
 [8,]   6.07  66.81 -33.79 -19.28  25.44
 [9,]  16.11   8.61 -58.17   7.76  -1.40
[10,]  -6.28 -13.34 -10.53  12.19     NA
[11,]  -1.50  46.81  -6.71  -4.79     NA
[12,]  -8.24   0.97  19.22  -4.74     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/1omet1289570535.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/2omet1289570535.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/3omet1289570535.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/4hdee1289570535.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]
[1,] 167.16 179.84 174.44 180.35 181.88 157.67 195.82 189.91 191.09 198.85
[2,] 196.07 194.54 191.07 192.82 193.17 195.16 202.43 210.37 195.98 212.09
[3,] 201.51 199.98 199.10 198.31 195.72 223.04 224.95 246.25 271.69 237.38
[4,] 239.36 224.69 230.98 233.47 256.70 253.41 238.41 259.73 295.55 270.29
[5,] 239.36 224.69 230.98 233.47 328.76 253.41 238.41 329.34 326.54 335.15
       [,11]   [,12]
[1,] 205.810 204.310
[2,] 208.425 205.280
[3,] 218.945 213.195
[4,] 274.330 294.380
[5,] 321.810 368.620

$n
 [1] 5 5 5 5 5 5 5 5 5 5 4 4

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 170.9214 178.6761 170.8997 169.5868 150.8299 181.8807 199.5266 211.3723
[2,] 232.0986 221.2839 227.3003 227.0332 240.6101 264.1993 250.3734 281.1277
         [,9]    [,10]    [,11]   [,12]
[1,] 201.3341 196.2560 166.8801 142.806
[2,] 342.0459 278.5040 271.0100 283.584

$out
[1] 369.59 425.00 439.72 362.23 348.55 328.18

$group
[1] 1 2 3 4 6 7

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12"

Warning message:
In bxp(list(stats = c(167.16, 196.07, 201.51, 239.36, 239.36, 179.84,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5hdee1289570535.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]
[1,] -14.67 -5.40 -0.79 -33.47 -24.21 -28.46 -14.58 -33.79 -1.40 -13.340 -6.710
[2,]  -6.97 -3.47 -0.79 -10.94  -3.29 -20.37 -12.52 -19.28 -1.40 -11.935 -5.750
[3,]   3.91 -0.88  1.75  -2.59   1.99   7.27   1.16   6.07  7.76  -8.405 -3.145
[4,]  12.68  6.29  2.49  12.82  19.79  15.37  21.32  25.44  8.61   2.955 22.655
[5,]  12.68 14.72  5.91  23.23  27.32  38.15  50.43  66.81 16.11  12.190 46.810
      [,12]
[1,] -8.240
[2,] -6.490
[3,] -1.885
[4,] 10.095
[5,] 19.220

$n
 [1] 5 5 5 5 5 5 5 5 5 4 4 4

$conf
         [,1]      [,2]       [,3]      [,4]      [,5]      [,6]      [,7]
[1,] -9.97464 -7.776391 -0.5676397 -19.37876 -14.31827 -17.98379 -22.75126
[2,] 17.79464  6.016391  4.0676397  14.19876  18.29827  32.52379  25.07126
          [,8]       [,9]    [,10]     [,11]     [,12]
[1,] -25.52904  0.6869592 -20.1681 -25.58495 -14.98715
[2,]  37.66904 14.8330408   3.3581  19.29495  11.21715

$out
[1]  55.41 -77.49 -58.17

$group
[1] 1 3 9

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12"

Warning message:
In bxp(list(stats = c(-14.6700000000000, -6.97, 3.91, 12.68, 12.68,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6rndz1289570535.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,] 167.160 195.720 220.140 191.090 157.67
[2,] 180.095 198.705 266.465 208.310 191.07
[3,] 194.165 230.725 329.050 224.820 195.18
[4,] 203.370 324.175 365.910 236.415 246.25
[5,] 212.090 368.620 439.720 256.700 271.69

$n
[1] 12 12 12 12 10

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]
[1,] 183.5491 173.4973 283.6925 212.0011 167.6099
[2,] 204.7809 287.9527 374.4075 237.6389 222.7501

$out
numeric(0)

$group
numeric(0)

$names
[1] "1" "2" "3" "4" NA 

Warning message:
In bxp(list(stats = c(167.16, 180.095, 194.165, 203.37, 212.09,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7rndz1289570535.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,] 231.2460 195.720 209.6150
[2,] 235.1520 199.540 215.4300
[3,] 243.0938 216.070 224.6100
[4,] 248.4750 231.165 233.2812
[5,] 256.1700 271.690 245.7650

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 237.0170 201.6456 216.4679
[2,] 249.1705 230.4944 232.7521

$out
numeric(0)

$group
numeric(0)

$names
[1] "mean"     "median"   "midrange"

Warning message:
In bxp(list(stats = c(231.246, 235.152, 243.09375, 248.475, 256.17,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1omet1289570535.ps tmp/1omet1289570535.png",intern=TRUE))
character(0)
> try(system("convert tmp/2omet1289570535.ps tmp/2omet1289570535.png",intern=TRUE))
character(0)
> try(system("convert tmp/3omet1289570535.ps tmp/3omet1289570535.png",intern=TRUE))
character(0)
> try(system("convert tmp/4hdee1289570535.ps tmp/4hdee1289570535.png",intern=TRUE))
character(0)
> try(system("convert tmp/5hdee1289570535.ps tmp/5hdee1289570535.png",intern=TRUE))
character(0)
> try(system("convert tmp/6rndz1289570535.ps tmp/6rndz1289570535.png",intern=TRUE))
character(0)
> try(system("convert tmp/7rndz1289570535.ps tmp/7rndz1289570535.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.924   1.454   2.139