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(8.5968,8.5114,8.3884,8.2671,8.2410,8.3177,8.4070,8.3917,8.4145,8.5245,8.6289,8.6622,8.9055,8.9770,9.1264,9.1120,9.0576,9.2106,9.2637,9.3107,9.6744,9.5780,9.4166,9.4359,9.2275,9.1828,9.0594,9.1358,9.2208,9.1137,9.2689,9.2489,9.1679,9.1051,9.0818,9.0961,9.1733,9.1455,9.2265,9.1541,9.1559,9.1182,9.1856,9.2378,9.0682,9.0105,8.9939,9.0228,9.1368,9.1763,9.2346,9.1653,9.1277,9.1430,9.1962,9.1861,9.0920,9.0620,8.9981,8.9819,9.0476,9.0852,9.0884,9.1670,9.1931,9.2628,9.4276,9.3398,9.3342,9.4223,9.5614,9.4316,9.3111,9.3414,9.4017,9.3346,9.3310,9.2349,9.2170,9.2098,9.2665,9.2533,9.1008,9.0377,9.0795,9.1896,9.2992,9.2372,9.2061,9.3290,9.1842,9.3231,9.2835,9.1735,9.2889,9.4319)
> 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,] 8.5968 8.9055 9.2275 9.1733 9.1368 9.0476 9.3111 9.0795   NA
 [2,] 8.5114 8.9770 9.1828 9.1455 9.1763 9.0852 9.3414 9.1896   NA
 [3,] 8.3884 9.1264 9.0594 9.2265 9.2346 9.0884 9.4017 9.2992   NA
 [4,] 8.2671 9.1120 9.1358 9.1541 9.1653 9.1670 9.3346 9.2372   NA
 [5,] 8.2410 9.0576 9.2208 9.1559 9.1277 9.1931 9.3310 9.2061   NA
 [6,] 8.3177 9.2106 9.1137 9.1182 9.1430 9.2628 9.2349 9.3290   NA
 [7,] 8.4070 9.2637 9.2689 9.1856 9.1962 9.4276 9.2170 9.1842   NA
 [8,] 8.3917 9.3107 9.2489 9.2378 9.1861 9.3398 9.2098 9.3231   NA
 [9,] 8.4145 9.6744 9.1679 9.0682 9.0920 9.3342 9.2665 9.2835   NA
[10,] 8.5245 9.5780 9.1051 9.0105 9.0620 9.4223 9.2533 9.1735   NA
[11,] 8.6289 9.4166 9.0818 8.9939 8.9981 9.5614 9.1008 9.2889   NA
[12,] 8.6622 9.4359 9.0961 9.0228 8.9819 9.4316 9.0377 9.4319   NA
> darr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8] [,9]
 [1,] -0.0854  0.0715 -0.0447 -0.0278  0.0395  0.0376  0.0303  0.1101   NA
 [2,] -0.1230  0.1494 -0.1234  0.0810  0.0583  0.0032  0.0603  0.1096   NA
 [3,] -0.1213 -0.0144  0.0764 -0.0724 -0.0693  0.0786 -0.0671 -0.0620   NA
 [4,] -0.0261 -0.0544  0.0850  0.0018 -0.0376  0.0261 -0.0036 -0.0311   NA
 [5,]  0.0767  0.1530 -0.1071 -0.0377  0.0153  0.0697 -0.0961  0.1229   NA
 [6,]  0.0893  0.0531  0.1552  0.0674  0.0532  0.1648 -0.0179 -0.1448   NA
 [7,] -0.0153  0.0470 -0.0200  0.0522 -0.0101 -0.0878 -0.0072  0.1389   NA
 [8,]  0.0228  0.3637 -0.0810 -0.1696 -0.0941 -0.0056  0.0567 -0.0396   NA
 [9,]  0.1100 -0.0964 -0.0628 -0.0577 -0.0300  0.0881 -0.0132 -0.1100   NA
[10,]  0.1044 -0.1614 -0.0233 -0.0166 -0.0639  0.1391 -0.1525  0.1154   NA
[11,]  0.0333  0.0193  0.0143  0.0289 -0.0162 -0.1298 -0.0631  0.1430   NA
[12,]  0.2433 -0.2084  0.0772  0.1140  0.0657 -0.1205  0.0418      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/16r5t1244024113.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/2pi6u1244024114.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/31oxe1244024114.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/45pnk1244024114.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,] 8.90550 8.9770 9.05940 9.1120 9.05760 9.11370 9.1842 9.18610 9.06820
[2,] 8.97655 9.0311 9.07390 9.1239 9.09265 9.11595 9.1849 9.19795 9.08010
[3,] 9.10815 9.1609 9.17645 9.1597 9.17450 9.17680 9.2066 9.24335 9.21720
[4,] 9.20040 9.1862 9.26690 9.2021 9.21345 9.24885 9.2663 9.31690 9.30885
[5,] 9.31110 9.3414 9.40170 9.2372 9.33100 9.32900 9.2689 9.33980 9.33420
       [,10]   [,11]   [,12]
