R version 2.7.2 (2008-08-25)
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.31,100.61,100.84,100.64,100.36,100.37,100.75,101.45,101.63,101.63,101.75,101.76,102.14,102.35,102.87,103.55,103.92,104.13,104.58,104.39,104.29,104.93,105.17,105.52,105.57,105.11,104.97,105.09,105.33,105.60,105.14,105.56,105.09,104.87,104.81,104.60,104.07,104.41,104.02,103.62,103.68,103.48,103.27,103.58,103.41,103.31,103.13,103.36,103.14,103.73,103.80,103.46,103.60,103.67,103.68,104.04,104.79,105.26,105.62,105.82,105.71,105.74,106.26,106.59,106.36,106.33,106.37,106.17,106.16,106.37,106.18,106.03,106.08,106.24,105.97,106.28,106.52,106.72,106.80,107.06,107.02,107.04,107.09,107.15,107.50,107.90,108.02,107.78,107.70)
> 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 101.76 105.52 104.60 103.36 105.82 106.03 107.15
 [2,] 100.31 102.14 105.57 104.07 103.14 105.71 106.08 107.50
 [3,] 100.61 102.35 105.11 104.41 103.73 105.74 106.24 107.90
 [4,] 100.84 102.87 104.97 104.02 103.80 106.26 105.97 108.02
 [5,] 100.64 103.55 105.09 103.62 103.46 106.59 106.28 107.78
 [6,] 100.36 103.92 105.33 103.68 103.60 106.36 106.52 107.70
 [7,] 100.37 104.13 105.60 103.48 103.67 106.33 106.72     NA
 [8,] 100.75 104.58 105.14 103.27 103.68 106.37 106.80     NA
 [9,] 101.45 104.39 105.56 103.58 104.04 106.17 107.06     NA
[10,] 101.63 104.29 105.09 103.41 104.79 106.16 107.02     NA
[11,] 101.63 104.93 104.87 103.31 105.26 106.37 107.04     NA
[12,] 101.75 105.17 104.81 103.13 105.62 106.18 107.09     NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,]  0.31  0.38  0.05 -0.53 -0.22 -0.11  0.05  0.35
 [2,]  0.30  0.21 -0.46  0.34  0.59  0.03  0.16  0.40
 [3,]  0.23  0.52 -0.14 -0.39  0.07  0.52 -0.27  0.12
 [4,] -0.20  0.68  0.12 -0.40 -0.34  0.33  0.31 -0.24
 [5,] -0.28  0.37  0.24  0.06  0.14 -0.23  0.24 -0.08
 [6,]  0.01  0.21  0.27 -0.20  0.07 -0.03  0.20    NA
 [7,]  0.38  0.45 -0.46 -0.21  0.01  0.04  0.08    NA
 [8,]  0.70 -0.19  0.42  0.31  0.36 -0.20  0.26    NA
 [9,]  0.18 -0.10 -0.47 -0.17  0.75 -0.01 -0.04    NA
[10,]  0.00  0.64 -0.22 -0.10  0.47  0.21  0.02    NA
[11,]  0.12  0.24 -0.06 -0.18  0.36 -0.19  0.05    NA
[12,]  0.01  0.35 -0.21  0.23  0.20 -0.15  0.06    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/17or41225460276.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/2gyes1225460276.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/3yqzl1225460276.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/4mo9l1225460276.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.310 100.61 100.840 100.640 100.360 100.370 100.750 101.450
[2,] 102.560 102.640 103.04 103.335 103.505 103.640 103.575 103.475 103.810
[3,] 105.060 104.820 104.76 104.495 104.355 104.625 104.130 104.580 104.390
[4,] 105.925 105.895 105.99 106.115 106.435 106.440 105.965 105.755 105.865
[5,] 107.150 107.500 107.90 108.020 107.780 107.700 106.720 106.800 107.060
       [,10]   [,11]  [,12]
