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(68897,38683,44720,39525,45315,50380,40600,36279,42438,38064,31879,11379,70249,39253,47060,41697,38708,49267,39018,32228,40870,39383,34571,12066,70938,34077,45409,40809,37013,44953,37848,32745,43412,34931,33008,8620,68906,39556,50669,36432,40891,48428,36222,33425,39401,37967,34801,12657,69116,41519,51321,38529,41547,52073,38401,40898,40439,41888,37898,8771,68184,50530,47221,41756,45633,48138,39486,39341,41117,41629,29722,7054,56676,34870,35117,30169,30936,35699,33228,27733,33666,35429,27438,8170,62557)
> par1 = '52'
> #'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] 85
> (np <- floor(n / par1))
[1] 1
> 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] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1
[39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1
> arr
       [,1]  [,2]
 [1,] 68897 41547
 [2,] 38683 52073
 [3,] 44720 38401
 [4,] 39525 40898
 [5,] 45315 40439
 [6,] 50380 41888
 [7,] 40600 37898
 [8,] 36279  8771
 [9,] 42438 68184
[10,] 38064 50530
[11,] 31879 47221
[12,] 11379 41756
[13,] 70249 45633
[14,] 39253 48138
[15,] 47060 39486
[16,] 41697 39341
[17,] 38708 41117
[18,] 49267 41629
[19,] 39018 29722
[20,] 32228  7054
[21,] 40870 56676
[22,] 39383 34870
[23,] 34571 35117
[24,] 12066 30169
[25,] 70938 30936
[26,] 34077 35699
[27,] 45409 33228
[28,] 40809 27733
[29,] 37013 33666
[30,] 44953 35429
[31,] 37848 27438
[32,] 32745  8170
[33,] 43412 62557
[34,] 34931    NA
[35,] 33008    NA
[36,]  8620    NA
[37,] 68906    NA
[38,] 39556    NA
[39,] 50669    NA
[40,] 36432    NA
[41,] 40891    NA
[42,] 48428    NA
[43,] 36222    NA
[44,] 33425    NA
[45,] 39401    NA
[46,] 37967    NA
[47,] 34801    NA
[48,] 12657    NA
[49,] 69116    NA
[50,] 41519    NA
[51,] 51321    NA
[52,] 38529    NA
> darr
        [,1]   [,2]
 [1,] -30214  10526
 [2,]   6037 -13672
 [3,]  -5195   2497
 [4,]   5790   -459
 [5,]   5065   1449
 [6,]  -9780  -3990
 [7,]  -4321 -29127
 [8,]   6159  59413
 [9,]  -4374 -17654
[10,]  -6185  -3309
[11,] -20500  -5465
[12,]  58870   3877
[13,] -30996   2505
[14,]   7807  -8652
[15,]  -5363   -145
[16,]  -2989   1776
[17,]  10559    512
[18,] -10249 -11907
[19,]  -6790 -22668
[20,]   8642  49622
[21,]  -1487 -21806
[22,]  -4812    247
[23,] -22505  -4948
[24,]  58872    767
[25,] -36861   4763
[26,]  11332  -2471
[27,]  -4600  -5495
[28,]  -3796   5933
[29,]   7940   1763
[30,]  -7105  -7991
[31,]  -5103 -19268
[32,]  10667  54387
[33,]  -8481     NA
[34,]  -1923     NA
[35,] -24388     NA
[36,]  60286     NA
[37,] -29350     NA
[38,]  11113     NA
[39,] -14237     NA
[40,]   4459     NA
[41,]   7537     NA
[42,] -12206     NA
[43,]  -2797     NA
[44,]   5976     NA
[45,]  -1434     NA
[46,]  -3166     NA
[47,] -22144     NA
[48,]  56459     NA
[49,] -27597     NA
[50,]   9802     NA
[51,] -12792     NA
[52,]   3018     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/1f6q11272115734.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/2f6q11272115734.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/38x7m1272115734.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/48x7m1272115734.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] [,10] [,11]
[1,] 41547 38683 38401.0 39525.0 40439 41888 37898  8771 42438 38064 31879
[2,] 41547 38683 38401.0 39525.0 40439 41888 37898  8771 42438 38064 31879
[3,] 55222 45378 41560.5 40211.5 42877 46134 39249 22525 55311 44297 39550
[4,] 68897 52073 44720.0 40898.0 45315 50380 40600 36279 68184 50530 47221
[5,] 68897 52073 44720.0 40898.0 45315 50380 40600 36279 68184 50530 47221
       [,12] [,13]   [,14] [,15] [,16]   [,17] [,18] [,19] [,20] [,21]   [,22]
