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(23.98,24.24,24.92,25.46,25.84,26.08,26.18,26.34,26.42,26.38,26.04,25.58,25.65,25.56,25.62,25.62,25.69,25.68,25.68,25.83,25.93,26.11,24.72,24.62,24.65,25.24,25.56,25.9,25.87,25.78,25.78,25.74,25.78,25.73,24.67,24.31,24.56,25,25.38,25.99,26.22,26.19,26.22,26.22,26.61,26.72,25.46,25.48,25.59,25.88,26,26.97,27.2,27.19,27.19,27.19,27.26,26.9,26.11,25.87,26.02,26.31,26.37,26.52,26.86,26.92,26.98,26.98,27.03,26.75,26.39,26.3,26.3,26.52,26.53,26.98,27.22,27.34,27.41,27.47,27.46,27.53,27.21,26.91,26.95,26.91,27.39,27.62,27.79,27.88,27.9,28.09,28.46,28.73,27.93,27.61,27.65,28.19,28.98,28.99,29.02,29,29.04,29.19,29.23,29.26,29.02,28.47,28.53,28.48,28.68,28.89,29.2,29.21,29.15,29.22,29.34,29.13,28.84,28.76) > par1 = '120' > #'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] 120 > (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] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [38] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [75] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [112] 1 1 1 1 1 1 1 1 1 > arr [,1] [,2] [1,] 23.98 NA [2,] 24.24 NA [3,] 24.92 NA [4,] 25.46 NA [5,] 25.84 NA [6,] 26.08 NA [7,] 26.18 NA [8,] 26.34 NA [9,] 26.42 NA [10,] 26.38 NA [11,] 26.04 NA [12,] 25.58 NA [13,] 25.65 NA [14,] 25.56 NA [15,] 25.62 NA [16,] 25.62 NA [17,] 25.69 NA [18,] 25.68 NA [19,] 25.68 NA [20,] 25.83 NA [21,] 25.93 NA [22,] 26.11 NA [23,] 24.72 NA [24,] 24.62 NA [25,] 24.65 NA [26,] 25.24 NA [27,] 25.56 NA [28,] 25.90 NA [29,] 25.87 NA [30,] 25.78 NA [31,] 25.78 NA [32,] 25.74 NA [33,] 25.78 NA [34,] 25.73 NA [35,] 24.67 NA [36,] 24.31 NA [37,] 24.56 NA [38,] 25.00 NA [39,] 25.38 NA [40,] 25.99 NA [41,] 26.22 NA [42,] 26.19 NA [43,] 26.22 NA [44,] 26.22 NA [45,] 26.61 NA [46,] 26.72 NA [47,] 25.46 NA [48,] 25.48 NA [49,] 25.59 NA [50,] 25.88 NA [51,] 26.00 NA [52,] 26.97 NA [53,] 27.20 NA [54,] 27.19 NA [55,] 27.19 NA [56,] 27.19 NA [57,] 27.26 NA [58,] 26.90 NA [59,] 26.11 NA [60,] 25.87 NA [61,] 26.02 NA [62,] 26.31 NA [63,] 26.37 NA [64,] 26.52 NA [65,] 26.86 NA [66,] 26.92 NA [67,] 26.98 NA [68,] 26.98 NA [69,] 27.03 NA [70,] 26.75 NA [71,] 26.39 NA [72,] 26.30 NA [73,] 26.30 NA [74,] 26.52 NA [75,] 26.53 NA [76,] 26.98 NA [77,] 27.22 NA [78,] 27.34 NA [79,] 27.41 NA [80,] 27.47 NA [81,] 27.46 NA [82,] 27.53 NA [83,] 27.21 NA [84,] 26.91 NA [85,] 26.95 NA [86,] 26.91 NA [87,] 27.39 NA [88,] 27.62 NA [89,] 27.79 NA [90,] 27.88 NA [91,] 27.90 NA [92,] 28.09 NA [93,] 28.46 NA [94,] 28.73 NA [95,] 27.93 NA [96,] 27.61 NA [97,] 27.65 NA [98,] 28.19 NA [99,] 28.98 NA [100,] 28.99 NA [101,] 29.02 NA [102,] 29.00 NA [103,] 29.04 NA [104,] 29.19 NA [105,] 29.23 NA [106,] 29.26 NA [107,] 29.02 NA [108,] 28.47 NA [109,] 28.53 NA [110,] 28.48 NA [111,] 28.68 NA [112,] 28.89 NA [113,] 29.20 NA [114,] 29.21 NA [115,] 29.15 NA [116,] 29.22 NA [117,] 29.34 NA [118,] 29.13 NA [119,] 28.84 NA [120,] 28.76 NA > darr [,1] [,2] [1,] 0.26 NA [2,] 0.68 NA [3,] 0.54 NA [4,] 0.38 NA [5,] 0.24 NA [6,] 0.10 NA [7,] 0.16 NA [8,] 0.08 NA [9,] -0.04 NA [10,] -0.34 NA [11,] -0.46 NA [12,] 0.07 NA [13,] -0.09 NA [14,] 0.06 NA [15,] 0.00 NA [16,] 0.07 NA [17,] -0.01 NA [18,] 0.00 NA [19,] 0.15 NA [20,] 0.10 NA [21,] 0.18 NA [22,] -1.39 NA [23,] -0.10 NA [24,] 0.03 NA [25,] 0.59 NA [26,] 0.32 NA [27,] 0.34 NA [28,] -0.03 NA [29,] -0.09 NA [30,] 0.00 NA [31,] -0.04 NA [32,] 0.04 NA [33,] -0.05 NA [34,] -1.06 NA [35,] -0.36 NA [36,] 0.25 NA [37,] 0.44 NA [38,] 0.38 NA [39,] 0.61 NA [40,] 0.23 NA [41,] -0.03 NA [42,] 0.03 NA [43,] 0.00 NA [44,] 0.39 NA [45,] 0.11 NA [46,] -1.26 NA [47,] 0.02 NA [48,] 0.11 NA [49,] 0.29 NA [50,] 0.12 NA [51,] 0.97 NA [52,] 0.23 NA [53,] -0.01 NA [54,] 0.00 NA [55,] 0.00 NA [56,] 0.07 NA [57,] -0.36 NA [58,] -0.79 NA [59,] -0.24 NA [60,] 0.15 NA [61,] 0.29 NA [62,] 0.06 NA [63,] 0.15 NA [64,] 0.34 NA [65,] 0.06 NA [66,] 0.06 NA [67,] 0.00 NA [68,] 0.05 NA [69,] -0.28 NA [70,] -0.36 NA [71,] -0.09 NA [72,] 0.00 NA [73,] 0.22 NA [74,] 0.01 NA [75,] 0.45 NA [76,] 0.24 NA [77,] 0.12 NA [78,] 0.07 NA [79,] 0.06 NA [80,] -0.01 NA [81,] 0.07 NA [82,] -0.32 NA [83,] -0.30 NA [84,] 0.04 NA [85,] -0.04 NA [86,] 0.48 NA [87,] 0.23 NA [88,] 0.17 NA [89,] 0.09 NA [90,] 0.02 NA [91,] 0.19 NA [92,] 0.37 NA [93,] 0.27 NA [94,] -0.80 NA [95,] -0.32 NA [96,] 0.04 NA [97,] 0.54 NA [98,] 0.79 NA [99,] 0.01 NA [100,] 0.03 NA [101,] -0.02 NA [102,] 0.04 NA [103,] 0.15 NA [104,] 0.04 NA [105,] 0.03 NA [106,] -0.24 NA [107,] -0.55 NA [108,] 0.06 NA [109,] -0.05 NA [110,] 0.20 NA [111,] 0.21 NA [112,] 0.31 NA [113,] 0.01 NA [114,] -0.06 NA [115,] 0.07 NA [116,] 0.12 NA [117,] -0.21 NA [118,] -0.29 NA [119,] -0.08 NA [120,] 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/13hh31259762378.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/2z7hn1259762378.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/3ltnl1259762378.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/4r71j1259762378.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] [,12] [1,] 23.98 24.24 24.92 25.46 25.84 26.08 26.18 26.34 26.42 26.38 26.04 25.58 [2,] 23.98 24.24 24.92 25.46 25.84 26.08 26.18 26.34 26.42 26.38 26.04 25.58 [3,] 23.98 24.24 24.92 25.46 25.84 26.08 26.18 26.34 26.42 26.38 26.04 25.58 [4,] 23.98 24.24 24.92 25.46 25.84 26.08 26.18 26.34 26.42 26.38 26.04 25.58 [5,] 23.98 24.24 24.92 25.46 25.84 26.08 26.18 26.34 26.42 26.38 26.04 25.58 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 25.65 25.56 25.62 25.62 25.69 25.68 25.68 25.83 25.93 26.11 24.72 24.62 [2,] 25.65 25.56 25.62 25.62 25.69 25.68 25.68 25.83 25.93 26.11 24.72 24.62 [3,] 25.65 25.56 25.62 25.62 25.69 25.68 25.68 25.83 25.93 26.11 24.72 24.62 [4,] 25.65 25.56 25.62 25.62 25.69 25.68 25.68 25.83 25.93 26.11 24.72 24.62 [5,] 25.65 25.56 25.62 25.62 25.69 25.68 25.68 25.83 25.93 26.11 24.72 24.62 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 24.65 25.24 25.56 25.9 25.87 25.78 25.78 25.74 25.78 25.73 24.67 24.31 [2,] 24.65 25.24 25.56 25.9 25.87 25.78 25.78 25.74 25.78 25.73 24.67 24.31 [3,] 24.65 25.24 25.56 25.9 25.87 25.78 25.78 25.74 25.78 25.73 24.67 24.31 [4,] 24.65 25.24 25.56 25.9 25.87 25.78 25.78 25.74 25.78 25.73 24.67 24.31 [5,] 24.65 25.24 25.56 25.9 25.87 25.78 25.78 25.74 25.78 25.73 24.67 24.31 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 24.56 25 25.38 25.99 26.22 26.19 26.22 26.22 26.61 26.72 25.46 25.48 [2,] 24.56 25 25.38 25.99 26.22 26.19 26.22 26.22 26.61 26.72 25.46 25.48 [3,] 24.56 25 25.38 25.99 26.22 26.19 26.22 26.22 26.61 26.72 25.46 25.48 [4,] 24.56 25 25.38 25.99 26.22 26.19 26.22 26.22 26.61 26.72 25.46 25.48 [5,] 24.56 25 25.38 25.99 26.22 26.19 26.22 26.22 26.61 26.72 25.46 25.48 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 25.59 25.88 26 26.97 27.2 27.19 27.19 27.19 27.26 26.9 26.11 25.87 [2,] 25.59 25.88 26 26.97 27.2 27.19 27.19 27.19 27.26 26.9 26.11 25.87 [3,] 25.59 25.88 26 26.97 27.2 27.19 27.19 27.19 27.26 26.9 26.11 25.87 [4,] 25.59 25.88 26 26.97 27.2 27.19 27.19 27.19 27.26 26.9 26.11 25.87 [5,] 25.59 25.88 26 26.97 27.2 27.19 27.19 27.19 27.26 26.9 26.11 25.87 [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 26.02 26.31 26.37 26.52 26.86 26.92 26.98 26.98 27.03 26.75 26.39 26.3 [2,] 26.02 26.31 26.37 26.52 26.86 26.92 26.98 26.98 27.03 26.75 26.39 26.3 [3,] 26.02 26.31 26.37 26.52 26.86 26.92 26.98 26.98 27.03 26.75 26.39 26.3 [4,] 26.02 26.31 26.37 26.52 26.86 26.92 26.98 26.98 27.03 26.75 26.39 26.3 [5,] 26.02 26.31 26.37 26.52 26.86 26.92 26.98 26.98 27.03 26.75 26.39 26.3 [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [1,] 26.3 26.52 26.53 26.98 27.22 27.34 27.41 27.47 27.46 27.53 27.21 26.91 [2,] 26.3 26.52 26.53 26.98 27.22 27.34 27.41 27.47 27.46 27.53 27.21 26.91 [3,] 26.3 26.52 26.53 26.98 27.22 27.34 27.41 27.47 27.46 27.53 27.21 26.91 [4,] 26.3 26.52 26.53 26.98 27.22 27.34 27.41 27.47 27.46 27.53 27.21 26.91 [5,] 26.3 26.52 26.53 26.98 27.22 27.34 27.41 27.47 27.46 27.53 27.21 26.91 [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [1,] 26.95 26.91 27.39 27.62 27.79 27.88 27.9 28.09 28.46 28.73 27.93 27.61 [2,] 26.95 26.91 27.39 27.62 27.79 27.88 27.9 28.09 28.46 28.73 27.93 27.61 [3,] 26.95 26.91 27.39 27.62 27.79 27.88 27.9 28.09 28.46 28.73 27.93 27.61 [4,] 26.95 26.91 27.39 27.62 27.79 27.88 27.9 28.09 28.46 28.73 27.93 27.61 [5,] 26.95 26.91 27.39 27.62 27.79 27.88 27.9 28.09 28.46 28.73 27.93 27.61 [,97] [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [1,] 27.65 28.19 28.98 28.99 29.02 29 29.04 29.19 29.23 29.26 29.02 [2,] 27.65 28.19 28.98 28.99 29.02 29 29.04 29.19 29.23 29.26 29.02 [3,] 27.65 28.19 28.98 28.99 29.02 29 29.04 29.19 29.23 29.26 29.02 [4,] 27.65 28.19 28.98 28.99 29.02 29 29.04 29.19 29.23 29.26 29.02 [5,] 27.65 28.19 28.98 28.99 29.02 29 29.04 29.19 29.23 29.26 29.02 [,108] [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [1,] 28.47 28.53 28.48 28.68 28.89 29.2 29.21 29.15 29.22 29.34 [2,] 28.47 28.53 28.48 28.68 28.89 29.2 29.21 29.15 29.22 29.34 [3,] 28.47 28.53 28.48 28.68 28.89 29.2 29.21 29.15 29.22 29.34 [4,] 28.47 28.53 28.48 28.68 28.89 29.2 29.21 29.15 29.22 29.34 [5,] 28.47 28.53 28.48 28.68 28.89 29.2 29.21 29.15 29.22 29.34 [,118] [,119] [,120] [1,] 29.13 28.84 28.76 [2,] 29.13 28.84 28.76 [3,] 29.13 28.84 28.76 [4,] 29.13 28.84 28.76 [5,] 29.13 28.84 28.76 $n [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [38] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [75] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [112] 1 1 1 1 1 1 1 1 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 23.98 24.24 24.92 25.46 25.84 26.08 26.18 26.34 26.42 26.38 26.04 25.58 [2,] 23.98 24.24 24.92 25.46 25.84 26.08 26.18 26.34 26.42 26.38 26.04 25.58 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 25.65 25.56 25.62 25.62 25.69 25.68 25.68 25.83 25.93 26.11 24.72 24.62 [2,] 25.65 25.56 25.62 25.62 25.69 25.68 25.68 25.83 25.93 26.11 24.72 24.62 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 24.65 25.24 25.56 25.9 25.87 25.78 25.78 25.74 25.78 25.73 24.67 24.31 [2,] 24.65 25.24 25.56 25.9 25.87 25.78 25.78 25.74 25.78 25.73 24.67 24.31 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 24.56 25 25.38 25.99 26.22 26.19 26.22 26.22 26.61 26.72 25.46 25.48 [2,] 24.56 25 25.38 25.99 26.22 26.19 26.22 26.22 26.61 26.72 25.46 25.48 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 25.59 25.88 26 26.97 27.2 27.19 27.19 27.19 27.26 26.9 26.11 25.87 [2,] 25.59 25.88 26 26.97 27.2 27.19 27.19 27.19 27.26 26.9 26.11 25.87 [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 26.02 26.31 26.37 26.52 26.86 26.92 26.98 26.98 27.03 26.75 26.39 26.3 [2,] 26.02 26.31 26.37 26.52 26.86 26.92 26.98 26.98 27.03 26.75 26.39 26.3 [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [1,] 26.3 26.52 26.53 26.98 27.22 27.34 27.41 27.47 27.46 27.53 27.21 26.91 [2,] 26.3 26.52 26.53 26.98 27.22 27.34 27.41 27.47 27.46 27.53 27.21 26.91 [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [1,] 26.95 26.91 27.39 27.62 27.79 27.88 27.9 28.09 28.46 28.73 27.93 27.61 [2,] 26.95 26.91 27.39 27.62 27.79 27.88 27.9 28.09 28.46 28.73 27.93 27.61 [,97] [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [1,] 27.65 28.19 28.98 28.99 29.02 29 29.04 29.19 29.23 29.26 29.02 [2,] 27.65 28.19 28.98 28.99 29.02 29 29.04 29.19 29.23 29.26 29.02 [,108] [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [1,] 28.47 28.53 28.48 28.68 28.89 29.2 29.21 29.15 29.22 29.34 [2,] 28.47 28.53 28.48 28.68 28.89 29.2 29.21 29.15 29.22 29.34 [,118] [,119] [,120] [1,] 29.13 28.84 28.76 [2,] 29.13 28.84 28.76 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" [13] "13" "14" "15" "16" "17" "18" "19" "20" "21" "22" "23" "24" [25] "25" "26" "27" "28" "29" "30" "31" "32" "33" "34" "35" "36" [37] "37" "38" "39" "40" "41" "42" "43" "44" "45" "46" "47" "48" [49] "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60" [61] "61" "62" "63" "64" "65" "66" "67" "68" "69" "70" "71" "72" [73] "73" "74" "75" "76" "77" "78" "79" "80" "81" "82" "83" "84" [85] "85" "86" "87" "88" "89" "90" "91" "92" "93" "94" "95" "96" [97] "97" "98" "99" "100" "101" "102" "103" "104" "105" "106" "107" "108" [109] "109" "110" "111" "112" "113" "114" "115" "116" "117" "118" "119" "120" > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/59qdt1259762378.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] [,12] [,13] [1,] 0.26 0.68 0.54 0.38 0.24 0.1 0.16 0.08 -0.04 -0.34 -0.46 0.07 -0.09 [2,] 0.26 0.68 0.54 0.38 0.24 0.1 0.16 0.08 -0.04 -0.34 -0.46 0.07 -0.09 [3,] 0.26 0.68 0.54 0.38 0.24 0.1 0.16 0.08 -0.04 -0.34 -0.46 0.07 -0.09 [4,] 0.26 0.68 0.54 0.38 0.24 0.1 0.16 0.08 -0.04 -0.34 -0.46 0.07 -0.09 [5,] 0.26 0.68 0.54 0.38 0.24 0.1 0.16 0.08 -0.04 -0.34 -0.46 0.07 -0.09 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 0.06 0 0.07 -0.01 0 0.15 0.1 0.18 -1.39 -0.1 0.03 0.59 [2,] 0.06 0 0.07 -0.01 0 0.15 0.1 0.18 -1.39 -0.1 0.03 0.59 [3,] 0.06 0 0.07 -0.01 0 0.15 0.1 0.18 -1.39 -0.1 0.03 0.59 [4,] 0.06 0 0.07 -0.01 0 0.15 0.1 0.18 -1.39 -0.1 0.03 0.59 [5,] 0.06 0 0.07 -0.01 0 0.15 0.1 0.18 -1.39 -0.1 0.03 0.59 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 0.32 0.34 -0.03 -0.09 0 -0.04 0.04 -0.05 -1.06 -0.36 0.25 0.44 [2,] 0.32 0.34 -0.03 -0.09 0 -0.04 0.04 -0.05 -1.06 -0.36 0.25 0.44 [3,] 0.32 0.34 -0.03 -0.09 0 -0.04 0.04 -0.05 -1.06 -0.36 0.25 0.44 [4,] 0.32 0.34 -0.03 -0.09 0 -0.04 0.04 -0.05 -1.06 -0.36 0.25 0.44 [5,] 0.32 0.34 -0.03 -0.09 0 -0.04 0.04 -0.05 -1.06 -0.36 0.25 0.44 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 0.38 0.61 0.23 -0.03 0.03 0 0.39 0.11 -1.26 0.02 0.11 0.29 [2,] 0.38 0.61 0.23 -0.03 0.03 0 0.39 0.11 -1.26 0.02 0.11 0.29 [3,] 0.38 0.61 0.23 -0.03 0.03 0 0.39 0.11 -1.26 0.02 0.11 0.29 [4,] 0.38 0.61 0.23 -0.03 0.03 0 0.39 0.11 -1.26 0.02 0.11 0.29 [5,] 0.38 0.61 0.23 -0.03 0.03 0 0.39 0.11 -1.26 0.02 0.11 0.29 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 0.12 0.97 0.23 -0.01 0 0 0.07 -0.36 -0.79 -0.24 0.15 0.29 [2,] 0.12 0.97 0.23 -0.01 0 0 0.07 -0.36 -0.79 -0.24 0.15 0.29 [3,] 0.12 0.97 0.23 -0.01 0 0 0.07 -0.36 -0.79 -0.24 0.15 0.29 [4,] 0.12 0.97 0.23 -0.01 0 0 0.07 -0.36 -0.79 -0.24 0.15 0.29 [5,] 0.12 0.97 0.23 -0.01 0 0 0.07 -0.36 -0.79 -0.24 0.15 0.29 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 0.06 0.15 0.34 0.06 0.06 0 0.05 -0.28 -0.36 -0.09 0 0.22 [2,] 0.06 0.15 0.34 0.06 0.06 0 0.05 -0.28 -0.36 -0.09 0 0.22 [3,] 0.06 0.15 0.34 0.06 0.06 0 0.05 -0.28 -0.36 -0.09 0 0.22 [4,] 0.06 0.15 0.34 0.06 0.06 0 0.05 -0.28 -0.36 -0.09 0 0.22 [5,] 0.06 0.15 0.34 0.06 0.06 0 0.05 -0.28 -0.36 -0.09 0 0.22 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 0.01 0.45 0.24 0.12 0.07 0.06 -0.01 0.07 -0.32 -0.3 0.04 -0.04 [2,] 0.01 0.45 0.24 0.12 0.07 0.06 -0.01 0.07 -0.32 -0.3 0.04 -0.04 [3,] 0.01 0.45 0.24 0.12 0.07 0.06 -0.01 0.07 -0.32 -0.3 0.04 -0.04 [4,] 0.01 0.45 0.24 0.12 0.07 0.06 -0.01 0.07 -0.32 -0.3 0.04 -0.04 [5,] 0.01 0.45 0.24 0.12 0.07 0.06 -0.01 0.07 -0.32 -0.3 0.04 -0.04 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] 0.48 0.23 0.17 0.09 0.02 0.19 0.37 0.27 -0.8 -0.32 0.04 0.54 [2,] 0.48 0.23 0.17 0.09 0.02 0.19 0.37 0.27 -0.8 -0.32 0.04 0.54 [3,] 0.48 0.23 0.17 0.09 0.02 0.19 0.37 0.27 -0.8 -0.32 0.04 0.54 [4,] 0.48 0.23 0.17 0.09 0.02 0.19 0.37 0.27 -0.8 -0.32 0.04 0.54 [5,] 0.48 0.23 0.17 0.09 0.02 0.19 0.37 0.27 -0.8 -0.32 0.04 0.54 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 0.79 0.01 0.03 -0.02 0.04 0.15 0.04 0.03 -0.24 -0.55 0.06 [2,] 0.79 0.01 0.03 -0.02 0.04 0.15 0.04 0.03 -0.24 -0.55 0.06 [3,] 0.79 0.01 0.03 -0.02 0.04 0.15 0.04 0.03 -0.24 -0.55 0.06 [4,] 0.79 0.01 0.03 -0.02 0.04 0.15 0.04 0.03 -0.24 -0.55 0.06 [5,] 0.79 0.01 0.03 -0.02 0.04 0.15 0.04 0.03 -0.24 -0.55 0.06 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] -0.05 0.2 0.21 0.31 0.01 -0.06 0.07 0.12 -0.21 -0.29 [2,] -0.05 0.2 0.21 0.31 0.01 -0.06 0.07 0.12 -0.21 -0.29 [3,] -0.05 0.2 0.21 0.31 0.01 -0.06 0.07 0.12 -0.21 -0.29 [4,] -0.05 0.2 0.21 0.31 0.01 -0.06 0.07 0.12 -0.21 -0.29 [5,] -0.05 0.2 0.21 0.31 0.01 -0.06 0.07 0.12 -0.21 -0.29 [,119] [,120] [1,] -0.08 NA [2,] -0.08 NA [3,] -0.08 NA [4,] -0.08 NA [5,] -0.08 NA $n [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [38] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [75] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [112] 1 1 1 1 1 1 1 1 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 0.26 0.68 0.54 0.38 0.24 0.1 0.16 0.08 -0.04 -0.34 -0.46 0.07 -0.09 [2,] 0.26 0.68 0.54 0.38 0.24 0.1 0.16 0.08 -0.04 -0.34 -0.46 0.07 -0.09 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 0.06 0 0.07 -0.01 0 0.15 0.1 0.18 -1.39 -0.1 0.03 0.59 [2,] 0.06 0 0.07 -0.01 0 0.15 0.1 0.18 -1.39 -0.1 0.03 0.59 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 0.32 0.34 -0.03 -0.09 0 -0.04 0.04 -0.05 -1.06 -0.36 0.25 0.44 [2,] 0.32 0.34 -0.03 -0.09 0 -0.04 0.04 -0.05 -1.06 -0.36 0.25 0.44 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 0.38 0.61 0.23 -0.03 0.03 0 0.39 0.11 -1.26 0.02 0.11 0.29 [2,] 0.38 0.61 0.23 -0.03 0.03 0 0.39 0.11 -1.26 0.02 0.11 0.29 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 0.12 0.97 0.23 -0.01 0 0 0.07 -0.36 -0.79 -0.24 0.15 0.29 [2,] 0.12 0.97 0.23 -0.01 0 0 0.07 -0.36 -0.79 -0.24 0.15 0.29 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 0.06 0.15 0.34 0.06 0.06 0 0.05 -0.28 -0.36 -0.09 0 0.22 [2,] 0.06 0.15 0.34 0.06 0.06 0 0.05 -0.28 -0.36 -0.09 0 0.22 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 0.01 0.45 0.24 0.12 0.07 0.06 -0.01 0.07 -0.32 -0.3 0.04 -0.04 [2,] 0.01 0.45 0.24 0.12 0.07 0.06 -0.01 0.07 -0.32 -0.3 0.04 -0.04 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] 0.48 0.23 0.17 0.09 0.02 0.19 0.37 0.27 -0.8 -0.32 0.04 0.54 [2,] 0.48 0.23 0.17 0.09 0.02 0.19 0.37 0.27 -0.8 -0.32 0.04 0.54 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 0.79 0.01 0.03 -0.02 0.04 0.15 0.04 0.03 -0.24 -0.55 0.06 [2,] 0.79 0.01 0.03 -0.02 0.04 0.15 0.04 0.03 -0.24 -0.55 0.06 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] -0.05 0.2 0.21 0.31 0.01 -0.06 0.07 0.12 -0.21 -0.29 [2,] -0.05 0.2 0.21 0.31 0.01 -0.06 0.07 0.12 -0.21 -0.29 [,119] [,120] [1,] -0.08 NA [2,] -0.08 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" [13] "13" "14" "15" "16" "17" "18" "19" "20" "21" "22" "23" "24" [25] "25" "26" "27" "28" "29" "30" "31" "32" "33" "34" "35" "36" [37] "37" "38" "39" "40" "41" "42" "43" "44" "45" "46" "47" "48" [49] "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60" [61] "61" "62" "63" "64" "65" "66" "67" "68" "69" "70" "71" "72" [73] "73" "74" "75" "76" "77" "78" "79" "80" "81" "82" "83" "84" [85] "85" "86" "87" "88" "89" "90" "91" "92" "93" "94" "95" "96" [97] "97" "98" "99" "100" "101" "102" "103" "104" "105" "106" "107" "108" [109] "109" "110" "111" "112" "113" "114" "115" "116" "117" "118" "119" "120" > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6bw7q1259762378.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,] 23.980 NA [2,] 25.805 NA [3,] 26.570 NA [4,] 27.720 NA [5,] 29.340 NA $n [1] 120 0 $conf [,1] [,2] [1,] 26.29379 NA [2,] 26.84621 NA $out numeric(0) $group numeric(0) $names [1] "1" NA > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7weou1259762378.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,] 23.980 23.980 23.980 [2,] 25.805 25.805 25.805 [3,] 26.570 26.570 26.570 [4,] 27.720 27.720 27.720 [5,] 29.340 29.340 29.340 $n [1] 120 120 120 $conf [,1] [,2] [,3] [1,] 26.29379 26.29379 26.29379 [2,] 26.84621 26.84621 26.84621 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > system("convert tmp/13hh31259762378.ps tmp/13hh31259762378.png") > system("convert tmp/2z7hn1259762378.ps tmp/2z7hn1259762378.png") > system("convert tmp/3ltnl1259762378.ps tmp/3ltnl1259762378.png") > system("convert tmp/4r71j1259762378.ps tmp/4r71j1259762378.png") > system("convert tmp/59qdt1259762378.ps tmp/59qdt1259762378.png") > system("convert tmp/6bw7q1259762378.ps tmp/6bw7q1259762378.png") > system("convert tmp/7weou1259762378.ps tmp/7weou1259762378.png") > > > proc.time() user system elapsed 1.568 1.024 1.965