[1,] 101.630 101.630 101.75
[2,] 103.850 104.090 103.97
[3,] 104.790 104.930 105.17
[4,] 105.625 105.815 105.90
[5,] 107.020 107.040 107.09

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 103.1803 103.0017 103.1121 102.9421 102.7183 103.0609 102.7027 103.2184
[2,] 106.9397 106.6383 106.4079 106.0479 105.9917 106.1891 105.5573 105.9416
         [,9]  [,10]    [,11]    [,12]
[1,] 103.1628 103.73 103.8999 104.0174
[2,] 105.6172 105.85 105.9601 106.3226

$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, 102.56, 105.06, 105.925, 107.15, 100.31,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5kqnc1225460276.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.530 0.030 -0.390 -0.40 -0.280 -0.200 -0.46 -0.200 -0.170 -0.22 -0.19
[2,] -0.165 0.095 -0.205 -0.29 -0.155 -0.010 -0.10  0.035 -0.135 -0.05 -0.12
[3,]  0.050 0.255  0.095 -0.04  0.100  0.070  0.04  0.310 -0.040  0.02  0.05
[4,]  0.330 0.370  0.375  0.32  0.240  0.205  0.23  0.390  0.085  0.34  0.18
[5,]  0.380 0.590  0.520  0.68  0.370  0.270  0.45  0.700  0.180  0.64  0.36
      [,12]
[1,] -0.210
[2,] -0.070
[3,]  0.060
[4,]  0.215
[5,]  0.350

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

$conf
           [,1]      [,2]       [,3]       [,4]       [,5]        [,6]
[1,] -0.2265141 0.1013811 -0.2289963 -0.3807548 -0.1206527 -0.05839453
[2,]  0.3265141 0.4086189  0.4189963  0.3007548  0.3206527  0.19839453
           [,7]       [,8]        [,9]      [,10]      [,11]      [,12]
[1,] -0.1570707 0.09799973 -0.17138045 -0.2129017 -0.1291552 -0.1101974
[2,]  0.2370707 0.52200027  0.09138045  0.2529017  0.2291552  0.2301974

$out
[1] -0.46 -0.47  0.75

$group
[1] 2 9 9

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

Warning message:
In bxp(list(stats = c(-0.530000000000001, -0.164999999999999, 0.0499999999999972,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6hh2j1225460276.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,] 100.000 101.760 104.810 103.130 103.140 105.710 105.97 107.15
[2,] 100.365 102.610 105.030 103.360 103.530 105.990 106.16 107.50
[3,] 100.695 104.025 105.125 103.600 103.705 106.220 106.62 107.74
[4,] 101.540 104.485 105.540 104.045 104.415 106.365 107.03 107.90
[5,] 101.750 105.170 105.600 104.600 105.620 106.590 107.09 108.02

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 100.1591 103.1698 104.8924 103.2876 103.3013 106.0490 106.2232 107.4820
[2,] 101.2309 104.8802 105.3576 103.9124 104.1087 106.3910 107.0168 107.9980

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(100, 100.365, 100.695, 101.54, 101.75, 101.76,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7i2jf1225460276.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.2800 104.1300 104.3462
[2,] 104.3493 104.4425 104.6200
[3,] 104.6004 104.6925 104.7875
[4,] 104.6554 104.8750 104.9387
[5,] 104.8214 105.1700 105.0300

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 104.4608 104.4952 104.6421
[2,] 104.7401 104.8898 104.9329

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(104.28, 104.349285714286, 104.600446428571,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/17or41225460276.ps tmp/17or41225460276.png")
> system("convert tmp/2gyes1225460276.ps tmp/2gyes1225460276.png")
> system("convert tmp/3yqzl1225460276.ps tmp/3yqzl1225460276.png")
> system("convert tmp/4mo9l1225460276.ps tmp/4mo9l1225460276.png")
> system("convert tmp/5kqnc1225460276.ps tmp/5kqnc1225460276.png")
> system("convert tmp/6hh2j1225460276.ps tmp/6hh2j1225460276.png")
> system("convert tmp/7i2jf1225460276.ps tmp/7i2jf1225460276.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.428   0.969   1.727