R version 2.7.0 (2008-04-22)
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.

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(100.00,100.35,100.38,100.52,100.34,99.97,99.88,100.26,100.93,100.88,100.73,100.74,100.25,100.29,100.57,101.24,101.69,101.89,102.07,102.43,102.51,102.71,103.22,103.79,103.99,103.83,103.55,103.24,103.77,104.37,104.61,104.21,104.77,104.33,104.14,104.37,104.20,103.58,103.51,103.39,103.11,103.28,102.83,102.56,102.62,102.66,102.72,102.92,103.26,103.02,103.33,103.57,103.61,103.85,104.22,104.15,104.52,105.27,105.60,105.99,106.23,106.40,106.25,106.74,106.96,106.74,106.59,106.65,106.56,106.69,106.66,106.40,105.99,105.99,106.38,106.41,106.75,106.90,107.29,107.24,107.56,107.45,107.35,107.63,107.88,108.21,108.78,108.94,108.66,108.38)
> 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] 90
> (np <- floor(n / par1))
[1] 7
> 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 7 7 7 7 7 7
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]
 [1,] 100.00 100.25 103.99 104.20 103.26 106.23 105.99 107.88
 [2,] 100.35 100.29 103.83 103.58 103.02 106.40 105.99 108.21
 [3,] 100.38 100.57 103.55 103.51 103.33 106.25 106.38 108.78
 [4,] 100.52 101.24 103.24 103.39 103.57 106.74 106.41 108.94
 [5,] 100.34 101.69 103.77 103.11 103.61 106.96 106.75 108.66
 [6,]  99.97 101.89 104.37 103.28 103.85 106.74 106.90 108.38
 [7,]  99.88 102.07 104.61 102.83 104.22 106.59 107.29     NA
 [8,] 100.26 102.43 104.21 102.56 104.15 106.65 107.24     NA
 [9,] 100.93 102.51 104.77 102.62 104.52 106.56 107.56     NA
[10,] 100.88 102.71 104.33 102.66 105.27 106.69 107.45     NA
[11,] 100.73 103.22 104.14 102.72 105.60 106.66 107.35     NA
[12,] 100.74 103.79 104.37 102.92 105.99 106.40 107.63     NA
> darr
       [,1] [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,]  0.35 0.04 -0.16 -0.62 -0.24  0.17  0.00  0.33
 [2,]  0.03 0.28 -0.28 -0.07  0.31 -0.15  0.39  0.57
 [3,]  0.14 0.67 -0.31 -0.12  0.24  0.49  0.03  0.16
 [4,] -0.18 0.45  0.53 -0.28  0.04  0.22  0.34 -0.28
 [5,] -0.37 0.20  0.60  0.17  0.24 -0.22  0.15 -0.28
 [6,] -0.09 0.18  0.24 -0.45  0.37 -0.15  0.39    NA
 [7,]  0.38 0.36 -0.40 -0.27 -0.07  0.06 -0.05    NA
 [8,]  0.67 0.08  0.56  0.06  0.37 -0.09  0.32    NA
 [9,] -0.05 0.20 -0.44  0.04  0.75  0.13 -0.11    NA
[10,] -0.15 0.51 -0.19  0.06  0.33 -0.03 -0.10    NA
[11,]  0.01 0.57  0.23  0.20  0.39 -0.26  0.28    NA
[12,] -0.49 0.20 -0.17  0.34  0.24 -0.41  0.25    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/1yb1g1225629509.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/2sd2p1225629510.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/3r3du1225629510.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/4brbl1225629510.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,] 100.000 100.290 100.380 100.520 100.340  99.970  99.88 100.260 100.930
[2,] 101.755 101.685 101.950 102.240 102.400 102.585 102.45 102.495 102.565
[3,] 104.095 103.705 103.530 103.480 103.690 104.110 104.22 104.150 104.520
[4,] 106.110 106.195 106.315 106.575 106.855 106.820 105.60 105.430 105.665
[5,] 107.880 108.210 108.780 108.940 108.660 108.380 107.29 107.240 107.560
       [,10]  [,11]   [,12]
