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(37,30,47,35,30,43,82,40,47,19,52,136,80,42,54,66,81,63,137,72,107,58,36,52,79,77,54,84,48,96,83,66,61,53,30,74,69,59,42,65,70,100,63,105,82,81,75,102,121,98,76,77,63,37,35,23,40,29,37,51,20,28,13,22,25,13,16,13,16,17,9,17,25,14,8,7,10,7,10,3)
> 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] 80
> (np <- floor(n / par1))
[1] 6
> 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] 7 7 7 7 7 7 7 7 6 6 6 6
> arr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,]   37   80   79   69  121   20   25
 [2,]   30   42   77   59   98   28   14
 [3,]   47   54   54   42   76   13    8
 [4,]   35   66   84   65   77   22    7
 [5,]   30   81   48   70   63   25   10
 [6,]   43   63   96  100   37   13    7
 [7,]   82  137   83   63   35   16   10
 [8,]   40   72   66  105   23   13    3
 [9,]   47  107   61   82   40   16   NA
[10,]   19   58   53   81   29   17   NA
[11,]   52   36   30   75   37    9   NA
[12,]  136   52   74  102   51   17   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,]   -7  -38   -2  -10  -23    8  -11
 [2,]   17   12  -23  -17  -22  -15   -6
 [3,]  -12   12   30   23    1    9   -1
 [4,]   -5   15  -36    5  -14    3    3
 [5,]   13  -18   48   30  -26  -12   -3
 [6,]   39   74  -13  -37   -2    3    3
 [7,]  -42  -65  -17   42  -12   -3   -7
 [8,]    7   35   -5  -23   17    3   NA
 [9,]  -28  -49   -8   -1  -11    1   NA
[10,]   33  -22  -23   -6    8   -8   NA
[11,]   84   16   44   27   14    8   NA
[12,]  -56   27   -5   19  -31    8   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/1slek1288699926.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/2slek1288699926.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/3slek1288699926.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/43uvn1288699926.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,]  20.0   14  8.0  7.0 10.0   7.0  10.0    3   16    17   9.0    17
[2,]  31.0   29 27.5 28.5 27.5  25.0  25.5   18   40    19  30.0    51
[3,]  69.0   42 47.0 65.0 48.0  43.0  63.0   40   54    41  36.5    63
[4,]  79.5   68 54.0 71.5 66.5  79.5  82.5   69   82    58  52.0   102
[5,] 121.0   98 76.0 84.0 81.0 100.0 137.0  105  107    81  75.0   136

$n
 [1] 7 7 7 7 7 7 7 7 6 6 6 6

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]      [,8]
[1,] 40.03658 18.70983 31.17463 39.32109 24.70983 10.45348 28.96052  9.543623
[2,] 97.96342 65.29017 62.82537 90.67891 71.29017 75.54652 97.03948 70.456377
         [,9]    [,10]    [,11]    [,12]
[1,] 26.90864 15.84374 22.30929 30.10335
[2,] 81.09136 66.15626 50.69071 95.89665

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(20, 31, 69, 79.5, 121, 14, 29, 42, 68, 98,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/53uvn1288699926.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,] -23.0 -23.0 -12.0 -14.0 -26.0 -37.0 -65.0  -23 -49.0   -23   8.0 -56.0
[2,] -17.0 -19.5   0.0  -9.5 -15.0  -7.5 -29.5   -5 -28.0   -22  14.0 -31.0
[3,] -10.0 -15.0   9.0   3.0  -3.0   3.0 -12.0    5  -9.5    -7  21.5   1.5
[4,]  -4.5   3.0  17.5   4.0  21.5  21.0  -5.0   17  -1.0     8  44.0  19.0
[5,]   8.0  17.0  30.0  15.0  48.0  39.0  -3.0   35   1.0    33  84.0  27.0

$n
 [1] 7 7 7 7 7 7 7 6 6 6 6 6

$conf
           [,1]       [,2]      [,3]      [,4]      [,5]      [,6]       [,7]
[1,] -17.464798 -28.436637 -1.450718 -5.061982 -24.79721 -14.01974 -26.631005
[2,]  -2.535202  -1.563363 19.450718 11.061982  18.79721  20.01974   2.631005
         [,8]       [,9]     [,10]     [,11]     [,12]
[1,] -9.19071 -26.915872 -26.35097  2.149031 -30.75161
[2,] 19.19071   7.915872  12.35097 40.850969  33.75161

$out
[1] -38 -36  74  42

$group
[1] 1 4 6 7

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

Warning message:
In bxp(list(stats = c(-23, -17, -10, -4.5, 8, -23, -19.5, -15, 3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6vlu81288699926.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]
[1,] 19.0  36.0 30.0  42.0  23.0  9.0    3
[2,] 32.5  53.0 53.5  64.0  36.0 13.0    7
[3,] 41.5  64.5 70.0  72.5  45.5 16.5    9
[4,] 49.5  80.5 81.0  91.0  76.5 21.0   12
[5,] 52.0 107.0 96.0 105.0 121.0 28.0   14

$n
[1] 12 12 12 12 12 12  8

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]      [,7]
[1,] 33.74619 51.95707 57.45707 60.18512 27.02768 12.85115  6.206928
[2,] 49.25381 77.04293 82.54293 84.81488 63.97232 20.14885 11.793072

$out
[1]  82 136 137  25

$group
[1] 1 1 2 7

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

Warning message:
In bxp(list(stats = c(19, 32.5, 41.5, 49.5, 52, 36, 53, 64.5, 80.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7vlu81288699926.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,] 39.83333 36.5 39.125
[2,] 44.41667 41.5 42.125
[3,] 50.28571 47.5 49.250
[4,] 59.84524 63.0 54.625
[5,] 72.00000 69.0 73.125

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 43.24864 37.69371 43.54867
[2,] 57.32279 57.30629 54.95133

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(39.8333333333333, 44.4166666666667, 50.2857142857143,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1slek1288699926.ps tmp/1slek1288699926.png",intern=TRUE))
character(0)
> try(system("convert tmp/2slek1288699926.ps tmp/2slek1288699926.png",intern=TRUE))
character(0)
> try(system("convert tmp/3slek1288699926.ps tmp/3slek1288699926.png",intern=TRUE))
character(0)
> try(system("convert tmp/43uvn1288699926.ps tmp/43uvn1288699926.png",intern=TRUE))
character(0)
> try(system("convert tmp/53uvn1288699926.ps tmp/53uvn1288699926.png",intern=TRUE))
character(0)
> try(system("convert tmp/6vlu81288699926.ps tmp/6vlu81288699926.png",intern=TRUE))
character(0)
> try(system("convert tmp/7vlu81288699926.ps tmp/7vlu81288699926.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.931   1.470   2.144