R version 2.8.0 (2008-10-20)
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(0.6348,0.634,0.62915,0.62168,0.1328,0.6089,0.60857,0.62672,0.62291,0.62393,0.61838,0.62012,0.61659,0.6116,0.61573,0.61407,0.62823,0.64405,0.6387,0.63633,0.63059,0.62994,0.63709,0.64217,0.65711,0.66977,0.68255,0.68902,0.71322,0.70224,0.70045,0.69919,0.69693,0.69763,0.69278,0.70196,0.69215,0.6769,0.67124,0.66532,0.67157,0.66428,0.66576,0.66942,0.6813,0.69144,0.69862,0.695,0.69867,0.68968,0.69233,0.68293,0.68399,0.66895,0.68756,0.68527,0.6776,0.68137,0.67933,0.67922)
> 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] 60
> (np <- floor(n / par1))
[1] 5
> 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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
         [,1]    [,2]    [,3]    [,4]    [,5] [,6]
 [1,] 0.63480 0.61659 0.65711 0.69215 0.69867   NA
 [2,] 0.63400 0.61160 0.66977 0.67690 0.68968   NA
 [3,] 0.62915 0.61573 0.68255 0.67124 0.69233   NA
 [4,] 0.62168 0.61407 0.68902 0.66532 0.68293   NA
 [5,] 0.13280 0.62823 0.71322 0.67157 0.68399   NA
 [6,] 0.60890 0.64405 0.70224 0.66428 0.66895   NA
 [7,] 0.60857 0.63870 0.70045 0.66576 0.68756   NA
 [8,] 0.62672 0.63633 0.69919 0.66942 0.68527   NA
 [9,] 0.62291 0.63059 0.69693 0.68130 0.67760   NA
[10,] 0.62393 0.62994 0.69763 0.69144 0.68137   NA
[11,] 0.61838 0.63709 0.69278 0.69862 0.67933   NA
[12,] 0.62012 0.64217 0.70196 0.69500 0.67922   NA
> darr
          [,1]     [,2]     [,3]     [,4]     [,5] [,6]
 [1,] -0.00080 -0.00499  0.01266 -0.01525 -0.00899   NA
 [2,] -0.00485  0.00413  0.01278 -0.00566  0.00265   NA
 [3,] -0.00747 -0.00166  0.00647 -0.00592 -0.00940   NA
 [4,] -0.48888  0.01416  0.02420  0.00625  0.00106   NA
 [5,]  0.47610  0.01582 -0.01098 -0.00729 -0.01504   NA
 [6,] -0.00033 -0.00535 -0.00179  0.00148  0.01861   NA
 [7,]  0.01815 -0.00237 -0.00126  0.00366 -0.00229   NA
 [8,] -0.00381 -0.00574 -0.00226  0.01188 -0.00767   NA
 [9,]  0.00102 -0.00065  0.00070  0.01014  0.00377   NA
[10,] -0.00555  0.00715 -0.00485  0.00718 -0.00204   NA
[11,]  0.00174  0.00508  0.00918 -0.00362 -0.00011   NA
[12,] -0.00353  0.01494 -0.00981  0.00367       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/freestat/rcomp/tmp/1xp191257771958.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/freestat/rcomp/tmp/2xeo61257771958.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/freestat/rcomp/tmp/3zu931257771958.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/freestat/rcomp/tmp/4rqbe1257771958.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,] 0.61659 0.61160 0.61573 0.61407 0.62823 0.60890 0.60857 0.62672 0.62291
[2,] 0.63480 0.63400 0.62915 0.62168 0.62823 0.64405 0.63870 0.63633 0.63059
[3,] 0.65711 0.66977 0.67124 0.66532 0.67157 0.66428 0.66576 0.66942 0.67760
[4,] 0.69215 0.67690 0.68255 0.68293 0.68399 0.66895 0.68756 0.68527 0.68130
[5,] 0.69867 0.68968 0.69233 0.68902 0.71322 0.70224 0.70045 0.69919 0.69693
       [,10]   [,11]   [,12]
[1,] 0.62393 0.61838 0.62012
[2,] 0.62994 0.63709 0.64217
[3,] 0.68137 0.67933 0.67922
[4,] 0.69144 0.69278 0.69500
[5,] 0.69763 0.69862 0.70196

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 5

