R version 2.10.1 (2009-12-14)
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(3893.9,3799.2,3769.6,3768.6,3854.9,3778.5,3779.7,3803.2,3900.3,3792.6,3767.4,3752.6,3829.6,3722.6,3692.9,3681,3762.9,3661.7,3633.1,3621.5,3710,3619.4,3595.2,3573.2,3650.1,3554.2,3537,3528.6,3597.1,3521.9,3516.5,3515.7,3600.2,3517.1,3513.7,3528.2,3608.3,3502.5,3502.5,3495.3,3543.8,3425.3,3418.4,3406.4,3446.1,3341.1,3347,3354.9,3399,3288.9,3279,3275.2,3314,3227.1,3225.3,3228.6,3287.1,3210.1,3213.1,3228,3287,3211,3199.8,3166.3,3164,3156.7,3156,3165.5,3179.2,3182.5,3179.5,3193.5,3219.6,3221.9,3210.1)
> 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] 75
> (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 6 6 6 6 6 6 6 6 6
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]
 [1,] 3893.9 3829.6 3650.1 3608.3 3399.0 3287.0 3219.6
 [2,] 3799.2 3722.6 3554.2 3502.5 3288.9 3211.0 3221.9
 [3,] 3769.6 3692.9 3537.0 3502.5 3279.0 3199.8 3210.1
 [4,] 3768.6 3681.0 3528.6 3495.3 3275.2 3166.3     NA
 [5,] 3854.9 3762.9 3597.1 3543.8 3314.0 3164.0     NA
 [6,] 3778.5 3661.7 3521.9 3425.3 3227.1 3156.7     NA
 [7,] 3779.7 3633.1 3516.5 3418.4 3225.3 3156.0     NA
 [8,] 3803.2 3621.5 3515.7 3406.4 3228.6 3165.5     NA
 [9,] 3900.3 3710.0 3600.2 3446.1 3287.1 3179.2     NA
[10,] 3792.6 3619.4 3517.1 3341.1 3210.1 3182.5     NA
[11,] 3767.4 3595.2 3513.7 3347.0 3213.1 3179.5     NA
[12,] 3752.6 3573.2 3528.2 3354.9 3228.0 3193.5     NA
> darr
        [,1]   [,2]  [,3]   [,4]   [,5]  [,6]  [,7]
 [1,]  -94.7 -107.0 -95.9 -105.8 -110.1 -76.0   2.3
 [2,]  -29.6  -29.7 -17.2    0.0   -9.9 -11.2 -11.8
 [3,]   -1.0  -11.9  -8.4   -7.2   -3.8 -33.5    NA
 [4,]   86.3   81.9  68.5   48.5   38.8  -2.3    NA
 [5,]  -76.4 -101.2 -75.2 -118.5  -86.9  -7.3    NA
 [6,]    1.2  -28.6  -5.4   -6.9   -1.8  -0.7    NA
 [7,]   23.5  -11.6  -0.8  -12.0    3.3   9.5    NA
 [8,]   97.1   88.5  84.5   39.7   58.5  13.7    NA
 [9,] -107.7  -90.6 -83.1 -105.0  -77.0   3.3    NA
[10,]  -25.2  -24.2  -3.4    5.9    3.0  -3.0    NA
[11,]  -14.8  -22.0  14.5    7.9   14.9  14.0    NA
[12,]   77.0   76.9  80.1   44.1   59.0  26.1    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/yougetitorg/rcomp/tmp/1p4a81305754790.ps",horizontal=F,onefile=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/yougetitorg/rcomp/tmp/26hv81305754790.ps",horizontal=F,onefile=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/yougetitorg/rcomp/tmp/3ne331305754790.ps",horizontal=F,onefile=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/yougetitorg/rcomp/tmp/4ionq1305754790.ps",horizontal=F,onefile=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,] 3219.60 3211.0 3199.80 3166.30 3164.00 3156.7 3156.00 3165.50 3179.20
[2,] 3343.00 3255.4 3244.55 3275.20 3314.00 3227.1 3225.30 3228.60 3287.10
[3,] 3608.30 3502.5 3502.50 3511.95 3570.45 3473.6 3467.45 3461.05 3523.15
[4,] 3739.85 3638.4 3614.95 3681.00 3762.90 3661.7 3633.10 3621.50 3710.00
[5,] 3893.90 3799.2 3769.60 3768.60 3854.90 3778.5 3779.70 3803.20 3900.30
      [,10]   [,11]   [,12]
