R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(225,243,270,289,273,274,220,271,288,276,266,239,201,242,239,273,280,294,212,264,272,262,238,227,250,245,270,288,298,281,218,284,281,277,276,222,255,267,261,263,264,278,248,320,305,301,274,220,235,252,272,280,305,299,246,307,325,302,274,251,272,253,292,288,258,295,231,250,268)
> 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] 69
> (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 6 6 6 6 6 6 6 6 5 5 5
> arr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,]  225  201  250  255  235  272
 [2,]  243  242  245  267  252  253
 [3,]  270  239  270  261  272  292
 [4,]  289  273  288  263  280  288
 [5,]  273  280  298  264  305  258
 [6,]  274  294  281  278  299  295
 [7,]  220  212  218  248  246  231
 [8,]  271  264  284  320  307  250
 [9,]  288  272  281  305  325  268
[10,]  276  262  277  301  302   NA
[11,]  266  238  276  274  274   NA
[12,]  239  227  222  220  251   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,]   18   41   -5   12   17  -19
 [2,]   27   -3   25   -6   20   39
 [3,]   19   34   18    2    8   -4
 [4,]  -16    7   10    1   25  -30
 [5,]    1   14  -17   14   -6   37
 [6,]  -54  -82  -63  -30  -53  -64
 [7,]   51   52   66   72   61   19
 [8,]   17    8   -3  -15   18   18
 [9,]  -12  -10   -4   -4  -23   NA
[10,]  -10  -24   -1  -27  -28   NA
[11,]  -27  -11  -54  -54  -23   NA
[12,]  -38   23   33   15   21   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/1ohvr1272470238.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/2ohvr1272470238.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/3g8uu1272470238.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/4g8uu1272470238.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,] 201.0 242.0  261  263 258.0 274.0 212.0 250.0 268.0   262   266   220
[2,] 225.0 243.0  261  273 264.0 278.0 218.0 264.0 272.0   276   266   222
[3,] 242.5 248.5  270  284 276.5 287.5 225.5 277.5 284.5   277   274   227
[4,] 255.0 253.0  272  288 298.0 295.0 246.0 307.0 305.0   301   274   239
[5,] 272.0 267.0  272  289 305.0 299.0 248.0 320.0 325.0   302   276   251

$n
 [1] 6 6 6 6 6 6 6 6 6 5 5 5

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 223.1490 242.0497 262.9046 274.3245 254.5689 276.5345 207.4391 249.7636
[2,] 261.8510 254.9503 277.0954 293.6755 298.4311 298.4655 243.5609 305.2364
         [,9]    [,10]    [,11]    [,12]
[1,] 263.2139 259.3351 268.3472 214.9878
[2,] 305.7861 294.6649 279.6528 239.0122

$out
[1] 239 292 238

$group
[1]  3  3 11

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

Warning message:
In bxp(list(stats = c(201, 225, 242.5, 255, 272, 242, 243, 248.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5g8uu1272470238.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,] -19.0 -6.0   -4  -30 -17.0 -64.0 51.0 -15.0  -23   -28   -54    15
[2,]  -5.0 -3.0    2  -16  -6.0 -64.0 51.0  -3.0  -12   -27   -54    15
[3,]  14.5 22.5   13    4   7.5 -58.5 56.5  12.5  -10   -24   -27    21
[4,]  18.0 27.0   19   10  14.0 -53.0 66.0  18.0   -4   -10   -23    23
[5,]  41.0 39.0   34   25  37.0 -53.0 72.0  18.0   -4    -1   -11    33

$n
 [1] 6 6 6 6 6 6 6 6 5 5 5 5

$conf
           [,1]      [,2]      [,3]      [,4]      [,5]      [,6]     [,7]
[1,] -0.3357429  3.149031  2.034451 -12.77084 -5.400646 -65.59536 46.82452
[2,] 29.3357429 41.850969 23.965549  20.77084 20.400646 -51.40464 66.17548
          [,8]      [,9]     [,10]      [,11]    [,12]
[1,] -1.045678 -15.65278 -36.01216 -48.904522 15.34722
[2,] 26.045678  -4.34722 -11.98784  -5.095478 26.65278

$out
[1] -82 -30  19 -38

$group
[1]  6  6  7 12

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

Warning message:
In bxp(list(stats = c(-19, -5, 14.5, 18, 41, -6, -3, 22.5, 27, 39,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6rzuf1272470238.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,] 220.0 201.0 218.0 220.0 235.0  231
[2,] 241.0 232.5 247.5 258.0 251.5  253
[3,] 270.5 252.0 276.5 265.5 277.0  268
[4,] 275.0 272.5 282.5 289.5 303.5  288
[5,] 289.0 294.0 298.0 320.0 325.0  295

$n
[1] 12 12 12 12 12  9

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]
[1,] 254.9924 233.7557 260.5363 251.1326 253.2825 249.5667
[2,] 286.0076 270.2443 292.4637 279.8674 300.7175 286.4333

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(220, 241, 270.5, 275, 289, 201, 232.5, 252,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7rzuf1272470238.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,] 229.1667 225.50 230.3750
[2,] 245.0000 245.50 244.3750
[3,] 273.5000 275.25 274.9375
[4,] 283.1333 280.75 285.1250
[5,] 289.8333 287.50 288.5000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 256.1071 259.1722 256.3512
[2,] 290.8929 291.3278 293.5238

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(229.166666666667, 245, 273.5, 283.133333333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1ohvr1272470238.ps tmp/1ohvr1272470238.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ohvr1272470238.ps tmp/2ohvr1272470238.png",intern=TRUE))
character(0)
> try(system("convert tmp/3g8uu1272470238.ps tmp/3g8uu1272470238.png",intern=TRUE))
character(0)
> try(system("convert tmp/4g8uu1272470238.ps tmp/4g8uu1272470238.png",intern=TRUE))
character(0)
> try(system("convert tmp/5g8uu1272470238.ps tmp/5g8uu1272470238.png",intern=TRUE))
character(0)
> try(system("convert tmp/6rzuf1272470238.ps tmp/6rzuf1272470238.png",intern=TRUE))
character(0)
> try(system("convert tmp/7rzuf1272470238.ps tmp/7rzuf1272470238.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.224   0.957   2.018