[1,] 100.880 100.73 100.740
[2,] 102.685 102.97 103.355
[3,] 104.330 104.14 104.370
[4,] 105.980 106.13 106.195
[5,] 107.450 107.35 107.630

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 101.6622 101.1856 101.0916 101.0584 101.2014 101.7443 102.3389 102.3973
[2,] 106.5278 106.2244 105.9684 105.9016 106.1786 106.4757 106.1011 105.9027
         [,9]    [,10]    [,11]   [,12]
[1,] 102.6687 102.3623 102.2529 102.674
[2,] 106.3713 106.2977 106.0271 106.066

$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(100, 101.755, 104.095, 106.11, 107.88, 100.29,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5u23q1225629510.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.62 -0.280 -0.310 -0.280 -0.37 -0.450 -0.40 -0.090 -0.440 -0.190 0.010
[2,] -0.20 -0.110 -0.045 -0.230 -0.25 -0.120 -0.17  0.070 -0.080 -0.125 0.105
[3,]  0.02  0.155  0.150  0.130  0.16  0.180 -0.05  0.320  0.040 -0.030 0.230
[4,]  0.25  0.350  0.365  0.395  0.22  0.305  0.21  0.465  0.165  0.195 0.335
[5,]  0.35  0.570  0.670  0.530  0.60  0.390  0.38  0.670  0.200  0.510 0.570
      [,12]
[1,] -0.490
[2,] -0.290
[3,]  0.200
[4,]  0.245
[5,]  0.340

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

$conf
           [,1]       [,2]        [,3]       [,4]       [,5]        [,6]
[1,] -0.2313765 -0.1019626 -0.07903189 -0.2191340 -0.1025487 -0.07380314
[2,]  0.2713765  0.4119626  0.37903189  0.4791340  0.4225487  0.43380314
           [,7]       [,8]       [,9]      [,10]      [,11]      [,12]
[1,] -0.2769299 0.08411237 -0.1063100 -0.2210988 0.09264771 -0.1194934
[2,]  0.1769299 0.55588763  0.1863100  0.1610988 0.36735229  0.5194934

$out
[1]  0.75 -0.26

$group
[1]  9 11

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

Warning message:
In bxp(list(stats = c(-0.620000000000005, -0.200000000000003, 0.0200000000000031,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6wa4o1225629510.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]
[1,]  99.880 100.250 103.240 102.560 103.020 106.230 105.990 107.88
[2,] 100.130 100.905 103.800 102.690 103.450 106.400 106.395 108.21
[3,] 100.365 101.980 104.175 103.015 104.000 106.620 107.070 108.52
[4,] 100.735 102.610 104.370 103.450 104.895 106.715 107.400 108.78
[5,] 100.930 103.790 104.770 104.200 105.990 106.960 107.630 108.94

$n
[1] 12 12 12 12 12 12 12  6

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 100.0891 101.2023 103.9150 102.6684 103.3409 106.4763 106.6116 108.1523
[2,] 100.6409 102.7577 104.4350 103.3616 104.6591 106.7637 107.5284 108.8877

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(99.88, 100.13, 100.365, 100.735, 100.93, 100.25,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7c2bm1225629510.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,] 103.9271 103.4800 103.9625
[2,] 103.9669 103.6975 104.1687
[3,] 104.2331 104.1250 104.3969
[4,] 104.3535 104.2750 104.6956
[5,] 104.5486 104.5200 104.8563

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 104.0568 103.8616 104.1566
[2,] 104.4095 104.3884 104.6372

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(103.927142857143, 103.966875, 104.233125,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1yb1g1225629509.ps tmp/1yb1g1225629509.png")
> system("convert tmp/2sd2p1225629510.ps tmp/2sd2p1225629510.png")
> system("convert tmp/3r3du1225629510.ps tmp/3r3du1225629510.png")
> system("convert tmp/4brbl1225629510.ps tmp/4brbl1225629510.png")
> system("convert tmp/5u23q1225629510.ps tmp/5u23q1225629510.png")
> system("convert tmp/6wa4o1225629510.ps tmp/6wa4o1225629510.png")
> system("convert tmp/7c2bm1225629510.ps tmp/7c2bm1225629510.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.589   1.581   2.826