R version 3.3.2 (2016-10-31) -- "Sincere Pumpkin Patch" 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(2341 + ,2115 + ,2402 + ,2180 + ,2453 + ,2507 + ,2679 + ,2622 + ,2618 + ,2648 + ,2523 + ,2473 + ,2513 + ,2466 + ,2544 + ,2537 + ,2564 + ,2582 + ,2716 + ,2904 + ,2851 + ,2932 + ,2772 + ,2811 + ,2935 + ,2783 + ,3003 + ,2995 + ,3127 + ,2985 + ,3287 + ,3236 + ,3252 + ,3228 + ,2856 + ,3176 + ,3362 + ,3036 + ,3330 + ,3251 + ,3318 + ,3238 + ,3597 + ,3708 + ,3902 + ,3745 + ,3426 + ,3526 + ,3483 + ,3458 + ,3824 + ,3696 + ,3518 + ,3814 + ,3996 + ,4136 + ,4037 + ,3915 + ,3760 + ,3955 + ,4160 + ,4115 + ,4202 + ,4018 + ,4233 + ,4029 + ,4401 + ,4645 + ,4491 + ,4379 + ,4394 + ,4472 + ,4614 + ,4160 + ,4328 + ,4202 + ,4635 + ,4542 + ,4920 + ,4774 + ,4698 + ,4916 + ,4703 + ,4616 + ,4873 + ,4375 + ,4801 + ,4427 + ,4684 + ,4648 + ,5225 + ,5174 + ,5181 + ,5266 + ,4839 + ,5032 + ,5221 + ,4658 + ,5014 + ,4980 + ,4952 + ,4946 + ,5365 + ,5456 + ,5397 + ,5436 + ,4995 + ,5019 + ,5249 + ,4799 + ,5137 + ,4979 + ,4951 + ,5265 + ,5612 + ,5572 + ,5403 + ,5373 + ,5252 + ,5437 + ,5296 + ,5011 + ,5294 + ,5335 + ,5398 + ,5396 + ,5724 + ,5898 + ,5718 + ,5625 + ,5380 + ,5488 + ,5678 + ,5224 + ,5596 + ,5184 + ,5620 + ,5531 + ,5816 + ,6086 + ,6175 + ,6112 + ,5813 + ,5740 + ,5821 + ,5294 + ,5881 + ,5589 + ,5845 + ,5706 + ,6355 + ,6404 + ,6426 + ,6375 + ,5869 + ,5994 + ,6105 + ,5792 + ,6011 + ,5968 + ,6255 + ,6208 + ,6897 + ,6814 + ,6897 + ,6596 + ,6188 + ,6406 + ,6548 + ,5842 + ,6555 + ,6424 + ,6596 + ,6645 + ,7203 + ,7128 + ,7133 + ,6778 + ,6593 + ,6591 + ,6120 + ,5612 + ,6070 + ,5983 + ,6145 + ,6303 + ,6588 + ,6640 + ,6719 + ,6575 + ,6487 + ,6510 + ,6365 + ,5844 + ,5974 + ,5880 + ,6279 + ,6342 + ,6598 + ,6801 + ,6529 + ,6369 + ,6028 + ,6187 + ,6164 + ,5866 + ,6198 + ,5898 + ,6462 + ,6063 + ,6496 + ,6678 + ,6554 + ,6513 + ,6210 + ,5928 + ,6268 + ,5582 + ,5869 + ,5764 + ,6082 + ,6062 + ,6810 + ,6727 + ,6537 + ,6175 + ,6014 + ,6109 + ,2341 + ,2115 + ,2402 + ,2180 + ,2453 + ,2507 + ,2679 + ,2622 + ,2618 + ,2648 + ,2523 + ,2473 + ,2513 + ,2466 + ,2544 + ,2537 + ,2564 + ,2582 + ,2716 + ,2904 + ,2851 + ,2932 + ,2772 + ,2811 + ,2935 + ,2783 + ,3003 + ,2995 + ,3127 + ,2985 + ,3287 + ,3236 + ,3252 + ,3228 + ,2856 + ,3176 + ,3362 + ,3036 + ,3330 + ,3251 + ,3318 + ,3238 + ,3597 + ,3708 + ,3902 + ,3745 + ,3426 + ,3526 + ,3483 + ,3458 + ,3824 + ,3696 + ,3518 + ,3814 + ,3996 + ,4136 + ,4037 + ,3915 + ,3760 + ,3955 + ,4160 + ,4115 + ,4202 + ,4018 + ,4233 + ,4029 + ,4401 + ,4645 + ,4491 + ,4379 + ,4394 + ,4472 + ,4614 + ,4160 + ,4328 + ,4202 + ,4635 + ,4542 + ,4920 + ,4774 + ,4698 + ,4916 + ,4703 + ,4616 + ,4873 + ,4375 + ,4801 + ,4427 + ,4684 + ,4648 + ,5225 + ,5174 + ,5181 + ,5266 + ,4839 + ,5032 + ,5221 + ,4658 + ,5014 + ,4980 + ,4952 + ,4946 + ,5365 + ,5456 + ,5397 + ,5436 + ,4995 + ,5019 + ,5249 + ,4799 + ,5137 + ,4979 + ,4951 + ,5265 + ,5612 + ,5572 + ,5403 + ,5373 + ,5252 + ,5437 + ,5296 + ,5011 + ,5294 + ,5335 + ,5398 + ,5396 + ,5724 + ,5898 + ,5718 + ,5625 + ,5380 + ,5488 + ,5678 + ,5224 + ,5596 + ,5184 + ,5620 + ,5531 + ,5816 + ,6086 + ,6175 + ,6112 + ,5813 + ,5740 + ,5821 + ,5294 + ,5881 + ,5589 + ,5845 + ,5706 + ,6355 + ,6404 + ,6426 + ,6375 + ,5869 + ,5994 + ,6105 + ,5792 + ,6011 + ,5968 + ,6255 + ,6208 + ,6897 + ,6814 + ,6897 + ,6596 + ,6188 + ,6406 + ,6548 + ,5842 + ,6555 + ,6424 + ,6596 + ,6645 + ,7203 + ,7128 + ,7133 + ,6778 + ,6593 + ,6591 + ,6120 + ,5612 + ,6070 + ,5983 + ,6145 + ,6303 + ,6588 + ,6640 + ,6719 + ,6575 + ,6487 + ,6510 + ,6365 + ,5844 + ,5974 + ,5880 + ,6279 + ,6342 + ,6598 + ,6801 + ,6529 + ,6369 + ,6028 + ,6187 + ,6164 + ,5866 + ,6198 + ,5898 + ,6462 + ,6063 + ,6496 + ,6678 + ,6554 + ,6513 + ,6210 + ,5928 + ,6268 + ,5582 + ,5869 + ,5764 + ,6082 + ,6062 + ,6810 + ,6727 + ,6537 + ,6175 + ,6014 + ,6109) > par1 = '12' > par1 <- '12' > #'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] 456 > (np <- floor(n / par1)) [1] 38 > 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] 38 38 38 38 38 38 38 38 38 38 38 38 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 2341 2513 2935 3362 3483 4160 4614 4873 5221 5249 5296 5678 5821 [2,] 2115 2466 2783 3036 3458 4115 4160 4375 4658 4799 5011 5224 5294 [3,] 2402 2544 3003 3330 3824 4202 4328 4801 5014 5137 5294 5596 5881 [4,] 2180 2537 2995 3251 3696 4018 4202 4427 4980 4979 5335 5184 5589 [5,] 2453 2564 3127 3318 3518 4233 4635 4684 4952 4951 5398 5620 5845 [6,] 2507 2582 2985 3238 3814 4029 4542 4648 4946 5265 5396 5531 5706 [7,] 2679 2716 3287 3597 3996 4401 4920 5225 5365 5612 5724 5816 6355 [8,] 2622 2904 3236 3708 4136 4645 4774 5174 5456 5572 5898 6086 6404 [9,] 2618 2851 3252 3902 4037 4491 4698 5181 5397 5403 5718 6175 6426 [10,] 2648 2932 3228 3745 3915 4379 4916 5266 5436 5373 5625 6112 6375 [11,] 2523 2772 2856 3426 3760 4394 4703 4839 4995 5252 5380 5813 5869 [12,] 2473 2811 3176 3526 3955 4472 4616 5032 5019 5437 5488 5740 5994 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 6105 6548 6120 6365 6164 6268 2341 2513 2935 3362 3483 4160 [2,] 5792 5842 5612 5844 5866 5582 2115 2466 2783 3036 3458 4115 [3,] 6011 6555 6070 5974 6198 5869 2402 2544 3003 3330 3824 4202 [4,] 5968 6424 5983 5880 5898 5764 2180 2537 2995 3251 3696 4018 [5,] 6255 6596 6145 6279 6462 6082 2453 2564 3127 3318 3518 4233 [6,] 6208 6645 6303 6342 6063 6062 2507 2582 2985 3238 3814 4029 [7,] 6897 7203 6588 6598 6496 6810 2679 2716 3287 3597 3996 4401 [8,] 6814 7128 6640 6801 6678 6727 2622 2904 3236 3708 4136 4645 [9,] 6897 7133 6719 6529 6554 6537 2618 2851 3252 3902 4037 4491 [10,] 6596 6778 6575 6369 6513 6175 2648 2932 3228 3745 3915 4379 [11,] 6188 6593 6487 6028 6210 6014 2523 2772 2856 3426 3760 4394 [12,] 6406 6591 6510 6187 5928 6109 2473 2811 3176 3526 3955 4472 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 4614 4873 5221 5249 5296 5678 5821 6105 6548 6120 6365 6164 [2,] 4160 4375 4658 4799 5011 5224 5294 5792 5842 5612 5844 5866 [3,] 4328 4801 5014 5137 5294 5596 5881 6011 6555 6070 5974 6198 [4,] 4202 4427 4980 4979 5335 5184 5589 5968 6424 5983 5880 5898 [5,] 4635 4684 4952 4951 5398 5620 5845 6255 6596 6145 6279 6462 [6,] 4542 4648 4946 5265 5396 5531 5706 6208 6645 6303 6342 6063 [7,] 4920 5225 5365 5612 5724 5816 6355 6897 7203 6588 6598 6496 [8,] 4774 5174 5456 5572 5898 6086 6404 6814 7128 6640 6801 6678 [9,] 4698 5181 5397 5403 5718 6175 6426 6897 7133 6719 6529 6554 [10,] 4916 5266 5436 5373 5625 6112 6375 6596 6778 6575 6369 6513 [11,] 4703 4839 4995 5252 5380 5813 5869 6188 6593 6487 6028 6210 [12,] 4616 5032 5019 5437 5488 5740 5994 6406 6591 6510 6187 5928 [,38] [,39] [1,] 6268 NA [2,] 5582 NA [3,] 5869 NA [4,] 5764 NA [5,] 6082 NA [6,] 6062 NA [7,] 6810 NA [8,] 6727 NA [9,] 6537 NA [10,] 6175 NA [11,] 6014 NA [12,] 6109 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] -226 -47 -152 -326 -25 -45 -454 -498 -563 -450 -285 -454 -527 [2,] 287 78 220 294 366 87 168 426 356 338 283 372 587 [3,] -222 -7 -8 -79 -128 -184 -126 -374 -34 -158 41 -412 -292 [4,] 273 27 132 67 -178 215 433 257 -28 -28 63 436 256 [5,] 54 18 -142 -80 296 -204 -93 -36 -6 314 -2 -89 -139 [6,] 172 134 302 359 182 372 378 577 419 347 328 285 649 [7,] -57 188 -51 111 140 244 -146 -51 91 -40 174 270 49 [8,] -4 -53 16 194 -99 -154 -76 7 -59 -169 -180 89 22 [9,] 30 81 -24 -157 -122 -112 218 85 39 -30 -93 -63 -51 [10,] -125 -160 -372 -319 -155 15 -213 -427 -441 -121 -245 -299 -506 [11,] -50 39 320 100 195 78 -87 193 24 185 108 -73 125 [12,] 40 124 186 -43 205 142 257 189 230 -141 190 81 111 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] -313 -706 -508 -521 -298 -686 -226 -47 -152 -326 -25 -45 [2,] 219 713 458 130 332 287 287 78 220 294 366 87 [3,] -43 -131 -87 -94 -300 -105 -222 -7 -8 -79 -128 -184 [4,] 287 172 162 399 564 318 273 27 132 67 -178 215 [5,] -47 49 158 63 -399 -20 54 18 -142 -80 296 -204 [6,] 689 558 285 256 433 748 172 134 302 359 182 372 [7,] -83 -75 52 203 182 -83 -57 188 -51 111 140 244 [8,] 83 5 79 -272 -124 -190 -4 -53 16 194 -99 -154 [9,] -301 -355 -144 -160 -41 -362 30 81 -24 -157 -122 -112 [10,] -408 -185 -88 -341 -303 -161 -125 -160 -372 -319 -155 15 [11,] 218 -2 23 159 -282 95 -50 39 320 100 195 78 [12,] 142 -471 -145 -23 340 -3768 40 124 186 -43 205 142 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] -454 -498 -563 -450 -285 -454 -527 -313 -706 -508 -521 -298 [2,] 168 426 356 338 283 372 587 219 713 458 130 332 [3,] -126 -374 -34 -158 41 -412 -292 -43 -131 -87 -94 -300 [4,] 433 257 -28 -28 63 436 256 287 172 162 399 564 [5,] -93 -36 -6 314 -2 -89 -139 -47 49 158 63 -399 [6,] 378 577 419 347 328 285 649 689 558 285 256 433 [7,] -146 -51 91 -40 174 270 49 -83 -75 52 203 182 [8,] -76 7 -59 -169 -180 89 22 83 5 79 -272 -124 [9,] 218 85 39 -30 -93 -63 -51 -301 -355 -144 -160 -41 [10,] -213 -427 -441 -121 -245 -299 -506 -408 -185 -88 -341 -303 [11,] -87 193 24 185 108 -73 125 218 -2 23 159 -282 [12,] 257 189 230 -141 190 81 111 142 -471 -145 -23 340 [,38] [,39] [1,] -686 NA [2,] 287 NA [3,] -105 NA [4,] 318 NA [5,] -20 NA [6,] 748 NA [7,] -83 NA [8,] -190 NA [9,] -362 NA [10,] -161 NA [11,] 95 NA [12,] 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/1bj0a1489659127.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/2unuk1489659127.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/3rvoo1489659127.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/4aqiv1489659127.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(arr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 2341 2115 2402 2180 2453 2507 2679 2622 2618 2648 2523 2473 [2,] 3483 3458 3824 3696 3518 3814 3996 4136 4037 3915 3760 3955 [3,] 5249 4799 5137 4980 4952 5265 5612 5572 5403 5436 5252 5437 [4,] 6120 5612 5974 5880 6145 6063 6588 6678 6537 6375 6028 6109 [5,] 6548 5866 6555 6424 6596 6645 7203 7128 7133 6778 6593 6591 $n [1] 38 38 38 38 38 38 38 38 38 38 38 38 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 4573.111 4246.909 4585.934 4420.219 4278.674 4688.559 4947.645 4920.46 [2,] 5924.889 5351.091 5688.066 5539.781 5625.326 5841.441 6276.355 6223.54 [,9] [,10] [,11] [,12] [1,] 4762.225 4805.478 4670.689 4884.909 [2,] 6043.775 6066.522 5833.311 5989.091 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5evvn1489659127.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(darr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] -706 78 -412 -178 -204 134 -146 -272 -362 -506 -282 -145 [2,] -521 219 -222 63 -93 285 -57 -154 -157 -372 -2 -23 [3,] -450 294 -126 215 -20 359 52 -53 -63 -245 95 124 [4,] -226 372 -43 318 54 558 182 22 30 -155 185 190 [5,] -25 587 41 564 158 748 270 194 218 15 320 340 $n [1] 38 38 38 38 38 38 38 38 38 38 38 37 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -525.6114 254.7846 -171.87946 149.641 -57.67755 289.0274 -9.258053 [2,] -374.3886 333.2154 -80.12054 280.359 17.67755 428.9726 113.258053 [,8] [,9] [,10] [,11] [,12] [1,] -98.110533 -110.92994 -300.6192 47.07006 68.67316 [2,] -7.889467 -15.07006 -189.3808 142.92994 179.32684 $out [1] 713 713 296 314 -399 296 314 -399 -471 -3768 -471 $group [1] 2 2 5 5 5 5 5 5 12 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-706, -521, -450, -226, -25, 78, 219, 294, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6f6to1489659127.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] [1,] 2115.0 2466.0 2783 3036.0 3458.0 4018.0 4160.0 4375.0 4658.0 4799.0 5011.0 [2,] 2371.5 2540.5 2960 3284.5 3607.0 4137.5 4435.0 4666.0 4966.0 5058.0 5315.5 [3,] 2490.0 2649.0 3065 3394.0 3819.0 4306.0 4625.5 4856.0 5016.5 5258.5 5397.0 [4,] 2620.0 2831.0 3232 3652.5 3975.5 4436.5 4738.5 5177.5 5381.0 5420.0 5671.5 [5,] 2679.0 2932.0 3287 3902.0 4136.0 4645.0 4920.0 5266.0 5456.0 5612.0 5898.0 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] 5184.0 5294.0 5792.0 6424.0 5612.0 5844.0 5866.0 5582.0 2115.0 2466.0 [2,] 5563.5 5763.5 6058.0 6551.5 6095.0 6001.0 5995.5 5941.5 2371.5 2540.5 [3,] 5709.0 5875.0 6231.5 6594.5 6395.0 6310.5 6204.0 6095.5 2490.0 2649.0 [4,] 5951.0 6365.0 6705.0 6953.0 6581.5 6449.0 6504.5 6402.5 2620.0 2831.0 [5,] 6175.0 6426.0 6897.0 7203.0 6719.0 6801.0 6678.0 6810.0 2679.0 2932.0 [,22] [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [1,] 2783 3036.0 3458.0 4018.0 4160.0 4375.0 4658.0 4799.0 5011.0 5184.0 [2,] 2960 3284.5 3607.0 4137.5 4435.0 4666.0 4966.0 5058.0 5315.5 5563.5 [3,] 3065 3394.0 3819.0 4306.0 4625.5 4856.0 5016.5 5258.5 5397.0 5709.0 [4,] 3232 3652.5 3975.5 4436.5 4738.5 5177.5 5381.0 5420.0 5671.5 5951.0 [5,] 3287 3902.0 4136.0 4645.0 4920.0 5266.0 5456.0 5612.0 5898.0 6175.0 [,32] [,33] [,34] [,35] [,36] [,37] [,38] [,39] [1,] 5294.0 5792.0 6424.0 5612.0 5844.0 5866.0 5582.0 NA [2,] 5763.5 6058.0 6551.5 6095.0 6001.0 5995.5 5941.5 NA [3,] 5875.0 6231.5 6594.5 6395.0 6310.5 6204.0 6095.5 NA [4,] 6365.0 6705.0 6953.0 6581.5 6449.0 6504.5 6402.5 NA [5,] 6426.0 6897.0 7203.0 6719.0 6801.0 6678.0 6810.0 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 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2376.657 2516.501 2940.939 3226.153 3650.925 4169.624 4487.072 4622.701 [2,] 2603.343 2781.499 3189.061 3561.847 3987.075 4442.376 4763.928 5089.299 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 4827.216 5093.389 5234.626 5532.259 5600.652 5936.399 6411.373 6173.104 [2,] 5205.784 5423.611 5559.374 5885.741 6149.348 6526.601 6777.627 6616.896 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 6106.164 5971.842 5885.235 2376.657 2516.501 2940.939 3226.153 3650.925 [2,] 6514.836 6436.158 6305.765 2603.343 2781.499 3189.061 3561.847 3987.075 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 4169.624 4487.072 4622.701 4827.216 5093.389 5234.626 5532.259 5600.652 [2,] 4442.376 4763.928 5089.299 5205.784 5423.611 5559.374 5885.741 6149.348 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [1,] 5936.399 6411.373 6173.104 6106.164 5971.842 5885.235 NA [2,] 6526.601 6777.627 6616.896 6514.836 6436.158 6305.765 NA $out [1] 5842 5842 $group [1] 15 34 $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" NA Warning message: In bxp(list(stats = c(2115, 2371.5, 2490, 2620, 2679, 2466, 2540.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7ga2h1489659127.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,] 4528.000 4799.0 4613.375 [2,] 4864.342 5058.5 4898.625 [3,] 4926.816 5258.5 4968.375 [4,] 5243.184 5436.5 5266.688 [5,] 5337.000 5612.0 5465.875 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 4754.023 5086.092 4800.499 [2,] 5099.608 5430.908 5136.251 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(4528, 4864.34210526316, 4926.81578947368, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1bj0a1489659127.ps tmp/1bj0a1489659127.png",intern=TRUE)) character(0) > try(system("convert tmp/2unuk1489659127.ps tmp/2unuk1489659127.png",intern=TRUE)) character(0) > try(system("convert tmp/3rvoo1489659127.ps tmp/3rvoo1489659127.png",intern=TRUE)) character(0) > try(system("convert tmp/4aqiv1489659127.ps tmp/4aqiv1489659127.png",intern=TRUE)) character(0) > try(system("convert tmp/5evvn1489659127.ps tmp/5evvn1489659127.png",intern=TRUE)) character(0) > try(system("convert tmp/6f6to1489659127.ps tmp/6f6to1489659127.png",intern=TRUE)) character(0) > try(system("convert tmp/7ga2h1489659127.ps tmp/7ga2h1489659127.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.785 0.215 3.095