R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: i686-pc-linux-gnu (32-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(221102 + ,220892 + ,225537 + ,219334 + ,216126 + ,224573 + ,213991 + ,221865 + ,215382 + ,213962 + ,217009 + ,318110 + ,219662 + ,212650 + ,209307 + ,210541 + ,210609 + ,208910 + ,207541 + ,207699 + ,205005 + ,205747 + ,203644 + ,229937 + ,214446 + ,210194 + ,206535 + ,216524 + ,198243 + ,208274 + ,207493 + ,215525 + ,207562 + ,213355 + ,209048 + ,220497 + ,214563 + ,211571 + ,216385 + ,211496 + ,209683 + ,206304 + ,213925 + ,204829 + ,205729 + ,200296 + ,207960 + ,207729 + ,208327 + ,207794 + ,213700 + ,213136 + ,210354 + ,202205 + ,220425 + ,210038 + ,205201 + ,197441 + ,200021 + ,203165 + ,201699 + ,201003 + ,232308 + ,211412 + ,209661 + ,213587 + ,193611 + ,192543 + ,196040 + ,191407 + ,189726 + ,191272 + ,187924 + ,198213 + ,199352 + ,199289 + ,195475 + ,198045 + ,197615 + ,189015 + ,189668 + ,189120 + ,194168 + ,192304 + ,185913 + ,197599 + ,186085 + ,190566 + ,187054 + ,193222 + ,189856 + ,190608 + ,190588 + ,186773 + ,183510 + ,180106 + ,180150 + ,178412 + ,182353 + ,203805 + ,186054 + ,184290 + ,187483 + ,187111 + ,189561 + ,184439 + ,182985 + ,183828 + ,184036 + ,183214 + ,183464 + ,173718 + ,180210 + ,171252 + ,172705 + ,174006 + ,172043 + ,169445 + ,169449 + ,177073 + ,170799 + ,171648 + ,172220 + ,165795 + ,167466 + ,165528 + ,162851 + ,165864 + ,162094 + ,162385 + ,164293 + ,165983 + ,159680 + ,161739 + ,159302 + ,167795 + ,164242 + ,159743 + ,160887 + ,163844 + ,161172 + ,159330 + ,155570 + ,156749 + ,155012 + ,163419 + ,153630 + ,154535 + ,151543 + ,152955 + ,150166 + ,151416 + ,150332 + ,152196 + ,153422 + ,147435) > par1 = '52' > par1 <- '52' > #'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] 3 > 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] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 [39] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 > arr [,1] [,2] [,3] [,4] [1,] 221102 210354 189561 NA [2,] 220892 202205 184439 NA [3,] 225537 220425 182985 NA [4,] 219334 210038 183828 NA [5,] 216126 205201 184036 NA [6,] 224573 197441 183214 NA [7,] 213991 200021 183464 NA [8,] 221865 203165 173718 NA [9,] 215382 201699 180210 NA [10,] 213962 201003 171252 NA [11,] 217009 232308 172705 NA [12,] 318110 211412 174006 NA [13,] 219662 209661 172043 NA [14,] 212650 213587 169445 NA [15,] 209307 193611 169449 NA [16,] 210541 192543 177073 NA [17,] 210609 196040 170799 NA [18,] 208910 191407 171648 NA [19,] 207541 189726 172220 NA [20,] 207699 191272 165795 NA [21,] 205005 187924 167466 NA [22,] 205747 198213 165528 NA [23,] 203644 199352 162851 NA [24,] 229937 199289 165864 NA [25,] 214446 195475 162094 NA [26,] 210194 198045 162385 NA [27,] 206535 197615 164293 NA [28,] 216524 189015 165983 NA [29,] 198243 189668 159680 NA [30,] 208274 189120 161739 NA [31,] 207493 194168 159302 NA [32,] 215525 192304 167795 NA [33,] 207562 185913 164242 NA [34,] 213355 197599 159743 NA [35,] 209048 186085 160887 NA [36,] 220497 190566 163844 NA [37,] 214563 187054 161172 NA [38,] 211571 193222 159330 NA [39,] 216385 189856 155570 NA [40,] 211496 190608 156749 NA [41,] 209683 190588 155012 NA [42,] 206304 186773 163419 NA [43,] 213925 183510 153630 NA [44,] 204829 180106 154535 NA [45,] 205729 180150 151543 NA [46,] 200296 178412 152955 NA [47,] 207960 182353 150166 NA [48,] 207729 203805 151416 NA [49,] 208327 186054 150332 NA [50,] 207794 184290 152196 NA [51,] 213700 187483 153422 NA [52,] 213136 187111 147435 NA > darr [,1] [,2] [,3] [,4] [1,] -210 -8149 -5122 NA [2,] 4645 18220 -1454 NA [3,] -6203 -10387 843 NA [4,] -3208 -4837 208 NA [5,] 8447 -7760 -822 NA [6,] -10582 2580 250 NA [7,] 7874 3144 -9746 NA [8,] -6483 -1466 6492 NA [9,] -1420 -696 -8958 NA [10,] 3047 31305 1453 NA [11,] 101101 -20896 1301 NA [12,] -98448 -1751 -1963 NA [13,] -7012 3926 -2598 NA [14,] -3343 -19976 4 NA [15,] 1234 -1068 7624 NA [16,] 68 3497 -6274 NA [17,] -1699 -4633 849 NA [18,] -1369 -1681 572 NA [19,] 158 1546 -6425 NA [20,] -2694 -3348 1671 NA [21,] 742 10289 -1938 NA [22,] -2103 1139 -2677 NA [23,] 26293 -63 3013 NA [24,] -15491 -3814 -3770 NA [25,] -4252 2570 291 NA [26,] -3659 -430 1908 NA [27,] 9989 -8600 1690 NA [28,] -18281 653 -6303 NA [29,] 10031 -548 2059 NA [30,] -781 5048 -2437 NA [31,] 8032 -1864 8493 NA [32,] -7963 -6391 -3553 NA [33,] 5793 11686 -4499 NA [34,] -4307 -11514 1144 NA [35,] 11449 4481 2957 NA [36,] -5934 -3512 -2672 NA [37,] -2992 6168 -1842 NA [38,] 4814 -3366 -3760 NA [39,] -4889 752 1179 NA [40,] -1813 -20 -1737 NA [41,] -3379 -3815 8407 NA [42,] 7621 -3263 -9789 NA [43,] -9096 -3404 905 NA [44,] 900 44 -2992 NA [45,] -5433 -1738 1412 NA [46,] 7664 3941 -2789 NA [47,] -231 21452 1250 NA [48,] 598 -17751 -1084 NA [49,] -533 -1764 1864 NA [50,] 5906 3193 1226 NA [51,] -564 -372 -5987 NA [52,] -2782 2450 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/1nlnt1389428935.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/29ptd1389428935.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/3igk71389428935.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/4com11389428935.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] [1,] 189561.0 184439.0 182985 183828 184036.0 183214.0 183464.0 173718.0 [2,] 199957.5 193322.0 201705 196933 194618.5 190327.5 191742.5 188441.5 [3,] 210354.0 202205.0 220425 210038 205201.0 197441.0 200021.0 203165.0 [4,] 215728.0 211548.5 222981 214686 210663.5 211007.0 207006.0 212515.0 [5,] 221102.0 220892.0 225537 219334 216126.0 224573.0 213991.0 221865.0 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [1,] 180210.0 171252.0 172705.0 174006 172043.0 169445.0 169449 177073 170799.0 [2,] 190954.5 186127.5 194857.0 192709 190852.0 191047.5 181530 184808 183419.5 [3,] 201699.0 201003.0 217009.0 211412 209661.0 212650.0 193611 192543 196040.0 [4,] 208540.5 207482.5 224658.5 264761 214661.5 213118.5 201459 201542 203324.5 [5,] 215382.0 213962.0 232308.0 318110 219662.0 213587.0 209307 210541 210609.0 [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 171648.0 172220.0 165795.0 167466.0 165528.0 162851.0 165864.0 162094.0 [2,] 181527.5 180973.0 178533.5 177695.0 181870.5 181101.5 182576.5 178784.5 [3,] 191407.0 189726.0 191272.0 187924.0 198213.0 199352.0 199289.0 195475.0 [4,] 200158.5 198633.5 199485.5 196464.5 201980.0 201498.0 214613.0 204960.5 [5,] 208910.0 207541.0 207699.0 205005.0 205747.0 203644.0 229937.0 214446.0 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [1,] 162385.0 164293 165983.0 159680.0 161739.0 159302.0 167795.0 164242.0 [2,] 180215.0 180954 177499.0 174674.0 175429.5 176735.0 180049.5 175077.5 [3,] 198045.0 197615 189015.0 189668.0 189120.0 194168.0 192304.0 185913.0 [4,] 204119.5 202075 202769.5 193955.5 198697.0 200830.5 203914.5 196737.5 [5,] 210194.0 206535 216524.0 198243.0 208274.0 207493.0 215525.0 207562.0 [,34] [,35] [,36] [,37] [,38] [,39] [,40] [,41] [1,] 159743 160887.0 163844.0 161172.0 159330.0 155570.0 156749.0 155012.0 [2,] 178671 173486.0 177205.0 174113.0 176276.0 172713.0 173678.5 172800.0 [3,] 197599 186085.0 190566.0 187054.0 193222.0 189856.0 190608.0 190588.0 [4,] 205477 197566.5 205531.5 200808.5 202396.5 203120.5 201052.0 200135.5 [5,] 213355 209048.0 220497.0 214563.0 211571.0 216385.0 211496.0 209683.0 [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 163419.0 153630.0 154535.0 151543.0 152955.0 150166.0 151416.0 150332.0 [2,] 175096.0 168570.0 167320.5 165846.5 165683.5 166259.5 177610.5 168193.0 [3,] 186773.0 183510.0 180106.0 180150.0 178412.0 182353.0 203805.0 186054.0 [4,] 196538.5 198717.5 192467.5 192939.5 189354.0 195156.5 205767.0 197190.5 [5,] 206304.0 213925.0 204829.0 205729.0 200296.0 207960.0 207729.0 208327.0 [,50] [,51] [,52] [1,] 152196 153422.0 147435.0 [2,] 168243 170452.5 167273.0 [3,] 184290 187483.0 187111.0 [4,] 196042 200591.5 200123.5 [5,] 207794 213700.0 213136.0 $n [1] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 [39] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 195967.9 185578.5 201016.7 193843.5 190564.5 178576.9 186097.4 181204.8 [2,] 224740.1 218831.5 239833.3 226232.5 219837.5 216305.1 213944.6 225125.2 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 185656.8 181522.7 189823.7 145685.2 187941.7 192516.5 175431.5 177278 [2,] 217741.2 220483.3 244194.3 277138.8 231380.3 232783.5 211790.5 207808 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 177882.4 174411.6 173615.9 172159.3 170802.2 179868.8 180746 170064.9 [2,] 214197.6 208402.4 205836.1 210384.7 205045.8 216557.2 217958 228513.1 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 171596.9 176239 178348.1 165962.9 172079.2 167895.1 172187.8 170534 [2,] 219353.1 219851 216881.9 212067.1 207256.8 210344.9 216148.2 214074 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 166154.5 173146.2 164118.4 164726.2 162702 169394.5 162117.9 165637.5 [2,] 205671.5 222051.8 208051.6 216405.8 211406 217049.5 217594.1 215578.5 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 165652.2 167212.9 156009 157166.6 155435.4 156819.5 155992.8 178120.3 [2,] 215523.8 206333.1 211011 203045.4 204864.6 200004.5 208713.2 229489.7 [,49] [,50] [,51] [,52] [1,] 159602.1 158931.4 159989.8 157144.3 [2,] 212505.9 209648.6 214976.2 217077.7 $out numeric(0) $group numeric(0) $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" Warning message: In bxp(list(stats = c(189561, 199957.5, 210354, 215728, 221102, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5822b1389428935.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] [1,] -8149.0 -1454.0 -10387 -4837.0 -7760.0 -10582 -9746 -6483.0 -8958 1453 [2,] -6635.5 1595.5 -8295 -4022.5 -4291.0 -5166 -3301 -3974.5 -5189 2250 [3,] -5122.0 4645.0 -6203 -3208.0 -822.0 250 3144 -1466.0 -1420 3047 [4,] -2666.0 11432.5 -2680 -1500.0 3812.5 1415 5509 2513.0 -1058 17176 [5,] -210.0 18220.0 843 208.0 8447.0 2580 7874 6492.0 -696 31305 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [1,] -20896.0 -98448.0 -7012 -19976.0 -1068 -6274.0 -4633 -1681.0 -6425.0 [2,] -9797.5 -50205.5 -4805 -11659.5 83 -3103.0 -3166 -1525.0 -3133.5 [3,] 1301.0 -1963.0 -2598 -3343.0 1234 68.0 -1699 -1369.0 158.0 [4,] 51201.0 -1857.0 664 -1669.5 4429 1782.5 -425 -398.5 852.0 [5,] 101101.0 -1751.0 3926 4.0 7624 3497.0 849 572.0 1546.0 [,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] -3348.0 -1938.0 -2677 -63 -15491.0 -4252.0 -3659.0 -8600.0 -18281 [2,] -3021.0 -598.0 -2390 1475 -9652.5 -1980.5 -2044.5 -3455.0 -12292 [3,] -2694.0 742.0 -2103 3013 -3814.0 291.0 -430.0 1690.0 -6303 [4,] -511.5 5515.5 -482 14653 -3792.0 1430.5 739.0 5839.5 -2825 [5,] 1671.0 10289.0 1139 26293 -3770.0 2570.0 1908.0 9989.0 653 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -548.0 -2437.0 -1864.0 -7963 -4499.0 -11514.0 2957 -5934 -2992 -3760 [2,] 755.5 -1609.0 3084.0 -7177 647.0 -7910.5 3719 -4723 -2417 -3563 [3,] 2059.0 -781.0 8032.0 -6391 5793.0 -4307.0 4481 -3512 -1842 -3366 [4,] 6045.0 2133.5 8262.5 -4972 8739.5 -1581.5 7965 -3092 2163 724 [5,] 10031.0 5048.0 8493.0 -3553 11686.0 1144.0 11449 -2672 6168 4814 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] -4889.0 -1813.0 -3815 -9789 -9096.0 -2992 -5433.0 -2789.0 -231.0 -17751.0 [2,] -2068.5 -1775.0 -3597 -6526 -6250.0 -1474 -3585.5 576.0 509.5 -9417.5 [3,] 752.0 -1737.0 -3379 -3263 -3404.0 44 -1738.0 3941.0 1250.0 -1084.0 [4,] 965.5 -878.5 2514 2179 -1249.5 472 -163.0 5802.5 11351.0 -243.0 [5,] 1179.0 -20.0 8407 7621 905.0 900 1412.0 7664.0 21452.0 598.0 [,49] [,50] [,51] [,52] [1,] -1764.0 1226.0 -5987.0 -2782 [2,] -1148.5 2209.5 -3275.5 -2782 [3,] -533.0 3193.0 -564.0 -166 [4,] 665.5 4549.5 -468.0 2450 [5,] 1864.0 5906.0 -372.0 2450 $n [1] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 [39] 3 3 3 3 3 3 3 3 3 3 3 3 3 2 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -8743.031 -4328.443 -11325.078 -5509.0584 -8214.121 -5753.277 -4892.6 [2,] -1500.969 13618.443 -1080.922 -906.9416 6570.121 6253.277 11180.6 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] -7383.985 -5188.354 -10568.7 -54342.65 -46067.15 -7586.895 -12456.012 [2,] 4451.985 2348.354 16662.7 56944.65 42141.15 2390.895 5770.012 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] -2730.48 -4388.619 -4199.377 -2396.6084 -3477.627 -4983.1996 -4834.817 [2,] 5198.48 4524.619 801.377 -341.3916 3793.627 -404.8004 6318.817 [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] -3843.5032 -9008.149 -9160.027 -2820.56 -2969.146 -6788.568 -14938.924 [2,] -362.4968 15034.149 1532.027 3402.56 2109.146 10168.568 2332.924 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [1,] -2766.153 -4194.959 3308.103 -8402.431 -1589.087 -10080.399 607.7418 [2,] 6884.153 2632.959 12755.897 -4379.569 13175.087 1466.399 8354.2582 [,36] [,37] [,38] [,39] [,40] [,41] [,42] [1,] -4999.82 -6019.937 -7276.659 -2015.656 -2554.7993 -8953.536 -11203.818 [2,] -2024.18 2335.937 544.659 3519.656 -919.2007 2195.536 4677.818 [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] -7965.523 -1731.167 -4860.05 -826.6835 -8639.762 -9453.102 -2187.755 [2,] 1157.523 1819.167 1384.05 8708.6835 11139.762 7285.102 1121.755 [,50] [,51] [,52] [1,] 1058.421 -3125.039 -6011.341 [2,] 5327.579 1997.039 5679.341 $out numeric(0) $group numeric(0) $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" Warning message: In bxp(list(stats = c(-8149, -6635.5, -5122, -2666, -210, -1454, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6db5p1389428935.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] [1,] 198243.0 178412.0 147435.0 NA [2,] 207714.0 187703.5 158025.5 NA [3,] 211533.5 192882.5 164267.5 NA [4,] 216255.5 200512.0 172131.5 NA [5,] 225537.0 213587.0 189561.0 NA $n [1] 52 52 52 0 $conf [,1] [,2] [,3] [,4] [1,] 209662 190076.1 161176.8 NA [2,] 213405 195688.9 167358.2 NA $out [1] 318110 229937 220425 232308 $group [1] 1 1 2 2 $names [1] "1" "2" "3" NA > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/73rpp1389428935.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,] 177221.0 178412.0 177518.8 [2,] 185702.2 187703.5 186187.6 [3,] 189829.2 192882.5 191334.0 [4,] 198462.0 200512.0 198984.5 [5,] 209649.0 217009.0 212343.0 $n [1] 52 52 52 $conf [,1] [,2] [,3] [1,] 187033.4 190076.1 188530.1 [2,] 192624.9 195688.9 194137.9 $out [1] 234509.3 220425.0 228735.0 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > try(system("convert tmp/1nlnt1389428935.ps tmp/1nlnt1389428935.png",intern=TRUE)) character(0) > try(system("convert tmp/29ptd1389428935.ps tmp/29ptd1389428935.png",intern=TRUE)) character(0) > try(system("convert tmp/3igk71389428935.ps tmp/3igk71389428935.png",intern=TRUE)) character(0) > try(system("convert tmp/4com11389428935.ps tmp/4com11389428935.png",intern=TRUE)) character(0) > try(system("convert tmp/5822b1389428935.ps tmp/5822b1389428935.png",intern=TRUE)) character(0) > try(system("convert tmp/6db5p1389428935.ps tmp/6db5p1389428935.png",intern=TRUE)) character(0) > try(system("convert tmp/73rpp1389428935.ps tmp/73rpp1389428935.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 6.589 1.519 8.061