[1,] 11379.0 45633 39253.0 39486 39341 38708.0 41629 29722  7054 40870 34870.0
[2,] 11379.0 45633 39253.0 39486 39341 38708.0 41629 29722  7054 40870 34870.0
[3,] 26567.5 57941 43695.5 43273 40519 39912.5 45448 34370 19641 48773 37126.5
[4,] 41756.0 70249 48138.0 47060 41697 41117.0 49267 39018 32228 56676 39383.0
[5,] 41756.0 70249 48138.0 47060 41697 41117.0 49267 39018 32228 56676 39383.0
     [,23]   [,24] [,25] [,26]   [,27] [,28]   [,29] [,30] [,31]   [,32]
[1,] 34571 12066.0 30936 34077 33228.0 27733 33666.0 35429 27438  8170.0
[2,] 34571 12066.0 30936 34077 33228.0 27733 33666.0 35429 27438  8170.0
[3,] 34844 21117.5 50937 34888 39318.5 34271 35339.5 40191 32643 20457.5
[4,] 35117 30169.0 70938 35699 45409.0 40809 37013.0 44953 37848 32745.0
[5,] 35117 30169.0 70938 35699 45409.0 40809 37013.0 44953 37848 32745.0
       [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44]
[1,] 43412.0 34931 33008  8620 68906 39556 50669 36432 40891 48428 36222 33425
[2,] 43412.0 34931 33008  8620 68906 39556 50669 36432 40891 48428 36222 33425
[3,] 52984.5 34931 33008  8620 68906 39556 50669 36432 40891 48428 36222 33425
[4,] 62557.0 34931 33008  8620 68906 39556 50669 36432 40891 48428 36222 33425
[5,] 62557.0 34931 33008  8620 68906 39556 50669 36432 40891 48428 36222 33425
     [,45] [,46] [,47] [,48] [,49] [,50] [,51] [,52]
[1,] 39401 37967 34801 12657 69116 41519 51321 38529
[2,] 39401 37967 34801 12657 69116 41519 51321 38529
[3,] 39401 37967 34801 12657 69116 41519 51321 38529
[4,] 39401 37967 34801 12657 69116 41519 51321 38529
[5,] 39401 37967 34801 12657 69116 41519 51321 38529

$n
 [1] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1
[39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]      [,8]
[1,] 24665.79 30418.31 34500.73 38677.54 37429.39 36646.49 36230.25 -8207.727
[2,] 85778.21 60337.69 48620.27 41745.46 48324.61 55621.51 42267.75 53257.727
         [,9]    [,10]    [,11]     [,12]   [,13]    [,14]    [,15]    [,16]
[1,] 26546.83 30369.63 22409.48 -7370.557 30439.3 33768.92 34811.11 37886.81
[2,] 84075.17 58224.37 56690.52 60505.557 85442.7 53622.08 51734.89 43151.19
        [,17]    [,18]    [,19]     [,20]    [,21]    [,22]    [,23]      [,24]
[1,] 37221.10 36914.61 23984.24 -8484.116 31114.08 32084.45 34233.99   892.3086
[2,] 42603.90 53981.39 44755.76 47766.116 66431.92 42168.55 35454.01 41342.6914
         [,25]    [,26]    [,27]    [,28]    [,29]    [,30]    [,31]     [,32]
[1,]  6245.617 33075.86 25709.54 19662.12 31600.14 29550.51 21012.65 -6998.396
[2,] 95628.383 36700.14 52927.46 48879.88 39078.86 50831.49 44273.35 47913.396
        [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44]
[1,] 31595.16 34931 33008  8620 68906 39556 50669 36432 40891 48428 36222 33425
[2,] 74373.84 34931 33008  8620 68906 39556 50669 36432 40891 48428 36222 33425
     [,45] [,46] [,47] [,48] [,49] [,50] [,51] [,52]
