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(0.102099945804887,0.737601051471779,0.0106618293088252,0.52239377142817,1.07142787304473,-0.914061367174038,-0.428462409725616,-0.0654398048209609,1.19700131954875,-0.415035627082705,0.790229966140711,-1.20057857422849,1.10379820587161,-0.273135171804026,-0.267822467388643,1.34265393731741,-1.54293043754615,-0.0119730328107786,0.222952413225187,-0.345660667035868,0.264305650805457,1.72740775878021,1.58805033299532,0.377632331984582,0.0174696730264222,-0.844095886303813,0.141125552395322,0.849526022031545,0.652168611346395,-0.919947768004273,0.433795873152471,-1.43446391109305,-0.649941295444581,0.290571683853168,-0.430213956986155,0.231554627911427,1.29492092986709,-0.39281211697063,-0.587917646775748,-0.806911791305758,-0.392543947667208,-1.0340885974776,-0.75641012970889,-0.301789001672712,-0.547361456817899)
> 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] 45
> (np <- floor(n / par1))
[1] 3
> 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] 4 4 4 4 4 4 4 4 4 3 3 3
> arr
             [,1]        [,2]        [,3]       [,4]
 [1,]  0.10209995  1.10379821  0.01746967  1.2949209
 [2,]  0.73760105 -0.27313517 -0.84409589 -0.3928121
 [3,]  0.01066183 -0.26782247  0.14112555 -0.5879176
 [4,]  0.52239377  1.34265394  0.84952602 -0.8069118
 [5,]  1.07142787 -1.54293044  0.65216861 -0.3925439
 [6,] -0.91406137 -0.01197303 -0.91994777 -1.0340886
 [7,] -0.42846241  0.22295241  0.43379587 -0.7564101
 [8,] -0.06543980 -0.34566067 -1.43446391 -0.3017890
 [9,]  1.19700132  0.26430565 -0.64994130 -0.5473615
[10,] -0.41503563  1.72740776  0.29057168         NA
[11,]  0.79022997  1.58805033 -0.43021396         NA
[12,] -1.20057857  0.37763233  0.23155463         NA
> darr
            [,1]         [,2]       [,3]       [,4]
 [1,]  0.6355011 -1.376933378 -0.8615656 -1.6877330
 [2,] -0.7269392  0.005312704  0.9852214 -0.1951055
 [3,]  0.5117319  1.610476405  0.7084005 -0.2189941
 [4,]  0.5490341 -2.885584375 -0.1973574  0.4143678
 [5,] -1.9854892  1.530957405 -1.5721164 -0.6415446
 [6,]  0.4855990  0.234925446  1.3537436  0.2776785
 [7,]  0.3630226 -0.568613080 -1.8682598  0.4546211
 [8,]  1.2624411  0.609966318  0.7845226 -0.2455725
 [9,] -1.6120369  1.463102108  0.9405130         NA
[10,]  1.2052656 -0.139357426 -0.7207856         NA
[11,] -1.9908085 -1.210418001  0.6617686         NA
[12,]  2.3043768 -0.360162659  1.0633663         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/1kfif1260810974.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/2u2a71260810974.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/3x5k81260810974.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/4aywp1260810974.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]
[1,] 0.01746967 -0.8440959 -0.58791765 -0.8069118 -1.5429304 -1.03408860
[2,] 0.05978481 -0.6184540 -0.42787006 -0.1422590 -0.9677372 -0.97701818
[3,] 0.60294908 -0.3329736 -0.12858032  0.6859599  0.1298123 -0.91700457
[4,] 1.19935957  0.2322329  0.07589369  1.0960900  0.8617982 -0.46301720
[5,] 1.29492093  0.7376011  0.14112555  1.3426539  1.0714279 -0.01197303
           [,7]       [,8]       [,9]       [,10]     [,11]      [,12]
[1,] -0.7564101 -1.4344639 -0.6499413 -0.41503563 -0.430214 -1.2005786
[2,] -0.5924363 -0.8900623 -0.5986514 -0.06223197  0.180008 -0.4845120
[3,] -0.1027550 -0.3237248 -0.1415279  0.29057168  0.790230  0.2315546
[4,]  0.3283741 -0.1836144  0.7306535  1.00898972  1.189140  0.3045935
[5,]  0.4337959 -0.0654398  1.1970013  1.72740776  1.588050  0.3776323

$n
 [1] 4 4 4 4 4 4 4 4 4 3 3 3

$conf
          [,1]       [,2]       [,3]       [,4]      [,5]       [,6]       [,7]
