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 = '144' > par1 <- '144' > #'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 2 2 2 2 2 2 2 2 2 2 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 [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 > arr [,1] [,2] [1,] 254 443 [2,] 200 165 [3,] 165 275 [4,] 123 304 [5,] 162 303 [6,] 145 342 [7,] 145 322 [8,] 161 291 [9,] 155 300 [10,] 173 491 [11,] 160 266 [12,] 47 176 [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 > darr [,1] [,2] [1,] -54 -278 [2,] -35 110 [3,] -42 29 [4,] 39 -1 [5,] -17 39 [6,] 0 -20 [7,] 16 -31 [8,] -6 9 [9,] 18 191 [10,] -13 -225 [11,] -113 -90 [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 > 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/1pb1q1445015189.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/260rd1445015189.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/3muog1445015189.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/49gqy1445015189.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] [1,] 254.0 165.0 165 123.0 162.0 145.0 145.0 161 155.0 173 160 47.0 [2,] 254.0 165.0 165 123.0 162.0 145.0 145.0 161 155.0 173 160 47.0 [3,] 348.5 182.5 220 213.5 232.5 243.5 233.5 226 227.5 332 213 111.5 [4,] 443.0 200.0 275 304.0 303.0 342.0 322.0 291 300.0 491 266 176.0 [5,] 443.0 200.0 275 304.0 303.0 342.0 322.0 291 300.0 491 266 176.0 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 232 143 161 159 243 192 157 143 221 227 132 41 [2,] 232 143 161 159 243 192 157 143 221 227 132 41 [3,] 232 143 161 159 243 192 157 143 221 227 132 41 [4,] 232 143 161 159 243 192 157 143 221 227 132 41 [5,] 232 143 161 159 243 192 157 143 221 227 132 41 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 273 182 188 162 140 186 178 236 202 184 119 16 [2,] 273 182 188 162 140 186 178 236 202 184 119 16 [3,] 273 182 188 162 140 186 178 236 202 184 119 16 [4,] 273 182 188 162 140 186 178 236 202 184 119 16 [5,] 273 182 188 162 140 186 178 236 202 184 119 16 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 340 151 240 235 174 309 174 207 209 171 117 10 [2,] 340 151 240 235 174 309 174 207 209 171 117 10 [3,] 340 151 240 235 174 309 174 207 209 171 117 10 [4,] 340 151 240 235 174 309 174 207 209 171 117 10 [5,] 340 151 240 235 174 309 174 207 209 171 117 10 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 339 139 186 155 153 222 102 107 188 162 185 24 [2,] 339 139 186 155 153 222 102 107 188 162 185 24 [3,] 339 139 186 155 153 222 102 107 188 162 185 24 [4,] 339 139 186 155 153 222 102 107 188 162 185 24 [5,] 339 139 186 155 153 222 102 107 188 162 185 24 [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 394 209 248 254 202 258 215 309 240 258 276 48 [2,] 394 209 248 254 202 258 215 309 240 258 276 48 [3,] 394 209 248 254 202 258 215 309 240 258 276 48 [4,] 394 209 248 254 202 258 215 309 240 258 276 48 [5,] 394 209 248 254 202 258 215 309 240 258 276 48 [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [1,] 455 345 311 346 310 297 300 274 292 304 186 14 [2,] 455 345 311 346 310 297 300 274 292 304 186 14 [3,] 455 345 311 346 310 297 300 274 292 304 186 14 [4,] 455 345 311 346 310 297 300 274 292 304 186 14 [5,] 455 345 311 346 310 297 300 274 292 304 186 14 [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [1,] 321 206 160 217 204 246 234 175 364 328 158 40 [2,] 321 206 160 217 204 246 234 175 364 328 158 40 [3,] 321 206 160 217 204 246 234 175 364 328 158 40 [4,] 321 206 160 217 204 246 234 175 364 328 158 40 [5,] 321 206 160 217 204 246 234 175 364 328 158 40 [,97] [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [1,] 556 193 221 278 230 253 240 252 228 306 206 [2,] 556 193 221 278 230 253 240 252 228 306 206 [3,] 556 193 221 278 230 253 240 252 228 306 206 [4,] 556 193 221 278 230 253 240 252 228 306 206 [5,] 556 193 221 278 230 253 240 252 228 306 206 [,108] [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [1,] 48 557 279 399 364 306 471 293 333 316 [2,] 48 557 279 399 364 306 471 293 333 316 [3,] 48 557 279 399 364 306 471 293 333 316 [4,] 48 557 279 399 364 306 471 293 333 316 [5,] 48 557 279 399 364 306 471 293 333 316 [,118] [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [1,] 329 265 61 679 428 394 352 387 590 177 [2,] 329 265 61 679 428 394 352 387 590 177 [3,] 329 265 61 679 428 394 352 387 590 177 [4,] 329 265 61 679 428 394 352 387 590 177 [5,] 329 265 61 679 428 394 352 387 590 177 [,128] [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [1,] 199 203 255 261 115 537 172 425 244 313 [2,] 199 203 255 261 115 537 172 425 244 313 [3,] 199 203 255 261 115 537 172 425 244 313 [4,] 199 203 255 261 115 537 172 425 244 313 [5,] 199 203 255 261 115 537 172 425 244 313 [,138] [,139] [,140] [,141] [,142] [,143] [,144] [1,] 335 222 223 179 335 286 154 [2,] 335 222 223 179 335 286 154 [3,] 335 222 223 179 335 286 154 [4,] 335 222 223 179 335 286 154 [5,] 335 222 223 179 335 286 154 $n [1] 2 2 2 2 2 2 2 2 2 2 2 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 [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 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 137.3438 143.397 97.10484 11.2816 74.97075 23.40594 35.75052 [2,] 559.6562 221.603 342.89516 415.7184 390.02925 463.59406 431.24948 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 80.76027 65.50184 -23.27873 94.57376 -32.6225 232 143 161 159 [2,] 371.23973 389.49816 687.27873 331.42624 255.6225 232 143 161 159 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] 243 192 157 143 221 227 132 41 273 182 188 162 [2,] 243 192 157 143 221 227 132 41 273 182 188 162 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 140 186 178 236 202 184 119 16 340 151 240 235 [2,] 140 186 178 236 202 184 119 16 340 151 240 235 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [,51] [,52] [1,] 174 309 174 207 209 171 117 10 339 139 186 155 [2,] 174 309 174 207 209 171 117 10 339 139 186 155 [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [1,] 153 222 102 107 188 162 185 24 394 209 248 254 [2,] 153 222 102 107 188 162 185 24 394 209 248 254 [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [,75] [,76] [1,] 202 258 215 309 240 258 276 48 455 345 311 346 [2,] 202 258 215 309 240 258 276 48 455 345 311 346 [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [,87] [,88] [1,] 310 297 300 274 292 304 186 14 321 206 160 217 [2,] 310 297 300 274 292 304 186 14 321 206 160 217 [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [,99] [,100] [1,] 204 246 234 175 364 328 158 40 556 193 221 278 [2,] 204 246 234 175 364 328 158 40 556 193 221 278 [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [,109] [,110] [1,] 230 253 240 252 228 306 206 48 557 279 [2,] 230 253 240 252 228 306 206 48 557 279 [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [,119] [,120] [1,] 399 364 306 471 293 333 316 329 265 61 [2,] 399 364 306 471 293 333 316 329 265 61 [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [,129] [,130] [1,] 679 428 394 352 387 590 177 199 203 255 [2,] 679 428 394 352 387 590 177 199 203 255 [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [,139] [,140] [1,] 261 115 537 172 425 244 313 335 222 223 [2,] 261 115 537 172 425 244 313 335 222 223 [,141] [,142] [,143] [,144] [1,] 179 335 286 154 [2,] 179 335 286 154 $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" Warning message: In bxp(list(stats = c(254, 254, 348.5, 443, 443, 165, 165, 182.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/54u3o1445015189.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,] -278 -35.0 -42.0 -1 -17 -20 -31.0 -6.0 18.0 -225 -113.0 185 -89 [2,] -278 -35.0 -42.0 -1 -17 -20 -31.0 -6.0 18.0 -225 -113.0 185 -89 [3,] -166 37.5 -6.5 19 11 -10 -7.5 1.5 104.5 -119 -101.5 185 -89 [4,] -54 110.0 29.0 39 39 0 16.0 9.0 191.0 -13 -90.0 185 -89 [5,] -54 110.0 29.0 39 39 0 16.0 9.0 191.0 -13 -90.0 185 -89 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 18 -2 84 -51 -35 -14 78 6 -95 -91 232 -91 [2,] 18 -2 84 -51 -35 -14 78 6 -95 -91 232 -91 [3,] 18 -2 84 -51 -35 -14 78 6 -95 -91 232 -91 [4,] 18 -2 84 -51 -35 -14 78 6 -95 -91 232 -91 [5,] 18 -2 84 -51 -35 -14 78 6 -95 -91 232 -91 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 6 -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 [2,] 6 -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 [3,] 6 -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 [4,] 6 -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 [5,] 6 -26 -22 46 -8 58 -34 -18 -65 -103 324 -189 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 89 -5 -61 135 -135 33 2 -38 -54 -107 329 -200 [2,] 89 -5 -61 135 -135 33 2 -38 -54 -107 329 -200 [3,] 89 -5 -61 135 -135 33 2 -38 -54 -107 329 -200 [4,] 89 -5 -61 135 -135 33 2 -38 -54 -107 329 -200 [5,] 89 -5 -61 135 -135 33 2 -38 -54 -107 329 -200 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 47 -31 -2 69 -120 5 81 -26 23 -161 370 -185 [2,] 47 -31 -2 69 -120 5 81 -26 23 -161 370 -185 [3,] 47 -31 -2 69 -120 5 81 -26 23 -161 370 -185 [4,] 47 -31 -2 69 -120 5 81 -26 23 -161 370 -185 [5,] 47 -31 -2 69 -120 5 81 -26 23 -161 370 -185 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 39 6 -52 56 -43 94 -69 18 18 -228 407 -110 [2,] 39 6 -52 56 -43 94 -69 18 18 -228 407 -110 [3,] 39 6 -52 56 -43 94 -69 18 18 -228 407 -110 [4,] 39 6 -52 56 -43 94 -69 18 18 -228 407 -110 [5,] 39 6 -52 56 -43 94 -69 18 18 -228 407 -110 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] -34 35 -36 -13 3 -26 18 12 -118 -172 307 -115 [2,] -34 35 -36 -13 3 -26 18 12 -118 -172 307 -115 [3,] -34 35 -36 -13 3 -26 18 12 -118 -172 307 -115 [4,] -34 35 -36 -13 3 -26 18 12 -118 -172 307 -115 [5,] -34 35 -36 -13 3 -26 18 12 -118 -172 307 -115 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] -46 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 [2,] -46 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 [3,] -46 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 [4,] -46 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 [5,] -46 57 -13 42 -12 -59 189 -36 -170 -118 516 -363 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 28 57 -48 23 -13 12 -24 78 -100 -158 509 [2,] 28 57 -48 23 -13 12 -24 78 -100 -158 509 [3,] 28 57 -48 23 -13 12 -24 78 -100 -158 509 [4,] 28 57 -48 23 -13 12 -24 78 -100 -158 509 [5,] 28 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] [1,] 1 -44 156 -49 -132 289 [2,] 1 -44 156 -49 -132 289 [3,] 1 -44 156 -49 -132 289 [4,] 1 -44 156 -49 -132 289 [5,] 1 -44 156 -49 -132 289 $n [1] 2 2 2 2 2 2 2 2 2 2 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 [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 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -416.25923 -124.4982 -85.82324 -25.68915 -51.56481 -32.34457 -60.00975 [2,] 84.25923 199.4982 72.82324 63.68915 73.56481 12.34457 45.00975 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] -15.25843 -88.78057 -355.8525 -127.19626 185 -89 18 -2 84 [2,] 18.25843 297.78057 117.8525 -75.80374 185 -89 18 -2 84 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] -51 -35 -14 78 6 -95 -91 232 -91 6 -26 -22 [2,] -51 -35 -14 78 6 -95 -91 232 -91 6 -26 -22 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 46 -8 58 -34 -18 -65 -103 324 -189 89 -5 -61 [2,] 46 -8 58 -34 -18 -65 -103 324 -189 89 -5 -61 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [,51] [,52] [1,] 135 -135 33 2 -38 -54 -107 329 -200 47 -31 -2 [2,] 135 -135 33 2 -38 -54 -107 329 -200 47 -31 -2 [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [1,] 69 -120 5 81 -26 23 -161 370 -185 39 6 -52 [2,] 69 -120 5 81 -26 23 -161 370 -185 39 6 -52 [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [,75] [,76] [1,] 56 -43 94 -69 18 18 -228 407 -110 -34 35 -36 [2,] 56 -43 94 -69 18 18 -228 407 -110 -34 35 -36 [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [,87] [,88] [1,] -13 3 -26 18 12 -118 -172 307 -115 -46 57 -13 [2,] -13 3 -26 18 12 -118 -172 307 -115 -46 57 -13 [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [,99] [,100] [1,] 42 -12 -59 189 -36 -170 -118 516 -363 28 57 -48 [2,] 42 -12 -59 189 -36 -170 -118 516 -363 28 57 -48 [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [,109] [,110] [1,] 23 -13 12 -24 78 -100 -158 509 -278 120 [2,] 23 -13 12 -24 78 -100 -158 509 -278 120 [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [,119] [,120] [1,] -35 -58 165 -178 40 -17 13 -64 -204 618 [2,] -35 -58 165 -178 40 -17 13 -64 -204 618 [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [,129] [,130] [1,] -251 -34 -42 35 203 -413 22 4 52 6 [2,] -251 -34 -42 35 203 -413 22 4 52 6 [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [,139] [,140] [1,] -146 422 -365 253 -181 69 22 -113 1 -44 [2,] -146 422 -365 253 -181 69 22 -113 1 -44 [,141] [,142] [,143] [,144] [1,] 156 -49 -132 289 [2,] 156 -49 -132 289 $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" Warning message: In bxp(list(stats = c(-278, -278, -166, -54, -54, -35, -35, 37.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6xotw1445015189.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,] 10.0 266.0 [2,] 162.0 270.5 [3,] 221.0 301.5 [4,] 298.5 332.0 [5,] 471.0 342.0 $n [1] 144 12 $conf [,1] [,2] [1,] 203.0275 273.4494 [2,] 238.9725 329.5506 $out [1] 556 557 679 590 537 443 165 491 176 $group [1] 1 1 1 1 1 2 2 2 2 $names [1] "1" NA > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7x0an1445015189.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,] 10.00 10.00 10.00 [2,] 176.00 176.00 176.00 [3,] 227.25 227.25 227.25 [4,] 305.00 305.00 305.00 [5,] 471.00 471.00 471.00 $n [1] 144 144 144 $conf [,1] [,2] [,3] [1,] 210.265 210.265 210.265 [2,] 244.235 244.235 244.235 $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/1pb1q1445015189.ps tmp/1pb1q1445015189.png",intern=TRUE)) character(0) > try(system("convert tmp/260rd1445015189.ps tmp/260rd1445015189.png",intern=TRUE)) character(0) > try(system("convert tmp/3muog1445015189.ps tmp/3muog1445015189.png",intern=TRUE)) character(0) > try(system("convert tmp/49gqy1445015189.ps tmp/49gqy1445015189.png",intern=TRUE)) character(0) > try(system("convert tmp/54u3o1445015189.ps tmp/54u3o1445015189.png",intern=TRUE)) character(0) > try(system("convert tmp/6xotw1445015189.ps tmp/6xotw1445015189.png",intern=TRUE)) character(0) > try(system("convert tmp/7x0an1445015189.ps tmp/7x0an1445015189.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.736 0.472 3.239