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(681,608,678,739,707,684,596,729,732,386,331,703,739,651,653,628,658,724,674,631,686,392,348,784,875,659,629,679,615,716,723,641,754,729,535,922,828,710,858,775,785,1006,789,734,906,536,394,991,841,892,782,813,793,978,775,797,946,594,438,1022,868)
> 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,]  681  739  875  828  841  868
 [2,]  608  651  659  710  892   NA
 [3,]  678  653  629  858  782   NA
 [4,]  739  628  679  775  813   NA
 [5,]  707  658  615  785  793   NA
 [6,]  684  724  716 1006  978   NA
 [7,]  596  674  723  789  775   NA
 [8,]  729  631  641  734  797   NA
 [9,]  732  686  754  906  946   NA
[10,]  386  392  729  536  594   NA
[11,]  331  348  535  394  438   NA
[12,]  703  784  922  991 1022   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,]  -73  -88 -216 -118   51   NA
 [2,]   70    2  -30  148 -110   NA
 [3,]   61  -25   50  -83   31   NA
 [4,]  -32   30  -64   10  -20   NA
 [5,]  -23   66  101  221  185   NA
 [6,]  -88  -50    7 -217 -203   NA
 [7,]  133  -43  -82  -55   22   NA
 [8,]    3   55  113  172  149   NA
 [9,] -346 -294  -25 -370 -352   NA
[10,]  -55  -44 -194 -142 -156   NA
[11,]  372  436  387  597  584   NA
[12,]   36   91  -94 -150 -154   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/15nwz1290817424.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/25nwz1290817424.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/35nwz1290817424.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/4yed31290817424.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,] 681.0  608  629  628  615  684  596  631  686   386   331   703
[2,] 739.0  651  653  679  658  716  674  641  732   392   348   784
[3,] 834.5  659  678  739  707  724  723  729  754   536   394   922
[4,] 868.0  710  782  775  785  978  775  734  906   594   438   991
[5,] 875.0  710  858  813  793 1006  789  797  946   729   535  1022

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 751.2908 617.3107 586.8489 671.1666 617.2621 538.8715 651.6337 663.2864
[2,] 917.7092 700.6893 769.1511 806.8334 796.7379 909.1285 794.3663 794.7136
        [,9]    [,10]    [,11]     [,12]
[1,] 631.052 393.2673 330.4062  775.7343
[2,] 876.948 678.7327 457.5938 1068.2657

$out
[1] 892

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(681, 739, 834.5, 868, 875, 608, 651, 659,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5yed31290817424.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,] -118 -110  -83  -64  -23 -217  -82    3 -370  -194   372  -154
[2,] -118  -30  -25  -32   66 -203  -55   55 -352  -156   387  -150
[3,]  -88    2   31  -20  101  -88  -43  113 -346  -142   436   -94
[4,]  -73   70   50   10  185  -50   22  149 -294   -55   584    36
[5,]  -73  148   61   30  221    7  133  172 -294   -44   597    91

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

$conf
           [,1]      [,2]      [,3]       [,4]     [,5]       [,6]      [,7]
[1,] -119.79689 -68.65975 -21.99481 -49.677094  16.9149 -196.10941 -97.40801
[2,]  -56.20311  72.65975  83.99481   9.677094 185.0851   20.10941  11.40801
          [,8]      [,9]      [,10]    [,11]      [,12]
[1,]  46.57984 -386.9827 -213.36635 296.8003 -225.42713
[2,] 179.42016 -305.0173  -70.63365 575.1997   37.42713

$out
[1] -216   51  -25

$group
[1] 1 1 9

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

Warning message:
In bxp(list(stats = c(-118, -118, -88, -73, -73, -110, -30, 2, 70,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/685c51290817424.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,] 596.0 628.0 535.0  536  594.0  868
[2,] 602.0 629.5 635.0  722  778.5  868
[3,] 682.5 655.5 697.5  787  805.0  868
[4,] 718.0 705.0 741.5  882  919.0  868
[5,] 739.0 784.0 875.0 1006 1022.0  868

$n
[1] 12 12 12 12 12  1

$conf
         [,1]    [,2]     [,3]     [,4]    [,5] [,6]
[1,] 629.5916 621.064 648.9246 714.0229 740.917  868
[2,] 735.4084 689.936 746.0754 859.9771 869.083  868

$out
[1] 386 331 392 348 922 394 438

$group
[1] 1 1 2 2 3 4 5

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

Warning message:
In bxp(list(stats = c(596, 602, 682.5, 718, 739, 628, 629.5, 655.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/785c51290817424.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,] 704.0000 659.0 493.000
[2,] 705.2000 668.5 684.000
[3,] 715.8000 723.5 723.000
[4,] 805.0667 746.5 815.125
[5,] 884.4000 834.5 887.500

$n
[1] 12 12 12

$conf
         [,1]     [,2]    [,3]
[1,] 670.2501 687.9237 663.193
[2,] 761.3499 759.0763 782.807

$out
[1] 527.4 409.2 536.0 394.0 922.0 393.0

$group
[1] 1 1 2 2 2 3

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

Warning message:
In bxp(list(stats = c(704, 705.2, 715.8, 805.066666666667, 884.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/15nwz1290817424.ps tmp/15nwz1290817424.png",intern=TRUE))
character(0)
> try(system("convert tmp/25nwz1290817424.ps tmp/25nwz1290817424.png",intern=TRUE))
character(0)
> try(system("convert tmp/35nwz1290817424.ps tmp/35nwz1290817424.png",intern=TRUE))
character(0)
> try(system("convert tmp/4yed31290817424.ps tmp/4yed31290817424.png",intern=TRUE))
character(0)
> try(system("convert tmp/5yed31290817424.ps tmp/5yed31290817424.png",intern=TRUE))
character(0)
> try(system("convert tmp/685c51290817424.ps tmp/685c51290817424.png",intern=TRUE))
character(0)
> try(system("convert tmp/785c51290817424.ps tmp/785c51290817424.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.941   1.475   2.130