[1,] -0.297315 -1.0050163 -0.5265537 -0.2923358 -1.315521 -1.3230653 -0.8301952
[2,]  1.503213  0.3390690  0.2693930  1.6642556  1.575145 -0.5109438  0.6246852
           [,8]      [,9]      [,10]      [,11]      [,12]
[1,] -0.8818187 -1.191679 -0.6866111 -0.1303139 -0.4882780
[2,]  0.2343690  0.908623  1.2677545  1.7107739  0.9513872

$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(0.0174696730264222, 0.0597848094156546, 0.602949075838249,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5j5ji1260810974.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]
[1,] -1.6877330 -0.72693922 -0.2189941 -2.8855844 -1.9854892 0.2349254
[2,] -1.5323332 -0.46102238  0.1463689 -1.5414709 -1.7788028 0.2563020
[3,] -1.1192495 -0.09489641  0.6100662  0.1085052 -1.1068305 0.3816387
[4,] -0.1130322  0.49526707  1.1594384  0.4817010  0.4447064 0.9196713
[5,]  0.6355011  0.98522144  1.6104764  0.5490341  1.5309574 1.3537436
           [,7]       [,8]      [,9]      [,10]      [,11]      [,12]
[1,] -1.8682598 -0.2455725 -1.612037 -0.7207856 -1.9908085 -0.3601627
[2,] -1.2184364  0.1821969 -0.335762 -0.4300715 -1.6006133  0.3516018
[3,] -0.1027952  0.6972445  0.940513 -0.1393574 -1.2104180  1.0633663
[4,]  0.4088219  1.0234819  1.201808  0.5329541 -0.2743247  1.6838715
[5,]  0.4546211  1.2624411  1.463102  1.2052656  0.6617686  2.3043768

$n
 [1] 4 4 4 4 4 4 4 4 3 3 3 3

$conf
            [,1]       [,2]       [,3]      [,4]       [,5]       [,6]
[1,] -2.24049725 -0.8503651 -0.1902587 -1.489801 -2.8634028 -0.1424231
[2,]  0.00199831  0.6605723  1.4103911  1.706811  0.6497417  0.9057005
          [,7]       [,8]       [,9]      [,10]         [,11]      [,12]
[1,] -1.388329 0.03262937 -0.4620786 -1.0178423 -2.4202762337 -0.1519480
[2,]  1.182739 1.36185957  2.3431045  0.7391275 -0.0005597683  2.2786806

$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(-1.68773304683772, -1.53233321225668, -1.11924946850294,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/66p441260810974.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]
[1,] -1.20057857 -1.5429304 -1.43446391 -1.0340886
[2,] -0.42174902 -0.2704788 -0.74701859 -0.7564101
[3,]  0.05638089  0.2436290  0.07929761 -0.5473615
[4,]  0.76391551  1.2232261  0.36218378 -0.3925439
[5,]  1.19700132  1.7274078  0.84952602 -0.3017890

$n
[1] 12 12 12  9

$conf
           [,1]       [,2]       [,3]       [,4]
[1,] -0.4844087 -0.4376598 -0.4266170 -0.7389976
[2,]  0.5971704  0.9249179  0.5852123 -0.3557253

$out
[1] 1.294921

$group
[1] 4

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

Warning message:
In bxp(list(stats = c(-1.20057857422849, -0.421749018404161, 0.0563808875568561,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/77zca1260810974.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.72001769 -0.91700457 -0.81851113
[2,] -0.19512053 -0.23262637 -0.20766317
[3,] -0.09250027  0.01352867 -0.06386134
[4,]  0.50561505  0.44676038  0.52740828
[5,]  0.64935545  0.79022997  0.68457408

$n
[1] 12 12 12

$conf
           [,1]       [,2]       [,3]
[1,] -0.4121105 -0.2963442 -0.3991324
[2,]  0.2271099  0.3234015  0.2714097

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(-0.720017691366672, -0.195120534506250, -0.0925002692348851,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> try(system("convert tmp/1kfif1260810974.ps tmp/1kfif1260810974.png",intern=TRUE))
character(0)
> try(system("convert tmp/2u2a71260810974.ps tmp/2u2a71260810974.png",intern=TRUE))
character(0)
> try(system("convert tmp/3x5k81260810974.ps tmp/3x5k81260810974.png",intern=TRUE))
character(0)
> try(system("convert tmp/4aywp1260810974.ps tmp/4aywp1260810974.png",intern=TRUE))
character(0)
> try(system("convert tmp/5j5ji1260810974.ps tmp/5j5ji1260810974.png",intern=TRUE))
character(0)
> try(system("convert tmp/66p441260810974.ps tmp/66p441260810974.png",intern=TRUE))
character(0)
> try(system("convert tmp/77zca1260810974.ps tmp/77zca1260810974.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.184   0.946   1.556