R version 2.8.0 (2008-10-20) 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. 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(1772.2 + ,1769.5 + ,1768 + ,1794.8 + ,1823.4 + ,1856.9 + ,1866.9 + ,1869.8 + ,1843.8 + ,1837.1 + ,1857.7 + ,1840.3 + ,1914.6 + ,1972.9 + ,2050.1 + ,2086.2 + ,2112.5 + ,2147.6 + ,2190.4 + ,2194.1 + ,2216.2 + ,2218.6 + ,2233.5 + ,2307.2 + ,2350.4 + ,2368.2 + ,2353.8 + ,2316.5 + ,2305.5 + ,2308.4 + ,2334.4 + ,2381.2 + ,2449.7 + ,2490.3 + ,2523.5 + ,2537.6 + ,2526.1 + ,2545.9 + ,2542.7 + ,2584.3 + ,2600.2 + ,2593.9 + ,2618.9 + ,2591.3 + ,2521.2 + ,2536.6 + ,2596.1 + ,2656.6 + ,2710.3 + ,2778.8 + ,2775.5 + ,2785.2 + ,2847.7 + ,2834.4 + ,2839 + ,2802.6 + ,2819.3 + ,2872 + ,2918.4 + ,2977.8 + ,3031.2 + ,3064.7 + ,3093 + ,3100.6 + ,3141.1 + ,3180.4 + ,3240.3 + ,3265 + ,3338.2 + ,3376.6 + ,3422.5 + ,3432 + ,3516.3 + ,3564 + ,3636.3 + ,3724 + ,3815.4 + ,3828.1 + ,3853.3 + ,3884.5 + ,3918.7 + ,3919.6 + ,3950.8 + ,3981 + ,4063 + ,4132 + ,4160.3 + ,4178.3 + ,4244.1 + ,4256.5 + ,4283.4 + ,4263.3 + ,4256.6 + ,4264.3 + ,4302.3 + ,4256.6 + ,4374 + ,4398.8 + ,4433.9 + ,4446.3 + ,4525.8 + ,4633.1 + ,4677.5 + ,4754.5 + ,4876.2 + ,4932.6 + ,4906.3 + ,4953.1 + ,4909.6 + ,4922.2 + ,4873.5 + ,4854.3 + ,4795.3 + ,4831.9 + ,4913.3 + ,4977.5 + ,5090.7 + ,5128.9 + ,5154.1 + ,5191.5 + ,5251.8 + ,5356.1 + ,5451.9 + ,5450.8 + ,5469.4 + ,5684.6 + ,5740.3 + ,5816.2 + ,5825.9 + ,5831.4 + ,5873.3 + ,5889.5 + ,5908.5 + ,5787.4 + ,5776.6 + ,5883.5 + ,6005.7 + ,5957.8 + ,6030.2 + ,5955.1 + ,5857.3 + ,5889.1 + ,5866.4 + ,5871 + ,5944 + ,6077.6 + ,6197.5 + ,6325.6 + ,6448.3 + ,6559.6 + ,6623.3 + ,6677.3 + ,6740.3 + ,6797.3 + ,6903.5 + ,6955.9 + ,7022.8 + ,7051 + ,7119 + ,7153.4 + ,7193 + ,7269.5 + ,7332.6 + ,7458 + ,7496.6 + ,7592.9 + ,7632.1 + ,7734 + ,7806.6 + ,7865 + ,7927.4 + ,7944.7 + ,8027.7 + ,8059.6 + ,8059.5 + ,7988.9 + ,7950.2 + ,8003.8 + ,8037.5 + ,8069 + ,8157.6 + ,8244.3 + ,8329.4 + ,8417 + ,8432.5 + ,8486.4 + ,8531.1 + ,8643.8 + ,8727.9 + ,8847.3 + ,8904.3 + ,9003.2 + ,9025.3 + ,9044.7 + ,9120.7 + ,9184.3 + ,9247.2 + ,9407.1 + ,9488.9 + ,9592.5 + ,9666.2 + ,9809.6 + ,9932.7 + ,10008.9 + ,10103.4 + ,10194.3 + ,10328.8 + ,10507.6 + ,10601.2 + ,10684 + ,10819.9 + ,11014.3 + ,11043 + ,11258.5 + ,11267.9 + ,11334.5 + ,11297.2 + ,11371.3 + ,11340.1 + ,11380.1 + ,11477.9 + ,11538.8 + ,11596.4 + ,11598.8 + ,11645.8 + ,11738.7 + ,11935.5 + ,12042.8 + ,12127.6 + ,12213.8 + ,12303.5 + ,12410.3 + ,12534.1 + ,12587.5 + ,12683.2 + ,12748.7 + ,12915.9 + ,12962.5 + ,12965.9 + ,13060.7 + ,13099.9 + ,13204 + ,13321.1 + ,13391.2 + ,13366.9 + ,13415.3 + ,13324.6 + ,13141.9 + ,12925.4 + ,12901.5 + ,12973 + ,13155) > par1 = '4' > #'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] 252 > (np <- floor(n / par1)) [1] 63 > 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] 63 63 63 63 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 1772.2 1823.4 1843.8 1914.6 2112.5 2216.2 2350.4 2305.5 2449.7 2526.1 [2,] 1769.5 1856.9 1837.1 1972.9 2147.6 2218.6 2368.2 2308.4 2490.3 2545.9 [3,] 1768.0 1866.9 1857.7 2050.1 2190.4 2233.5 2353.8 2334.4 2523.5 2542.7 [4,] 1794.8 1869.8 1840.3 2086.2 2194.1 2307.2 2316.5 2381.2 2537.6 2584.3 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [1,] 2600.2 2521.2 2710.3 2847.7 2819.3 3031.2 3141.1 3338.2 3516.3 3815.4 [2,] 2593.9 2536.6 2778.8 2834.4 2872.0 3064.7 3180.4 3376.6 3564.0 3828.1 [3,] 2618.9 2596.1 2775.5 2839.0 2918.4 3093.0 3240.3 3422.5 3636.3 3853.3 [4,] 2591.3 2656.6 2785.2 2802.6 2977.8 3100.6 3265.0 3432.0 3724.0 3884.5 [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [1,] 3918.7 4063.0 4244.1 4256.6 4374.0 4525.8 4876.2 4909.6 4795.3 5090.7 [2,] 3919.6 4132.0 4256.5 4264.3 4398.8 4633.1 4932.6 4922.2 4831.9 5128.9 [3,] 3950.8 4160.3 4283.4 4302.3 4433.9 4677.5 4906.3 4873.5 4913.3 5154.1 [4,] 3981.0 4178.3 4263.3 4256.6 4446.3 4754.5 4953.1 4854.3 4977.5 5191.5 [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 5251.8 5469.4 5825.9 5908.5 6005.7 5857.3 5944.0 6448.3 6740.3 7022.8 [2,] 5356.1 5684.6 5831.4 5787.4 5957.8 5889.1 6077.6 6559.6 6797.3 7051.0 [3,] 5451.9 5740.3 5873.3 5776.6 6030.2 5866.4 6197.5 6623.3 6903.5 7119.0 [4,] 5450.8 5816.2 5889.5 5883.5 5955.1 5871.0 6325.6 6677.3 6955.9 7153.4 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 7193.0 7496.6 7806.6 8027.7 7950.2 8157.6 8432.5 8727.9 9025.3 9247.2 [2,] 7269.5 7592.9 7865.0 8059.6 8003.8 8244.3 8486.4 8847.3 9044.7 9407.1 [3,] 7332.6 7632.1 7927.4 8059.5 8037.5 8329.4 8531.1 8904.3 9120.7 9488.9 [4,] 7458.0 7734.0 7944.7 7988.9 8069.0 8417.0 8643.8 9003.2 9184.3 9592.5 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [1,] 9666.2 10103.4 10601.2 11043.0 11297.2 11477.9 11645.8 12127.6 12534.1 [2,] 9809.6 10194.3 10684.0 11258.5 11371.3 11538.8 11738.7 12213.8 12587.5 [3,] 9932.7 10328.8 10819.9 11267.9 11340.1 11596.4 11935.5 12303.5 12683.2 [4,] 10008.9 10507.6 11014.3 11334.5 11380.1 11598.8 12042.8 12410.3 12748.7 [,60] [,61] [,62] [,63] [,64] [1,] 12915.9 13099.9 13366.9 12925.4 NA [2,] 12962.5 13204.0 13415.3 12901.5 NA [3,] 12965.9 13321.1 13324.6 12973.0 NA [4,] 13060.7 13391.2 13141.9 13155.0 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] -2.7 33.5 -6.7 58.3 35.1 2.4 17.8 2.9 40.6 19.8 -6.3 15.4 68.5 [2,] -1.5 10.0 20.6 77.2 42.8 14.9 -14.4 26.0 33.2 -3.2 25.0 59.5 -3.3 [3,] 26.8 2.9 -17.4 36.1 3.7 73.7 -37.3 46.8 14.1 41.6 -27.6 60.5 9.7 [4,] 28.6 -26.0 74.3 26.3 22.1 43.2 -11.0 68.5 -11.5 15.9 -70.1 53.7 62.5 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] -13.3 52.7 33.5 39.3 38.4 47.7 12.7 0.9 69.0 12.4 7.7 24.8 [2,] 4.6 46.4 28.3 59.9 45.9 72.3 25.2 31.2 28.3 26.9 38.0 35.1 [3,] -36.4 59.4 7.6 24.7 9.5 87.7 31.2 30.2 18.0 -20.1 -45.7 12.4 [4,] 16.7 53.4 40.5 73.2 84.3 91.4 34.2 82.0 65.8 -6.7 117.4 79.5 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 107.3 56.4 12.6 36.6 38.2 104.3 215.2 5.5 -121.1 -47.9 31.8 133.6 [2,] 44.4 -26.3 -48.7 81.4 25.2 95.8 55.7 41.9 -10.8 72.4 -22.7 119.9 [3,] 77.0 46.8 -19.2 64.2 37.4 -1.1 75.9 16.2 106.9 -75.1 4.6 128.1 [4,] 121.7 -43.5 -59.0 113.2 60.3 18.6 9.7 19.0 122.2 -97.8 73.0 122.7 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 111.3 57.0 28.2 76.5 96.3 58.4 31.9 53.6 86.7 53.9 119.4 19.4 [2,] 63.7 106.2 68.0 63.1 39.2 62.4 -0.1 33.7 85.1 44.7 57.0 76.0 [3,] 54.0 52.4 34.4 125.4 101.9 17.3 -70.6 31.5 87.6 112.7 98.9 63.6 [4,] 63.0 66.9 39.6 38.6 72.6 83.0 -38.7 88.6 15.5 84.1 22.1 62.9 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 159.9 143.4 90.9 82.8 215.5 74.1 60.9 92.9 86.2 53.4 46.6 104.1 [2,] 81.8 123.1 134.5 135.9 9.4 -31.2 57.6 196.8 89.7 95.7 3.4 117.1 [3,] 103.6 76.2 178.8 194.4 66.6 40.0 2.4 107.3 106.8 65.5 94.8 70.1 [4,] 73.7 94.5 93.6 28.7 -37.3 97.8 47.0 84.8 123.8 167.2 39.2 -24.3 [,62] [,63] [,64] [1,] 48.4 -23.9 NA [2,] -90.7 71.5 NA [3,] -182.7 182.0 NA [4,] -216.5 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/freestat/rcomp/tmp/1j24c1272277047.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/2tclf1272277047.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/3tclf1272277047.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/443li1272277047.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] [1,] 1772.20 1769.50 1768.00 1794.8 [2,] 3086.15 3122.55 3166.65 3182.8 [3,] 5469.40 5684.60 5740.30 5816.2 [4,] 8580.20 8666.85 8717.70 8823.5 [5,] 13366.90 13415.30 13324.60 13391.2 $n [1] 63 63 63 63 $conf [,1] [,2] [,3] [,4] [1,] 4375.747 4580.944 4635.301 4693.355 [2,] 6563.053 6788.256 6845.299 6939.045 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/543li1272277047.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] [1,] -47.90 -48.70 -75.10 -70.10 [2,] 16.60 17.75 8.55 16.70 [3,] 46.60 42.80 40.00 53.55 [4,] 79.65 72.35 76.60 83.00 [5,] 159.90 135.90 128.10 167.20 $n [1] 63 63 63 62 $conf [,1] [,2] [,3] [,4] [1,] 34.04919 31.93125 26.45388 40.24623 [2,] 59.15081 53.66875 53.54612 66.85377 $out [1] 215.2 -121.1 215.5 196.8 -90.7 178.8 194.4 -182.7 182.0 -97.8 [11] -216.5 $group [1] 1 1 1 2 2 3 3 3 3 4 4 $names [1] "1" "2" "3" "4" > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6xu231272277047.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,] 1768.00 1823.40 1837.10 1914.60 2112.50 2216.20 2316.50 2305.50 2449.70 [2,] 1768.75 1840.15 1838.70 1943.75 2130.05 2217.40 2333.45 2306.95 2470.00 [3,] 1770.85 1861.90 1842.05 2011.50 2169.00 2226.05 2352.10 2321.40 2506.90 [4,] 1783.50 1868.35 1850.75 2068.15 2192.25 2270.35 2361.00 2357.80 2530.55 [5,] 1794.80 1869.80 1857.70 2086.20 2194.10 2307.20 2368.20 2381.20 2537.60 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] 2526.1 2591.30 2521.20 2710.30 2802.60 2819.30 3031.20 3141.10 3338.20 [2,] 2534.4 2592.60 2528.90 2742.90 2818.50 2845.65 3047.95 3160.75 3357.40 [3,] 2544.3 2597.05 2566.35 2777.15 2836.70 2895.20 3078.85 3210.35 3399.55 [4,] 2565.1 2609.55 2626.35 2782.00 2843.35 2948.10 3096.80 3252.65 3427.25 [5,] 2584.3 2618.90 2656.60 2785.20 2847.70 2977.80 3100.60 3265.00 3432.00 [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27] [1,] 3516.30 3815.40 3918.70 4063.00 4244.10 4256.60 4374.00 4525.80 4876.20 [2,] 3540.15 3821.75 3919.15 4097.50 4250.30 4256.60 4386.40 4579.45 4891.25 [3,] 3600.15 3840.70 3935.20 4146.15 4259.90 4260.45 4416.35 4655.30 4919.45 [4,] 3680.15 3868.90 3965.90 4169.30 4273.35 4283.30 4440.10 4716.00 4942.85 [5,] 3724.00 3884.50 3981.00 4178.30 4283.40 4302.30 4446.30 4754.50 4953.10 [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 4854.30 4795.3 5090.7 5251.80 5469.40 5825.90 5776.60 5955.10 5857.30 [2,] 4863.90 4813.6 5109.8 5303.95 5577.00 5828.65 5782.00 5956.45 5861.85 [3,] 4891.55 4872.6 5141.5 5403.45 5712.45 5852.35 5835.45 5981.75 5868.70 [4,] 4915.90 4945.4 5172.8 5451.35 5778.25 5881.40 5896.00 6017.95 5880.05 [5,] 4922.20 4977.5 5191.5 5451.90 5816.20 5889.50 5908.50 6030.20 5889.10 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [1,] 5944.00 6448.30 6740.3 7022.8 7193.00 7496.60 7806.60 7988.90 7950.20 [2,] 6010.80 6503.95 6768.8 7036.9 7231.25 7544.75 7835.80 8008.30 7977.00 [3,] 6137.55 6591.45 6850.4 7085.0 7301.05 7612.50 7896.20 8043.60 8020.65 [4,] 6261.55 6650.30 6929.7 7136.2 7395.30 7683.05 7936.05 8059.55 8053.25 [5,] 6325.60 6677.30 6955.9 7153.4 7458.00 7734.00 7944.70 8059.60 8069.00 [,46] [,47] [,48] [,49] [,50] [,51] [,52] [,53] [,54] [1,] 8157.60 8432.50 8727.90 9025.3 9247.20 9666.20 10103.40 10601.20 11043.00 [2,] 8200.95 8459.45 8787.60 9035.0 9327.15 9737.90 10148.85 10642.60 11150.75 [3,] 8286.85 8508.75 8875.80 9082.7 9448.00 9871.15 10261.55 10751.95 11263.20 [4,] 8373.20 8587.45 8953.75 9152.5 9540.70 9970.80 10418.20 10917.10 11301.20 [5,] 8417.00 8643.80 9003.20 9184.3 9592.50 10008.90 10507.60 11014.30 11334.50 [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 11297.20 11477.90 11645.80 12127.60 12534.10 12915.9 13099.90 13141.90 [2,] 11318.65 11508.35 11692.25 12170.70 12560.80 12939.2 13151.95 13233.25 [3,] 11355.70 11567.60 11837.10 12258.65 12635.35 12964.2 13262.55 13345.75 [4,] 11375.70 11597.60 11989.15 12356.90 12715.95 13013.3 13356.15 13391.10 [5,] 11380.10 11598.80 12042.80 12410.30 12748.70 13060.7 13391.20 13415.30 [,63] [,64] [1,] 12901.50 NA [2,] 12913.45 NA [3,] 12949.20 NA [4,] 13064.00 NA [5,] 13155.00 NA $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 [39] 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1759.197 1839.622 1832.530 1913.224 2119.862 2184.220 2330.336 2281.228 [2,] 1782.502 1884.178 1851.570 2109.776 2218.138 2267.881 2373.865 2361.572 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 2459.066 2520.047 2583.660 2489.364 2746.261 2817.068 2814.264 3040.258 [2,] 2554.735 2568.553 2610.441 2643.336 2808.039 2856.331 2976.135 3117.442 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 3137.749 3344.369 3489.55 3803.451 3898.267 4089.428 4241.690 4239.357 [2,] 3282.951 3454.732 3710.75 3877.948 3972.133 4202.872 4278.109 4281.543 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 4373.927 4547.426 4878.686 4850.47 4768.478 5091.73 5287.004 5553.463 [2,] 4458.773 4763.174 4960.214 4932.63 4976.722 5191.27 5519.896 5871.438 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 5810.678 5745.39 5933.165 5854.322 5939.458 6475.834 6723.289 7006.553 [2,] 5894.023 5925.51 6030.335 5883.078 6335.642 6707.066 6977.511 7163.447 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 7171.45 7503.243 7817.003 8003.112 7960.412 8150.772 8407.63 8744.541 [2,] 7430.65 7721.757 7975.397 8084.088 8080.887 8422.927 8609.87 9007.059 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [1,] 8989.875 9279.296 9687.159 10048.76 10535.10 11144.34 11310.63 11497.09 [2,] 9175.525 9616.704 10055.141 10474.34 10968.80 11382.06 11400.77 11638.11 [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [1,] 11602.55 12111.55 12512.78 12905.66 13101.23 13221.05 12830.27 NA [2,] 12071.65 12405.75 12757.92 13022.74 13423.87 13470.45 13068.13 NA $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" "53" "54" "55" "56" "57" "58" "59" "60" [61] "61" "62" "63" NA Warning message: In bxp(list(stats = c(1768, 1768.75, 1770.85, 1783.5, 1794.8, 1823.4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7xu231272277047.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,] 6223.054 5469.40 5833.175 [2,] 6248.488 5577.00 5863.938 [3,] 6296.547 5712.45 5918.438 [4,] 6340.510 5778.25 5972.663 [5,] 6361.849 5816.20 6003.150 $n [1] 4 4 4 $conf [,1] [,2] [,3] [1,] 6223.849 5553.463 5832.545 [2,] 6369.244 5871.438 6004.330 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(6223.05396825397, 6248.4880952381, 6296.54682539683, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1j24c1272277047.ps tmp/1j24c1272277047.png",intern=TRUE)) character(0) > try(system("convert tmp/2tclf1272277047.ps tmp/2tclf1272277047.png",intern=TRUE)) character(0) > try(system("convert tmp/3tclf1272277047.ps tmp/3tclf1272277047.png",intern=TRUE)) character(0) > try(system("convert tmp/443li1272277047.ps tmp/443li1272277047.png",intern=TRUE)) character(0) > try(system("convert tmp/543li1272277047.ps tmp/543li1272277047.png",intern=TRUE)) character(0) > try(system("convert tmp/6xu231272277047.ps tmp/6xu231272277047.png",intern=TRUE)) character(0) > try(system("convert tmp/7xu231272277047.ps tmp/7xu231272277047.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.937 1.477 2.277