R version 3.2.2 (2015-08-14) -- "Fire Safety" Copyright (C) 2015 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(254,200,165,123,162,145,145,161,155,173,160,47,232,143,161,159,243,192,157,143,221,227,132,41,273,182,188,162,140,186,178,236,202,184,119,16,340,151,240,235,174,309,174,207,209,171,117,10,339,139,186,155,153,222,102,107,188,162,185,24,394,209,248,254,202,258,215,309,240,258,276,48,455,345,311,346,310,297,300,274,292,304,186,14,321,206,160,217,204,246,234,175,364,328,158,40,556,193,221,278,230,253,240,252,228,306,206,48,557,279,399,364,306,471,293,333,316,329,265,61,679,428,394,352,387,590,177,199,203,255,261,115,537,172,425,244,313,335,222,223,179,335,286,154,443,165,275,304,303,342,322,291,300,491,266,176) > par1 = '155' > par1 <- '155' > #'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] 156 > (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 1 1 1 1 1 1 1 1 1 1 1 1 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 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [149] 1 1 1 1 1 1 1 > arr [,1] [,2] [1,] 254 176 [2,] 200 NA [3,] 165 NA [4,] 123 NA [5,] 162 NA [6,] 145 NA [7,] 145 NA [8,] 161 NA [9,] 155 NA [10,] 173 NA [11,] 160 NA [12,] 47 NA [13,] 232 NA [14,] 143 NA [15,] 161 NA [16,] 159 NA [17,] 243 NA [18,] 192 NA [19,] 157 NA [20,] 143 NA [21,] 221 NA [22,] 227 NA [23,] 132 NA [24,] 41 NA [25,] 273 NA [26,] 182 NA [27,] 188 NA [28,] 162 NA [29,] 140 NA [30,] 186 NA [31,] 178 NA [32,] 236 NA [33,] 202 NA [34,] 184 NA [35,] 119 NA [36,] 16 NA [37,] 340 NA [38,] 151 NA [39,] 240 NA [40,] 235 NA [41,] 174 NA [42,] 309 NA [43,] 174 NA [44,] 207 NA [45,] 209 NA [46,] 171 NA [47,] 117 NA [48,] 10 NA [49,] 339 NA [50,] 139 NA [51,] 186 NA [52,] 155 NA [53,] 153 NA [54,] 222 NA [55,] 102 NA [56,] 107 NA [57,] 188 NA [58,] 162 NA [59,] 185 NA [60,] 24 NA [61,] 394 NA [62,] 209 NA [63,] 248 NA [64,] 254 NA [65,] 202 NA [66,] 258 NA [67,] 215 NA [68,] 309 NA [69,] 240 NA [70,] 258 NA [71,] 276 NA [72,] 48 NA [73,] 455 NA [74,] 345 NA [75,] 311 NA [76,] 346 NA [77,] 310 NA [78,] 297 NA [79,] 300 NA [80,] 274 NA [81,] 292 NA [82,] 304 NA [83,] 186 NA [84,] 14 NA [85,] 321 NA [86,] 206 NA [87,] 160 NA [88,] 217 NA [89,] 204 NA [90,] 246 NA [91,] 234 NA [92,] 175 NA [93,] 364 NA [94,] 328 NA [95,] 158 NA [96,] 40 NA [97,] 556 NA [98,] 193 NA [99,] 221 NA [100,] 278 NA [101,] 230 NA [102,] 253 NA [103,] 240 NA [104,] 252 NA [105,] 228 NA [106,] 306 NA [107,] 206 NA [108,] 48 NA [109,] 557 NA [110,] 279 NA [111,] 399 NA [112,] 364 NA [113,] 306 NA [114,] 471 NA [115,] 293 NA [116,] 333 NA [117,] 316 NA [118,] 329 NA [119,] 265 NA [120,] 61 NA [121,] 679 NA [122,] 428 NA [123,] 394 NA [124,] 352 NA [125,] 387 NA [126,] 590 NA [127,] 177 NA [128,] 199 NA [129,] 203 NA [130,] 255 NA [131,] 261 NA [132,] 115 NA [133,] 537 NA [134,] 172 NA [135,] 425 NA [136,] 244 NA [137,] 313 NA [138,] 335 NA [139,] 222 NA [140,] 223 NA [141,] 179 NA [142,] 335 NA [143,] 286 NA [144,] 154 NA [145,] 443 NA [146,] 165 NA [147,] 275 NA [148,] 304 NA [149,] 303 NA [150,] 342 NA [151,] 322 NA [152,] 291 NA [153,] 300 NA [154,] 491 NA [155,] 266 NA > darr [,1] [,2] [1,] -54 NA [2,] -35 NA [3,] -42 NA [4,] 39 NA [5,] -17 NA [6,] 0 NA [7,] 16 NA [8,] -6 NA [9,] 18 NA [10,] -13 NA [11,] -113 NA [12,] 185 NA [13,] -89 NA [14,] 18 NA [15,] -2 NA [16,] 84 NA [17,] -51 NA [18,] -35 NA [19,] -14 NA [20,] 78 NA [21,] 6 NA [22,] -95 NA [23,] -91 NA [24,] 232 NA [25,] -91 NA [26,] 6 NA [27,] -26 NA [28,] -22 NA [29,] 46 NA [30,] -8 NA [31,] 58 NA [32,] -34 NA [33,] -18 NA [34,] -65 NA [35,] -103 NA [36,] 324 NA [37,] -189 NA [38,] 89 NA [39,] -5 NA [40,] -61 NA [41,] 135 NA [42,] -135 NA [43,] 33 NA [44,] 2 NA [45,] -38 NA [46,] -54 NA [47,] -107 NA [48,] 329 NA [49,] -200 NA [50,] 47 NA [51,] -31 NA [52,] -2 NA [53,] 69 NA [54,] -120 NA [55,] 5 NA [56,] 81 NA [57,] -26 NA [58,] 23 NA [59,] -161 NA [60,] 370 NA [61,] -185 NA [62,] 39 NA [63,] 6 NA [64,] -52 NA [65,] 56 NA [66,] -43 NA [67,] 94 NA [68,] -69 NA [69,] 18 NA [70,] 18 NA [71,] -228 NA [72,] 407 NA [73,] -110 NA [74,] -34 NA [75,] 35 NA [76,] -36 NA [77,] -13 NA [78,] 3 NA [79,] -26 NA [80,] 18 NA [81,] 12 NA [82,] -118 NA [83,] -172 NA [84,] 307 NA [85,] -115 NA [86,] -46 NA [87,] 57 NA [88,] -13 NA [89,] 42 NA [90,] -12 NA [91,] -59 NA [92,] 189 NA [93,] -36 NA [94,] -170 NA [95,] -118 NA [96,] 516 NA [97,] -363 NA [98,] 28 NA [99,] 57 NA [100,] -48 NA [101,] 23 NA [102,] -13 NA [103,] 12 NA [104,] -24 NA [105,] 78 NA [106,] -100 NA [107,] -158 NA [108,] 509 NA [109,] -278 NA [110,] 120 NA [111,] -35 NA [112,] -58 NA [113,] 165 NA [114,] -178 NA [115,] 40 NA [116,] -17 NA [117,] 13 NA [118,] -64 NA [119,] -204 NA [120,] 618 NA [121,] -251 NA [122,] -34 NA [123,] -42 NA [124,] 35 NA [125,] 203 NA [126,] -413 NA [127,] 22 NA [128,] 4 NA [129,] 52 NA [130,] 6 NA [131,] -146 NA [132,] 422 NA [133,] -365 NA [134,] 253 NA [135,] -181 NA [136,] 69 NA [137,] 22 NA [138,] -113 NA [139,] 1 NA [140,] -44 NA [141,] 156 NA [142,] -49 NA [143,] -132 NA [144,] 289 NA [145,] -278 NA [146,] 110 NA [147,] 29 NA [148,] -1 NA [149,] 39 NA [150,] -20 NA [151,] -31 NA [152,] 9 NA [153,] 191 NA [154,] -225 NA [155,] -90 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/14mfj1445024274.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=13.888888888889,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/2jj6d1445024274.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=13.888888888889,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/3bdct1445024274.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=13.888888888889,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/4dtx51445024274.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=13.888888888889,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] [,14] [1,] 176 200 165 123 162 145 145 161 155 173 160 47 232 143 [2,] 176 200 165 123 162 145 145 161 155 173 160 47 232 143 [3,] 215 200 165 123 162 145 145 161 155 173 160 47 232 143 [4,] 254 200 165 123 162 145 145 161 155 173 160 47 232 143 [5,] 254 200 165 123 162 145 145 161 155 173 160 47 232 143 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 161 159 243 192 157 143 221 227 132 41 273 182 [2,] 161 159 243 192 157 143 221 227 132 41 273 182 [3,] 161 159 243 192 157 143 221 227 132 41 273 182 [4,] 161 159 243 192 157 143 221 227 132 41 273 182 [5,] 161 159 243 192 157 143 221 227 132 41 273 182 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] 188 162 140 186 178 236 202 184 119 16 340 151 [2,] 188 162 140 186 178 236 202 184 119 16 340 151 [3,] 188 162 140 186 178 236 202 184 119 16 340 151 [4,] 188 162 140 186 178 236 202 184 119 16 340 151 [5,] 188 162 140 186 178 236 202 184 119 16 340 151 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 240 235 174 309 174 207 209 171 117 10 339 139 [2,] 240 235 174 309 174 207 209 171 117 10 339 139 [3,] 240 235 174 309 174 207 209 171 117 10 339 139 [4,] 240 235 174 309 174 207 209 171 117 10 339 139 [5,] 240 235 174 309 174 207 209 171 117 10 339 139 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 186 155 153 222 102 107 188 162 185 24 394 209 [2,] 186 155 153 222 102 107 188 162 185 24 394 209 [3,] 186 155 153 222 102 107 188 162 185 24 394 209 [4,] 186 155 153 222 102 107 188 162 185 24 394 209 [5,] 186 155 153 222 102 107 188 162 185 24 394 209 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] 248 254 202 258 215 309 240 258 276 48 455 345 [2,] 248 254 202 258 215 309 240 258 276 48 455 345 [3,] 248 254 202 258 215 309 240 258 276 48 455 345 [4,] 248 254 202 258 215 309 240 258 276 48 455 345 [5,] 248 254 202 258 215 309 240 258 276 48 455 345 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] 311 346 310 297 300 274 292 304 186 14 321 206 [2,] 311 346 310 297 300 274 292 304 186 14 321 206 [3,] 311 346 310 297 300 274 292 304 186 14 321 206 [4,] 311 346 310 297 300 274 292 304 186 14 321 206 [5,] 311 346 310 297 300 274 292 304 186 14 321 206 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 160 217 204 246 234 175 364 328 158 40 556 193 [2,] 160 217 204 246 234 175 364 328 158 40 556 193 [3,] 160 217 204 246 234 175 364 328 158 40 556 193 [4,] 160 217 204 246 234 175 364 328 158 40 556 193 [5,] 160 217 204 246 234 175 364 328 158 40 556 193 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 221 278 230 253 240 252 228 306 206 48 [2,] 221 278 230 253 240 252 228 306 206 48 [3,] 221 278 230 253 240 252 228 306 206 48 [4,] 221 278 230 253 240 252 228 306 206 48 [5,] 221 278 230 253 240 252 228 306 206 48 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 557 279 399 364 306 471 293 333 316 329 [2,] 557 279 399 364 306 471 293 333 316 329 [3,] 557 279 399 364 306 471 293 333 316 329 [4,] 557 279 399 364 306 471 293 333 316 329 [5,] 557 279 399 364 306 471 293 333 316 329 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] 265 61 679 428 394 352 387 590 177 199 [2,] 265 61 679 428 394 352 387 590 177 199 [3,] 265 61 679 428 394 352 387 590 177 199 [4,] 265 61 679 428 394 352 387 590 177 199 [5,] 265 61 679 428 394 352 387 590 177 199 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] 203 255 261 115 537 172 425 244 313 335 [2,] 203 255 261 115 537 172 425 244 313 335 [3,] 203 255 261 115 537 172 425 244 313 335 [4,] 203 255 261 115 537 172 425 244 313 335 [5,] 203 255 261 115 537 172 425 244 313 335 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] 222 223 179 335 286 154 443 165 275 304 [2,] 222 223 179 335 286 154 443 165 275 304 [3,] 222 223 179 335 286 154 443 165 275 304 [4,] 222 223 179 335 286 154 443 165 275 304 [5,] 222 223 179 335 286 154 443 165 275 304 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [1,] 303 342 322 291 300 491 266 [2,] 303 342 322 291 300 491 266 [3,] 303 342 322 291 300 491 266 [4,] 303 342 322 291 300 491 266 [5,] 303 342 322 291 300 491 266 $n [1] 2 1 1 1 1 1 1 1 1 1 1 1 1 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 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [149] 1 1 1 1 1 1 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 127.8562 200 165 123 162 145 145 161 155 173 160 47 232 [2,] 302.1438 200 165 123 162 145 145 161 155 173 160 47 232 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 143 161 159 243 192 157 143 221 227 132 41 273 [2,] 143 161 159 243 192 157 143 221 227 132 41 273 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 182 188 162 140 186 178 236 202 184 119 16 340 [2,] 182 188 162 140 186 178 236 202 184 119 16 340 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 151 240 235 174 309 174 207 209 171 117 10 339 [2,] 151 240 235 174 309 174 207 209 171 117 10 339 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 139 186 155 153 222 102 107 188 162 185 24 394 [2,] 139 186 155 153 222 102 107 188 162 185 24 394 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 209 248 254 202 258 215 309 240 258 276 48 455 [2,] 209 248 254 202 258 215 309 240 258 276 48 455 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 345 311 346 310 297 300 274 292 304 186 14 321 [2,] 345 311 346 310 297 300 274 292 304 186 14 321 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] 206 160 217 204 246 234 175 364 328 158 40 556 [2,] 206 160 217 204 246 234 175 364 328 158 40 556 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 193 221 278 230 253 240 252 228 306 206 48 [2,] 193 221 278 230 253 240 252 228 306 206 48 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 557 279 399 364 306 471 293 333 316 329 [2,] 557 279 399 364 306 471 293 333 316 329 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] 265 61 679 428 394 352 387 590 177 199 [2,] 265 61 679 428 394 352 387 590 177 199 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] 203 255 261 115 537 172 425 244 313 335 [2,] 203 255 261 115 537 172 425 244 313 335 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] 222 223 179 335 286 154 443 165 275 304 [2,] 222 223 179 335 286 154 443 165 275 304 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [1,] 303 342 322 291 300 491 266 [2,] 303 342 322 291 300 491 266 $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" [121] "121" "122" "123" "124" "125" "126" "127" "128" "129" "130" "131" "132" [133] "133" "134" "135" "136" "137" "138" "139" "140" "141" "142" "143" "144" [145] "145" "146" "147" "148" "149" "150" "151" "152" "153" "154" "155" Warning message: In bxp(list(stats = c(176, 176, 215, 254, 254, 200, 200, 200, 200, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5bts21445024274.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=13.888888888889,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,] -54 -35 -42 39 -17 0 16 -6 18 -13 -113 185 -89 18 [2,] -54 -35 -42 39 -17 0 16 -6 18 -13 -113 185 -89 18 [3,] -54 -35 -42 39 -17 0 16 -6 18 -13 -113 185 -89 18 [4,] -54 -35 -42 39 -17 0 16 -6 18 -13 -113 185 -89 18 [5,] -54 -35 -42 39 -17 0 16 -6 18 -13 -113 185 -89 18 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] -2 84 -51 -35 -14 78 6 -95 -91 232 -91 6 [2,] -2 84 -51 -35 -14 78 6 -95 -91 232 -91 6 [3,] -2 84 -51 -35 -14 78 6 -95 -91 232 -91 6 [4,] -2 84 -51 -35 -14 78 6 -95 -91 232 -91 6 [5,] -2 84 -51 -35 -14 78 6 -95 -91 232 -91 6 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 89 [2,] -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 89 [3,] -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 89 [4,] -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 89 [5,] -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 89 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] -5 -61 135 -135 33 2 -38 -54 -107 329 -200 47 [2,] -5 -61 135 -135 33 2 -38 -54 -107 329 -200 47 [3,] -5 -61 135 -135 33 2 -38 -54 -107 329 -200 47 [4,] -5 -61 135 -135 33 2 -38 -54 -107 329 -200 47 [5,] -5 -61 135 -135 33 2 -38 -54 -107 329 -200 47 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] -31 -2 69 -120 5 81 -26 23 -161 370 -185 39 [2,] -31 -2 69 -120 5 81 -26 23 -161 370 -185 39 [3,] -31 -2 69 -120 5 81 -26 23 -161 370 -185 39 [4,] -31 -2 69 -120 5 81 -26 23 -161 370 -185 39 [5,] -31 -2 69 -120 5 81 -26 23 -161 370 -185 39 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] 6 -52 56 -43 94 -69 18 18 -228 407 -110 -34 [2,] 6 -52 56 -43 94 -69 18 18 -228 407 -110 -34 [3,] 6 -52 56 -43 94 -69 18 18 -228 407 -110 -34 [4,] 6 -52 56 -43 94 -69 18 18 -228 407 -110 -34 [5,] 6 -52 56 -43 94 -69 18 18 -228 407 -110 -34 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] 35 -36 -13 3 -26 18 12 -118 -172 307 -115 -46 [2,] 35 -36 -13 3 -26 18 12 -118 -172 307 -115 -46 [3,] 35 -36 -13 3 -26 18 12 -118 -172 307 -115 -46 [4,] 35 -36 -13 3 -26 18 12 -118 -172 307 -115 -46 [5,] 35 -36 -13 3 -26 18 12 -118 -172 307 -115 -46 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 28 [2,] 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 28 [3,] 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 28 [4,] 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 28 [5,] 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 28 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 57 -48 23 -13 12 -24 78 -100 -158 509 [2,] 57 -48 23 -13 12 -24 78 -100 -158 509 [3,] 57 -48 23 -13 12 -24 78 -100 -158 509 [4,] 57 -48 23 -13 12 -24 78 -100 -158 509 [5,] 57 -48 23 -13 12 -24 78 -100 -158 509 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] -278 120 -35 -58 165 -178 40 -17 13 -64 [2,] -278 120 -35 -58 165 -178 40 -17 13 -64 [3,] -278 120 -35 -58 165 -178 40 -17 13 -64 [4,] -278 120 -35 -58 165 -178 40 -17 13 -64 [5,] -278 120 -35 -58 165 -178 40 -17 13 -64 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] -204 618 -251 -34 -42 35 203 -413 22 4 [2,] -204 618 -251 -34 -42 35 203 -413 22 4 [3,] -204 618 -251 -34 -42 35 203 -413 22 4 [4,] -204 618 -251 -34 -42 35 203 -413 22 4 [5,] -204 618 -251 -34 -42 35 203 -413 22 4 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] 52 6 -146 422 -365 253 -181 69 22 -113 [2,] 52 6 -146 422 -365 253 -181 69 22 -113 [3,] 52 6 -146 422 -365 253 -181 69 22 -113 [4,] 52 6 -146 422 -365 253 -181 69 22 -113 [5,] 52 6 -146 422 -365 253 -181 69 22 -113 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] 1 -44 156 -49 -132 289 -278 110 29 -1 [2,] 1 -44 156 -49 -132 289 -278 110 29 -1 [3,] 1 -44 156 -49 -132 289 -278 110 29 -1 [4,] 1 -44 156 -49 -132 289 -278 110 29 -1 [5,] 1 -44 156 -49 -132 289 -278 110 29 -1 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [1,] 39 -20 -31 9 191 -225 -90 [2,] 39 -20 -31 9 191 -225 -90 [3,] 39 -20 -31 9 191 -225 -90 [4,] 39 -20 -31 9 191 -225 -90 [5,] 39 -20 -31 9 191 -225 -90 $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 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [149] 1 1 1 1 1 1 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] -54 -35 -42 39 -17 0 16 -6 18 -13 -113 185 -89 18 [2,] -54 -35 -42 39 -17 0 16 -6 18 -13 -113 185 -89 18 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] -2 84 -51 -35 -14 78 6 -95 -91 232 -91 6 [2,] -2 84 -51 -35 -14 78 6 -95 -91 232 -91 6 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 89 [2,] -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 89 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] -5 -61 135 -135 33 2 -38 -54 -107 329 -200 47 [2,] -5 -61 135 -135 33 2 -38 -54 -107 329 -200 47 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] -31 -2 69 -120 5 81 -26 23 -161 370 -185 39 [2,] -31 -2 69 -120 5 81 -26 23 -161 370 -185 39 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] 6 -52 56 -43 94 -69 18 18 -228 407 -110 -34 [2,] 6 -52 56 -43 94 -69 18 18 -228 407 -110 -34 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] 35 -36 -13 3 -26 18 12 -118 -172 307 -115 -46 [2,] 35 -36 -13 3 -26 18 12 -118 -172 307 -115 -46 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 28 [2,] 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 28 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 57 -48 23 -13 12 -24 78 -100 -158 509 [2,] 57 -48 23 -13 12 -24 78 -100 -158 509 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] -278 120 -35 -58 165 -178 40 -17 13 -64 [2,] -278 120 -35 -58 165 -178 40 -17 13 -64 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] -204 618 -251 -34 -42 35 203 -413 22 4 [2,] -204 618 -251 -34 -42 35 203 -413 22 4 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] 52 6 -146 422 -365 253 -181 69 22 -113 [2,] 52 6 -146 422 -365 253 -181 69 22 -113 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] 1 -44 156 -49 -132 289 -278 110 29 -1 [2,] 1 -44 156 -49 -132 289 -278 110 29 -1 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [1,] 39 -20 -31 9 191 -225 -90 [2,] 39 -20 -31 9 191 -225 -90 $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" [121] "121" "122" "123" "124" "125" "126" "127" "128" "129" "130" "131" "132" [133] "133" "134" "135" "136" "137" "138" "139" "140" "141" "142" "143" "144" [145] "145" "146" "147" "148" "149" "150" "151" "152" "153" "154" "155" > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6otmu1445024274.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=13.888888888889,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,] 10 176 [2,] 165 176 [3,] 227 176 [4,] 304 176 [5,] 491 176 $n [1] 155 1 $conf [,1] [,2] [1,] 209.3597 176 [2,] 244.6403 176 $out [1] 556 557 679 590 537 $group [1] 1 1 1 1 1 $names [1] "1" NA > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7rf7y1445024274.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=13.888888888889,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,] 10 10 10 [2,] 165 165 165 [3,] 223 223 223 [4,] 304 304 304 [5,] 491 491 491 $n [1] 155 155 155 $conf [,1] [,2] [,3] [1,] 205.3597 205.3597 205.3597 [2,] 240.6403 240.6403 240.6403 $out [1] 556 557 679 590 537 556 557 679 590 537 556 557 679 590 537 $group [1] 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > try(system("convert tmp/14mfj1445024274.ps tmp/14mfj1445024274.png",intern=TRUE)) character(0) > try(system("convert tmp/2jj6d1445024274.ps tmp/2jj6d1445024274.png",intern=TRUE)) character(0) > try(system("convert tmp/3bdct1445024274.ps tmp/3bdct1445024274.png",intern=TRUE)) character(0) > try(system("convert tmp/4dtx51445024274.ps tmp/4dtx51445024274.png",intern=TRUE)) character(0) > try(system("convert tmp/5bts21445024274.ps tmp/5bts21445024274.png",intern=TRUE)) character(0) > try(system("convert tmp/6otmu1445024274.ps tmp/6otmu1445024274.png",intern=TRUE)) character(0) > try(system("convert tmp/7rf7y1445024274.ps tmp/7rf7y1445024274.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 3.008 0.548 3.578