[1,] 3182.5 3179.50 3193.50
[2,] 3210.1 3213.10 3228.00
[3,] 3429.1 3430.35 3441.55
[4,] 3619.4 3595.20 3573.20
[5,] 3792.6 3767.40 3752.60

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 3371.308 3273.779 3281.303 3250.196 3280.895 3193.269 3204.406 3207.617
[2,] 3845.292 3731.221 3723.697 3773.704 3860.005 3753.931 3730.494 3714.483
         [,9]    [,10]    [,11]    [,12]
[1,] 3250.366 3165.088 3183.883 3218.885
[2,] 3795.934 3693.112 3676.817 3664.215

$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(3219.6, 3343, 3608.3, 3739.85, 3893.9, 3211,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/yougetitorg/rcomp/tmp/5oxg41305754790.ps",horizontal=F,onefile=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,] -110.10 -29.70 -11.9 -2.3 -118.50 -6.9 -12.00 13.7 -107.70 -25.2 -22.00
[2,] -106.40 -23.40 -11.9 38.8 -101.20 -6.9 -11.60 39.7 -105.00 -24.2 -14.80
[3,]  -95.90 -11.80  -7.8 58.5  -81.65 -3.6   1.25 71.5  -86.85  -3.2  10.95
[4,]  -85.35 -10.55  -3.8 81.9  -75.20 -0.7   9.50 88.5  -77.00   3.0  14.50
[5,]  -76.00   0.00  -1.0 86.3  -75.20  1.2  23.50 97.1  -77.00   5.9  14.90
     [,12]
[1,] 26.10
[2,] 44.10
[3,] 67.95
[4,] 77.00
[5,] 80.10

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

$conf
           [,1]       [,2]       [,3]     [,4]      [,5]       [,6]      [,7]
[1,] -108.47072 -19.473813 -13.024762 30.69911 -98.42084 -7.5992003 -12.36018
[2,]  -83.32928  -4.126187  -2.575238 86.30089 -64.87916  0.3992003  14.86018
          [,8]      [,9]     [,10]     [,11]    [,12]
[1,]  40.02242 -104.9109 -20.74488 -7.949446 46.72844
[2,] 102.97758  -68.7891  14.34488 29.849446 89.17156

$out
[1]   2.3 -33.5  -7.3 -28.6   3.3

$group
[1] 1 3 5 6 9

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

Warning message:
In bxp(list(stats = c(-110.1, -106.4, -95.9, -85.3500000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/yougetitorg/rcomp/tmp/6f8vd1305754790.ps",horizontal=F,onefile=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,] 3752.60 3573.20 3513.70 3341.10 3210.1 3156.00 3210.10
[2,] 3769.10 3620.45 3516.80 3380.65 3226.2 3164.75 3214.85
[3,] 3786.15 3671.35 3528.40 3435.70 3251.9 3179.35 3219.60
[4,] 3829.05 3716.30 3575.65 3502.50 3288.0 3196.65 3220.75
[5,] 3900.30 3829.60 3650.10 3608.30 3314.0 3211.00 3221.90

$n
[1] 12 12 12 12 12 12  3

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]   [,6]     [,7]
[1,] 3758.806 3627.632 3501.558 3380.123 3223.713 3164.8 3214.218
[2,] 3813.494 3715.068 3555.242 3491.277 3280.087 3193.9 3224.982

$out
[1] 3399 3287

$group
[1] 5 6

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

Warning message:
In bxp(list(stats = c(3752.6, 3769.1, 3786.15, 3829.05, 3900.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/yougetitorg/rcomp/tmp/726gw1305754790.ps",horizontal=F,onefile=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,] 3435.983 3429.10 3410.700
[2,] 3449.317 3451.30 3424.044
[3,] 3459.342 3488.05 3442.831
[4,] 3503.158 3517.55 3495.631
[5,] 3555.357 3608.30 3546.450

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 3434.784 3457.833 3410.180
[2,] 3483.899 3518.267 3475.483

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(3435.98333333333, 3449.31666666667, 3459.34166666667,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1p4a81305754790.ps tmp/1p4a81305754790.png",intern=TRUE))
character(0)
> try(system("convert tmp/26hv81305754790.ps tmp/26hv81305754790.png",intern=TRUE))
character(0)
> try(system("convert tmp/3ne331305754790.ps tmp/3ne331305754790.png",intern=TRUE))
character(0)
> try(system("convert tmp/4ionq1305754790.ps tmp/4ionq1305754790.png",intern=TRUE))
character(0)
> try(system("convert tmp/5oxg41305754790.ps tmp/5oxg41305754790.png",intern=TRUE))
character(0)
> try(system("convert tmp/6f8vd1305754790.ps tmp/6f8vd1305754790.png",intern=TRUE))
character(0)
> try(system("convert tmp/726gw1305754790.ps tmp/726gw1305754790.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.590   1.300   2.073