[1,] 39401 37967 34801 12657 69116 41519 51321 38529
[2,] 39401 37967 34801 12657 69116 41519 51321 38529

$out
numeric(0)

$group
numeric(0)

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30"
[31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45"
[46] "46" "47" "48" "49" "50" "51" "52"

Warning message:
In bxp(list(stats = c(41547, 41547, 55222, 68897, 68897, 38683,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5jop71272115734.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,] -30214 -13672.0 -5195 -459.0 1449 -9780 -29127  6159 -17654 -6185 -20500.0
[2,] -30214 -13672.0 -5195 -459.0 1449 -9780 -29127  6159 -17654 -6185 -20500.0
[3,]  -9844  -3817.5 -1349 2665.5 3257 -6885 -16724 32786 -11014 -4747 -12982.5
[4,]  10526   6037.0  2497 5790.0 5065 -3990  -4321 59413  -4374 -3309  -5465.0
[5,]  10526   6037.0  2497 5790.0 5065 -3990  -4321 59413  -4374 -3309  -5465.0
       [,12]    [,13]   [,14] [,15]   [,16]   [,17]  [,18]  [,19] [,20]
[1,]  3877.0 -30996.0 -8652.0 -5363 -2989.0   512.0 -11907 -22668  8642
[2,]  3877.0 -30996.0 -8652.0 -5363 -2989.0   512.0 -11907 -22668  8642
[3,] 31373.5 -14245.5  -422.5 -2754  -606.5  5535.5 -11078 -14729 29132
[4,] 58870.0   2505.0  7807.0  -145  1776.0 10559.0 -10249  -6790 49622
[5,] 58870.0   2505.0  7807.0  -145  1776.0 10559.0 -10249  -6790 49622
        [,21]   [,22]    [,23]   [,24]  [,25]   [,26]   [,27]   [,28]  [,29]
[1,] -21806.0 -4812.0 -22505.0   767.0 -36861 -2471.0 -5495.0 -3796.0 1763.0
[2,] -21806.0 -4812.0 -22505.0   767.0 -36861 -2471.0 -5495.0 -3796.0 1763.0
[3,] -11646.5 -2282.5 -13726.5 29819.5 -16049  4430.5 -5047.5  1068.5 4851.5
[4,]  -1487.0   247.0  -4948.0 58872.0   4763 11332.0 -4600.0  5933.0 7940.0
[5,]  -1487.0   247.0  -4948.0 58872.0   4763 11332.0 -4600.0  5933.0 7940.0
     [,30]    [,31] [,32] [,33] [,34]  [,35] [,36]  [,37] [,38]  [,39] [,40]
[1,] -7991 -19268.0 10667 -8481 -1923 -24388 60286 -29350 11113 -14237  4459
[2,] -7991 -19268.0 10667 -8481 -1923 -24388 60286 -29350 11113 -14237  4459
[3,] -7548 -12185.5 32527 -8481 -1923 -24388 60286 -29350 11113 -14237  4459
[4,] -7105  -5103.0 54387 -8481 -1923 -24388 60286 -29350 11113 -14237  4459
[5,] -7105  -5103.0 54387 -8481 -1923 -24388 60286 -29350 11113 -14237  4459
     [,41]  [,42] [,43] [,44] [,45] [,46]  [,47] [,48]  [,49] [,50]  [,51]
[1,]  7537 -12206 -2797  5976 -1434 -3166 -22144 56459 -27597  9802 -12792
[2,]  7537 -12206 -2797  5976 -1434 -3166 -22144 56459 -27597  9802 -12792
[3,]  7537 -12206 -2797  5976 -1434 -3166 -22144 56459 -27597  9802 -12792
[4,]  7537 -12206 -2797  5976 -1434 -3166 -22144 56459 -27597  9802 -12792
[5,]  7537 -12206 -2797  5976 -1434 -3166 -22144 56459 -27597  9802 -12792
     [,52]
[1,]  3018
[2,]  3018
[3,]  3018
[4,]  3018
[5,]  3018

$n
 [1] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1
[39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1

$conf
         [,1]      [,2]      [,3]      [,4]     [,5]        [,6]      [,7]
[1,] -55359.9 -25836.96 -9942.723 -4316.062 -782.899 -13353.7543 -44437.98
[2,]  35671.9  18201.96  7244.723  9647.062 7296.899   -416.2457  10989.98
         [,8]       [,9]     [,10]      [,11]     [,12]     [,13]     [,14]
[1,] -26710.9 -25850.797 -7960.150 -29780.034 -30066.26 -51673.78 -18810.97
[2,]  92282.9   3822.797 -1533.850   3815.034  92813.26  23182.78  17965.97
         [,15]     [,16]     [,17]      [,18]      [,19]     [,20]     [,21]
[1,] -8583.699 -5930.095 -5689.297 -12930.365 -32468.358 -16652.03 -34347.47
[2,]  3075.699  4717.095 16760.297  -9225.635   3010.358  74916.03  11054.47
        [,22]      [,23]     [,24]     [,25]     [,26]    [,27]     [,28]
[1,] -7934.56 -33341.685 -35097.07 -62552.53 -10990.61 -6047.42 -9801.018
[2,]  3369.56   5888.685  94736.07  30454.53  19851.61 -4047.58 11938.018
         [,29]     [,30]      [,31]     [,32] [,33] [,34]  [,35] [,36]  [,37]
[1,] -2049.622 -8537.865 -28011.045 -16318.24 -8481 -1923 -24388 60286 -29350
[2,] 11752.622 -6558.135   3640.045  81372.24 -8481 -1923 -24388 60286 -29350
     [,38]  [,39] [,40] [,41]  [,42] [,43] [,44] [,45] [,46]  [,47] [,48]
[1,] 11113 -14237  4459  7537 -12206 -2797  5976 -1434 -3166 -22144 56459
[2,] 11113 -14237  4459  7537 -12206 -2797  5976 -1434 -3166 -22144 56459
      [,49] [,50]  [,51] [,52]
[1,] -27597  9802 -12792  3018
[2,] -27597  9802 -12792  3018

$out
numeric(0)

$group
numeric(0)

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30"
[31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45"
[46] "46" "47" "48" "49" "50" "51" "52"

Warning message:
In bxp(list(stats = c(-30214, -30214, -9844, 10526, 10526, -13672,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6jop71272115734.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]
[1,] 31879.0 27438
[2,] 35576.5 33228
[3,] 39392.0 39341
[4,] 45134.0 41888
[5,] 51321.0 52073

$n
[1] 52 33

$conf
         [,1]     [,2]
[1,] 37297.89 36959.13
[2,] 41486.11 41722.87

$out
 [1] 68897 11379 70249 12066 70938  8620 68906 12657 69116  8771 68184  7054
[13] 56676  8170 62557

$group
 [1] 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2

$names
[1] "1" NA 

> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7jop71272115734.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,] 19641.0 19641.0 19641.0
[2,] 34822.5 34822.5 34822.5
[3,] 39553.0 39553.0 39553.0
[4,] 45413.0 45413.0 45413.0
[5,] 57941.0 57941.0 57941.0

$n
[1] 52 52 52

$conf
         [,1]     [,2]     [,3]
[1,] 37232.55 37232.55 37232.55
[2,] 41873.45 41873.45 41873.45

$out
 [1]  8620 68906 12657 69116  8620 68906 12657 69116  8620 68906 12657 69116

$group
 [1] 1 1 1 1 2 2 2 2 3 3 3 3

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

> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1f6q11272115734.ps tmp/1f6q11272115734.png",intern=TRUE))
character(0)
> try(system("convert tmp/2f6q11272115734.ps tmp/2f6q11272115734.png",intern=TRUE))
character(0)
> try(system("convert tmp/38x7m1272115734.ps tmp/38x7m1272115734.png",intern=TRUE))
character(0)
> try(system("convert tmp/48x7m1272115734.ps tmp/48x7m1272115734.png",intern=TRUE))
character(0)
> try(system("convert tmp/5jop71272115734.ps tmp/5jop71272115734.png",intern=TRUE))
character(0)
> try(system("convert tmp/6jop71272115734.ps tmp/6jop71272115734.png",intern=TRUE))
character(0)
> try(system("convert tmp/7jop71272115734.ps tmp/7jop71272115734.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.336   0.958   3.386