R version 2.6.2 (2008-02-08) Copyright (C) 2008 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. Natural language support but running in an English locale 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(1.327 + ,1.3231 + ,1.3106 + ,1.3084 + ,1.3006 + ,1.3018 + ,1.2988 + ,1.2984 + ,1.2893 + ,1.2941 + ,1.2953 + ,1.2909 + ,1.2922 + ,1.2958 + ,1.2936 + ,1.304 + ,1.3005 + ,1.2978 + ,1.2901 + ,1.2921 + ,1.2972 + ,1.2954 + ,1.302 + ,1.302 + ,1.2925 + ,1.2955 + ,1.2987 + ,1.2991 + ,1.3007 + ,1.2956 + ,1.3022 + ,1.3082 + ,1.3137 + ,1.3119 + ,1.3132 + ,1.3145 + ,1.3145 + ,1.3106 + ,1.3134 + ,1.316 + ,1.323 + ,1.3211 + ,1.3225 + ,1.3163 + ,1.3083 + ,1.31 + ,1.3135 + ,1.3152 + ,1.3155 + ,1.3156 + ,1.3218 + ,1.3183 + ,1.3226 + ,1.3325 + ,1.3303 + ,1.3296 + ,1.3293 + ,1.3351 + ,1.3327 + ,1.3265 + ,1.3347 + ,1.3348 + ,1.3352 + ,1.3318 + ,1.3366 + ,1.3358 + ,1.3352 + ,1.3373 + ,1.3426 + ,1.3418 + ,1.3467 + ,1.3532 + ,1.355 + ,1.3549 + ,1.3577 + ,1.3601 + ,1.3606 + ,1.3557 + ,1.3582 + ,1.3649 + ,1.3596 + ,1.3643 + ,1.3605 + ,1.3588 + ,1.3613 + ,1.3561 + ,1.3615 + ,1.3558 + ,1.3535 + ,1.3527 + ,1.3486 + ,1.3549 + ,1.3538 + ,1.3574 + ,1.3516 + ,1.3477 + ,1.3444 + ,1.3454 + ,1.349 + ,1.3448 + ,1.3441 + ,1.3453 + ,1.3509 + ,1.342 + ,1.3453 + ,1.3436 + ,1.3482 + ,1.3532 + ,1.3513 + ,1.347 + ,1.3349 + ,1.3355 + ,1.3345 + ,1.3287 + ,1.3304 + ,1.3314 + ,1.3404 + ,1.3403 + ,1.3427 + ,1.3397 + ,1.3441 + ,1.3461 + ,1.346 + ,1.3438 + ,1.3467 + ,1.3505 + ,1.3588 + ,1.3601 + ,1.3618 + ,1.364 + ,1.3596 + ,1.3621 + ,1.3666 + ,1.3753 + ,1.3788 + ,1.3782 + ,1.3781 + ,1.3771 + ,1.3779 + ,1.382 + ,1.3803 + ,1.3821 + ,1.3833 + ,1.3743 + ,1.3722 + ,1.3651 + ,1.3659 + ,1.3707 + ,1.3663 + ,1.3664 + ,1.3694 + ,1.3818 + ,1.3794 + ,1.3794 + ,1.3729 + ,1.365 + ,1.3651 + ,1.3591 + ,1.3476 + ,1.3405 + ,1.3454 + ,1.3476 + ,1.3508 + ,1.3493 + ,1.3574 + ,1.3615 + ,1.3658 + ,1.3664 + ,1.3631 + ,1.361 + ,1.3705 + ,1.3632 + ,1.358 + ,1.3588 + ,1.3669 + ,1.3696 + ,1.3795 + ,1.3824 + ,1.3885 + ,1.3897 + ,1.386 + ,1.3877 + ,1.3867 + ,1.3975 + ,1.403 + ,1.4049 + ,1.4113 + ,1.4106 + ,1.4127 + ,1.418 + ,1.4179 + ,1.4232 + ,1.4165 + ,1.4195 + ,1.4109 + ,1.4136 + ,1.4089 + ,1.4037 + ,1.4146 + ,1.4199 + ,1.4173 + ,1.4226 + ,1.415 + ,1.42 + ,1.4299 + ,1.4288 + ,1.4166 + ,1.4254 + ,1.423 + ,1.4309 + ,1.4384 + ,1.4391 + ,1.4407 + ,1.4447 + ,1.4423 + ,1.4479 + ,1.4488 + ,1.4547 + ,1.4722 + ,1.4666 + ,1.4683 + ,1.4579 + ,1.4607 + ,1.47 + ,1.4639 + ,1.4651 + ,1.4654 + ,1.4785 + ,1.4814 + ,1.4829 + ,1.4809 + ,1.4845 + ,1.4874 + ,1.4747 + ,1.4738 + ,1.4761 + ,1.4666 + ,1.4741 + ,1.472 + ,1.4554 + ,1.4649 + ,1.4718 + ,1.4672 + ,1.4675 + ,1.4683 + ,1.4509 + ,1.4393 + ,1.4416 + ,1.4385 + ,1.4349 + ,1.438 + ,1.4398 + ,1.4516 + ,1.4692 + ,1.4721) > 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] 255 > (np <- floor(n / par1)) [1] 21 > 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] 22 22 22 21 21 21 21 21 21 21 21 21 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 1.3270 1.2922 1.2925 1.3145 1.3155 1.3347 1.3550 1.3613 1.3444 1.3513 [2,] 1.3231 1.2958 1.2955 1.3106 1.3156 1.3348 1.3549 1.3561 1.3454 1.3470 [3,] 1.3106 1.2936 1.2987 1.3134 1.3218 1.3352 1.3577 1.3615 1.3490 1.3349 [4,] 1.3084 1.3040 1.2991 1.3160 1.3183 1.3318 1.3601 1.3558 1.3448 1.3355 [5,] 1.3006 1.3005 1.3007 1.3230 1.3226 1.3366 1.3606 1.3535 1.3441 1.3345 [6,] 1.3018 1.2978 1.2956 1.3211 1.3325 1.3358 1.3557 1.3527 1.3453 1.3287 [7,] 1.2988 1.2901 1.3022 1.3225 1.3303 1.3352 1.3582 1.3486 1.3509 1.3304 [8,] 1.2984 1.2921 1.3082 1.3163 1.3296 1.3373 1.3649 1.3549 1.3420 1.3314 [9,] 1.2893 1.2972 1.3137 1.3083 1.3293 1.3426 1.3596 1.3538 1.3453 1.3404 [10,] 1.2941 1.2954 1.3119 1.3100 1.3351 1.3418 1.3643 1.3574 1.3436 1.3403 [11,] 1.2953 1.3020 1.3132 1.3135 1.3327 1.3467 1.3605 1.3516 1.3482 1.3427 [12,] 1.2909 1.3020 1.3145 1.3152 1.3265 1.3532 1.3588 1.3477 1.3532 1.3397 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [1,] 1.3441 1.3666 1.3722 1.3651 1.3631 1.3860 1.4165 1.4299 1.4488 1.4814 [2,] 1.3461 1.3753 1.3651 1.3591 1.3610 1.3877 1.4195 1.4288 1.4547 1.4829 [3,] 1.3460 1.3788 1.3659 1.3476 1.3705 1.3867 1.4109 1.4166 1.4722 1.4809 [4,] 1.3438 1.3782 1.3707 1.3405 1.3632 1.3975 1.4136 1.4254 1.4666 1.4845 [5,] 1.3467 1.3781 1.3663 1.3454 1.3580 1.4030 1.4089 1.4230 1.4683 1.4874 [6,] 1.3505 1.3771 1.3664 1.3476 1.3588 1.4049 1.4037 1.4309 1.4579 1.4747 [7,] 1.3588 1.3779 1.3694 1.3508 1.3669 1.4113 1.4146 1.4384 1.4607 1.4738 [8,] 1.3601 1.3820 1.3818 1.3493 1.3696 1.4106 1.4199 1.4391 1.4700 1.4761 [9,] 1.3618 1.3803 1.3794 1.3574 1.3795 1.4127 1.4173 1.4407 1.4639 1.4666 [10,] 1.3640 1.3821 1.3794 1.3615 1.3824 1.4180 1.4226 1.4447 1.4651 1.4741 [11,] 1.3596 1.3833 1.3729 1.3658 1.3885 1.4179 1.4150 1.4423 1.4654 1.4720 [12,] 1.3621 1.3743 1.3650 1.3664 1.3897 1.4232 1.4200 1.4479 1.4785 1.4554 [,21] [,22] [1,] 1.4649 1.4516 [2,] 1.4718 1.4692 [3,] 1.4672 1.4721 [4,] 1.4675 NA [5,] 1.4683 NA [6,] 1.4509 NA [7,] 1.4393 NA [8,] 1.4416 NA [9,] 1.4385 NA [10,] 1.4349 NA [11,] 1.4380 NA [12,] 1.4398 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -0.0039 0.0036 0.0030 -0.0039 0.0001 0.0001 -0.0001 -0.0052 0.0010 [2,] -0.0125 -0.0022 0.0032 0.0028 0.0062 0.0004 0.0028 0.0054 0.0036 [3,] -0.0022 0.0104 0.0004 0.0026 -0.0035 -0.0034 0.0024 -0.0057 -0.0042 [4,] -0.0078 -0.0035 0.0016 0.0070 0.0043 0.0048 0.0005 -0.0023 -0.0007 [5,] 0.0012 -0.0027 -0.0051 -0.0019 0.0099 -0.0008 -0.0049 -0.0008 0.0012 [6,] -0.0030 -0.0077 0.0066 0.0014 -0.0022 -0.0006 0.0025 -0.0041 0.0056 [7,] -0.0004 0.0020 0.0060 -0.0062 -0.0007 0.0021 0.0067 0.0063 -0.0089 [8,] -0.0091 0.0051 0.0055 -0.0080 -0.0003 0.0053 -0.0053 -0.0011 0.0033 [9,] 0.0048 -0.0018 -0.0018 0.0017 0.0058 -0.0008 0.0047 0.0036 -0.0017 [10,] 0.0012 0.0066 0.0013 0.0035 -0.0024 0.0049 -0.0038 -0.0058 0.0046 [11,] -0.0044 0.0000 0.0013 0.0017 -0.0062 0.0065 -0.0017 -0.0039 0.0050 [12,] 0.0013 -0.0095 0.0000 0.0003 0.0082 0.0018 0.0025 -0.0033 -0.0019 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] -0.0043 0.0020 0.0087 -0.0071 -0.0060 -0.0021 0.0017 0.0030 -0.0011 [2,] -0.0121 -0.0001 0.0035 0.0008 -0.0115 0.0095 -0.0010 -0.0086 -0.0122 [3,] 0.0006 -0.0022 -0.0006 0.0048 -0.0071 -0.0073 0.0108 0.0027 0.0088 [4,] -0.0010 0.0029 -0.0001 -0.0044 0.0049 -0.0052 0.0055 -0.0047 -0.0024 [5,] -0.0058 0.0038 -0.0010 0.0001 0.0022 0.0008 0.0019 -0.0052 0.0079 [6,] 0.0017 0.0083 0.0008 0.0030 0.0032 0.0081 0.0064 0.0109 0.0075 [7,] 0.0010 0.0013 0.0041 0.0124 -0.0015 0.0027 -0.0007 0.0053 0.0007 [8,] 0.0090 0.0017 -0.0017 -0.0024 0.0081 0.0099 0.0021 -0.0026 0.0016 [9,] -0.0001 0.0022 0.0018 0.0000 0.0041 0.0029 0.0053 0.0053 0.0040 [10,] 0.0024 -0.0044 0.0012 -0.0065 0.0043 0.0061 -0.0001 -0.0076 -0.0024 [11,] -0.0030 0.0025 -0.0090 -0.0079 0.0006 0.0012 0.0053 0.0050 0.0056 [12,] 0.0044 0.0045 -0.0021 0.0001 -0.0033 -0.0037 -0.0067 0.0099 0.0009 [,19] [,20] [,21] [,22] [1,] 0.0059 0.0015 0.0069 0.0176 [2,] 0.0175 -0.0020 -0.0046 0.0029 [3,] -0.0056 0.0036 0.0003 NA [4,] 0.0017 0.0029 0.0008 NA [5,] -0.0104 -0.0127 -0.0174 NA [6,] 0.0028 -0.0009 -0.0116 NA [7,] 0.0093 0.0023 0.0023 NA [8,] -0.0061 -0.0095 -0.0031 NA [9,] 0.0012 0.0075 -0.0036 NA [10,] 0.0003 -0.0021 0.0031 NA [11,] 0.0131 -0.0166 0.0018 NA [12,] 0.0029 0.0095 0.0118 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/freestat/rcomp/tmp/1oh681209291775.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/freestat/rcomp/tmp/2sc861209291775.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/freestat/rcomp/tmp/3udxx1209291775.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/freestat/rcomp/tmp/42fcr1209291775.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] [1,] 1.2922 1.2955 1.2936 1.2991 1.3005 1.2956 1.2901 1.2921 1.2893 1.2941 [2,] 1.3347 1.3348 1.3349 1.3318 1.3345 1.3325 1.3304 1.3314 1.3404 1.3403 [3,] 1.3622 1.3576 1.3596 1.3558 1.3535 1.3527 1.3582 1.3601 1.3596 1.3640 [4,] 1.4165 1.4195 1.4109 1.3975 1.4030 1.4037 1.4113 1.4106 1.4127 1.4180 [5,] 1.4814 1.4829 1.4809 1.4845 1.4874 1.4747 1.4738 1.4761 1.4666 1.4741 [,11] [,12] [1,] 1.2953 1.2909 [2,] 1.3427 1.3397 [3,] 1.3605 1.3621 [4,] 1.4150 1.4200 [5,] 1.4720 1.4785 $n [1] 22 22 22 21 21 21 21 21 21 21 21 21 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1.334645 1.329068 1.333999 1.333148 1.329882 1.328151 1.330307 1.332793 [2,] 1.389755 1.386132 1.385201 1.378452 1.377118 1.377249 1.386093 1.387407 [,9] [,10] [,11] [,12] [1,] 1.334672 1.337210 1.335572 1.334414 [2,] 1.384528 1.390790 1.385428 1.389786 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(1.2922, 1.3347, 1.3622, 1.4165, 1.4814, 1.2955, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/56tw71209291775.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] [1,] -0.00710 -0.0125 -0.0073 -0.0078 -0.0127 -0.0116 -0.0089 -0.0095 -0.0036 [2,] -0.00390 -0.0046 -0.0035 -0.0024 -0.0051 -0.0009 -0.0004 -0.0031 -0.0001 [3,] 0.00055 0.0006 0.0003 0.0005 -0.0008 0.0025 0.0021 -0.0003 0.0022 [4,] 0.00300 0.0035 0.0027 0.0029 0.0012 0.0064 0.0053 0.0051 0.0047 [5,] 0.00870 0.0095 0.0108 0.0070 0.0099 0.0109 0.0124 0.0099 0.0075 [,10] [,11] [,12] [1,] -0.0076 -0.0166 -0.0095 [2,] -0.0024 -0.0039 -0.0021 [3,] 0.0012 0.0012 0.0009 [4,] 0.0035 0.0050 0.0044 [5,] 0.0066 0.0131 0.0118 $n [1] 22 22 21 21 21 21 21 21 21 21 21 21 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.001774314 -0.002128543 -0.001837662 -0.001327357 -0.002972141 [2,] 0.002874314 0.003328543 0.002437662 0.002327357 0.001372141 [,6] [,7] [,8] [,9] [,10] [1,] -1.692515e-05 0.0001347297 -0.003127231 0.0005450355 -0.0008342272 [2,] 5.016925e-03 0.0040652703 0.002527231 0.0038549645 0.0032342272 [,11] [,12] [1,] -0.00186858 -0.001341098 [2,] 0.00426858 0.003141098 $out [1] 0.0176 0.0175 -0.0174 $group [1] 1 2 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.0071000000000001, -0.00390000000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6jjt51209291775.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] [1,] 1.2893 1.29010 1.29250 1.30830 1.31550 1.3318 1.35490 1.34770 1.34200 [2,] 1.2947 1.29290 1.29715 1.31200 1.32005 1.3350 1.35670 1.35215 1.34425 [3,] 1.2997 1.29650 1.30145 1.31485 1.32790 1.3362 1.35920 1.35435 1.34530 [4,] 1.3095 1.30125 1.31255 1.31870 1.33140 1.3422 1.36055 1.35675 1.34860 [5,] 1.3270 1.30400 1.31450 1.32300 1.33510 1.3467 1.36490 1.36150 1.35320 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] 1.32870 1.34380 1.37430 1.36500 1.3405 1.35800 1.38600 1.40370 1.41660 [2,] 1.33295 1.34605 1.37620 1.36610 1.3476 1.36205 1.39260 1.41225 1.42710 [3,] 1.33760 1.35465 1.37815 1.37005 1.3541 1.36825 1.40775 1.41575 1.43465 [4,] 1.34155 1.36095 1.38115 1.37615 1.3633 1.38095 1.41530 1.41970 1.44150 [5,] 1.35130 1.36400 1.38330 1.38180 1.3664 1.38970 1.42320 1.42260 1.44790 [,19] [,20] [,21] [,22] [1,] 1.44880 1.46660 1.43490 1.45160 [2,] 1.45930 1.47290 1.43890 1.46040 [3,] 1.46525 1.47540 1.44625 1.46920 [4,] 1.46915 1.48215 1.46735 1.47065 [5,] 1.47850 1.48740 1.47180 1.47210 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1.292950 1.292692 1.294426 1.311794 1.322723 1.332916 1.357444 1.352252 [2,] 1.306450 1.300308 1.308474 1.317906 1.333077 1.339484 1.360956 1.356448 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 1.343316 1.333677 1.347854 1.375892 1.365466 1.346939 1.359630 1.397396 [2,] 1.347284 1.341523 1.361446 1.380408 1.374634 1.361261 1.376870 1.418104 [,17] [,18] [,19] [,20] [,21] [,22] [1,] 1.412352 1.428082 1.460757 1.471181 1.433274 1.459850 [2,] 1.419148 1.441218 1.469743 1.479619 1.459226 1.478550 $out [1] 1.3532 1.3666 1.4554 $group [1] 6 12 20 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" NA Warning message: In bxp(list(stats = c(1.2893, 1.2947, 1.2997, 1.3095, 1.327, 1.2901, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7kqhq1209291775.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,] 1.366210 1.3527 1.364650 [2,] 1.368076 1.3567 1.369331 [3,] 1.371058 1.3596 1.371981 [4,] 1.372540 1.3613 1.377700 [5,] 1.372727 1.3640 1.379850 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1.369022 1.357502 1.368164 [2,] 1.373094 1.361698 1.375798 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1.36620952380952, 1.36807619047619, 1.37105822510823, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1oh681209291775.ps tmp/1oh681209291775.png") > system("convert tmp/2sc861209291775.ps tmp/2sc861209291775.png") > system("convert tmp/3udxx1209291775.ps tmp/3udxx1209291775.png") > system("convert tmp/42fcr1209291775.ps tmp/42fcr1209291775.png") > system("convert tmp/56tw71209291775.ps tmp/56tw71209291775.png") > system("convert tmp/6jjt51209291775.ps tmp/6jjt51209291775.png") > system("convert tmp/7kqhq1209291775.ps tmp/7kqhq1209291775.png") > > > proc.time() user system elapsed 2.172 1.391 2.370