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(101.02,101.15,101.51,101.75,101.8,101.8,101.8,101.82,101.99,102.25,102.34,102.35,102.35,102.39,102.49,102.67,102.68,102.7,102.71,102.72,102.83,102.92,103.04,103.08,103.09,103.11,103.18,103.18,103.22,103.25,103.25,103.25,103.47,103.57,103.66,103.7,103.7,103.75,103.85,104.02,104.13,104.17,104.18,104.2,104.5,104.78,104.88,104.89,104.9,104.95,105.24,105.35,105.44,105.46,105.47,105.48,105.75,106.1,106.19,106.23,106.24,106.25,106.35,106.48,106.52,106.55,106.55,106.56,106.89,107.09,107.24,107.28,107.3,107.31,107.47,107.35,107.31,107.32,107.32,107.34,107.53,107.72,107.75,107.79,107.81,107.9,107.8,107.86,107.8,107.74,107.75,107.83,107.8,107.81,107.86,107.83)
> 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] 96
> (np <- floor(n / par1))
[1] 8
> 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] 8 8 8 8 8 8 8 8 8 8 8 8
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8] [,9]
 [1,] 101.02 102.35 103.09 103.70 104.90 106.24 107.30 107.81   NA
 [2,] 101.15 102.39 103.11 103.75 104.95 106.25 107.31 107.90   NA
 [3,] 101.51 102.49 103.18 103.85 105.24 106.35 107.47 107.80   NA
 [4,] 101.75 102.67 103.18 104.02 105.35 106.48 107.35 107.86   NA
 [5,] 101.80 102.68 103.22 104.13 105.44 106.52 107.31 107.80   NA
 [6,] 101.80 102.70 103.25 104.17 105.46 106.55 107.32 107.74   NA
 [7,] 101.80 102.71 103.25 104.18 105.47 106.55 107.32 107.75   NA
 [8,] 101.82 102.72 103.25 104.20 105.48 106.56 107.34 107.83   NA
 [9,] 101.99 102.83 103.47 104.50 105.75 106.89 107.53 107.80   NA
[10,] 102.25 102.92 103.57 104.78 106.10 107.09 107.72 107.81   NA
[11,] 102.34 103.04 103.66 104.88 106.19 107.24 107.75 107.86   NA
[12,] 102.35 103.08 103.70 104.89 106.23 107.28 107.79 107.83   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6]  [,7]  [,8] [,9]
 [1,] 0.13 0.04 0.02 0.05 0.05 0.01  0.01  0.09   NA
 [2,] 0.36 0.10 0.07 0.10 0.29 0.10  0.16 -0.10   NA
 [3,] 0.24 0.18 0.00 0.17 0.11 0.13 -0.12  0.06   NA
 [4,] 0.05 0.01 0.04 0.11 0.09 0.04 -0.04 -0.06   NA
 [5,] 0.00 0.02 0.03 0.04 0.02 0.03  0.01 -0.06   NA
 [6,] 0.00 0.01 0.00 0.01 0.01 0.00  0.00  0.01   NA
 [7,] 0.02 0.01 0.00 0.02 0.01 0.01  0.02  0.08   NA
 [8,] 0.17 0.11 0.22 0.30 0.27 0.33  0.19 -0.03   NA
 [9,] 0.26 0.09 0.10 0.28 0.35 0.20  0.19  0.01   NA
[10,] 0.09 0.12 0.09 0.10 0.09 0.15  0.03  0.05   NA
[11,] 0.01 0.04 0.04 0.01 0.04 0.04  0.04 -0.03   NA
[12,] 0.00 0.01 0.00 0.01 0.01 0.02  0.02    NA   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/1spwz1290108399.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/2spwz1290108399.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/3spwz1290108399.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/4qk2r1290108399.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]
[1,] 101.02 101.15 101.510 101.750 101.800 101.800 101.800 101.820 101.990
[2,] 102.72 102.75 102.835 102.925 102.950 102.975 102.980 102.985 103.150
[3,] 104.30 104.35 104.545 104.685 104.785 104.815 104.825 104.840 105.125
[4,] 106.77 106.78 106.910 106.915 106.915 106.935 106.935 106.950 107.210
[5,] 107.81 107.90 107.800 107.860 107.800 107.740 107.750 107.830 107.800
       [,10]   [,11]   [,12]
