R version 3.3.1 (2016-06-21) -- "Bug in Your Hair" Copyright (C) 2016 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(20 + ,25 + ,15 + ,15 + ,25 + ,25 + ,25 + ,21 + ,30 + ,25 + ,20 + ,40 + ,13 + ,30 + ,25 + ,20 + ,25 + ,20 + ,25 + ,20 + ,20 + ,15 + ,15 + ,12 + ,20 + ,5 + ,20 + ,15 + ,25 + ,22 + ,20 + ,22 + ,25 + ,20 + ,20 + ,35 + ,30 + ,25 + ,20 + ,20 + ,20 + ,25 + ,25 + ,15 + ,20 + ,35 + ,25 + ,25 + ,30 + ,23 + ,10 + ,22 + ,25 + ,25 + ,22 + ,30 + ,20 + ,25 + ,25 + ,22 + ,25 + ,25 + ,25 + ,22 + ,25 + ,12 + ,18 + ,20 + ,20 + ,22 + ,30 + ,25 + ,22 + ,20 + ,50 + ,30 + ,25 + ,20 + ,30 + ,22 + ,25 + ,30 + ,22 + ,25 + ,22 + ,22 + ,25 + ,25 + ,25 + ,20 + ,22 + ,15 + ,20 + ,30 + ,20 + ,25 + ,30 + ,35 + ,22 + ,12 + ,30 + ,15 + ,10 + ,30 + ,9 + ,25 + ,20 + ,20 + ,35 + ,25 + ,35 + ,30 + ,12 + ,25 + ,15 + ,25 + ,25 + ,20 + ,20 + ,6 + ,15 + ,40 + ,20 + ,40 + ,25 + ,25 + ,20 + ,15 + ,15 + ,22 + ,24 + ,22 + ,20 + ,25 + ,25 + ,25 + ,35 + ,40 + ,20 + ,22 + ,22 + ,20 + ,25 + ,25 + ,18 + ,25 + ,20 + ,25 + ,30 + ,20 + ,22 + ,35 + ,22 + ,25 + ,25 + ,25 + ,25 + ,22 + ,23 + ,35 + ,15 + ,25 + ,18 + ,22 + ,25 + ,25 + ,28 + ,30 + ,20 + ,25 + ,25 + ,30 + ,22 + ,30 + ,10 + ,10 + ,25 + ,20 + ,22 + ,25 + ,25 + ,15 + ,22 + ,25 + ,25 + ,28 + ,22 + ,30 + ,25 + ,20 + ,25 + ,25 + ,20 + ,30 + ,20 + ,30 + ,50 + ,19 + ,20 + ,28 + ,20 + ,25 + ,35 + ,25 + ,25 + ,15 + ,16 + ,20 + ,20 + ,25 + ,30 + ,20 + ,25 + ,25 + ,25 + ,20 + ,20 + ,25 + ,25 + ,30 + ,22 + ,20 + ,25 + ,25 + ,18 + ,18 + ,20 + ,25 + ,25 + ,30 + ,25 + ,20 + ,25 + ,20 + ,20 + ,20 + ,22 + ,18 + ,22 + ,20 + ,15 + ,25 + ,25 + ,20 + ,25 + ,15 + ,22 + ,25 + ,25 + ,15 + ,12 + ,25 + ,30 + ,22 + ,15 + ,22 + ,25 + ,12 + ,18 + ,30 + ,25 + ,25 + ,40 + ,24 + ,25 + ,15 + ,25 + ,20 + ,25 + ,25 + ,25 + ,20 + ,30 + ,20 + ,25 + ,30 + ,22 + ,25 + ,25 + ,25 + ,50 + ,19 + ,50 + ,25 + ,35 + ,20 + ,20 + ,20 + ,20 + ,20 + ,25 + ,25 + ,25 + ,20 + ,20 + ,20 + ,20 + ,25 + ,18 + ,25 + ,22 + ,22 + ,30 + ,30 + ,8 + ,20 + ,25 + ,30 + ,50 + ,22 + ,20 + ,10 + ,25 + ,25 + ,25 + ,25 + ,18 + ,25 + ,20 + ,25 + ,30 + ,18 + ,20 + ,25 + ,22 + ,22 + ,20 + ,20 + ,25 + ,20 + ,20 + ,20 + ,20 + ,25 + ,20 + ,10 + ,20 + ,25 + ,30 + ,25 + ,50 + ,30 + ,30 + ,50 + ,15 + ,25 + ,25 + ,22 + ,20 + ,22 + ,30 + ,25 + ,18 + ,22 + ,22 + ,30 + ,40 + ,25 + ,20 + ,10 + ,20 + ,9 + ,15 + ,20 + ,15 + ,20 + ,30 + ,12 + ,15 + ,12 + ,20 + ,15 + ,12 + ,25 + ,20 + ,25 + ,25 + ,25 + ,30 + ,20 + ,25 + ,15 + ,15 + ,22 + ,10 + ,15 + ,10 + ,20 + ,25 + ,20 + ,20 + ,38 + ,20 + ,20 + ,20 + ,40 + ,25 + ,25 + ,30 + ,25 + ,10 + ,20 + ,25 + ,12 + ,15 + ,25 + ,20 + ,22 + ,22 + ,20 + ,25 + ,25 + ,25 + ,15 + ,40 + ,20 + ,20 + ,16 + ,25 + ,15 + ,20 + ,25 + ,20 + ,30 + ,50 + ,20 + ,25 + ,20 + ,30 + ,30 + ,25 + ,25 + ,12 + ,25 + ,25 + ,25 + ,20 + ,20 + ,20 + ,15 + ,20 + ,25 + ,15 + ,25 + ,50 + ,30 + ,20 + ,20 + ,25 + ,12 + ,15 + ,20 + ,20 + ,35 + ,22 + ,15 + ,18 + ,30 + ,22 + ,12 + ,12 + ,20 + ,20 + ,15 + ,25 + ,15 + ,20 + ,20 + ,25 + ,18 + ,30 + ,20 + ,25 + ,25 + ,25 + ,20 + ,20 + ,25 + ,20 + ,22 + ,15 + ,15 + ,22 + ,20 + ,10 + ,25 + ,20 + ,20 + ,15 + ,12 + ,20 + ,5 + ,20 + ,15 + ,15 + ,25 + ,25 + ,25 + ,15 + ,25 + ,22 + ,25 + ,20 + ,18 + ,22 + ,25 + ,35 + ,25 + ,25 + ,25 + ,35 + ,30 + ,22 + ,30 + ,50 + ,15 + ,25 + ,24 + ,20 + ,25 + ,25 + ,25 + ,12 + ,15 + ,22 + ,25 + ,25 + ,25 + ,25 + ,15 + ,20 + ,20 + ,15 + ,35 + ,30 + ,20 + ,22 + ,65 + ,20 + ,25 + ,22 + ,20 + ,25 + ,25 + ,20 + ,25 + ,15 + ,20 + ,12 + ,15 + ,10 + ,25 + ,15 + ,30 + ,35 + ,25 + ,25 + ,25 + ,25 + ,25 + ,40 + ,40 + ,25 + ,25 + ,20 + ,25 + ,25 + ,22 + ,25 + ,30 + ,25 + ,25 + ,30 + ,25 + ,25 + ,30 + ,25 + ,25 + ,20 + ,22 + ,22 + ,20 + ,25 + ,22 + ,25 + ,22 + ,40 + ,25 + ,25 + ,25 + ,22 + ,20 + ,35 + ,20 + ,35 + ,25 + ,22 + ,25 + ,25 + ,25 + ,25 + ,25 + ,40 + ,25 + ,30 + ,25 + ,20 + ,25 + ,25 + ,30 + ,22 + ,22 + ,20 + ,15 + ,15 + ,25 + ,25 + ,20 + ,20 + ,15 + ,25 + ,15 + ,20 + ,22 + ,25 + ,15 + ,15 + ,18 + ,5 + ,15 + ,25 + ,18 + ,40 + ,25 + ,25 + ,20 + ,30 + ,20 + ,25 + ,25 + ,25 + ,22 + ,22 + ,25 + ,25 + ,30 + ,25 + ,25 + ,25 + ,25 + ,20 + ,20 + ,25 + ,25 + ,25 + ,25 + ,20 + ,30 + ,25 + ,22 + ,30 + ,20 + ,20 + ,30 + ,25 + ,25 + ,30 + ,20 + ,25 + ,25 + ,24 + ,25 + ,30 + ,18 + ,15 + ,22 + ,22 + ,25 + ,22 + ,22 + ,25 + ,15 + ,20 + ,22 + ,18 + ,35 + ,20 + ,20 + ,20 + ,25 + ,25 + ,30 + ,15 + ,25 + ,22 + ,26 + ,25 + ,20 + ,25 + ,25 + ,25 + ,22 + ,25 + ,25 + ,20 + ,22 + ,30 + ,15 + ,30 + ,25 + ,20 + ,25 + ,25 + ,35 + ,22 + ,20 + ,25 + ,20 + ,20 + ,18 + ,20 + ,22 + ,25 + ,10 + ,20 + ,25 + ,20 + ,20 + ,30 + ,25 + ,20 + ,15 + ,20 + ,25 + ,10 + ,20 + ,25 + ,22 + ,22 + ,25 + ,25 + ,15 + ,25 + ,20 + ,10 + ,25 + ,16 + ,25 + ,35 + ,25 + ,15 + ,25 + ,25 + ,30 + ,25 + ,10 + ,22 + ,20 + ,25 + ,20 + ,20 + ,25 + ,22 + ,18 + ,30 + ,19 + ,25 + ,20 + ,25 + ,20 + ,25 + ,20 + ,22 + ,12 + ,30 + ,12 + ,22 + ,25 + ,25 + ,25 + ,25 + ,30 + ,30 + ,10 + ,22 + ,22 + ,25 + ,20 + ,22 + ,20 + ,25 + ,20 + ,15 + ,25 + ,20 + ,25 + ,20 + ,30 + ,15 + ,40 + ,25 + ,20 + ,22 + ,22 + ,30 + ,20 + ,40 + ,20 + ,25 + ,20 + ,25 + ,20 + ,50 + ,50 + ,25 + ,25 + ,40 + ,30 + ,22 + ,30 + ,20 + ,25 + ,25 + ,30 + ,25 + ,25 + ,20 + ,18 + ,18 + ,28 + ,25 + ,22 + ,15 + ,40 + ,40 + ,12 + ,12 + ,18 + ,12 + ,25 + ,26 + ,18 + ,25 + ,22 + ,15 + ,25 + ,15 + ,15 + ,15 + ,25 + ,15 + ,12 + ,22 + ,20 + ,20 + ,25 + ,20 + ,12 + ,9 + ,15 + ,12 + ,15 + ,25 + ,20 + ,20 + ,15 + ,15 + ,30 + ,21 + ,25 + ,22 + ,22 + ,50 + ,15 + ,25 + ,15 + ,25 + ,22 + ,18 + ,50 + ,20 + ,50 + ,20 + ,20 + ,30 + ,25 + ,20 + ,22 + ,25 + ,50 + ,40 + ,25 + ,25 + ,25 + ,25 + ,30 + ,40 + ,25 + ,30 + ,20) > par1 = '5' > par1 <- '5' > #'GNU S' R Code compiled by R2WASP v. 1.2.327 (Mon, 30 Nov 2015 07:01:18 +0000) > #Author: root > #To cite this work: Wessa P., (2015), Mean Plot (v1.0.5) 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) > x <- na.omit(x) > (n <- length(x)) [1] 900 > (np <- floor(n / par1)) [1] 180 > 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] 180 180 180 180 180 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 20 25 20 20 20 5 20 35 20 35 10 30 25 12 [2,] 25 25 40 25 15 20 22 30 25 25 22 20 25 18 [3,] 15 21 13 20 15 15 25 25 25 25 25 25 25 20 [4,] 15 30 30 25 12 25 20 20 15 30 25 25 22 20 [5,] 25 25 25 20 20 22 20 20 20 23 22 22 25 22 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 30 30 25 22 22 25 30 25 35 25 15 25 [2,] 25 25 30 25 15 30 15 20 30 25 40 20 [3,] 22 20 22 25 20 35 10 20 12 20 20 15 [4,] 20 30 25 25 30 22 30 35 25 20 40 15 [5,] 50 22 22 20 20 12 9 25 15 6 25 22 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] 24 25 22 25 22 25 15 25 25 10 25 28 [2,] 22 35 20 20 35 25 25 28 30 25 15 22 [3,] 20 40 25 25 22 22 18 30 22 20 22 30 [4,] 25 20 25 30 25 23 22 20 30 22 25 25 [5,] 25 22 18 20 25 35 25 25 10 25 25 20 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 25 30 20 15 30 20 22 18 25 20 15 15 [2,] 25 50 25 16 20 20 20 20 20 22 25 22 [3,] 20 19 35 20 25 25 25 25 25 18 25 25 [4,] 30 20 25 20 25 25 25 25 20 22 20 25 [5,] 20 28 25 25 25 30 18 30 20 20 25 15 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 12 22 25 15 25 30 50 20 25 20 22 20 [2,] 25 25 25 25 20 22 19 20 25 20 22 25 [3,] 30 12 40 20 30 25 50 20 25 25 30 30 [4,] 22 18 24 25 20 25 25 20 20 18 30 50 [5,] 15 30 25 25 25 25 35 20 20 25 8 22 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] 20 25 30 22 20 10 50 25 30 30 20 20 [2,] 10 18 18 20 20 20 30 25 25 40 9 30 [3,] 25 25 20 20 20 25 30 22 18 25 15 12 [4,] 25 20 25 25 25 30 50 20 22 20 20 15 [5,] 25 25 22 20 20 25 15 22 22 10 15 12 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] 20 25 25 15 20 40 10 25 25 20 20 20 [2,] 15 25 15 10 38 25 20 20 25 20 25 25 [3,] 12 25 15 20 20 25 25 22 25 16 20 20 [4,] 25 30 22 25 20 30 12 22 15 25 30 30 [5,] 20 20 10 20 20 25 15 20 40 15 50 30 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 25 25 20 30 15 15 12 15 30 20 15 25 [2,] 25 20 25 20 20 18 20 20 20 20 15 20 [3,] 12 20 15 20 20 30 20 20 25 25 22 20 [4,] 25 20 25 25 35 22 15 25 25 20 20 15 [5,] 25 15 50 12 22 12 25 18 25 22 10 12 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 20 25 22 25 35 15 25 25 20 20 [2,] 5 25 25 35 30 25 25 25 20 22 [3,] 20 25 20 25 22 24 12 25 15 65 [4,] 15 15 18 25 30 20 15 25 35 20 [5,] 15 25 22 25 50 25 22 15 30 25 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 22 25 10 25 40 25 25 30 22 22 [2,] 20 15 25 25 40 25 25 25 20 40 [3,] 25 20 15 25 25 22 30 25 25 25 [4,] 25 12 30 25 25 25 25 20 22 25 [5,] 20 15 35 25 20 30 25 22 25 25 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] 22 25 25 25 22 25 25 15 25 20 [2,] 20 22 25 20 22 25 15 15 18 30 [3,] 35 25 40 25 20 20 20 18 40 20 [4,] 20 25 25 25 15 20 22 5 25 25 [5,] 35 25 30 30 15 15 25 15 25 25 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] 25 30 20 25 30 25 24 22 25 35 [2,] 22 25 20 20 20 30 25 22 15 20 [3,] 22 25 25 30 20 20 30 25 20 20 [4,] 25 25 25 25 30 25 18 22 22 20 [5,] 25 25 25 22 25 25 15 22 18 25 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] 25 26 25 22 20 20 20 25 20 20 [2,] 30 25 22 30 25 25 22 20 15 25 [3,] 15 20 25 15 25 20 25 20 20 22 [4,] 25 25 25 30 35 20 10 30 25 22 [5,] 22 25 20 25 22 18 20 25 10 25 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [,156] [,157] [,158] [1,] 25 25 15 10 20 19 25 12 25 22 [2,] 15 16 25 22 25 25 20 22 30 25 [3,] 25 25 25 20 22 20 22 25 30 20 [4,] 20 35 30 25 18 25 12 25 10 22 [5,] 10 25 25 20 30 20 30 25 22 20 [,159] [,160] [,161] [,162] [,163] [,164] [,165] [,166] [,167] [,168] [1,] 25 25 25 20 25 25 20 25 25 12 [2,] 20 20 20 40 20 40 25 20 22 12 [3,] 15 30 22 20 50 30 25 18 15 18 [4,] 25 15 22 25 50 22 30 18 40 12 [5,] 20 40 30 20 25 30 25 28 40 25 [,169] [,170] [,171] [,172] [,173] [,174] [,175] [,176] [,177] [,178] [1,] 26 25 15 25 12 15 22 15 20 25 [2,] 18 15 12 20 15 15 22 25 50 20 [3,] 25 15 22 12 25 30 50 22 20 22 [4,] 22 15 20 9 20 21 15 18 20 25 [5,] 15 25 20 15 20 25 25 50 30 50 [,179] [,180] [,181] [1,] 40 30 NA [2,] 25 40 NA [3,] 25 25 NA [4,] 25 30 NA [5,] 25 20 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 5 0 20 5 -5 15 2 -5 5 -10 12 -10 0 6 [2,] -10 -4 -27 -5 0 -5 3 -5 0 0 3 5 0 2 [3,] 0 9 17 5 -3 10 -5 -5 -10 5 0 0 -3 0 [4,] 10 -5 -5 -5 8 -3 0 0 5 -7 -3 -3 3 2 [5,] 0 -5 -5 0 -15 -2 15 0 15 -13 8 3 -13 8 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] -5 -5 5 3 -7 5 -15 -5 -5 0 25 -5 [2,] -3 -5 -8 0 5 5 -5 0 -18 -5 -20 -5 [3,] -2 10 3 0 10 -13 20 15 13 0 20 0 [4,] 30 -8 -3 -5 -10 -10 -21 -10 -10 -14 -15 7 [5,] -20 3 0 2 5 18 16 10 10 9 0 2 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -2 10 -2 -5 13 0 10 3 5 15 -10 -6 [2,] -2 5 5 5 -13 -3 -7 2 -8 -5 7 8 [3,] 5 -20 0 5 3 1 4 -10 8 2 3 -5 [4,] 0 2 -7 -10 0 12 3 5 -20 3 0 -5 [5,] 0 0 7 2 0 -20 0 0 0 0 3 5 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 0 20 5 1 -10 0 -2 2 -5 2 10 7 [2,] -5 -31 10 4 5 5 5 5 5 -4 0 3 [3,] 10 1 -10 0 0 0 0 0 -5 4 -5 0 [4,] -10 8 0 5 0 5 -7 5 0 -2 5 -10 [5,] 10 -8 -10 5 -5 -8 0 -5 0 -5 -10 -3 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 13 3 0 10 -5 -8 -31 0 0 0 0 5 [2,] 5 -13 15 -5 10 3 31 0 0 5 8 5 [3,] -8 6 -16 5 -10 0 -25 0 -5 -7 0 20 [4,] -7 12 1 0 5 0 10 0 0 7 -22 -28 [5,] 7 -5 -10 0 5 25 -15 5 0 -3 12 -2 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] -10 -7 -12 -2 0 10 -20 0 -5 10 -11 10 [2,] 15 7 2 0 0 5 0 -3 -7 -15 6 -18 [3,] 0 -5 5 5 5 5 20 -2 4 -5 5 3 [4,] 0 5 -3 -5 -5 -5 -35 2 0 -10 -5 -3 [5,] 0 5 0 0 -10 25 10 8 8 10 5 8 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] -5 0 -10 -5 18 -15 10 -5 0 0 5 5 [2,] -3 0 0 10 -18 0 5 2 0 -4 -5 -5 [3,] 13 5 7 5 0 5 -13 0 -10 9 10 10 [4,] -5 -10 -12 -5 0 -5 3 -2 25 -10 20 0 [5,] 5 5 5 0 20 -15 10 5 -20 5 -30 -5 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 0 -5 5 -10 5 3 8 5 -10 0 0 -5 [2,] -13 0 -10 0 0 12 0 0 5 5 7 0 [3,] 13 0 10 5 15 -8 -5 5 0 -5 -2 -5 [4,] 0 -5 25 -13 -13 -10 10 -7 0 2 -10 -3 [5,] 0 5 -20 3 -7 0 -10 12 -5 -7 15 8 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] -15 0 3 10 -5 10 0 0 0 2 [2,] 15 0 -5 -10 -8 -1 -13 0 -5 43 [3,] -5 -10 -2 0 8 -4 3 0 20 -45 [4,] 0 10 4 0 20 5 7 -10 -5 5 [5,] 10 -3 3 10 -35 0 3 5 -10 -3 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] -2 -10 15 0 0 0 0 -5 -2 18 [2,] 5 5 -10 0 -15 -3 5 0 5 -15 [3,] 0 -8 15 0 0 3 -5 -5 -3 0 [4,] -5 3 5 0 -5 5 0 2 3 0 [5,] 5 -5 -10 15 5 -5 5 0 -3 -3 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] -2 -3 0 -5 0 0 -10 0 -7 10 [2,] 15 3 15 5 -2 -5 5 3 22 -10 [3,] -15 0 -15 0 -5 0 2 -13 -15 5 [4,] 15 0 5 5 0 -5 3 10 0 0 [5,] -10 0 -5 -8 10 10 -10 10 -5 0 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] -3 -5 0 -5 -10 5 1 0 -10 -15 [2,] 0 0 5 10 0 -10 5 3 5 0 [3,] 3 0 0 -5 10 5 -12 -3 2 0 [4,] 0 0 0 -3 -5 0 -3 0 -4 5 [5,] 5 -5 0 8 0 -1 7 3 17 0 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] 5 -1 -3 8 5 5 2 -5 -5 5 [2,] -15 -5 3 -15 0 -5 3 0 5 -3 [3,] 10 5 0 15 10 0 -15 10 5 0 [4,] -3 0 -5 -5 -13 -2 10 -5 -15 3 [5,] 4 0 2 -5 -2 2 5 -5 10 0 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [,156] [,157] [,158] [1,] -10 -9 10 12 5 6 -5 10 5 3 [2,] 10 9 0 -2 -3 -5 2 3 0 -5 [3,] -5 10 5 5 -4 5 -10 0 -20 2 [4,] -10 -10 -5 -5 12 -5 18 0 12 -2 [5,] 15 -10 -15 0 -11 5 -18 0 0 5 [,159] [,160] [,161] [,162] [,163] [,164] [,165] [,166] [,167] [,168] [1,] -5 -5 -5 20 -5 15 5 -5 -3 0 [2,] -5 10 2 -20 30 -10 0 -2 -7 6 [3,] 10 -15 0 5 0 -8 5 0 25 -6 [4,] -5 25 8 -5 -25 8 -5 10 0 13 [5,] 5 -15 -10 5 0 -10 0 -3 -28 1 [,169] [,170] [,171] [,172] [,173] [,174] [,175] [,176] [,177] [,178] [1,] -8 -10 -3 -5 3 0 0 10 30 -5 [2,] 7 0 10 -8 10 15 28 -3 -30 2 [3,] -3 0 -2 -3 -5 -9 -35 -4 0 3 [4,] -7 10 0 6 0 4 10 32 10 25 [5,] 10 -10 5 -3 -5 -3 -10 -30 -5 -10 [,179] [,180] [,181] [1,] -15 10 NA [2,] 0 -15 NA [3,] 0 5 NA [4,] 0 -10 NA [5,] 5 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/18bjm1476394262.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/2pzrc1476394262.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/3msnf1476394262.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/4t2z41476394262.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] [1,] 15 15 13 15 15 [2,] 20 20 20 20 20 [3,] 25 22 22 25 22 [4,] 25 25 25 25 25 [5,] 30 30 30 30 30 $n [1] 180 180 180 180 180 $conf [,1] [,2] [,3] [,4] [,5] [1,] 24.41117 21.41117 21.41117 24.41117 21.41117 [2,] 25.58883 22.58883 22.58883 25.58883 22.58883 $out [1] 5 35 35 10 12 35 10 12 50 10 50 40 10 12 35 10 40 35 10 12 12 12 40 40 40 [26] 35 35 50 10 40 9 10 38 5 35 40 40 40 40 12 12 50 40 35 10 12 40 35 12 40 [51] 50 12 12 12 12 65 35 40 40 50 12 50 12 35 40 50 50 12 35 35 12 5 35 10 35 [76] 12 10 50 40 12 9 50 12 9 6 35 10 35 8 10 12 10 40 50 50 12 12 10 12 50 [101] 35 35 10 10 40 40 50 50 $group [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 [38] 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 [75] 4 4 4 4 4 4 4 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 $names [1] "1" "2" "3" "4" "5" Warning message: In bxp(list(stats = c(15, 20, 25, 25, 30, 15, 20, 22, 25, 30, 13, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5z4fx1476394262.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] [1,] -20 -20 -20 -20 -20 [2,] -5 -5 -5 -5 -5 [3,] 0 0 0 0 0 [4,] 5 5 5 5 5 [5,] 20 15 20 20 20 $n [1] 180 180 180 180 179 $conf [,1] [,2] [,3] [,4] [,5] [1,] -1.177662 -1.177662 -1.177662 -1.177662 -1.180947 [2,] 1.177662 1.177662 1.177662 1.177662 1.180947 $out [1] 25 -31 30 -27 -31 31 43 22 30 28 -30 -25 -45 25 -35 30 -21 -22 -28 [20] -35 25 25 25 -25 32 25 25 25 -30 -35 -28 -30 $group [1] 1 1 1 2 2 2 2 2 2 2 2 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 5 5 5 5 5 5 $names [1] "1" "2" "3" "4" "5" > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6acgh1476394262.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] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 15 25 13 20 12 5 20 20 15 23 22 20 25 18 [2,] 15 25 20 20 15 15 20 20 20 25 22 22 25 18 [3,] 20 25 25 20 15 20 20 25 20 25 22 25 25 20 [4,] 25 25 30 25 20 22 22 30 25 30 25 25 25 20 [5,] 25 25 40 25 20 25 25 35 25 35 25 25 25 22 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 20 20 22 20 20 12 9 20 12 20 15 15 [2,] 22 22 22 22 20 22 10 20 15 20 20 15 [3,] 25 25 25 25 20 25 15 25 25 20 25 20 [4,] 30 30 25 25 22 30 30 25 30 25 40 22 [5,] 30 30 25 25 22 35 30 25 35 25 40 25 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] 20 20 18 20 22 22 15 25 10 20 22 20 [2,] 22 22 20 20 22 23 18 25 22 20 22 22 [3,] 24 25 22 25 25 25 22 25 25 22 25 25 [4,] 25 35 25 25 25 25 25 28 30 25 25 28 [5,] 25 40 25 30 25 25 25 30 30 25 25 30 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 20 19 25 15 25 20 18 18 20 18 15 15 [2,] 20 20 25 16 25 20 20 20 20 20 20 15 [3,] 25 28 25 20 25 25 22 25 20 20 25 22 [4,] 25 30 25 20 25 25 25 25 25 22 25 25 [5,] 30 30 25 25 25 30 25 30 25 22 25 25 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 12 12 25 15 20 25 19 20 20 18 22 20 [2,] 15 18 25 20 20 25 25 20 20 20 22 22 [3,] 22 22 25 25 25 25 35 20 25 20 22 25 [4,] 25 25 25 25 25 25 50 20 25 25 30 30 [5,] 30 30 25 25 30 25 50 20 25 25 30 30 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] 20 18 18 20 20 20 15 20 18 10 9 12 [2,] 20 20 20 20 20 20 30 22 22 20 15 12 [3,] 25 25 22 20 20 25 30 22 22 25 15 15 [4,] 25 25 25 22 20 25 50 25 25 30 20 20 [5,] 25 25 30 25 20 30 50 25 25 40 20 30 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] 12 25 10 10 20 25 10 20 25 15 20 20 [2,] 15 25 15 15 20 25 12 20 25 16 20 20 [3,] 20 25 15 20 20 25 15 22 25 20 25 25 [4,] 20 25 22 20 20 30 20 22 25 20 30 30 [5,] 25 25 25 25 20 30 25 25 25 25 30 30 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 25 20 15 20 20 12 12 15 25 20 10 12 [2,] 25 20 20 20 20 15 15 18 25 20 15 15 [3,] 25 20 25 20 20 18 20 20 25 20 15 20 [4,] 25 20 25 25 22 22 20 20 25 22 20 20 [5,] 25 20 25 30 22 30 25 20 25 25 22 25 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 15 25 18 25 30 15 12 25 15 20 [2,] 15 25 20 25 30 20 15 25 20 20 [3,] 15 25 22 25 30 24 22 25 20 22 [4,] 20 25 22 25 35 25 25 25 30 25 [5,] 20 25 25 25 35 25 25 25 35 25 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 20 12 10 25 20 25 25 20 20 25 [2,] 20 15 15 25 25 25 25 22 22 25 [3,] 22 15 25 25 25 25 25 25 22 25 [4,] 25 20 30 25 40 25 25 25 25 25 [5,] 25 25 35 25 40 25 25 25 25 25 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] 20 25 25 25 15 15 15 15 25 20 [2,] 20 25 25 25 15 20 20 15 25 20 [3,] 22 25 25 25 20 20 22 15 25 25 [4,] 35 25 30 25 22 25 25 15 25 25 [5,] 35 25 30 25 22 25 25 15 25 30 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] 22 25 20 20 20 25 15 22 15 20 [2,] 22 25 20 22 20 25 18 22 18 20 [3,] 25 25 25 25 25 25 24 22 20 20 [4,] 25 25 25 25 30 25 25 22 22 25 [5,] 25 25 25 25 30 25 30 22 25 25 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] 22 25 20 15 20 20 20 20 10 20 [2,] 22 25 22 22 22 20 20 20 15 22 [3,] 25 25 25 25 25 20 20 25 20 22 [4,] 25 25 25 30 25 20 22 25 20 25 [5,] 25 25 25 30 25 20 25 30 25 25 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [,156] [,157] [,158] [1,] 10 25 25 20 18 19 20 22 10 20 [2,] 15 25 25 20 20 20 20 22 22 20 [3,] 20 25 25 20 22 20 22 25 25 22 [4,] 25 25 25 22 25 25 25 25 30 22 [5,] 25 25 25 25 30 25 30 25 30 25 [,159] [,160] [,161] [,162] [,163] [,164] [,165] [,166] [,167] [,168] [1,] 15 15 20 20 20 22 25 18 15 12 [2,] 20 20 22 20 25 25 25 18 22 12 [3,] 20 25 22 20 25 30 25 20 25 12 [4,] 25 30 25 25 50 30 25 25 40 18 [5,] 25 40 25 25 50 30 25 28 40 25 [,169] [,170] [,171] [,172] [,173] [,174] [,175] [,176] [,177] [,178] [1,] 15 15 12 9 12 15 22 15 20 20 [2,] 18 15 15 12 15 15 22 18 20 22 [3,] 22 15 20 15 20 21 22 22 20 25 [4,] 25 25 20 20 20 25 25 25 30 25 [5,] 26 25 22 25 25 30 25 25 30 25 [,179] [,180] [,181] [1,] 25 20 NA [2,] 25 25 NA [3,] 25 30 NA [4,] 25 30 NA [5,] 25 30 NA $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 [38] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 [75] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 [112] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 [149] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 12.93403 25 17.93403 16.46701 11.46701 15.05382 18.58681 17.93403 [2,] 27.06597 25 32.06597 23.53299 18.53299 24.94618 21.41319 32.06597 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 16.46701 21.46701 19.88021 22.88021 25 18.58681 19.34722 19.34722 [2,] 23.53299 28.53299 24.11979 27.11979 25 21.41319 30.65278 30.65278 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 22.88021 22.88021 18.58681 19.34722 0.8680504 21.46701 14.40104 16.46701 [2,] 27.11979 27.11979 21.41319 30.65278 29.1319496 28.53299 35.59896 23.53299 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 10.86805 15.05382 21.88021 15.81423 18.46701 21.46701 22.88021 23.58681 [2,] 39.13195 24.94618 26.11979 34.18577 25.53299 28.53299 27.11979 26.41319 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 17.05382 22.88021 19.34722 18.46701 22.88021 20.76042 21.46701 20.93403 [2,] 26.94618 27.11979 30.65278 25.53299 27.11979 29.23958 28.53299 35.06597 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 25 17.17361 25 21.46701 18.46701 21.46701 16.46701 18.58681 21.46701 [2,] 25 22.82639 25 28.53299 25.53299 28.53299 23.53299 21.41319 28.53299 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [1,] 14.93403 14.93403 17.05382 25 21.46701 21.46701 25 17.33506 20 [2,] 29.06597 29.06597 26.94618 25 28.53299 28.53299 25 52.66494 20 [,59] [,60] [,61] [,62] [,63] [,64] [,65] [,66] [1,] 21.46701 16.46701 16.34722 19.34722 21.46701 21.46701 18.46701 18.58681 [2,] 28.53299 23.53299 27.65278 30.65278 28.53299 28.53299 25.53299 21.41319 [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] 20 21.46701 15.86805 19.88021 19.88021 17.93403 11.46701 9.34722 [2,] 20 28.53299 44.13195 24.11979 24.11979 32.06597 18.53299 20.65278 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [1,] 16.46701 25 10.05382 16.46701 20 21.46701 9.34722 20.58681 25 [2,] 23.53299 25 19.94618 23.53299 20 28.53299 20.65278 23.41319 25 [,84] [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [1,] 17.17361 17.93403 17.93403 25 20 21.46701 16.46701 18.58681 13.05382 [2,] 22.82639 32.06597 32.06597 25 20 28.53299 23.53299 21.41319 22.94618 [,93] [,94] [,95] [,96] [,97] [,98] [,99] [,100] [1,] 16.46701 18.58681 25 18.58681 11.46701 16.46701 11.46701 25 [2,] 23.53299 21.41319 25 21.41319 18.53299 23.53299 18.53299 25 [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 20.58681 25 26.46701 20.46701 14.93403 25 12.93403 18.46701 [2,] 23.41319 25 33.53299 27.53299 29.06597 25 27.06597 25.53299 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [1,] 18.46701 11.46701 14.40104 25 14.40104 25 25 22.88021 19.88021 [2,] 25.53299 18.53299 35.59896 25 35.59896 25 25 27.11979 24.11979 [,118] [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [1,] 25 11.40104 25 21.46701 25 15.05382 16.46701 18.46701 15 [2,] 25 32.59896 25 28.53299 25 24.94618 23.53299 25.53299 15 [,127] [,128] [,129] [,130] [,131] [,132] [,133] [,134] [,135] [1,] 25 21.46701 22.88021 25 21.46701 22.88021 17.93403 25 19.05382 [2,] 25 28.53299 27.11979 25 28.53299 27.11979 32.06597 25 28.94618 [,136] [,137] [,138] [,139] [,140] [,141] [,142] [,143] [,144] [1,] 22 17.17361 16.46701 22.88021 25 22.88021 19.34722 22.88021 20 [2,] 22 22.82639 23.53299 27.11979 25 27.11979 30.65278 27.11979 20 [,145] [,146] [,147] [,148] [,149] [,150] [,151] [,152] [1,] 18.58681 21.46701 16.46701 19.88021 12.93403 25 25 18.58681 [2,] 21.41319 28.53299 23.53299 24.11979 27.06597 25 25 21.41319 [,153] [,154] [,155] [,156] [,157] [,158] [,159] [,160] [1,] 18.46701 16.46701 18.46701 22.88021 19.34722 20.58681 16.46701 17.93403 [2,] 25.53299 23.53299 25.53299 27.11979 30.65278 23.41319 23.53299 32.06597 [,161] [,162] [,163] [,164] [,165] [,166] [,167] [,168] [1,] 19.88021 16.46701 7.335063 26.46701 25 15.05382 12.28125 7.760415 [2,] 24.11979 23.53299 42.664937 33.53299 25 24.94618 37.71875 16.239585 [,169] [,170] [,171] [,172] [,173] [,174] [,175] [,176] [1,] 17.05382 7.934025 16.46701 9.34722 16.46701 13.93403 19.88021 17.05382 [2,] 26.94618 22.065975 23.53299 20.65278 23.53299 28.06597 24.11979 26.94618 [,177] [,178] [,179] [,180] [,181] [1,] 12.93403 22.88021 25 26.46701 NA [2,] 27.06597 27.11979 25 33.53299 NA $out [1] 21 30 10 30 22 12 50 30 15 30 35 6 35 35 20 10 15 50 20 35 30 20 40 24 30 [26] 22 8 50 10 25 10 30 30 20 38 40 15 40 50 12 25 15 50 12 15 35 25 30 20 5 [51] 15 35 22 50 15 65 22 30 30 30 22 40 22 40 20 30 18 5 18 40 30 30 30 20 25 [76] 35 30 15 26 20 35 25 18 10 16 35 15 30 10 12 12 30 40 40 20 30 50 15 50 50 [101] 50 40 40 $group [1] 2 2 11 12 13 14 15 17 19 19 22 24 31 32 34 36 37 40 [19] 41 41 43 43 53 53 56 56 61 62 63 67 68 71 76 76 79 80 [37] 83 83 85 87 88 88 89 90 91 91 94 95 95 99 100 102 103 103 [55] 106 108 114 114 115 116 118 118 120 121 122 122 126 126 127 127 130 132 [73] 134 134 136 138 139 139 140 140 143 144 144 145 150 150 151 151 152 155 [91] 156 161 162 164 165 165 175 175 176 177 178 179 180 $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" "156" [157] "157" "158" "159" "160" "161" "162" "163" "164" "165" "166" "167" "168" [169] "169" "170" "171" "172" "173" "174" "175" "176" "177" "178" "179" "180" [181] NA Warning message: In bxp(list(stats = c(15, 15, 20, 25, 25, 25, 25, 25, 25, 25, 13, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7h6jb1476394262.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,] 23.11667 22 22.5 [2,] 23.11667 22 22.5 [3,] 23.18333 22 22.5 [4,] 23.20000 25 22.5 [5,] 23.20000 25 22.5 $n [1] 5 5 5 $conf [,1] [,2] [,3] [1,] 23.12445 19.88021 22.5 [2,] 23.24222 24.11979 22.5 $out [1] 22.87778 23.45000 $group [1] 1 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(23.1166666666667, 23.1166666666667, 23.1833333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/18bjm1476394262.ps tmp/18bjm1476394262.png",intern=TRUE)) character(0) > try(system("convert tmp/2pzrc1476394262.ps tmp/2pzrc1476394262.png",intern=TRUE)) character(0) > try(system("convert tmp/3msnf1476394262.ps tmp/3msnf1476394262.png",intern=TRUE)) character(0) > try(system("convert tmp/4t2z41476394262.ps tmp/4t2z41476394262.png",intern=TRUE)) character(0) > try(system("convert tmp/5z4fx1476394262.ps tmp/5z4fx1476394262.png",intern=TRUE)) character(0) > try(system("convert tmp/6acgh1476394262.ps tmp/6acgh1476394262.png",intern=TRUE)) character(0) > try(system("convert tmp/7h6jb1476394262.ps tmp/7h6jb1476394262.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.855 0.248 3.166