$conf
          [,1]     [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.6165866 0.639457 0.6335077 0.6220409 0.6321701 0.6466857 0.6312356
[2,] 0.6976334 0.700083 0.7089723 0.7085991 0.7109699 0.6818743 0.7002844
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,] 0.6348391 0.6417684 0.6379143 0.6399796 0.6418905
[2,] 0.7040009 0.7134316 0.7248257 0.7186804 0.7165495

$out
[1] 0.1328

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(0.61659, 0.6348, 0.65711, 0.69215, 0.69867,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5rdi91257771958.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.01525 -0.00566 -0.00940 0.00106 -0.01504 -0.00535 -0.00237 -0.00767
[2,] -0.00899 -0.00485 -0.00747 0.00106 -0.01098 -0.00179 -0.00229 -0.00574
[3,] -0.00499  0.00265 -0.00592 0.00625 -0.00729 -0.00033 -0.00126 -0.00381
[4,] -0.00080  0.00413 -0.00166 0.01416  0.01582  0.00148  0.00366 -0.00226
[5,] -0.00080  0.01278  0.00647 0.02420  0.01582  0.00148  0.00366 -0.00226
         [,9]    [,10]    [,11]     [,12]
[1,] -0.00065 -0.00555 -0.00362 -0.009810
[2,]  0.00070 -0.00485 -0.00011 -0.006670
[3,]  0.00102 -0.00204  0.00174  0.000070
[4,]  0.00377  0.00715  0.00508  0.009305
[5,]  0.00377  0.00718  0.00918  0.014940

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 4

$conf
              [,1]         [,2]         [,3]         [,4]        [,5]
[1,] -0.0107770334 -0.003695245 -0.010025331 -0.003006427 -0.02622681
[2,]  0.0007970334  0.008995245 -0.001814669  0.015506427  0.01164681
             [,6]         [,7]         [,8]         [,9]       [,10]
[1,] -0.002640574 -0.005464255 -0.006268959 -0.001149254 -0.01051917
[2,]  0.001980574  0.002944255 -0.001351041  0.003189254  0.00643917
            [,11]       [,12]
[1,] -0.001927241 -0.01255025
[2,]  0.005407241  0.01269025

$out
[1]  0.01266 -0.48888  0.47610  0.01861  0.01815  0.01188  0.01014

$group
[1] 1 4 5 6 7 8 9

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

Warning message:
In bxp(list(stats = c(-0.0152500000000001, -0.00899000000000005,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6ov0b1257771958.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]
[1,] 0.608570 0.611600 0.669770 0.664280 0.668950   NA
[2,] 0.613640 0.616160 0.685785 0.667590 0.679275   NA
[3,] 0.622295 0.630265 0.697280 0.674235 0.683460   NA
[4,] 0.627935 0.637895 0.701205 0.691795 0.688620   NA
[5,] 0.634800 0.644050 0.713220 0.698620 0.698670   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]      [,2]      [,3]     [,4]      [,5] [,6]
[1,] 0.615775 0.6203515 0.6902468 0.663195 0.6791977   NA
[2,] 0.628815 0.6401785 0.7043132 0.685275 0.6877223   NA

$out
[1] 0.13280 0.65711

$group
[1] 1 3

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

Warning message:
In bxp(list(stats = c(0.60857, 0.61364, 0.622295, 0.627935, 0.6348,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7c4im1257771958.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,] 0.654604 0.657110 0.6523050
[2,] 0.657037 0.665540 0.6558975
[3,] 0.660036 0.670505 0.6585950
[4,] 0.664124 0.678410 0.6633025
[5,] 0.667694 0.681370 0.6685850

$n
[1] 12 12 12

$conf
          [,1]      [,2]      [,3]
[1,] 0.6568036 0.6646349 0.6552175
[2,] 0.6632684 0.6763751 0.6619725

$out
[1] 0.565962

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(0.654604, 0.657037, 0.660036, 0.664124, 0.667694,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1xp191257771958.ps tmp/1xp191257771958.png")
> system("convert tmp/2xeo61257771958.ps tmp/2xeo61257771958.png")
> system("convert tmp/3zu931257771958.ps tmp/3zu931257771958.png")
> system("convert tmp/4rqbe1257771958.ps tmp/4rqbe1257771958.png")
> system("convert tmp/5rdi91257771958.ps tmp/5rdi91257771958.png")
> system("convert tmp/6ov0b1257771958.ps tmp/6ov0b1257771958.png")
> system("convert tmp/7c4im1257771958.ps tmp/7c4im1257771958.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.856   1.457   2.223