R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- c(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 = '12' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > par1 <- as.numeric(par1) > (n <- length(x)) [1] 900 > (np <- floor(n / par1)) [1] 75 > 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] 75 75 75 75 75 75 75 75 75 75 75 75 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 20 13 20 30 30 25 22 22 30 35 15 20 18 [2,] 25 30 5 25 23 25 20 22 35 25 40 25 25 [3,] 15 25 20 20 10 25 50 25 22 35 20 25 20 [4,] 15 20 15 20 22 22 30 25 12 30 40 25 25 [5,] 25 25 25 20 25 25 25 25 30 12 25 35 30 [6,] 25 20 22 25 25 12 20 20 15 25 25 40 20 [7,] 25 25 20 25 22 18 30 22 10 15 20 20 22 [8,] 21 20 22 15 30 20 22 15 30 25 15 22 35 [9,] 30 20 25 20 20 20 25 20 9 25 15 22 22 [10,] 25 15 20 35 25 22 30 30 25 20 22 20 25 [11,] 20 15 20 25 25 30 22 20 20 20 24 25 25 [12,] 40 12 35 25 22 25 25 25 20 6 22 25 25 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 25 20 25 20 25 20 25 15 30 25 22 20 [2,] 22 25 15 30 15 25 30 25 22 15 25 20 [3,] 23 25 22 20 16 25 25 25 15 25 25 25 [4,] 35 30 25 30 20 30 20 20 22 20 25 25 [5,] 15 22 25 50 20 22 25 25 25 25 50 25 [6,] 25 30 28 19 25 20 20 15 12 25 19 20 [7,] 18 10 22 20 30 25 20 22 18 25 50 20 [8,] 22 10 30 28 20 25 20 25 30 20 25 20 [9,] 25 25 25 20 25 18 22 25 25 30 35 20 [10,] 25 20 20 25 25 18 18 15 25 20 20 25 [11,] 28 22 25 35 25 20 22 12 40 25 20 18 [12,] 30 25 25 25 20 25 20 25 24 30 20 25 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 22 25 22 20 20 20 12 10 25 22 20 12 [2,] 22 25 22 25 22 9 25 15 25 20 25 25 [3,] 30 25 20 30 30 15 20 10 30 25 20 25 [4,] 30 25 20 25 25 20 25 20 25 25 30 25 [5,] 8 18 25 50 18 15 25 25 10 25 50 20 [6,] 20 25 20 30 22 20 25 20 20 15 20 20 [7,] 25 20 20 30 22 30 30 20 25 40 25 20 [8,] 30 25 20 50 30 12 20 38 12 20 20 15 [9,] 50 30 20 15 40 15 25 20 15 20 30 20 [10,] 22 18 25 25 25 12 15 20 25 16 30 25 [11,] 20 20 20 25 20 20 15 20 20 25 25 15 [12,] 10 25 10 22 10 15 22 40 22 15 25 25 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 50 18 25 15 20 22 25 25 22 15 20 25 [2,] 30 30 18 15 15 25 24 15 20 30 25 25 [3,] 20 22 30 22 15 35 20 20 25 35 25 20 [4,] 20 12 20 20 25 25 25 20 25 25 22 22 [5,] 25 12 25 10 25 25 25 15 20 25 25 22 [6,] 12 20 25 25 25 25 25 35 25 25 30 20 [7,] 15 20 25 20 15 35 12 30 15 25 25 25 [8,] 20 15 20 20 25 30 15 20 20 25 25 22 [9,] 20 25 20 15 22 22 22 22 12 40 30 25 [10,] 35 15 25 12 25 30 25 65 15 40 25 22 [11,] 22 20 20 20 20 50 25 20 10 25 25 40 [12,] 15 20 22 5 18 15 25 25 25 25 30 25 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 25 25 20 25 30 25 30 30 18 25 30 20 [2,] 25 25 15 15 20 25 20 18 35 20 25 22 [3,] 22 40 15 15 25 20 20 15 20 25 20 25 [4,] 20 25 25 18 25 20 30 22 20 25 25 10 [5,] 35 30 25 5 25 25 25 22 20 25 25 20 [6,] 20 25 20 15 22 25 25 25 25 22 35 25 [7,] 35 20 20 25 22 25 30 22 25 25 22 20 [8,] 25 25 15 18 25 25 20 22 30 25 20 20 [9,] 22 25 25 40 25 20 25 25 15 20 25 30 [10,] 25 30 15 25 30 30 25 15 25 22 20 25 [11,] 25 22 20 25 25 25 24 20 22 30 20 20 [12,] 25 22 22 20 25 22 25 22 26 15 18 15 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 20 10 22 25 25 15 30 40 18 26 22 20 [2,] 25 25 20 20 30 25 20 30 28 18 20 15 [3,] 10 16 25 25 30 20 40 22 25 25 20 15 [4,] 20 25 20 20 10 25 20 30 22 22 25 30 [5,] 25 35 20 22 22 20 25 20 15 15 20 21 [6,] 22 25 25 12 22 30 20 25 40 25 12 25 [7,] 22 15 22 30 25 15 25 25 40 15 9 22 [8,] 25 25 18 12 20 40 20 30 12 15 15 22 [9,] 25 25 30 22 22 25 50 25 12 15 12 50 [10,] 15 30 19 25 20 20 50 25 18 25 15 15 [11,] 25 25 25 25 25 22 25 20 12 15 25 25 [12,] 20 10 20 25 20 22 25 18 25 12 20 15 [,74] [,75] [,76] [1,] 25 25 NA [2,] 22 50 NA [3,] 18 40 NA [4,] 50 25 NA [5,] 20 25 NA [6,] 50 25 NA [7,] 20 25 NA [8,] 20 30 NA [9,] 30 40 NA [10,] 25 25 NA [11,] 20 30 NA [12,] 22 20 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 5 17 -15 -5 -7 0 -2 0 5 -10 25 5 7 [2,] -10 -5 15 -5 -13 0 30 3 -13 10 -20 0 -5 [3,] 0 -5 -5 0 12 -3 -20 0 -10 -5 20 0 5 [4,] 10 5 10 0 3 3 -5 0 18 -18 -15 10 5 [5,] 0 -5 -3 5 0 -13 -5 -5 -15 13 0 5 -10 [6,] 0 5 -2 0 -3 6 10 2 -5 -10 -5 -20 2 [7,] -4 -5 2 -10 8 2 -8 -7 20 10 -5 2 13 [8,] 9 0 3 5 -10 0 3 5 -21 0 0 0 -13 [9,] -5 -5 -5 15 5 2 5 10 16 -5 7 -2 3 [10,] -5 0 0 -10 0 8 -8 -10 -5 0 2 5 0 [11,] 20 -3 15 0 -3 -5 3 5 0 -14 -2 0 0 [12,] -27 8 -5 5 3 -3 -3 5 15 9 -2 -7 0 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] -3 5 -10 10 -10 5 5 10 -8 -10 3 0 [2,] 1 0 7 -10 1 0 -5 0 -7 10 0 5 [3,] 12 5 3 10 4 5 -5 -5 7 -5 0 0 [4,] -20 -8 0 20 0 -8 5 5 3 5 25 0 [5,] 10 8 3 -31 5 -2 -5 -10 -13 0 -31 -5 [6,] -7 -20 -6 1 5 5 0 7 6 0 31 0 [7,] 4 0 8 8 -10 0 0 3 12 -5 -25 0 [8,] 3 15 -5 -8 5 -7 2 0 -5 10 10 0 [9,] 0 -5 -5 5 0 0 -4 -10 0 -10 -15 5 [10,] 3 2 5 10 0 2 4 -3 15 5 0 -7 [11,] 2 3 0 -10 -5 5 -2 13 -16 5 0 7 [12,] -10 0 -5 0 0 0 -5 5 1 -8 0 -3 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 0 0 0 5 2 -11 13 5 0 -2 5 13 [2,] 8 0 -2 5 8 6 -5 -5 5 5 -5 0 [3,] 0 0 0 -5 -5 5 5 10 -5 0 10 0 [4,] -22 -7 5 25 -7 -5 0 5 -15 0 20 -5 [5,] 12 7 -5 -20 4 5 0 -5 10 -10 -30 0 [6,] 5 -5 0 0 0 10 5 0 5 25 5 0 [7,] 5 5 0 20 8 -18 -10 18 -13 -20 -5 -5 [8,] 20 5 0 -35 10 3 5 -18 3 0 10 5 [9,] -28 -12 5 10 -15 -3 -10 0 10 -4 0 5 [10,] -2 2 -5 0 -5 8 0 0 -5 9 -5 -10 [11,] -10 5 -10 -3 -10 -5 7 20 2 -10 0 10 [12,] 15 -3 10 -2 10 -3 -12 -15 0 5 -13 25 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] -20 12 -7 0 -5 3 -1 -10 -2 15 5 0 [2,] -10 -8 12 7 0 10 -4 5 5 5 0 -5 [3,] 0 -10 -10 -2 10 -10 5 0 0 -10 -3 2 [4,] 5 0 5 -10 0 0 0 -5 -5 0 3 0 [5,] -13 8 0 15 0 0 0 20 5 0 5 -2 [6,] 3 0 0 -5 -10 10 -13 -5 -10 0 -5 5 [7,] 5 -5 -5 0 10 -5 3 -10 5 0 0 -3 [8,] 0 10 0 -5 -3 -8 7 2 -8 15 5 3 [9,] 15 -10 5 -3 3 8 3 43 3 0 -5 -3 [10,] -13 5 -5 8 -5 20 0 -45 -5 -15 0 18 [11,] -7 0 2 -15 -2 -35 0 5 15 0 5 -15 [12,] 3 5 -7 15 4 10 0 -3 -10 -5 -5 0 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 0 0 -5 -10 -10 0 -10 -12 17 -5 -5 2 [2,] -3 15 0 0 5 -5 0 -3 -15 5 -5 3 [3,] -2 -15 10 3 0 0 10 7 0 0 5 -15 [4,] 15 5 0 -13 0 5 -5 0 0 0 0 10 [5,] -15 -5 -5 10 -3 0 0 3 5 -3 10 5 [6,] 15 -5 0 10 0 0 5 -3 0 3 -13 -5 [7,] -10 5 -5 -7 3 0 -10 0 5 0 -2 0 [8,] -3 0 10 22 0 -5 5 3 -15 -5 5 10 [9,] 3 5 -10 -15 5 10 0 -10 10 2 -5 -5 [10,] 0 -8 5 0 -5 -5 -1 5 -3 8 0 -5 [11,] 0 0 2 -5 0 -3 1 2 4 -15 -2 -5 [12,] 0 -2 3 10 0 8 5 -4 -1 15 2 5 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 5 15 -2 -5 5 10 -10 -10 10 -8 -2 -5 [2,] -15 -9 5 5 0 -5 20 -8 -3 7 0 0 [3,] 10 9 -5 -5 -20 5 -20 8 -3 -3 5 15 [4,] 5 10 0 2 12 -5 5 -10 -7 -7 -5 -9 [5,] -3 -10 5 -10 0 10 -5 5 25 10 -8 4 [6,] 0 -10 -3 18 3 -15 5 0 0 -10 -3 -3 [7,] 3 10 -4 -18 -5 25 -5 5 -28 0 6 0 [8,] 0 0 12 10 2 -15 30 -5 0 0 -3 28 [9,] -10 5 -11 3 -2 -5 0 0 6 10 3 -35 [10,] 10 -5 6 0 5 2 -25 -5 -6 -10 10 10 [11,] -5 -15 -5 0 -5 0 0 -2 13 -3 -5 -10 [12,] -10 12 5 0 -5 8 15 0 1 10 0 10 [,74] [,75] [,76] [1,] -3 25 NA [2,] -4 -10 NA [3,] 32 -15 NA [4,] -30 0 NA [5,] 30 0 NA [6,] -30 0 NA [7,] 0 5 NA [8,] 10 10 NA [9,] -5 -15 NA [10,] -5 5 NA [11,] 2 -10 NA [12,] 3 NA NA > arr.mean <- array(NA,dim=par1) > arr.median <- array(NA,dim=par1) > arr.midrange <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.mean[j] <- mean(arr[j,],na.rm=TRUE) + arr.median[j] <- median(arr[j,],na.rm=TRUE) + arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2 + } > overall.mean <- mean(x) > overall.median <- median(x) > overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2 > postscript(file="/var/www/html/rcomp/tmp/1xru71289987105.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.mean,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/270cs1289987105.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.median,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/370cs1289987105.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.midrange,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/470cs1289987105.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(arr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 13 15 15 15 15 15 15 15 15 15 15 15 [2,] 20 20 20 20 20 20 20 20 20 20 20 20 [3,] 22 25 22 25 25 25 22 21 25 25 22 22 [4,] 25 25 25 25 25 25 25 25 25 25 25 25 [5,] 30 30 30 30 30 30 30 30 30 30 30 30 $n [1] 75 75 75 75 75 75 75 75 75 75 75 75 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 21.08779 24.08779 21.08779 24.08779 24.08779 24.08779 21.08779 20.08779 [2,] 22.91221 25.91221 22.91221 25.91221 25.91221 25.91221 22.91221 21.91221 [,9] [,10] [,11] [,12] [1,] 24.08779 24.08779 21.08779 21.08779 [2,] 25.91221 25.91221 22.91221 22.91221 $out [1] 35 12 10 12 50 10 40 5 35 40 9 35 50 10 50 35 10 35 35 40 10 40 40 12 40 [26] 35 12 10 10 50 12 35 50 50 8 50 10 50 12 10 35 5 35 12 40 12 12 35 35 12 [51] 40 12 50 10 10 50 40 35 12 35 40 9 35 10 50 12 38 12 12 40 12 9 35 50 40 [76] 12 40 40 50 12 12 50 40 35 12 35 12 65 40 50 35 12 40 50 10 40 12 40 12 35 [101] 6 10 10 10 40 5 10 12 $group [1] 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 4 4 [26] 4 4 4 4 4 5 5 5 5 5 5 5 5 5 5 5 5 5 6 6 6 6 6 6 6 [51] 6 6 6 7 7 7 7 7 7 7 7 7 8 8 8 8 8 8 8 8 8 9 9 9 9 [76] 9 9 9 9 9 9 9 9 10 10 10 10 10 10 10 11 11 11 11 11 11 11 12 12 12 [101] 12 12 12 12 12 12 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(13, 20, 22, 25, 30, 15, 20, 25, 25, 30, 15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5irbv1289987105.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(darr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] -20 -20 -20 -20 -20 -20 -20 -15 -15 -15 -15 -15 [2,] -5 -5 -5 -5 -5 -5 -5 -3 -5 -5 -5 -3 [3,] 0 0 0 0 0 0 0 2 0 0 0 0 [4,] 5 5 5 5 5 5 5 6 5 5 2 5 [5,] 17 20 20 20 20 18 20 15 16 20 10 15 $n [1] 75 75 75 75 75 75 75 75 75 75 75 74 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.824427 -1.824427 -1.824427 -1.824427 -1.824427 -1.824427 -1.824427 [2,] 1.824427 1.824427 1.824427 1.824427 1.824427 1.824427 1.824427 [,8] [,9] [,10] [,11] [,12] [1,] 0.3580158 -1.824427 -1.824427 -1.277099 -1.469370 [2,] 3.6419842 1.824427 1.824427 1.277099 1.469370 $out [1] 25 25 30 32 25 -22 25 -30 -31 -31 -30 25 30 31 25 -30 -25 25 -28 [20] -21 20 -35 -18 22 30 28 -28 43 -35 -45 -25 20 15 13 -16 20 -35 15 [39] 13 -27 25 $group [1] 1 1 2 3 4 4 4 4 5 5 5 5 5 6 6 6 7 7 7 8 8 8 8 8 8 [26] 8 9 9 9 10 10 11 11 11 11 11 11 11 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6irbv1289987105.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(arr) > names(z) <- c(1:np) > (boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 15 12 15.0 15 20 18.0 20 15 9.0 6.0 15.0 20 18 18.0 [2,] 20 15 20.0 20 22 20.0 22 20 13.5 17.5 17.5 21 21 22.0 [3,] 25 20 20.0 25 24 23.5 25 22 21.0 25.0 22.0 25 25 25.0 [4,] 25 25 23.5 25 25 25.0 30 25 30.0 27.5 25.0 25 25 26.5 [5,] 30 30 25.0 30 25 30.0 30 30 35.0 35.0 25.0 25 30 30.0 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 20.0 20 19 15.0 18.0 18 12.0 12.0 15 19 18 8 [2,] 20.0 22 20 20.0 20.0 20 15.0 20.0 20 20 20 20 [3,] 23.5 25 25 22.5 23.5 21 23.5 24.5 25 25 20 22 [4,] 25.0 25 30 25.0 25.0 25 25.0 27.5 25 30 25 30 [5,] 30.0 28 35 30.0 30.0 30 25.0 30.0 30 35 25 30 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] 18 20 15.0 10.0 9.0 12.0 10.0 10.0 15 20 12.0 12.0 [2,] 20 20 23.5 20.0 13.5 17.5 17.5 17.5 18 20 17.5 17.5 [3,] 25 20 25.0 22.0 15.0 23.5 20.0 23.5 21 25 20.0 20.0 [4,] 25 22 30.0 27.5 20.0 25.0 22.5 25.0 25 30 25.0 27.5 [5,] 30 25 30.0 30.0 20.0 30.0 25.0 30.0 25 30 25.0 35.0 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 12 18.0 5.0 15.0 15.0 15 15.0 10 15.0 22.0 20.0 20 [2,] 15 20.0 13.5 16.5 23.5 21 20.0 15 25.0 25.0 22.0 22 [3,] 20 23.5 17.5 21.0 25.0 25 21.0 20 25.0 25.0 23.5 25 [4,] 21 25.0 20.0 25.0 32.5 25 27.5 25 32.5 27.5 25.0 25 [5,] 30 30.0 25.0 25.0 35.0 25 35.0 25 40.0 30.0 25.0 25 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 20.0 15.0 5 22.0 20 20.0 15.0 15.0 15 18.0 15 15 [2,] 23.5 15.0 15 23.5 21 22.0 19.0 20.0 21 20.0 20 20 [3,] 25.0 20.0 19 25.0 25 25.0 22.0 23.5 25 23.5 20 22 [4,] 27.5 23.5 25 25.0 25 27.5 23.5 25.5 25 25.0 25 25 [5,] 30.0 25.0 40 25.0 30 30.0 30.0 30.0 30 30.0 30 25 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] 10.0 18 20.0 20 15 20 18 12.0 12.0 9.0 15.0 18.0 [2,] 15.5 20 20.0 20 20 20 21 13.5 15.0 13.5 15.0 20.0 [3,] 25.0 21 23.5 22 22 25 25 20.0 16.5 20.0 21.5 22.0 [4,] 25.0 25 25.0 25 25 35 30 26.5 25.0 21.0 25.0 27.5 [5,] 35.0 30 30.0 30 30 50 40 40.0 26.0 25.0 30.0 30.0 [,75] [,76] [1,] 20 NA [2,] 25 NA [3,] 25 NA [4,] 35 NA [5,] 50 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 [26] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 [51] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 [76] 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 22.71947 15.43893 18.40363 22.71947 22.63168 21.21947 21.35115 19.71947 [2,] 27.28053 24.56107 21.59637 27.28053 25.36832 25.78053 28.64885 24.28053 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 13.47424 20.43893 18.5792 23.17557 23.17557 22.94752 21.21947 23.63168 [2,] 28.52576 29.56107 25.4208 26.82443 26.82443 27.05248 25.78053 26.36832 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 20.43893 20.21947 21.21947 18.71947 18.93893 21.0792 22.71947 20.43893 [2,] 29.56107 24.78053 25.78053 23.28053 28.06107 27.9208 27.28053 29.56107 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 17.71947 17.43893 22.71947 19.08779 22.03531 18.5792 12.03531 20.0792 [2,] 22.28053 26.56107 27.28053 20.91221 27.96469 25.4208 17.96469 26.9208 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 17.71947 20.0792 17.80725 20.43893 16.5792 15.43893 17.26336 21.21947 [2,] 22.28053 26.9208 24.19275 29.56107 23.4208 24.56107 22.73664 25.78053 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 14.53531 17.12309 20.89504 23.17557 17.5792 15.43893 21.5792 23.85973 [2,] 20.46469 24.87691 29.10496 26.82443 24.4208 24.56107 28.4208 26.14027 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [1,] 22.13168 23.63168 23.17557 16.12309 14.43893 24.31584 23.17557 22.49141 [2,] 24.86832 26.36832 26.82443 23.87691 23.56107 25.68416 26.82443 27.50859 [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [1,] 19.94752 20.99141 23.17557 21.21947 17.71947 19.71947 20.66699 18.71947 [2,] 24.05248 26.00859 26.82443 25.78053 22.28053 24.28053 29.33301 23.28053 [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 21.21947 19.71947 19.71947 18.1584 20.89504 14.07061 11.93893 16.5792 [2,] 25.78053 24.28053 24.28053 31.8416 29.10496 25.92939 21.06107 23.4208 [,73] [,74] [,75] [,76] [1,] 16.93893 18.5792 20.43893 NA [2,] 26.06107 25.4208 29.56107 NA $out [1] 40 5 35 35 30 10 30 12 50 40 40 35 40 35 35 15 10 10 15 30 50 40 50 50 50 [26] 10 50 50 40 30 38 40 40 50 50 50 12 65 20 40 35 35 40 30 20 30 35 35 10 10 [51] 12 12 10 40 50 50 50 $group [1] 1 3 3 4 5 5 5 6 7 11 11 12 12 13 14 14 15 15 16 16 17 22 24 24 26 [26] 28 29 29 30 31 33 33 35 36 38 43 44 45 48 49 50 50 51 54 54 54 58 60 61 62 [51] 65 65 66 67 73 74 74 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45" [46] "46" "47" "48" "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" "74" "75" [76] NA Warning message: In bxp(list(stats = c(15, 20, 25, 25, 30, 12, 15, 20, 25, 30, 15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7tjsy1289987105.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.mean,arr.median,arr.midrange)) > names(z) <- list('mean','median','midrange') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 22.53333 21.0 22.5 [2,] 22.88000 22.0 22.5 [3,] 23.21333 23.5 22.5 [4,] 23.62667 25.0 22.5 [5,] 24.30667 25.0 22.5 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 22.87277 22.13168 22.5 [2,] 23.55389 24.86832 22.5 $out [1] 21.64 $group [1] 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(22.5333333333333, 22.88, 23.2133333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1xru71289987105.ps tmp/1xru71289987105.png",intern=TRUE)) character(0) > try(system("convert tmp/270cs1289987105.ps tmp/270cs1289987105.png",intern=TRUE)) character(0) > try(system("convert tmp/370cs1289987105.ps tmp/370cs1289987105.png",intern=TRUE)) character(0) > try(system("convert tmp/470cs1289987105.ps tmp/470cs1289987105.png",intern=TRUE)) character(0) > try(system("convert tmp/5irbv1289987105.ps tmp/5irbv1289987105.png",intern=TRUE)) character(0) > try(system("convert tmp/6irbv1289987105.ps tmp/6irbv1289987105.png",intern=TRUE)) character(0) > try(system("convert tmp/7tjsy1289987105.ps tmp/7tjsy1289987105.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.467 1.005 4.925