[1,] 102.250 102.340 102.350
[2,] 103.245 103.350 103.390
[3,] 105.440 105.535 105.560
[4,] 107.405 107.495 107.535
[5,] 107.810 107.860 107.830

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 102.0376 102.0988 102.2686 102.4561 102.5701 102.6029 102.6157 102.6251
[2,] 106.5624 106.6012 106.8214 106.9139 106.9999 107.0271 107.0343 107.0549
         [,9]    [,10]    [,11]    [,12]
[1,] 102.8570 103.1162 103.2195 103.2445
[2,] 107.3930 107.7638 107.8505 107.8755

$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(101.02, 102.72, 104.3, 106.77, 107.81, 101.15,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5qk2r1290108399.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,] 0.010 -0.100 -0.120 -0.060 0.000 0.000 0.000 -0.030 0.010  0.03 -0.03
[2,] 0.015  0.085  0.030 -0.015 0.005 0.000 0.010  0.140 0.095  0.07  0.01
[3,] 0.045  0.100  0.120  0.040 0.020 0.005 0.015  0.205 0.195  0.09  0.04
[4,] 0.070  0.225  0.175  0.070 0.030 0.010 0.020  0.285 0.270  0.11  0.04
[5,] 0.130  0.360  0.240  0.110 0.040 0.010 0.020  0.330 0.350  0.15  0.04
     [,12]
[1,] 0.000
[2,] 0.005
[3,] 0.010
[4,] 0.015
[5,] 0.020

$n
 [1] 8 8 8 8 8 8 8 8 8 8 8 7

$conf
           [,1]       [,2]       [,3]        [,4]        [,5]          [,6]
[1,] 0.01427621 0.02179399 0.03900092 -0.00748222 0.006034641 -0.0005861436
[2,] 0.07572379 0.17820601 0.20099908  0.08748222 0.033965359  0.0105861436
            [,7]      [,8]       [,9]      [,10]      [,11]       [,12]
[1,] 0.009413856 0.1240009 0.09724249 0.06765543 0.02324157 0.004028161
[2,] 0.020586144 0.2859991 0.29275751 0.11234457 0.05675843 0.015971839

$out
[1] -0.06  0.08

$group
[1] 5 7

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

Warning message:
In bxp(list(stats = c(0.0100000000000051, 0.0150000000000006, 0.0450000000000017,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6qk2r1290108399.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]
[1,] 101.02 102.350 103.09 103.700 104.900 106.240 107.300 107.740   NA
[2,] 101.63 102.580 103.18 103.935 105.295 106.415 107.315 107.800   NA
[3,] 101.80 102.705 103.25 104.175 105.465 106.550 107.345 107.810   NA
[4,] 102.12 102.875 103.52 104.640 105.925 106.990 107.625 107.845   NA
[5,] 102.35 103.080 103.70 104.890 106.230 107.280 107.790 107.900   NA

$n
[1] 12 12 12 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 101.5765 102.5704 103.0949 103.8534 105.1777 106.2877 107.2036 107.7895
[2,] 102.0235 102.8396 103.4051 104.4966 105.7523 106.8123 107.4864 107.8305
     [,9]
[1,]   NA
[2,]   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(101.02, 101.63, 101.8, 102.12, 102.35, 102.35,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7v8d61290108399.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,] 104.5512 104.3000 104.7050
[2,] 104.7844 104.6150 104.8469
[3,] 104.8762 104.8200 104.9281
[4,] 105.1875 105.2825 105.2537
[5,] 105.3937 105.5600 105.4763

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 104.6924 104.5155 104.7425
[2,] 105.0601 105.1245 105.1137

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(104.55125, 104.784375, 104.87625, 105.1875,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1spwz1290108399.ps tmp/1spwz1290108399.png",intern=TRUE))
character(0)
> try(system("convert tmp/2spwz1290108399.ps tmp/2spwz1290108399.png",intern=TRUE))
character(0)
> try(system("convert tmp/3spwz1290108399.ps tmp/3spwz1290108399.png",intern=TRUE))
character(0)
> try(system("convert tmp/4qk2r1290108399.ps tmp/4qk2r1290108399.png",intern=TRUE))
character(0)
> try(system("convert tmp/5qk2r1290108399.ps tmp/5qk2r1290108399.png",intern=TRUE))
character(0)
> try(system("convert tmp/6qk2r1290108399.ps tmp/6qk2r1290108399.png",intern=TRUE))
character(0)
> try(system("convert tmp/7v8d61290108399.ps tmp/7v8d61290108399.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.254   0.969   3.383