[1,] 9.01050 8.62890 8.66220
[2,] 9.03625 8.99600 9.00235
[3,] 9.13930 9.09130 9.06690
[4,] 9.33780 9.35275 9.43175
[5,] 9.57800 9.56140 9.43590

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]    [,6]     [,7]     [,8]
[1,] 8.983104 9.074259 9.068637 9.116016 9.10702 9.10256 9.161129 9.176903
[2,] 9.233196 9.247541 9.284263 9.203384 9.24198 9.25104 9.252071 9.309797
         [,9]   [,10]    [,11]    [,12]
[1,] 9.089417 8.97085 8.892014 8.827031
[2,] 9.344983 9.30775 9.290586 9.306769

$out
 [1] 8.5968 8.5114 8.3884 8.2671 9.3346 8.2410 8.3177 8.4070 9.4276 8.3917
[11] 8.4145 9.6744 8.5245

$group
 [1]  1  2  3  4  4  5  6  7  7  8  9  9 10

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

Warning message:
In bxp(list(stats = c(8.9055, 8.97655, 9.10815, 9.2004, 9.3111,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/50lg91244024114.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]
[1,] -0.08540 -0.1234 -0.12130 -0.05440 -0.1071 -0.01790 -0.08780 -0.16960
[2,] -0.03625 -0.0599 -0.07085 -0.03435 -0.0669  0.01760 -0.01765 -0.08755
[3,]  0.03395  0.0593 -0.06455 -0.01485  0.0425  0.06030 -0.00865 -0.02260
[4,]  0.05550  0.0953  0.03100  0.01395  0.0998  0.12225  0.04960  0.03975
[5,]  0.11010  0.1494  0.07860  0.08500  0.1530  0.16480  0.13890  0.05670
         [,9]    [,10]    [,11]    [,12]
[1,] -0.11000 -0.16140 -0.12980 -0.20840
[2,] -0.07960 -0.10820 -0.03965 -0.03935
[3,] -0.04385 -0.01995  0.01680  0.06570
[4,]  0.03745  0.10990  0.03110  0.09560
[5,]  0.11000  0.13910  0.03330  0.24330

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

$conf
            [,1]        [,2]         [,3]        [,4]        [,5]        [,6]
[1,] -0.01730287 -0.02739695 -0.121444872 -0.04183107 -0.05062101 0.001841008
[2,]  0.08520287  0.14599695 -0.007655128  0.01213107  0.13562101 0.118758992
            [,7]        [,8]        [,9]      [,10]       [,11]       [,12]
[1,] -0.04621682 -0.09371161 -0.10923581 -0.1417838 -0.02272197 -0.01488996
[2,]  0.02891682  0.04851161  0.02153581  0.1018838  0.05632197  0.14628996

$out
[1] -0.1448  0.3637  0.1430

$group
[1]  6  8 11

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

Warning message:
In bxp(list(stats = c(-0.0854, -0.036249999999999, 0.0339500000000008,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6aauj1244024114.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,] 8.24100 8.90550 9.05940 8.99390 8.9819 9.04760 9.0377 9.07950   NA
[2,] 8.35305 9.08480 9.10060 9.04550 9.0770 9.12770 9.2134 9.18690   NA
[3,] 8.41075 9.23715 9.15185 9.14980 9.1399 9.29850 9.2599 9.26035   NA
[4,] 8.56065 9.42625 9.22415 9.17945 9.1812 9.42495 9.3328 9.31115   NA
[5,] 8.66220 9.67440 9.26890 9.23780 9.2346 9.56140 9.4017 9.43190   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]    [,7]     [,8]
[1,] 8.316062 9.081412 9.095498 9.088705 9.092374 9.162922 9.20544 9.203679
[2,] 8.505438 9.392888 9.208202 9.210895 9.187426 9.434078 9.31436 9.317021
     [,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(8.241, 8.35305, 8.41075, 8.56065, 8.6622,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7qk0b1244024114.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,] 9.059762 9.066900 9.099462
[2,] 9.073894 9.123725 9.158069
[3,] 9.118438 9.167700 9.169144
[4,] 9.142463 9.191700 9.206619
[5,] 9.162650 9.243350 9.258838

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 9.087163 9.136696 9.147000
[2,] 9.149712 9.198704 9.191288

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(9.0597625, 9.07389375, 9.1184375, 9.1424625,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/16r5t1244024113.ps tmp/16r5t1244024113.png")
> system("convert tmp/2pi6u1244024114.ps tmp/2pi6u1244024114.png")
> system("convert tmp/31oxe1244024114.ps tmp/31oxe1244024114.png")
> system("convert tmp/45pnk1244024114.ps tmp/45pnk1244024114.png")
> system("convert tmp/50lg91244024114.ps tmp/50lg91244024114.png")
> system("convert tmp/6aauj1244024114.ps tmp/6aauj1244024114.png")
> system("convert tmp/7qk0b1244024114.ps tmp/7qk0b1244024114.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.253   0.998   1.870