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(157.4,157.5,157.5,157.8,158.1,158.4,158.3,158.5,158.7,159.6,159.5,160,160.2,160.5,161.1,161.6,161.8,162.1,162.5,162.5,162.6,163.3,163.2,162.9,163.2,163.4,163.3,162.8,162.8,163.1,164,164.6,166.1,166.9,166.4,166.9,168.8,170.3,170,170.1,170.8,170.7,170.4,171,171.1,170.9,171.7,172.1,172.7,173.1,173.2,173.7,174.9,174.9,174.5,174.3,175.5,175.7,175.9,175.8,176.2,175.9,176.5,175.9,176.2,176.7,177,177,177,177.6,178.6,179,180.1,180.7,181,181.7,182.8,183.4,184.7,184.1,184.9,187.6,189.3,189.7,191,192,193,193.5,193.4,194.6,195.4,196.5,200.6,204,206.2,210,216.6,224.2,228.3,228.5,230.2,232.2,233.1,233,230.3,229.2,227.7,226,226.8,225.8,224.3,223.4,223.5,224.5,221.3,220.4) > par1 = '116' > #'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] 116 > (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 > arr [,1] [,2] [1,] 157.4 NA [2,] 157.5 NA [3,] 157.5 NA [4,] 157.8 NA [5,] 158.1 NA [6,] 158.4 NA [7,] 158.3 NA [8,] 158.5 NA [9,] 158.7 NA [10,] 159.6 NA [11,] 159.5 NA [12,] 160.0 NA [13,] 160.2 NA [14,] 160.5 NA [15,] 161.1 NA [16,] 161.6 NA [17,] 161.8 NA [18,] 162.1 NA [19,] 162.5 NA [20,] 162.5 NA [21,] 162.6 NA [22,] 163.3 NA [23,] 163.2 NA [24,] 162.9 NA [25,] 163.2 NA [26,] 163.4 NA [27,] 163.3 NA [28,] 162.8 NA [29,] 162.8 NA [30,] 163.1 NA [31,] 164.0 NA [32,] 164.6 NA [33,] 166.1 NA [34,] 166.9 NA [35,] 166.4 NA [36,] 166.9 NA [37,] 168.8 NA [38,] 170.3 NA [39,] 170.0 NA [40,] 170.1 NA [41,] 170.8 NA [42,] 170.7 NA [43,] 170.4 NA [44,] 171.0 NA [45,] 171.1 NA [46,] 170.9 NA [47,] 171.7 NA [48,] 172.1 NA [49,] 172.7 NA [50,] 173.1 NA [51,] 173.2 NA [52,] 173.7 NA [53,] 174.9 NA [54,] 174.9 NA [55,] 174.5 NA [56,] 174.3 NA [57,] 175.5 NA [58,] 175.7 NA [59,] 175.9 NA [60,] 175.8 NA [61,] 176.2 NA [62,] 175.9 NA [63,] 176.5 NA [64,] 175.9 NA [65,] 176.2 NA [66,] 176.7 NA [67,] 177.0 NA [68,] 177.0 NA [69,] 177.0 NA [70,] 177.6 NA [71,] 178.6 NA [72,] 179.0 NA [73,] 180.1 NA [74,] 180.7 NA [75,] 181.0 NA [76,] 181.7 NA [77,] 182.8 NA [78,] 183.4 NA [79,] 184.7 NA [80,] 184.1 NA [81,] 184.9 NA [82,] 187.6 NA [83,] 189.3 NA [84,] 189.7 NA [85,] 191.0 NA [86,] 192.0 NA [87,] 193.0 NA [88,] 193.5 NA [89,] 193.4 NA [90,] 194.6 NA [91,] 195.4 NA [92,] 196.5 NA [93,] 200.6 NA [94,] 204.0 NA [95,] 206.2 NA [96,] 210.0 NA [97,] 216.6 NA [98,] 224.2 NA [99,] 228.3 NA [100,] 228.5 NA [101,] 230.2 NA [102,] 232.2 NA [103,] 233.1 NA [104,] 233.0 NA [105,] 230.3 NA [106,] 229.2 NA [107,] 227.7 NA [108,] 226.0 NA [109,] 226.8 NA [110,] 225.8 NA [111,] 224.3 NA [112,] 223.4 NA [113,] 223.5 NA [114,] 224.5 NA [115,] 221.3 NA [116,] 220.4 NA > darr [,1] [,2] [1,] 0.1 NA [2,] 0.0 NA [3,] 0.3 NA [4,] 0.3 NA [5,] 0.3 NA [6,] -0.1 NA [7,] 0.2 NA [8,] 0.2 NA [9,] 0.9 NA [10,] -0.1 NA [11,] 0.5 NA [12,] 0.2 NA [13,] 0.3 NA [14,] 0.6 NA [15,] 0.5 NA [16,] 0.2 NA [17,] 0.3 NA [18,] 0.4 NA [19,] 0.0 NA [20,] 0.1 NA [21,] 0.7 NA [22,] -0.1 NA [23,] -0.3 NA [24,] 0.3 NA [25,] 0.2 NA [26,] -0.1 NA [27,] -0.5 NA [28,] 0.0 NA [29,] 0.3 NA [30,] 0.9 NA [31,] 0.6 NA [32,] 1.5 NA [33,] 0.8 NA [34,] -0.5 NA [35,] 0.5 NA [36,] 1.9 NA [37,] 1.5 NA [38,] -0.3 NA [39,] 0.1 NA [40,] 0.7 NA [41,] -0.1 NA [42,] -0.3 NA [43,] 0.6 NA [44,] 0.1 NA [45,] -0.2 NA [46,] 0.8 NA [47,] 0.4 NA [48,] 0.6 NA [49,] 0.4 NA [50,] 0.1 NA [51,] 0.5 NA [52,] 1.2 NA [53,] 0.0 NA [54,] -0.4 NA [55,] -0.2 NA [56,] 1.2 NA [57,] 0.2 NA [58,] 0.2 NA [59,] -0.1 NA [60,] 0.4 NA [61,] -0.3 NA [62,] 0.6 NA [63,] -0.6 NA [64,] 0.3 NA [65,] 0.5 NA [66,] 0.3 NA [67,] 0.0 NA [68,] 0.0 NA [69,] 0.6 NA [70,] 1.0 NA [71,] 0.4 NA [72,] 1.1 NA [73,] 0.6 NA [74,] 0.3 NA [75,] 0.7 NA [76,] 1.1 NA [77,] 0.6 NA [78,] 1.3 NA [79,] -0.6 NA [80,] 0.8 NA [81,] 2.7 NA [82,] 1.7 NA [83,] 0.4 NA [84,] 1.3 NA [85,] 1.0 NA [86,] 1.0 NA [87,] 0.5 NA [88,] -0.1 NA [89,] 1.2 NA [90,] 0.8 NA [91,] 1.1 NA [92,] 4.1 NA [93,] 3.4 NA [94,] 2.2 NA [95,] 3.8 NA [96,] 6.6 NA [97,] 7.6 NA [98,] 4.1 NA [99,] 0.2 NA [100,] 1.7 NA [101,] 2.0 NA [102,] 0.9 NA [103,] -0.1 NA [104,] -2.7 NA [105,] -1.1 NA [106,] -1.5 NA [107,] -1.7 NA [108,] 0.8 NA [109,] -1.0 NA [110,] -1.5 NA [111,] -0.9 NA [112,] 0.1 NA [113,] 1.0 NA [114,] -3.2 NA [115,] -0.9 NA [116,] 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/171ux1259793495.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/2067z1259793495.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/3ihw81259793495.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/4zy8e1259793495.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,] 157.4 157.5 157.5 157.8 158.1 158.4 158.3 158.5 158.7 159.6 159.5 160 [2,] 157.4 157.5 157.5 157.8 158.1 158.4 158.3 158.5 158.7 159.6 159.5 160 [3,] 157.4 157.5 157.5 157.8 158.1 158.4 158.3 158.5 158.7 159.6 159.5 160 [4,] 157.4 157.5 157.5 157.8 158.1 158.4 158.3 158.5 158.7 159.6 159.5 160 [5,] 157.4 157.5 157.5 157.8 158.1 158.4 158.3 158.5 158.7 159.6 159.5 160 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 160.2 160.5 161.1 161.6 161.8 162.1 162.5 162.5 162.6 163.3 163.2 162.9 [2,] 160.2 160.5 161.1 161.6 161.8 162.1 162.5 162.5 162.6 163.3 163.2 162.9 [3,] 160.2 160.5 161.1 161.6 161.8 162.1 162.5 162.5 162.6 163.3 163.2 162.9 [4,] 160.2 160.5 161.1 161.6 161.8 162.1 162.5 162.5 162.6 163.3 163.2 162.9 [5,] 160.2 160.5 161.1 161.6 161.8 162.1 162.5 162.5 162.6 163.3 163.2 162.9 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 163.2 163.4 163.3 162.8 162.8 163.1 164 164.6 166.1 166.9 166.4 166.9 [2,] 163.2 163.4 163.3 162.8 162.8 163.1 164 164.6 166.1 166.9 166.4 166.9 [3,] 163.2 163.4 163.3 162.8 162.8 163.1 164 164.6 166.1 166.9 166.4 166.9 [4,] 163.2 163.4 163.3 162.8 162.8 163.1 164 164.6 166.1 166.9 166.4 166.9 [5,] 163.2 163.4 163.3 162.8 162.8 163.1 164 164.6 166.1 166.9 166.4 166.9 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 168.8 170.3 170 170.1 170.8 170.7 170.4 171 171.1 170.9 171.7 172.1 [2,] 168.8 170.3 170 170.1 170.8 170.7 170.4 171 171.1 170.9 171.7 172.1 [3,] 168.8 170.3 170 170.1 170.8 170.7 170.4 171 171.1 170.9 171.7 172.1 [4,] 168.8 170.3 170 170.1 170.8 170.7 170.4 171 171.1 170.9 171.7 172.1 [5,] 168.8 170.3 170 170.1 170.8 170.7 170.4 171 171.1 170.9 171.7 172.1 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 172.7 173.1 173.2 173.7 174.9 174.9 174.5 174.3 175.5 175.7 175.9 175.8 [2,] 172.7 173.1 173.2 173.7 174.9 174.9 174.5 174.3 175.5 175.7 175.9 175.8 [3,] 172.7 173.1 173.2 173.7 174.9 174.9 174.5 174.3 175.5 175.7 175.9 175.8 [4,] 172.7 173.1 173.2 173.7 174.9 174.9 174.5 174.3 175.5 175.7 175.9 175.8 [5,] 172.7 173.1 173.2 173.7 174.9 174.9 174.5 174.3 175.5 175.7 175.9 175.8 [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 176.2 175.9 176.5 175.9 176.2 176.7 177 177 177 177.6 178.6 179 [2,] 176.2 175.9 176.5 175.9 176.2 176.7 177 177 177 177.6 178.6 179 [3,] 176.2 175.9 176.5 175.9 176.2 176.7 177 177 177 177.6 178.6 179 [4,] 176.2 175.9 176.5 175.9 176.2 176.7 177 177 177 177.6 178.6 179 [5,] 176.2 175.9 176.5 175.9 176.2 176.7 177 177 177 177.6 178.6 179 [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [1,] 180.1 180.7 181 181.7 182.8 183.4 184.7 184.1 184.9 187.6 189.3 189.7 [2,] 180.1 180.7 181 181.7 182.8 183.4 184.7 184.1 184.9 187.6 189.3 189.7 [3,] 180.1 180.7 181 181.7 182.8 183.4 184.7 184.1 184.9 187.6 189.3 189.7 [4,] 180.1 180.7 181 181.7 182.8 183.4 184.7 184.1 184.9 187.6 189.3 189.7 [5,] 180.1 180.7 181 181.7 182.8 183.4 184.7 184.1 184.9 187.6 189.3 189.7 [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [1,] 191 192 193 193.5 193.4 194.6 195.4 196.5 200.6 204 206.2 210 [2,] 191 192 193 193.5 193.4 194.6 195.4 196.5 200.6 204 206.2 210 [3,] 191 192 193 193.5 193.4 194.6 195.4 196.5 200.6 204 206.2 210 [4,] 191 192 193 193.5 193.4 194.6 195.4 196.5 200.6 204 206.2 210 [5,] 191 192 193 193.5 193.4 194.6 195.4 196.5 200.6 204 206.2 210 [,97] [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [1,] 216.6 224.2 228.3 228.5 230.2 232.2 233.1 233 230.3 229.2 227.7 [2,] 216.6 224.2 228.3 228.5 230.2 232.2 233.1 233 230.3 229.2 227.7 [3,] 216.6 224.2 228.3 228.5 230.2 232.2 233.1 233 230.3 229.2 227.7 [4,] 216.6 224.2 228.3 228.5 230.2 232.2 233.1 233 230.3 229.2 227.7 [5,] 216.6 224.2 228.3 228.5 230.2 232.2 233.1 233 230.3 229.2 227.7 [,108] [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [1,] 226 226.8 225.8 224.3 223.4 223.5 224.5 221.3 220.4 [2,] 226 226.8 225.8 224.3 223.4 223.5 224.5 221.3 220.4 [3,] 226 226.8 225.8 224.3 223.4 223.5 224.5 221.3 220.4 [4,] 226 226.8 225.8 224.3 223.4 223.5 224.5 221.3 220.4 [5,] 226 226.8 225.8 224.3 223.4 223.5 224.5 221.3 220.4 $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 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 157.4 157.5 157.5 157.8 158.1 158.4 158.3 158.5 158.7 159.6 159.5 160 [2,] 157.4 157.5 157.5 157.8 158.1 158.4 158.3 158.5 158.7 159.6 159.5 160 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 160.2 160.5 161.1 161.6 161.8 162.1 162.5 162.5 162.6 163.3 163.2 162.9 [2,] 160.2 160.5 161.1 161.6 161.8 162.1 162.5 162.5 162.6 163.3 163.2 162.9 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 163.2 163.4 163.3 162.8 162.8 163.1 164 164.6 166.1 166.9 166.4 166.9 [2,] 163.2 163.4 163.3 162.8 162.8 163.1 164 164.6 166.1 166.9 166.4 166.9 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 168.8 170.3 170 170.1 170.8 170.7 170.4 171 171.1 170.9 171.7 172.1 [2,] 168.8 170.3 170 170.1 170.8 170.7 170.4 171 171.1 170.9 171.7 172.1 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 172.7 173.1 173.2 173.7 174.9 174.9 174.5 174.3 175.5 175.7 175.9 175.8 [2,] 172.7 173.1 173.2 173.7 174.9 174.9 174.5 174.3 175.5 175.7 175.9 175.8 [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 176.2 175.9 176.5 175.9 176.2 176.7 177 177 177 177.6 178.6 179 [2,] 176.2 175.9 176.5 175.9 176.2 176.7 177 177 177 177.6 178.6 179 [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [1,] 180.1 180.7 181 181.7 182.8 183.4 184.7 184.1 184.9 187.6 189.3 189.7 [2,] 180.1 180.7 181 181.7 182.8 183.4 184.7 184.1 184.9 187.6 189.3 189.7 [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [1,] 191 192 193 193.5 193.4 194.6 195.4 196.5 200.6 204 206.2 210 [2,] 191 192 193 193.5 193.4 194.6 195.4 196.5 200.6 204 206.2 210 [,97] [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [1,] 216.6 224.2 228.3 228.5 230.2 232.2 233.1 233 230.3 229.2 227.7 [2,] 216.6 224.2 228.3 228.5 230.2 232.2 233.1 233 230.3 229.2 227.7 [,108] [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [1,] 226 226.8 225.8 224.3 223.4 223.5 224.5 221.3 220.4 [2,] 226 226.8 225.8 224.3 223.4 223.5 224.5 221.3 220.4 $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" > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5zohm1259793495.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] [,14] [1,] 0.1 0 0.3 0.3 0.3 -0.1 0.2 0.2 0.9 -0.1 0.5 0.2 0.3 0.6 [2,] 0.1 0 0.3 0.3 0.3 -0.1 0.2 0.2 0.9 -0.1 0.5 0.2 0.3 0.6 [3,] 0.1 0 0.3 0.3 0.3 -0.1 0.2 0.2 0.9 -0.1 0.5 0.2 0.3 0.6 [4,] 0.1 0 0.3 0.3 0.3 -0.1 0.2 0.2 0.9 -0.1 0.5 0.2 0.3 0.6 [5,] 0.1 0 0.3 0.3 0.3 -0.1 0.2 0.2 0.9 -0.1 0.5 0.2 0.3 0.6 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 0.5 0.2 0.3 0.4 0 0.1 0.7 -0.1 -0.3 0.3 0.2 -0.1 [2,] 0.5 0.2 0.3 0.4 0 0.1 0.7 -0.1 -0.3 0.3 0.2 -0.1 [3,] 0.5 0.2 0.3 0.4 0 0.1 0.7 -0.1 -0.3 0.3 0.2 -0.1 [4,] 0.5 0.2 0.3 0.4 0 0.1 0.7 -0.1 -0.3 0.3 0.2 -0.1 [5,] 0.5 0.2 0.3 0.4 0 0.1 0.7 -0.1 -0.3 0.3 0.2 -0.1 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -0.5 0 0.3 0.9 0.6 1.5 0.8 -0.5 0.5 1.9 1.5 -0.3 [2,] -0.5 0 0.3 0.9 0.6 1.5 0.8 -0.5 0.5 1.9 1.5 -0.3 [3,] -0.5 0 0.3 0.9 0.6 1.5 0.8 -0.5 0.5 1.9 1.5 -0.3 [4,] -0.5 0 0.3 0.9 0.6 1.5 0.8 -0.5 0.5 1.9 1.5 -0.3 [5,] -0.5 0 0.3 0.9 0.6 1.5 0.8 -0.5 0.5 1.9 1.5 -0.3 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 0.1 0.7 -0.1 -0.3 0.6 0.1 -0.2 0.8 0.4 0.6 0.4 0.1 [2,] 0.1 0.7 -0.1 -0.3 0.6 0.1 -0.2 0.8 0.4 0.6 0.4 0.1 [3,] 0.1 0.7 -0.1 -0.3 0.6 0.1 -0.2 0.8 0.4 0.6 0.4 0.1 [4,] 0.1 0.7 -0.1 -0.3 0.6 0.1 -0.2 0.8 0.4 0.6 0.4 0.1 [5,] 0.1 0.7 -0.1 -0.3 0.6 0.1 -0.2 0.8 0.4 0.6 0.4 0.1 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 0.5 1.2 0 -0.4 -0.2 1.2 0.2 0.2 -0.1 0.4 -0.3 0.6 [2,] 0.5 1.2 0 -0.4 -0.2 1.2 0.2 0.2 -0.1 0.4 -0.3 0.6 [3,] 0.5 1.2 0 -0.4 -0.2 1.2 0.2 0.2 -0.1 0.4 -0.3 0.6 [4,] 0.5 1.2 0 -0.4 -0.2 1.2 0.2 0.2 -0.1 0.4 -0.3 0.6 [5,] 0.5 1.2 0 -0.4 -0.2 1.2 0.2 0.2 -0.1 0.4 -0.3 0.6 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] -0.6 0.3 0.5 0.3 0 0 0.6 1 0.4 1.1 0.6 0.3 [2,] -0.6 0.3 0.5 0.3 0 0 0.6 1 0.4 1.1 0.6 0.3 [3,] -0.6 0.3 0.5 0.3 0 0 0.6 1 0.4 1.1 0.6 0.3 [4,] -0.6 0.3 0.5 0.3 0 0 0.6 1 0.4 1.1 0.6 0.3 [5,] -0.6 0.3 0.5 0.3 0 0 0.6 1 0.4 1.1 0.6 0.3 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] 0.7 1.1 0.6 1.3 -0.6 0.8 2.7 1.7 0.4 1.3 1 1 [2,] 0.7 1.1 0.6 1.3 -0.6 0.8 2.7 1.7 0.4 1.3 1 1 [3,] 0.7 1.1 0.6 1.3 -0.6 0.8 2.7 1.7 0.4 1.3 1 1 [4,] 0.7 1.1 0.6 1.3 -0.6 0.8 2.7 1.7 0.4 1.3 1 1 [5,] 0.7 1.1 0.6 1.3 -0.6 0.8 2.7 1.7 0.4 1.3 1 1 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 0.5 -0.1 1.2 0.8 1.1 4.1 3.4 2.2 3.8 6.6 7.6 4.1 [2,] 0.5 -0.1 1.2 0.8 1.1 4.1 3.4 2.2 3.8 6.6 7.6 4.1 [3,] 0.5 -0.1 1.2 0.8 1.1 4.1 3.4 2.2 3.8 6.6 7.6 4.1 [4,] 0.5 -0.1 1.2 0.8 1.1 4.1 3.4 2.2 3.8 6.6 7.6 4.1 [5,] 0.5 -0.1 1.2 0.8 1.1 4.1 3.4 2.2 3.8 6.6 7.6 4.1 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 0.2 1.7 2 0.9 -0.1 -2.7 -1.1 -1.5 -1.7 0.8 [2,] 0.2 1.7 2 0.9 -0.1 -2.7 -1.1 -1.5 -1.7 0.8 [3,] 0.2 1.7 2 0.9 -0.1 -2.7 -1.1 -1.5 -1.7 0.8 [4,] 0.2 1.7 2 0.9 -0.1 -2.7 -1.1 -1.5 -1.7 0.8 [5,] 0.2 1.7 2 0.9 -0.1 -2.7 -1.1 -1.5 -1.7 0.8 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [1,] -1 -1.5 -0.9 0.1 1 -3.2 -0.9 NA [2,] -1 -1.5 -0.9 0.1 1 -3.2 -0.9 NA [3,] -1 -1.5 -0.9 0.1 1 -3.2 -0.9 NA [4,] -1 -1.5 -0.9 0.1 1 -3.2 -0.9 NA [5,] -1 -1.5 -0.9 0.1 1 -3.2 -0.9 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 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 0.1 0 0.3 0.3 0.3 -0.1 0.2 0.2 0.9 -0.1 0.5 0.2 0.3 0.6 [2,] 0.1 0 0.3 0.3 0.3 -0.1 0.2 0.2 0.9 -0.1 0.5 0.2 0.3 0.6 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 0.5 0.2 0.3 0.4 0 0.1 0.7 -0.1 -0.3 0.3 0.2 -0.1 [2,] 0.5 0.2 0.3 0.4 0 0.1 0.7 -0.1 -0.3 0.3 0.2 -0.1 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -0.5 0 0.3 0.9 0.6 1.5 0.8 -0.5 0.5 1.9 1.5 -0.3 [2,] -0.5 0 0.3 0.9 0.6 1.5 0.8 -0.5 0.5 1.9 1.5 -0.3 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 0.1 0.7 -0.1 -0.3 0.6 0.1 -0.2 0.8 0.4 0.6 0.4 0.1 [2,] 0.1 0.7 -0.1 -0.3 0.6 0.1 -0.2 0.8 0.4 0.6 0.4 0.1 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 0.5 1.2 0 -0.4 -0.2 1.2 0.2 0.2 -0.1 0.4 -0.3 0.6 [2,] 0.5 1.2 0 -0.4 -0.2 1.2 0.2 0.2 -0.1 0.4 -0.3 0.6 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] -0.6 0.3 0.5 0.3 0 0 0.6 1 0.4 1.1 0.6 0.3 [2,] -0.6 0.3 0.5 0.3 0 0 0.6 1 0.4 1.1 0.6 0.3 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] 0.7 1.1 0.6 1.3 -0.6 0.8 2.7 1.7 0.4 1.3 1 1 [2,] 0.7 1.1 0.6 1.3 -0.6 0.8 2.7 1.7 0.4 1.3 1 1 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 0.5 -0.1 1.2 0.8 1.1 4.1 3.4 2.2 3.8 6.6 7.6 4.1 [2,] 0.5 -0.1 1.2 0.8 1.1 4.1 3.4 2.2 3.8 6.6 7.6 4.1 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 0.2 1.7 2 0.9 -0.1 -2.7 -1.1 -1.5 -1.7 0.8 [2,] 0.2 1.7 2 0.9 -0.1 -2.7 -1.1 -1.5 -1.7 0.8 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [1,] -1 -1.5 -0.9 0.1 1 -3.2 -0.9 NA [2,] -1 -1.5 -0.9 0.1 1 -3.2 -0.9 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" > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6t9v31259793495.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,] 157.40 NA [2,] 163.35 NA [3,] 175.75 NA [4,] 193.20 NA [5,] 233.10 NA $n [1] 116 0 $conf [,1] [,2] [1,] 171.3710 NA [2,] 180.1290 NA $out numeric(0) $group numeric(0) $names [1] "1" NA > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7hnlx1259793496.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,] 157.40 157.40 157.40 [2,] 163.35 163.35 163.35 [3,] 175.75 175.75 175.75 [4,] 193.20 193.20 193.20 [5,] 233.10 233.10 233.10 $n [1] 116 116 116 $conf [,1] [,2] [,3] [1,] 171.3710 171.3710 171.3710 [2,] 180.1290 180.1290 180.1290 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > system("convert tmp/171ux1259793495.ps tmp/171ux1259793495.png") > system("convert tmp/2067z1259793495.ps tmp/2067z1259793495.png") > system("convert tmp/3ihw81259793495.ps tmp/3ihw81259793495.png") > system("convert tmp/4zy8e1259793495.ps tmp/4zy8e1259793495.png") > system("convert tmp/5zohm1259793495.ps tmp/5zohm1259793495.png") > system("convert tmp/6t9v31259793495.ps tmp/6t9v31259793495.png") > system("convert tmp/7hnlx1259793496.ps tmp/7hnlx1259793496.png") > > > proc.time() user system elapsed 1.523 1.016 2.167