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(989,1215,2911,2372,2013,2050,1580,1407,903,709,490,206,1101,1189,2877,2489,2145,1837,1613,1296,849,642,475,224,920,1263,2999,2988,2163,2391,1556,1089,976,626,392,203,1052,1034,2353,3075,2309,2009,1464,1099,1035,792,406,187,862,822,2128,2264,1987,1728,1311,1152,945,704,526,361,1035,869,2698,2367,1926,1843,1404,1314,1007,865,587,339,1143,1807,2380,2337,2117,1789,1569,1305,952,810,473,278,993,1038,2257,2284,1747,1515,1233,882,1029,707,391,239,592,692,2127,1854,1468,1535,1203,880,821,604,315,139,528,654,1895,1598,1519,1242,1027,762,735,485,281,131,651,611,1898,1385,1047,1008,843,833,711,444,315,204,473,566,1611,1301,1154,1158,862,801,559,404,223,158,548,647,1757,1326,1308,1175,992,808,758,553,310,146) > 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,] 989 146 [2,] 1215 NA [3,] 2911 NA [4,] 2372 NA [5,] 2013 NA [6,] 2050 NA [7,] 1580 NA [8,] 1407 NA [9,] 903 NA [10,] 709 NA [11,] 490 NA [12,] 206 NA [13,] 1101 NA [14,] 1189 NA [15,] 2877 NA [16,] 2489 NA [17,] 2145 NA [18,] 1837 NA [19,] 1613 NA [20,] 1296 NA [21,] 849 NA [22,] 642 NA [23,] 475 NA [24,] 224 NA [25,] 920 NA [26,] 1263 NA [27,] 2999 NA [28,] 2988 NA [29,] 2163 NA [30,] 2391 NA [31,] 1556 NA [32,] 1089 NA [33,] 976 NA [34,] 626 NA [35,] 392 NA [36,] 203 NA [37,] 1052 NA [38,] 1034 NA [39,] 2353 NA [40,] 3075 NA [41,] 2309 NA [42,] 2009 NA [43,] 1464 NA [44,] 1099 NA [45,] 1035 NA [46,] 792 NA [47,] 406 NA [48,] 187 NA [49,] 862 NA [50,] 822 NA [51,] 2128 NA [52,] 2264 NA [53,] 1987 NA [54,] 1728 NA [55,] 1311 NA [56,] 1152 NA [57,] 945 NA [58,] 704 NA [59,] 526 NA [60,] 361 NA [61,] 1035 NA [62,] 869 NA [63,] 2698 NA [64,] 2367 NA [65,] 1926 NA [66,] 1843 NA [67,] 1404 NA [68,] 1314 NA [69,] 1007 NA [70,] 865 NA [71,] 587 NA [72,] 339 NA [73,] 1143 NA [74,] 1807 NA [75,] 2380 NA [76,] 2337 NA [77,] 2117 NA [78,] 1789 NA [79,] 1569 NA [80,] 1305 NA [81,] 952 NA [82,] 810 NA [83,] 473 NA [84,] 278 NA [85,] 993 NA [86,] 1038 NA [87,] 2257 NA [88,] 2284 NA [89,] 1747 NA [90,] 1515 NA [91,] 1233 NA [92,] 882 NA [93,] 1029 NA [94,] 707 NA [95,] 391 NA [96,] 239 NA [97,] 592 NA [98,] 692 NA [99,] 2127 NA [100,] 1854 NA [101,] 1468 NA [102,] 1535 NA [103,] 1203 NA [104,] 880 NA [105,] 821 NA [106,] 604 NA [107,] 315 NA [108,] 139 NA [109,] 528 NA [110,] 654 NA [111,] 1895 NA [112,] 1598 NA [113,] 1519 NA [114,] 1242 NA [115,] 1027 NA [116,] 762 NA [117,] 735 NA [118,] 485 NA [119,] 281 NA [120,] 131 NA [121,] 651 NA [122,] 611 NA [123,] 1898 NA [124,] 1385 NA [125,] 1047 NA [126,] 1008 NA [127,] 843 NA [128,] 833 NA [129,] 711 NA [130,] 444 NA [131,] 315 NA [132,] 204 NA [133,] 473 NA [134,] 566 NA [135,] 1611 NA [136,] 1301 NA [137,] 1154 NA [138,] 1158 NA [139,] 862 NA [140,] 801 NA [141,] 559 NA [142,] 404 NA [143,] 223 NA [144,] 158 NA [145,] 548 NA [146,] 647 NA [147,] 1757 NA [148,] 1326 NA [149,] 1308 NA [150,] 1175 NA [151,] 992 NA [152,] 808 NA [153,] 758 NA [154,] 553 NA [155,] 310 NA > darr [,1] [,2] [1,] 226 NA [2,] 1696 NA [3,] -539 NA [4,] -359 NA [5,] 37 NA [6,] -470 NA [7,] -173 NA [8,] -504 NA [9,] -194 NA [10,] -219 NA [11,] -284 NA [12,] 895 NA [13,] 88 NA [14,] 1688 NA [15,] -388 NA [16,] -344 NA [17,] -308 NA [18,] -224 NA [19,] -317 NA [20,] -447 NA [21,] -207 NA [22,] -167 NA [23,] -251 NA [24,] 696 NA [25,] 343 NA [26,] 1736 NA [27,] -11 NA [28,] -825 NA [29,] 228 NA [30,] -835 NA [31,] -467 NA [32,] -113 NA [33,] -350 NA [34,] -234 NA [35,] -189 NA [36,] 849 NA [37,] -18 NA [38,] 1319 NA [39,] 722 NA [40,] -766 NA [41,] -300 NA [42,] -545 NA [43,] -365 NA [44,] -64 NA [45,] -243 NA [46,] -386 NA [47,] -219 NA [48,] 675 NA [49,] -40 NA [50,] 1306 NA [51,] 136 NA [52,] -277 NA [53,] -259 NA [54,] -417 NA [55,] -159 NA [56,] -207 NA [57,] -241 NA [58,] -178 NA [59,] -165 NA [60,] 674 NA [61,] -166 NA [62,] 1829 NA [63,] -331 NA [64,] -441 NA [65,] -83 NA [66,] -439 NA [67,] -90 NA [68,] -307 NA [69,] -142 NA [70,] -278 NA [71,] -248 NA [72,] 804 NA [73,] 664 NA [74,] 573 NA [75,] -43 NA [76,] -220 NA [77,] -328 NA [78,] -220 NA [79,] -264 NA [80,] -353 NA [81,] -142 NA [82,] -337 NA [83,] -195 NA [84,] 715 NA [85,] 45 NA [86,] 1219 NA [87,] 27 NA [88,] -537 NA [89,] -232 NA [90,] -282 NA [91,] -351 NA [92,] 147 NA [93,] -322 NA [94,] -316 NA [95,] -152 NA [96,] 353 NA [97,] 100 NA [98,] 1435 NA [99,] -273 NA [100,] -386 NA [101,] 67 NA [102,] -332 NA [103,] -323 NA [104,] -59 NA [105,] -217 NA [106,] -289 NA [107,] -176 NA [108,] 389 NA [109,] 126 NA [110,] 1241 NA [111,] -297 NA [112,] -79 NA [113,] -277 NA [114,] -215 NA [115,] -265 NA [116,] -27 NA [117,] -250 NA [118,] -204 NA [119,] -150 NA [120,] 520 NA [121,] -40 NA [122,] 1287 NA [123,] -513 NA [124,] -338 NA [125,] -39 NA [126,] -165 NA [127,] -10 NA [128,] -122 NA [129,] -267 NA [130,] -129 NA [131,] -111 NA [132,] 269 NA [133,] 93 NA [134,] 1045 NA [135,] -310 NA [136,] -147 NA [137,] 4 NA [138,] -296 NA [139,] -61 NA [140,] -242 NA [141,] -155 NA [142,] -181 NA [143,] -65 NA [144,] 390 NA [145,] 99 NA [146,] 1110 NA [147,] -431 NA [148,] -18 NA [149,] -133 NA [150,] -183 NA [151,] -184 NA [152,] -50 NA [153,] -205 NA [154,] -243 NA [155,] -164 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/1j6ib1445024423.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/2ti681445024423.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/3r8ey1445024423.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/4u93w1445024423.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] [1,] 146.0 1215 2911 2372 2013 2050 1580 1407 903 709 490 206 1101 [2,] 146.0 1215 2911 2372 2013 2050 1580 1407 903 709 490 206 1101 [3,] 567.5 1215 2911 2372 2013 2050 1580 1407 903 709 490 206 1101 [4,] 989.0 1215 2911 2372 2013 2050 1580 1407 903 709 490 206 1101 [5,] 989.0 1215 2911 2372 2013 2050 1580 1407 903 709 490 206 1101 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 1189 2877 2489 2145 1837 1613 1296 849 642 475 224 920 [2,] 1189 2877 2489 2145 1837 1613 1296 849 642 475 224 920 [3,] 1189 2877 2489 2145 1837 1613 1296 849 642 475 224 920 [4,] 1189 2877 2489 2145 1837 1613 1296 849 642 475 224 920 [5,] 1189 2877 2489 2145 1837 1613 1296 849 642 475 224 920 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 1263 2999 2988 2163 2391 1556 1089 976 626 392 203 1052 [2,] 1263 2999 2988 2163 2391 1556 1089 976 626 392 203 1052 [3,] 1263 2999 2988 2163 2391 1556 1089 976 626 392 203 1052 [4,] 1263 2999 2988 2163 2391 1556 1089 976 626 392 203 1052 [5,] 1263 2999 2988 2163 2391 1556 1089 976 626 392 203 1052 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 1034 2353 3075 2309 2009 1464 1099 1035 792 406 187 862 [2,] 1034 2353 3075 2309 2009 1464 1099 1035 792 406 187 862 [3,] 1034 2353 3075 2309 2009 1464 1099 1035 792 406 187 862 [4,] 1034 2353 3075 2309 2009 1464 1099 1035 792 406 187 862 [5,] 1034 2353 3075 2309 2009 1464 1099 1035 792 406 187 862 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 822 2128 2264 1987 1728 1311 1152 945 704 526 361 1035 [2,] 822 2128 2264 1987 1728 1311 1152 945 704 526 361 1035 [3,] 822 2128 2264 1987 1728 1311 1152 945 704 526 361 1035 [4,] 822 2128 2264 1987 1728 1311 1152 945 704 526 361 1035 [5,] 822 2128 2264 1987 1728 1311 1152 945 704 526 361 1035 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 869 2698 2367 1926 1843 1404 1314 1007 865 587 339 1143 [2,] 869 2698 2367 1926 1843 1404 1314 1007 865 587 339 1143 [3,] 869 2698 2367 1926 1843 1404 1314 1007 865 587 339 1143 [4,] 869 2698 2367 1926 1843 1404 1314 1007 865 587 339 1143 [5,] 869 2698 2367 1926 1843 1404 1314 1007 865 587 339 1143 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 1807 2380 2337 2117 1789 1569 1305 952 810 473 278 993 [2,] 1807 2380 2337 2117 1789 1569 1305 952 810 473 278 993 [3,] 1807 2380 2337 2117 1789 1569 1305 952 810 473 278 993 [4,] 1807 2380 2337 2117 1789 1569 1305 952 810 473 278 993 [5,] 1807 2380 2337 2117 1789 1569 1305 952 810 473 278 993 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] 1038 2257 2284 1747 1515 1233 882 1029 707 391 239 592 [2,] 1038 2257 2284 1747 1515 1233 882 1029 707 391 239 592 [3,] 1038 2257 2284 1747 1515 1233 882 1029 707 391 239 592 [4,] 1038 2257 2284 1747 1515 1233 882 1029 707 391 239 592 [5,] 1038 2257 2284 1747 1515 1233 882 1029 707 391 239 592 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 692 2127 1854 1468 1535 1203 880 821 604 315 139 [2,] 692 2127 1854 1468 1535 1203 880 821 604 315 139 [3,] 692 2127 1854 1468 1535 1203 880 821 604 315 139 [4,] 692 2127 1854 1468 1535 1203 880 821 604 315 139 [5,] 692 2127 1854 1468 1535 1203 880 821 604 315 139 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 528 654 1895 1598 1519 1242 1027 762 735 485 [2,] 528 654 1895 1598 1519 1242 1027 762 735 485 [3,] 528 654 1895 1598 1519 1242 1027 762 735 485 [4,] 528 654 1895 1598 1519 1242 1027 762 735 485 [5,] 528 654 1895 1598 1519 1242 1027 762 735 485 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] 281 131 651 611 1898 1385 1047 1008 843 833 [2,] 281 131 651 611 1898 1385 1047 1008 843 833 [3,] 281 131 651 611 1898 1385 1047 1008 843 833 [4,] 281 131 651 611 1898 1385 1047 1008 843 833 [5,] 281 131 651 611 1898 1385 1047 1008 843 833 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] 711 444 315 204 473 566 1611 1301 1154 1158 [2,] 711 444 315 204 473 566 1611 1301 1154 1158 [3,] 711 444 315 204 473 566 1611 1301 1154 1158 [4,] 711 444 315 204 473 566 1611 1301 1154 1158 [5,] 711 444 315 204 473 566 1611 1301 1154 1158 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] 862 801 559 404 223 158 548 647 1757 1326 [2,] 862 801 559 404 223 158 548 647 1757 1326 [3,] 862 801 559 404 223 158 548 647 1757 1326 [4,] 862 801 559 404 223 158 548 647 1757 1326 [5,] 862 801 559 404 223 158 548 647 1757 1326 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [1,] 1308 1175 992 808 758 553 310 [2,] 1308 1175 992 808 758 553 310 [3,] 1308 1175 992 808 758 553 310 [4,] 1308 1175 992 808 758 553 310 [5,] 1308 1175 992 808 758 553 310 $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,] -374.3238 1215 2911 2372 2013 2050 1580 1407 903 709 490 206 1101 [2,] 1509.3238 1215 2911 2372 2013 2050 1580 1407 903 709 490 206 1101 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 1189 2877 2489 2145 1837 1613 1296 849 642 475 224 920 [2,] 1189 2877 2489 2145 1837 1613 1296 849 642 475 224 920 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 1263 2999 2988 2163 2391 1556 1089 976 626 392 203 1052 [2,] 1263 2999 2988 2163 2391 1556 1089 976 626 392 203 1052 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 1034 2353 3075 2309 2009 1464 1099 1035 792 406 187 862 [2,] 1034 2353 3075 2309 2009 1464 1099 1035 792 406 187 862 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 822 2128 2264 1987 1728 1311 1152 945 704 526 361 1035 [2,] 822 2128 2264 1987 1728 1311 1152 945 704 526 361 1035 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 869 2698 2367 1926 1843 1404 1314 1007 865 587 339 1143 [2,] 869 2698 2367 1926 1843 1404 1314 1007 865 587 339 1143 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 1807 2380 2337 2117 1789 1569 1305 952 810 473 278 993 [2,] 1807 2380 2337 2117 1789 1569 1305 952 810 473 278 993 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [1,] 1038 2257 2284 1747 1515 1233 882 1029 707 391 239 592 [2,] 1038 2257 2284 1747 1515 1233 882 1029 707 391 239 592 [,98] [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] 692 2127 1854 1468 1535 1203 880 821 604 315 139 [2,] 692 2127 1854 1468 1535 1203 880 821 604 315 139 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 528 654 1895 1598 1519 1242 1027 762 735 485 [2,] 528 654 1895 1598 1519 1242 1027 762 735 485 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] 281 131 651 611 1898 1385 1047 1008 843 833 [2,] 281 131 651 611 1898 1385 1047 1008 843 833 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] 711 444 315 204 473 566 1611 1301 1154 1158 [2,] 711 444 315 204 473 566 1611 1301 1154 1158 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] 862 801 559 404 223 158 548 647 1757 1326 [2,] 862 801 559 404 223 158 548 647 1757 1326 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [1,] 1308 1175 992 808 758 553 310 [2,] 1308 1175 992 808 758 553 310 $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(146, 146, 567.5, 989, 989, 1215, 1215, 1215, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5xll21445024423.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,] 226 1696 -539 -359 37 -470 -173 -504 -194 -219 -284 895 88 1688 [2,] 226 1696 -539 -359 37 -470 -173 -504 -194 -219 -284 895 88 1688 [3,] 226 1696 -539 -359 37 -470 -173 -504 -194 -219 -284 895 88 1688 [4,] 226 1696 -539 -359 37 -470 -173 -504 -194 -219 -284 895 88 1688 [5,] 226 1696 -539 -359 37 -470 -173 -504 -194 -219 -284 895 88 1688 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] -388 -344 -308 -224 -317 -447 -207 -167 -251 696 343 1736 [2,] -388 -344 -308 -224 -317 -447 -207 -167 -251 696 343 1736 [3,] -388 -344 -308 -224 -317 -447 -207 -167 -251 696 343 1736 [4,] -388 -344 -308 -224 -317 -447 -207 -167 -251 696 343 1736 [5,] -388 -344 -308 -224 -317 -447 -207 -167 -251 696 343 1736 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -11 -825 228 -835 -467 -113 -350 -234 -189 849 -18 1319 [2,] -11 -825 228 -835 -467 -113 -350 -234 -189 849 -18 1319 [3,] -11 -825 228 -835 -467 -113 -350 -234 -189 849 -18 1319 [4,] -11 -825 228 -835 -467 -113 -350 -234 -189 849 -18 1319 [5,] -11 -825 228 -835 -467 -113 -350 -234 -189 849 -18 1319 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 722 -766 -300 -545 -365 -64 -243 -386 -219 675 -40 1306 [2,] 722 -766 -300 -545 -365 -64 -243 -386 -219 675 -40 1306 [3,] 722 -766 -300 -545 -365 -64 -243 -386 -219 675 -40 1306 [4,] 722 -766 -300 -545 -365 -64 -243 -386 -219 675 -40 1306 [5,] 722 -766 -300 -545 -365 -64 -243 -386 -219 675 -40 1306 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 136 -277 -259 -417 -159 -207 -241 -178 -165 674 -166 1829 [2,] 136 -277 -259 -417 -159 -207 -241 -178 -165 674 -166 1829 [3,] 136 -277 -259 -417 -159 -207 -241 -178 -165 674 -166 1829 [4,] 136 -277 -259 -417 -159 -207 -241 -178 -165 674 -166 1829 [5,] 136 -277 -259 -417 -159 -207 -241 -178 -165 674 -166 1829 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] -331 -441 -83 -439 -90 -307 -142 -278 -248 804 664 573 [2,] -331 -441 -83 -439 -90 -307 -142 -278 -248 804 664 573 [3,] -331 -441 -83 -439 -90 -307 -142 -278 -248 804 664 573 [4,] -331 -441 -83 -439 -90 -307 -142 -278 -248 804 664 573 [5,] -331 -441 -83 -439 -90 -307 -142 -278 -248 804 664 573 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] -43 -220 -328 -220 -264 -353 -142 -337 -195 715 45 1219 [2,] -43 -220 -328 -220 -264 -353 -142 -337 -195 715 45 1219 [3,] -43 -220 -328 -220 -264 -353 -142 -337 -195 715 45 1219 [4,] -43 -220 -328 -220 -264 -353 -142 -337 -195 715 45 1219 [5,] -43 -220 -328 -220 -264 -353 -142 -337 -195 715 45 1219 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 27 -537 -232 -282 -351 147 -322 -316 -152 353 100 1435 [2,] 27 -537 -232 -282 -351 147 -322 -316 -152 353 100 1435 [3,] 27 -537 -232 -282 -351 147 -322 -316 -152 353 100 1435 [4,] 27 -537 -232 -282 -351 147 -322 -316 -152 353 100 1435 [5,] 27 -537 -232 -282 -351 147 -322 -316 -152 353 100 1435 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] -273 -386 67 -332 -323 -59 -217 -289 -176 389 [2,] -273 -386 67 -332 -323 -59 -217 -289 -176 389 [3,] -273 -386 67 -332 -323 -59 -217 -289 -176 389 [4,] -273 -386 67 -332 -323 -59 -217 -289 -176 389 [5,] -273 -386 67 -332 -323 -59 -217 -289 -176 389 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 126 1241 -297 -79 -277 -215 -265 -27 -250 -204 [2,] 126 1241 -297 -79 -277 -215 -265 -27 -250 -204 [3,] 126 1241 -297 -79 -277 -215 -265 -27 -250 -204 [4,] 126 1241 -297 -79 -277 -215 -265 -27 -250 -204 [5,] 126 1241 -297 -79 -277 -215 -265 -27 -250 -204 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] -150 520 -40 1287 -513 -338 -39 -165 -10 -122 [2,] -150 520 -40 1287 -513 -338 -39 -165 -10 -122 [3,] -150 520 -40 1287 -513 -338 -39 -165 -10 -122 [4,] -150 520 -40 1287 -513 -338 -39 -165 -10 -122 [5,] -150 520 -40 1287 -513 -338 -39 -165 -10 -122 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] -267 -129 -111 269 93 1045 -310 -147 4 -296 [2,] -267 -129 -111 269 93 1045 -310 -147 4 -296 [3,] -267 -129 -111 269 93 1045 -310 -147 4 -296 [4,] -267 -129 -111 269 93 1045 -310 -147 4 -296 [5,] -267 -129 -111 269 93 1045 -310 -147 4 -296 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] -61 -242 -155 -181 -65 390 99 1110 -431 -18 [2,] -61 -242 -155 -181 -65 390 99 1110 -431 -18 [3,] -61 -242 -155 -181 -65 390 99 1110 -431 -18 [4,] -61 -242 -155 -181 -65 390 99 1110 -431 -18 [5,] -61 -242 -155 -181 -65 390 99 1110 -431 -18 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [1,] -133 -183 -184 -50 -205 -243 -164 [2,] -133 -183 -184 -50 -205 -243 -164 [3,] -133 -183 -184 -50 -205 -243 -164 [4,] -133 -183 -184 -50 -205 -243 -164 [5,] -133 -183 -184 -50 -205 -243 -164 $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,] 226 1696 -539 -359 37 -470 -173 -504 -194 -219 -284 895 88 1688 [2,] 226 1696 -539 -359 37 -470 -173 -504 -194 -219 -284 895 88 1688 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] -388 -344 -308 -224 -317 -447 -207 -167 -251 696 343 1736 [2,] -388 -344 -308 -224 -317 -447 -207 -167 -251 696 343 1736 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -11 -825 228 -835 -467 -113 -350 -234 -189 849 -18 1319 [2,] -11 -825 228 -835 -467 -113 -350 -234 -189 849 -18 1319 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 722 -766 -300 -545 -365 -64 -243 -386 -219 675 -40 1306 [2,] 722 -766 -300 -545 -365 -64 -243 -386 -219 675 -40 1306 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 136 -277 -259 -417 -159 -207 -241 -178 -165 674 -166 1829 [2,] 136 -277 -259 -417 -159 -207 -241 -178 -165 674 -166 1829 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] -331 -441 -83 -439 -90 -307 -142 -278 -248 804 664 573 [2,] -331 -441 -83 -439 -90 -307 -142 -278 -248 804 664 573 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] -43 -220 -328 -220 -264 -353 -142 -337 -195 715 45 1219 [2,] -43 -220 -328 -220 -264 -353 -142 -337 -195 715 45 1219 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [,94] [,95] [,96] [,97] [,98] [1,] 27 -537 -232 -282 -351 147 -322 -316 -152 353 100 1435 [2,] 27 -537 -232 -282 -351 147 -322 -316 -152 353 100 1435 [,99] [,100] [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [1,] -273 -386 67 -332 -323 -59 -217 -289 -176 389 [2,] -273 -386 67 -332 -323 -59 -217 -289 -176 389 [,109] [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 126 1241 -297 -79 -277 -215 -265 -27 -250 -204 [2,] 126 1241 -297 -79 -277 -215 -265 -27 -250 -204 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [,128] [1,] -150 520 -40 1287 -513 -338 -39 -165 -10 -122 [2,] -150 520 -40 1287 -513 -338 -39 -165 -10 -122 [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [,137] [,138] [1,] -267 -129 -111 269 93 1045 -310 -147 4 -296 [2,] -267 -129 -111 269 93 1045 -310 -147 4 -296 [,139] [,140] [,141] [,142] [,143] [,144] [,145] [,146] [,147] [,148] [1,] -61 -242 -155 -181 -65 390 99 1110 -431 -18 [2,] -61 -242 -155 -181 -65 390 99 1110 -431 -18 [,149] [,150] [,151] [,152] [,153] [,154] [,155] [1,] -133 -183 -184 -50 -205 -243 -164 [2,] -133 -183 -184 -50 -205 -243 -164 $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/66sgm1445024423.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,] 131.0 146 [2,] 644.5 146 [3,] 1034.0 146 [4,] 1604.5 146 [5,] 2999.0 146 $n [1] 155 1 $conf [,1] [,2] [1,] 912.1677 146 [2,] 1155.8323 146 $out [1] 3075 $group [1] 1 $names [1] "1" NA > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7ug9i1445024423.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,] 131.0 131.0 131.0 [2,] 634.0 634.0 634.0 [3,] 1034.0 1034.0 1034.0 [4,] 1604.5 1604.5 1604.5 [5,] 2999.0 2999.0 2999.0 $n [1] 155 155 155 $conf [,1] [,2] [,3] [1,] 910.8352 910.8352 910.8352 [2,] 1157.1648 1157.1648 1157.1648 $out [1] 3075 3075 3075 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > try(system("convert tmp/1j6ib1445024423.ps tmp/1j6ib1445024423.png",intern=TRUE)) character(0) > try(system("convert tmp/2ti681445024423.ps tmp/2ti681445024423.png",intern=TRUE)) character(0) > try(system("convert tmp/3r8ey1445024423.ps tmp/3r8ey1445024423.png",intern=TRUE)) character(0) > try(system("convert tmp/4u93w1445024423.ps tmp/4u93w1445024423.png",intern=TRUE)) character(0) > try(system("convert tmp/5xll21445024423.ps tmp/5xll21445024423.png",intern=TRUE)) character(0) > try(system("convert tmp/66sgm1445024423.ps tmp/66sgm1445024423.png",intern=TRUE)) character(0) > try(system("convert tmp/7ug9i1445024423.ps tmp/7ug9i1445024423.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 3.038 0.472 3.541