R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) 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(93.6,103.5,127,117.5,111.5,137.6,103.2,86.9,124.4,113.6,101.6,148.5,108.3,117.2,128.7,116.5,131.7,139.9,107.4,96.1,126.5,116.4,109.8,148,111.4,117,141.7,120,132.1,146.7,122.5,99.6,122.7,139,117.8,125.5,134.5,121.3,126.7,117.7,123,132.1,113.1,89.2,121.7,105.3,85.3,105.3,72.2,92.1,97.2,78.6,78.1,93,81,65.9,88.6,85.7,76.3,96.8,76.8,85.6,119.2,91.4,95.7,112.3,95.2,82.8,111.3,108.2,97,124.4,99.3,117.6,131.5,114.2,116.8,116.5,105.4,89.2,115.8,111.4,106.4,128.4,107.7,111,129.8,130.5,142.9,159.9,84.1,75,100.7,106.8,97.4,113,76.9,87.3,103.7,92.1,92.9,112.2,88.7,74.6,101.5,119.7,120.7,153.5) > par1 = '108' > par1 <- '108' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Wessa P., (2012), Mean Plot (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_meanplot.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > par1 <- as.numeric(par1) > (n <- length(x)) [1] 108 > (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 > arr [,1] [,2] [1,] 93.6 NA [2,] 103.5 NA [3,] 127.0 NA [4,] 117.5 NA [5,] 111.5 NA [6,] 137.6 NA [7,] 103.2 NA [8,] 86.9 NA [9,] 124.4 NA [10,] 113.6 NA [11,] 101.6 NA [12,] 148.5 NA [13,] 108.3 NA [14,] 117.2 NA [15,] 128.7 NA [16,] 116.5 NA [17,] 131.7 NA [18,] 139.9 NA [19,] 107.4 NA [20,] 96.1 NA [21,] 126.5 NA [22,] 116.4 NA [23,] 109.8 NA [24,] 148.0 NA [25,] 111.4 NA [26,] 117.0 NA [27,] 141.7 NA [28,] 120.0 NA [29,] 132.1 NA [30,] 146.7 NA [31,] 122.5 NA [32,] 99.6 NA [33,] 122.7 NA [34,] 139.0 NA [35,] 117.8 NA [36,] 125.5 NA [37,] 134.5 NA [38,] 121.3 NA [39,] 126.7 NA [40,] 117.7 NA [41,] 123.0 NA [42,] 132.1 NA [43,] 113.1 NA [44,] 89.2 NA [45,] 121.7 NA [46,] 105.3 NA [47,] 85.3 NA [48,] 105.3 NA [49,] 72.2 NA [50,] 92.1 NA [51,] 97.2 NA [52,] 78.6 NA [53,] 78.1 NA [54,] 93.0 NA [55,] 81.0 NA [56,] 65.9 NA [57,] 88.6 NA [58,] 85.7 NA [59,] 76.3 NA [60,] 96.8 NA [61,] 76.8 NA [62,] 85.6 NA [63,] 119.2 NA [64,] 91.4 NA [65,] 95.7 NA [66,] 112.3 NA [67,] 95.2 NA [68,] 82.8 NA [69,] 111.3 NA [70,] 108.2 NA [71,] 97.0 NA [72,] 124.4 NA [73,] 99.3 NA [74,] 117.6 NA [75,] 131.5 NA [76,] 114.2 NA [77,] 116.8 NA [78,] 116.5 NA [79,] 105.4 NA [80,] 89.2 NA [81,] 115.8 NA [82,] 111.4 NA [83,] 106.4 NA [84,] 128.4 NA [85,] 107.7 NA [86,] 111.0 NA [87,] 129.8 NA [88,] 130.5 NA [89,] 142.9 NA [90,] 159.9 NA [91,] 84.1 NA [92,] 75.0 NA [93,] 100.7 NA [94,] 106.8 NA [95,] 97.4 NA [96,] 113.0 NA [97,] 76.9 NA [98,] 87.3 NA [99,] 103.7 NA [100,] 92.1 NA [101,] 92.9 NA [102,] 112.2 NA [103,] 88.7 NA [104,] 74.6 NA [105,] 101.5 NA [106,] 119.7 NA [107,] 120.7 NA [108,] 153.5 NA > darr [,1] [,2] [1,] 9.9 NA [2,] 23.5 NA [3,] -9.5 NA [4,] -6.0 NA [5,] 26.1 NA [6,] -34.4 NA [7,] -16.3 NA [8,] 37.5 NA [9,] -10.8 NA [10,] -12.0 NA [11,] 46.9 NA [12,] -40.2 NA [13,] 8.9 NA [14,] 11.5 NA [15,] -12.2 NA [16,] 15.2 NA [17,] 8.2 NA [18,] -32.5 NA [19,] -11.3 NA [20,] 30.4 NA [21,] -10.1 NA [22,] -6.6 NA [23,] 38.2 NA [24,] -36.6 NA [25,] 5.6 NA [26,] 24.7 NA [27,] -21.7 NA [28,] 12.1 NA [29,] 14.6 NA [30,] -24.2 NA [31,] -22.9 NA [32,] 23.1 NA [33,] 16.3 NA [34,] -21.2 NA [35,] 7.7 NA [36,] 9.0 NA [37,] -13.2 NA [38,] 5.4 NA [39,] -9.0 NA [40,] 5.3 NA [41,] 9.1 NA [42,] -19.0 NA [43,] -23.9 NA [44,] 32.5 NA [45,] -16.4 NA [46,] -20.0 NA [47,] 20.0 NA [48,] -33.1 NA [49,] 19.9 NA [50,] 5.1 NA [51,] -18.6 NA [52,] -0.5 NA [53,] 14.9 NA [54,] -12.0 NA [55,] -15.1 NA [56,] 22.7 NA [57,] -2.9 NA [58,] -9.4 NA [59,] 20.5 NA [60,] -20.0 NA [61,] 8.8 NA [62,] 33.6 NA [63,] -27.8 NA [64,] 4.3 NA [65,] 16.6 NA [66,] -17.1 NA [67,] -12.4 NA [68,] 28.5 NA [69,] -3.1 NA [70,] -11.2 NA [71,] 27.4 NA [72,] -25.1 NA [73,] 18.3 NA [74,] 13.9 NA [75,] -17.3 NA [76,] 2.6 NA [77,] -0.3 NA [78,] -11.1 NA [79,] -16.2 NA [80,] 26.6 NA [81,] -4.4 NA [82,] -5.0 NA [83,] 22.0 NA [84,] -20.7 NA [85,] 3.3 NA [86,] 18.8 NA [87,] 0.7 NA [88,] 12.4 NA [89,] 17.0 NA [90,] -75.8 NA [91,] -9.1 NA [92,] 25.7 NA [93,] 6.1 NA [94,] -9.4 NA [95,] 15.6 NA [96,] -36.1 NA [97,] 10.4 NA [98,] 16.4 NA [99,] -11.6 NA [100,] 0.8 NA [101,] 19.3 NA [102,] -23.5 NA [103,] -14.1 NA [104,] 26.9 NA [105,] 18.2 NA [106,] 1.0 NA [107,] 32.8 NA [108,] 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/wessaorg/rcomp/tmp/186f01424972777.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/22m0d1424972777.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3ums71424972777.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4yrmv1424972777.ps",horizontal=F,onefile=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] [,13] [1,] 93.6 103.5 127 117.5 111.5 137.6 103.2 86.9 124.4 113.6 101.6 148.5 108.3 [2,] 93.6 103.5 127 117.5 111.5 137.6 103.2 86.9 124.4 113.6 101.6 148.5 108.3 [3,] 93.6 103.5 127 117.5 111.5 137.6 103.2 86.9 124.4 113.6 101.6 148.5 108.3 [4,] 93.6 103.5 127 117.5 111.5 137.6 103.2 86.9 124.4 113.6 101.6 148.5 108.3 [5,] 93.6 103.5 127 117.5 111.5 137.6 103.2 86.9 124.4 113.6 101.6 148.5 108.3 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 117.2 128.7 116.5 131.7 139.9 107.4 96.1 126.5 116.4 109.8 148 111.4 [2,] 117.2 128.7 116.5 131.7 139.9 107.4 96.1 126.5 116.4 109.8 148 111.4 [3,] 117.2 128.7 116.5 131.7 139.9 107.4 96.1 126.5 116.4 109.8 148 111.4 [4,] 117.2 128.7 116.5 131.7 139.9 107.4 96.1 126.5 116.4 109.8 148 111.4 [5,] 117.2 128.7 116.5 131.7 139.9 107.4 96.1 126.5 116.4 109.8 148 111.4 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 117 141.7 120 132.1 146.7 122.5 99.6 122.7 139 117.8 125.5 134.5 [2,] 117 141.7 120 132.1 146.7 122.5 99.6 122.7 139 117.8 125.5 134.5 [3,] 117 141.7 120 132.1 146.7 122.5 99.6 122.7 139 117.8 125.5 134.5 [4,] 117 141.7 120 132.1 146.7 122.5 99.6 122.7 139 117.8 125.5 134.5 [5,] 117 141.7 120 132.1 146.7 122.5 99.6 122.7 139 117.8 125.5 134.5 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 121.3 126.7 117.7 123 132.1 113.1 89.2 121.7 105.3 85.3 105.3 72.2 [2,] 121.3 126.7 117.7 123 132.1 113.1 89.2 121.7 105.3 85.3 105.3 72.2 [3,] 121.3 126.7 117.7 123 132.1 113.1 89.2 121.7 105.3 85.3 105.3 72.2 [4,] 121.3 126.7 117.7 123 132.1 113.1 89.2 121.7 105.3 85.3 105.3 72.2 [5,] 121.3 126.7 117.7 123 132.1 113.1 89.2 121.7 105.3 85.3 105.3 72.2 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 92.1 97.2 78.6 78.1 93 81 65.9 88.6 85.7 76.3 96.8 76.8 [2,] 92.1 97.2 78.6 78.1 93 81 65.9 88.6 85.7 76.3 96.8 76.8 [3,] 92.1 97.2 78.6 78.1 93 81 65.9 88.6 85.7 76.3 96.8 76.8 [4,] 92.1 97.2 78.6 78.1 93 81 65.9 88.6 85.7 76.3 96.8 76.8 [5,] 92.1 97.2 78.6 78.1 93 81 65.9 88.6 85.7 76.3 96.8 76.8 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 85.6 119.2 91.4 95.7 112.3 95.2 82.8 111.3 108.2 97 124.4 99.3 [2,] 85.6 119.2 91.4 95.7 112.3 95.2 82.8 111.3 108.2 97 124.4 99.3 [3,] 85.6 119.2 91.4 95.7 112.3 95.2 82.8 111.3 108.2 97 124.4 99.3 [4,] 85.6 119.2 91.4 95.7 112.3 95.2 82.8 111.3 108.2 97 124.4 99.3 [5,] 85.6 119.2 91.4 95.7 112.3 95.2 82.8 111.3 108.2 97 124.4 99.3 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 117.6 131.5 114.2 116.8 116.5 105.4 89.2 115.8 111.4 106.4 128.4 107.7 [2,] 117.6 131.5 114.2 116.8 116.5 105.4 89.2 115.8 111.4 106.4 128.4 107.7 [3,] 117.6 131.5 114.2 116.8 116.5 105.4 89.2 115.8 111.4 106.4 128.4 107.7 [4,] 117.6 131.5 114.2 116.8 116.5 105.4 89.2 115.8 111.4 106.4 128.4 107.7 [5,] 117.6 131.5 114.2 116.8 116.5 105.4 89.2 115.8 111.4 106.4 128.4 107.7 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] 111 129.8 130.5 142.9 159.9 84.1 75 100.7 106.8 97.4 113 76.9 [2,] 111 129.8 130.5 142.9 159.9 84.1 75 100.7 106.8 97.4 113 76.9 [3,] 111 129.8 130.5 142.9 159.9 84.1 75 100.7 106.8 97.4 113 76.9 [4,] 111 129.8 130.5 142.9 159.9 84.1 75 100.7 106.8 97.4 113 76.9 [5,] 111 129.8 130.5 142.9 159.9 84.1 75 100.7 106.8 97.4 113 76.9 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 87.3 103.7 92.1 92.9 112.2 88.7 74.6 101.5 119.7 120.7 153.5 [2,] 87.3 103.7 92.1 92.9 112.2 88.7 74.6 101.5 119.7 120.7 153.5 [3,] 87.3 103.7 92.1 92.9 112.2 88.7 74.6 101.5 119.7 120.7 153.5 [4,] 87.3 103.7 92.1 92.9 112.2 88.7 74.6 101.5 119.7 120.7 153.5 [5,] 87.3 103.7 92.1 92.9 112.2 88.7 74.6 101.5 119.7 120.7 153.5 $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 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 93.6 103.5 127 117.5 111.5 137.6 103.2 86.9 124.4 113.6 101.6 148.5 108.3 [2,] 93.6 103.5 127 117.5 111.5 137.6 103.2 86.9 124.4 113.6 101.6 148.5 108.3 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 117.2 128.7 116.5 131.7 139.9 107.4 96.1 126.5 116.4 109.8 148 111.4 [2,] 117.2 128.7 116.5 131.7 139.9 107.4 96.1 126.5 116.4 109.8 148 111.4 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 117 141.7 120 132.1 146.7 122.5 99.6 122.7 139 117.8 125.5 134.5 [2,] 117 141.7 120 132.1 146.7 122.5 99.6 122.7 139 117.8 125.5 134.5 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 121.3 126.7 117.7 123 132.1 113.1 89.2 121.7 105.3 85.3 105.3 72.2 [2,] 121.3 126.7 117.7 123 132.1 113.1 89.2 121.7 105.3 85.3 105.3 72.2 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 92.1 97.2 78.6 78.1 93 81 65.9 88.6 85.7 76.3 96.8 76.8 [2,] 92.1 97.2 78.6 78.1 93 81 65.9 88.6 85.7 76.3 96.8 76.8 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 85.6 119.2 91.4 95.7 112.3 95.2 82.8 111.3 108.2 97 124.4 99.3 [2,] 85.6 119.2 91.4 95.7 112.3 95.2 82.8 111.3 108.2 97 124.4 99.3 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 117.6 131.5 114.2 116.8 116.5 105.4 89.2 115.8 111.4 106.4 128.4 107.7 [2,] 117.6 131.5 114.2 116.8 116.5 105.4 89.2 115.8 111.4 106.4 128.4 107.7 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] 111 129.8 130.5 142.9 159.9 84.1 75 100.7 106.8 97.4 113 76.9 [2,] 111 129.8 130.5 142.9 159.9 84.1 75 100.7 106.8 97.4 113 76.9 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 87.3 103.7 92.1 92.9 112.2 88.7 74.6 101.5 119.7 120.7 153.5 [2,] 87.3 103.7 92.1 92.9 112.2 88.7 74.6 101.5 119.7 120.7 153.5 $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" > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5ddlw1424972777.ps",horizontal=F,onefile=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,] 9.9 23.5 -9.5 -6 26.1 -34.4 -16.3 37.5 -10.8 -12 46.9 -40.2 8.9 [2,] 9.9 23.5 -9.5 -6 26.1 -34.4 -16.3 37.5 -10.8 -12 46.9 -40.2 8.9 [3,] 9.9 23.5 -9.5 -6 26.1 -34.4 -16.3 37.5 -10.8 -12 46.9 -40.2 8.9 [4,] 9.9 23.5 -9.5 -6 26.1 -34.4 -16.3 37.5 -10.8 -12 46.9 -40.2 8.9 [5,] 9.9 23.5 -9.5 -6 26.1 -34.4 -16.3 37.5 -10.8 -12 46.9 -40.2 8.9 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 11.5 -12.2 15.2 8.2 -32.5 -11.3 30.4 -10.1 -6.6 38.2 -36.6 5.6 [2,] 11.5 -12.2 15.2 8.2 -32.5 -11.3 30.4 -10.1 -6.6 38.2 -36.6 5.6 [3,] 11.5 -12.2 15.2 8.2 -32.5 -11.3 30.4 -10.1 -6.6 38.2 -36.6 5.6 [4,] 11.5 -12.2 15.2 8.2 -32.5 -11.3 30.4 -10.1 -6.6 38.2 -36.6 5.6 [5,] 11.5 -12.2 15.2 8.2 -32.5 -11.3 30.4 -10.1 -6.6 38.2 -36.6 5.6 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 24.7 -21.7 12.1 14.6 -24.2 -22.9 23.1 16.3 -21.2 7.7 9 -13.2 [2,] 24.7 -21.7 12.1 14.6 -24.2 -22.9 23.1 16.3 -21.2 7.7 9 -13.2 [3,] 24.7 -21.7 12.1 14.6 -24.2 -22.9 23.1 16.3 -21.2 7.7 9 -13.2 [4,] 24.7 -21.7 12.1 14.6 -24.2 -22.9 23.1 16.3 -21.2 7.7 9 -13.2 [5,] 24.7 -21.7 12.1 14.6 -24.2 -22.9 23.1 16.3 -21.2 7.7 9 -13.2 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 5.4 -9 5.3 9.1 -19 -23.9 32.5 -16.4 -20 20 -33.1 19.9 [2,] 5.4 -9 5.3 9.1 -19 -23.9 32.5 -16.4 -20 20 -33.1 19.9 [3,] 5.4 -9 5.3 9.1 -19 -23.9 32.5 -16.4 -20 20 -33.1 19.9 [4,] 5.4 -9 5.3 9.1 -19 -23.9 32.5 -16.4 -20 20 -33.1 19.9 [5,] 5.4 -9 5.3 9.1 -19 -23.9 32.5 -16.4 -20 20 -33.1 19.9 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 5.1 -18.6 -0.5 14.9 -12 -15.1 22.7 -2.9 -9.4 20.5 -20 8.8 [2,] 5.1 -18.6 -0.5 14.9 -12 -15.1 22.7 -2.9 -9.4 20.5 -20 8.8 [3,] 5.1 -18.6 -0.5 14.9 -12 -15.1 22.7 -2.9 -9.4 20.5 -20 8.8 [4,] 5.1 -18.6 -0.5 14.9 -12 -15.1 22.7 -2.9 -9.4 20.5 -20 8.8 [5,] 5.1 -18.6 -0.5 14.9 -12 -15.1 22.7 -2.9 -9.4 20.5 -20 8.8 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 33.6 -27.8 4.3 16.6 -17.1 -12.4 28.5 -3.1 -11.2 27.4 -25.1 18.3 [2,] 33.6 -27.8 4.3 16.6 -17.1 -12.4 28.5 -3.1 -11.2 27.4 -25.1 18.3 [3,] 33.6 -27.8 4.3 16.6 -17.1 -12.4 28.5 -3.1 -11.2 27.4 -25.1 18.3 [4,] 33.6 -27.8 4.3 16.6 -17.1 -12.4 28.5 -3.1 -11.2 27.4 -25.1 18.3 [5,] 33.6 -27.8 4.3 16.6 -17.1 -12.4 28.5 -3.1 -11.2 27.4 -25.1 18.3 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 13.9 -17.3 2.6 -0.3 -11.1 -16.2 26.6 -4.4 -5 22 -20.7 3.3 [2,] 13.9 -17.3 2.6 -0.3 -11.1 -16.2 26.6 -4.4 -5 22 -20.7 3.3 [3,] 13.9 -17.3 2.6 -0.3 -11.1 -16.2 26.6 -4.4 -5 22 -20.7 3.3 [4,] 13.9 -17.3 2.6 -0.3 -11.1 -16.2 26.6 -4.4 -5 22 -20.7 3.3 [5,] 13.9 -17.3 2.6 -0.3 -11.1 -16.2 26.6 -4.4 -5 22 -20.7 3.3 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] 18.8 0.7 12.4 17 -75.8 -9.1 25.7 6.1 -9.4 15.6 -36.1 10.4 [2,] 18.8 0.7 12.4 17 -75.8 -9.1 25.7 6.1 -9.4 15.6 -36.1 10.4 [3,] 18.8 0.7 12.4 17 -75.8 -9.1 25.7 6.1 -9.4 15.6 -36.1 10.4 [4,] 18.8 0.7 12.4 17 -75.8 -9.1 25.7 6.1 -9.4 15.6 -36.1 10.4 [5,] 18.8 0.7 12.4 17 -75.8 -9.1 25.7 6.1 -9.4 15.6 -36.1 10.4 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 16.4 -11.6 0.8 19.3 -23.5 -14.1 26.9 18.2 1 32.8 NA [2,] 16.4 -11.6 0.8 19.3 -23.5 -14.1 26.9 18.2 1 32.8 NA [3,] 16.4 -11.6 0.8 19.3 -23.5 -14.1 26.9 18.2 1 32.8 NA [4,] 16.4 -11.6 0.8 19.3 -23.5 -14.1 26.9 18.2 1 32.8 NA [5,] 16.4 -11.6 0.8 19.3 -23.5 -14.1 26.9 18.2 1 32.8 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 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 9.9 23.5 -9.5 -6 26.1 -34.4 -16.3 37.5 -10.8 -12 46.9 -40.2 8.9 [2,] 9.9 23.5 -9.5 -6 26.1 -34.4 -16.3 37.5 -10.8 -12 46.9 -40.2 8.9 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 11.5 -12.2 15.2 8.2 -32.5 -11.3 30.4 -10.1 -6.6 38.2 -36.6 5.6 [2,] 11.5 -12.2 15.2 8.2 -32.5 -11.3 30.4 -10.1 -6.6 38.2 -36.6 5.6 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 24.7 -21.7 12.1 14.6 -24.2 -22.9 23.1 16.3 -21.2 7.7 9 -13.2 [2,] 24.7 -21.7 12.1 14.6 -24.2 -22.9 23.1 16.3 -21.2 7.7 9 -13.2 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 5.4 -9 5.3 9.1 -19 -23.9 32.5 -16.4 -20 20 -33.1 19.9 [2,] 5.4 -9 5.3 9.1 -19 -23.9 32.5 -16.4 -20 20 -33.1 19.9 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 5.1 -18.6 -0.5 14.9 -12 -15.1 22.7 -2.9 -9.4 20.5 -20 8.8 [2,] 5.1 -18.6 -0.5 14.9 -12 -15.1 22.7 -2.9 -9.4 20.5 -20 8.8 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 33.6 -27.8 4.3 16.6 -17.1 -12.4 28.5 -3.1 -11.2 27.4 -25.1 18.3 [2,] 33.6 -27.8 4.3 16.6 -17.1 -12.4 28.5 -3.1 -11.2 27.4 -25.1 18.3 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 13.9 -17.3 2.6 -0.3 -11.1 -16.2 26.6 -4.4 -5 22 -20.7 3.3 [2,] 13.9 -17.3 2.6 -0.3 -11.1 -16.2 26.6 -4.4 -5 22 -20.7 3.3 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] 18.8 0.7 12.4 17 -75.8 -9.1 25.7 6.1 -9.4 15.6 -36.1 10.4 [2,] 18.8 0.7 12.4 17 -75.8 -9.1 25.7 6.1 -9.4 15.6 -36.1 10.4 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 16.4 -11.6 0.8 19.3 -23.5 -14.1 26.9 18.2 1 32.8 NA [2,] 16.4 -11.6 0.8 19.3 -23.5 -14.1 26.9 18.2 1 32.8 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" > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/67agp1424972777.ps",horizontal=F,onefile=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,] 65.90 NA [2,] 94.40 NA [3,] 111.35 NA [4,] 122.60 NA [5,] 159.90 NA $n [1] 108 0 $conf [,1] [,2] [1,] 107.0626 NA [2,] 115.6374 NA $out numeric(0) $group numeric(0) $names [1] "1" NA > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7fd9y1424972777.ps",horizontal=F,onefile=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,] 65.90 65.90 65.90 [2,] 94.40 94.40 94.40 [3,] 111.35 111.35 111.35 [4,] 122.60 122.60 122.60 [5,] 159.90 159.90 159.90 $n [1] 108 108 108 $conf [,1] [,2] [,3] [1,] 107.0626 107.0626 107.0626 [2,] 115.6374 115.6374 115.6374 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > try(system("convert tmp/186f01424972777.ps tmp/186f01424972777.png",intern=TRUE)) character(0) > try(system("convert tmp/22m0d1424972777.ps tmp/22m0d1424972777.png",intern=TRUE)) character(0) > try(system("convert tmp/3ums71424972777.ps tmp/3ums71424972777.png",intern=TRUE)) character(0) > try(system("convert tmp/4yrmv1424972777.ps tmp/4yrmv1424972777.png",intern=TRUE)) character(0) > try(system("convert tmp/5ddlw1424972777.ps tmp/5ddlw1424972777.png",intern=TRUE)) character(0) > try(system("convert tmp/67agp1424972777.ps tmp/67agp1424972777.png",intern=TRUE)) character(0) > try(system("convert tmp/7fd9y1424972777.ps tmp/7fd9y1424972777.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.798 0.422 3.242