R version 2.6.1 (2007-11-26) Copyright (C) 2007 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(87.28 + ,87.28 + ,87.09 + ,86.92 + ,87.59 + ,90.72 + ,90.69 + ,90.3 + ,89.55 + ,88.94 + ,88.41 + ,87.82 + ,87.07 + ,86.82 + ,86.4 + ,86.02 + ,85.66 + ,85.32 + ,85 + ,84.67 + ,83.94 + ,82.83 + ,81.95 + ,81.19 + ,80.48 + ,78.86 + ,69.47 + ,68.77 + ,70.06 + ,73.95 + ,75.8 + ,77.79 + ,81.57 + ,83.07 + ,84.34 + ,85.1 + ,85.25 + ,84.26 + ,83.63 + ,86.44 + ,85.3 + ,84.1 + ,83.36 + ,82.48 + ,81.58 + ,80.47 + ,79.34 + ,82.13 + ,81.69 + ,80.7 + ,79.88 + ,79.16 + ,78.38 + ,77.42 + ,76.47 + ,75.46 + ,74.48 + ,78.27 + ,80.7 + ,79.91 + ,78.75 + ,77.78 + ,81.14 + ,81.08 + ,80.03 + ,78.91 + ,78.01 + ,76.9 + ,75.97 + ,81.93 + ,80.27 + ,78.67 + ,77.42 + ,76.16 + ,74.7 + ,76.39 + ,76.04 + ,74.65 + ,73.29 + ,71.79 + ,74.39 + ,74.91 + ,74.54 + ,73.08 + ,72.75 + ,71.32 + ,70.38 + ,70.35 + ,70.01 + ,69.36 + ,67.77 + ,69.26 + ,69.8 + ,68.38 + ,67.62 + ,68.39 + ,66.95 + ,65.21 + ,66.64 + ,63.45 + ,60.66 + ,62.34 + ,60.32 + ,58.64 + ,60.46 + ,58.59 + ,61.87 + ,61.85 + ,67.44 + ,77.06 + ,91.74 + ,93.15 + ,94.15 + ,93.11 + ,91.51 + ,89.96 + ,88.16 + ,86.98 + ,88.03 + ,86.24 + ,84.65 + ,83.23 + ,81.7 + ,80.25 + ,78.8 + ,77.51 + ,76.2 + ,75.04 + ,74 + ,75.49 + ,77.14 + ,76.15 + ,76.27 + ,78.19 + ,76.49 + ,77.31 + ,76.65 + ,74.99 + ,73.51 + ,72.07 + ,70.59 + ,71.96 + ,76.29 + ,74.86 + ,74.93 + ,71.9 + ,71.01 + ,77.47 + ,75.78 + ,76.6 + ,76.07 + ,74.57 + ,73.02 + ,72.65 + ,73.16 + ,71.53 + ,69.78 + ,67.98 + ,69.96 + ,72.16 + ,70.47 + ,68.86 + ,67.37 + ,65.87 + ,72.16 + ,71.34 + ,69.93 + ,68.44 + ,67.16 + ,66.01 + ,67.25 + ,70.91 + ,69.75 + ,68.59 + ,67.48 + ,66.31 + ,64.81 + ,66.58 + ,65.97 + ,64.7 + ,64.7 + ,60.94 + ,59.08 + ,58.42 + ,57.77 + ,57.11 + ,53.31 + ,49.96 + ,49.4 + ,48.84 + ,48.3 + ,47.74 + ,47.24 + ,46.76 + ,46.29 + ,48.9 + ,49.23 + ,48.53 + ,48.03 + ,54.34 + ,53.79 + ,53.24 + ,52.96 + ,52.17 + ,51.7 + ,58.55 + ,78.2 + ,77.03 + ,76.19 + ,77.15 + ,75.87 + ,95.47 + ,109.67 + ,112.28 + ,112.01 + ,107.93 + ,105.96 + ,105.06 + ,102.98 + ,102.2 + ,105.23 + ,101.85 + ,99.89 + ,96.23 + ,94.76 + ,91.51 + ,91.63 + ,91.54 + ,85.23 + ,87.83 + ,87.38 + ,84.44 + ,85.19 + ,84.03 + ,86.73 + ,102.52 + ,104.45 + ,106.98 + ,107.02 + ,99.26 + ,94.45 + ,113.44 + ,157.33 + ,147.38 + ,171.89 + ,171.95 + ,132.71 + ,126.02 + ,121.18 + ,115.45 + ,110.48 + ,117.85 + ,117.63 + ,124.65 + ,109.59 + ,111.27 + ,99.78 + ,98.21 + ,99.2 + ,97.97 + ,89.55 + ,87.91 + ,93.34 + ,94.42 + ,93.2 + ,90.29 + ,91.46 + ,89.98 + ,88.35 + ,88.41 + ,82.44 + ,79.89 + ,75.69 + ,75.66 + ,84.5 + ,96.73 + ,87.48 + ,82.39 + ,83.48 + ,79.31 + ,78.16 + ,72.77 + ,72.45 + ,68.46 + ,67.62 + ,68.76 + ,70.07 + ,68.55 + ,65.3 + ,58.96 + ,59.17 + ,62.37 + ,66.28 + ,55.62 + ,55.23 + ,55.85 + ,56.75 + ,50.89 + ,53.88 + ,52.95 + ,55.08 + ,53.61 + ,58.78 + ,61.85 + ,55.91 + ,53.32 + ,46.41 + ,44.57 + ,50 + ,50 + ,53.36 + ,46.23 + ,50.45 + ,49.07 + ,45.85 + ,48.45 + ,49.96 + ,46.53 + ,50.51 + ,47.58 + ,48.05 + ,46.84 + ,47.67 + ,49.16 + ,55.54 + ,55.82 + ,58.22 + ,56.19 + ,57.77 + ,63.19 + ,54.76 + ,55.74 + ,62.54 + ,61.39 + ,69.6 + ,79.23 + ,80 + ,93.68 + ,107.63 + ,100.18 + ,97.3 + ,90.45 + ,80.64 + ,80.58 + ,75.82 + ,85.59 + ,89.35 + ,89.42 + ,104.73 + ,95.32 + ,89.27 + ,90.44 + ,86.97 + ,79.98 + ,81.22 + ,87.35 + ,83.64 + ,82.22 + ,94.4 + ,102.18) > 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] 360 > (np <- floor(n / par1)) [1] 30 > 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] 30 30 30 30 30 30 30 30 30 30 30 30 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 87.28 87.07 80.48 85.25 81.69 78.75 77.42 72.75 66.95 67.44 84.65 76.27 [2,] 87.28 86.82 78.86 84.26 80.70 77.78 76.16 71.32 65.21 77.06 83.23 78.19 [3,] 87.09 86.40 69.47 83.63 79.88 81.14 74.70 70.38 66.64 91.74 81.70 76.49 [4,] 86.92 86.02 68.77 86.44 79.16 81.08 76.39 70.35 63.45 93.15 80.25 77.31 [5,] 87.59 85.66 70.06 85.30 78.38 80.03 76.04 70.01 60.66 94.15 78.80 76.65 [6,] 90.72 85.32 73.95 84.10 77.42 78.91 74.65 69.36 62.34 93.11 77.51 74.99 [7,] 90.69 85.00 75.80 83.36 76.47 78.01 73.29 67.77 60.32 91.51 76.20 73.51 [8,] 90.30 84.67 77.79 82.48 75.46 76.90 71.79 69.26 58.64 89.96 75.04 72.07 [9,] 89.55 83.94 81.57 81.58 74.48 75.97 74.39 69.80 60.46 88.16 74.00 70.59 [10,] 88.94 82.83 83.07 80.47 78.27 81.93 74.91 68.38 58.59 86.98 75.49 71.96 [11,] 88.41 81.95 84.34 79.34 80.70 80.27 74.54 67.62 61.87 88.03 77.14 76.29 [12,] 87.82 81.19 85.10 82.13 79.91 78.67 73.08 68.39 61.85 86.24 76.15 74.86 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] 74.93 69.78 67.16 64.70 47.24 51.70 105.96 85.23 94.45 117.63 93.20 [2,] 71.90 67.98 66.01 60.94 46.76 58.55 105.06 87.83 113.44 124.65 90.29 [3,] 71.01 69.96 67.25 59.08 46.29 78.20 102.98 87.38 157.33 109.59 91.46 [4,] 77.47 72.16 70.91 58.42 48.90 77.03 102.20 84.44 147.38 111.27 89.98 [5,] 75.78 70.47 69.75 57.77 49.23 76.19 105.23 85.19 171.89 99.78 88.35 [6,] 76.60 68.86 68.59 57.11 48.53 77.15 101.85 84.03 171.95 98.21 88.41 [7,] 76.07 67.37 67.48 53.31 48.03 75.87 99.89 86.73 132.71 99.20 82.44 [8,] 74.57 65.87 66.31 49.96 54.34 95.47 96.23 102.52 126.02 97.97 79.89 [9,] 73.02 72.16 64.81 49.40 53.79 109.67 94.76 104.45 121.18 89.55 75.69 [10,] 72.65 71.34 66.58 48.84 53.24 112.28 91.51 106.98 115.45 87.91 75.66 [11,] 73.16 69.93 65.97 48.30 52.96 112.01 91.63 107.02 110.48 93.34 84.50 [12,] 71.53 68.44 64.70 47.74 52.17 107.93 91.54 99.26 117.85 94.42 96.73 [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [1,] 87.48 65.30 55.08 50.45 55.54 80.00 104.73 NA [2,] 82.39 58.96 53.61 49.07 55.82 93.68 95.32 NA [3,] 83.48 59.17 58.78 45.85 58.22 107.63 89.27 NA [4,] 79.31 62.37 61.85 48.45 56.19 100.18 90.44 NA [5,] 78.16 66.28 55.91 49.96 57.77 97.30 86.97 NA [6,] 72.77 55.62 53.32 46.53 63.19 90.45 79.98 NA [7,] 72.45 55.23 46.41 50.51 54.76 80.64 81.22 NA [8,] 68.46 55.85 44.57 47.58 55.74 80.58 87.35 NA [9,] 67.62 56.75 50.00 48.05 62.54 75.82 83.64 NA [10,] 68.76 50.89 50.00 46.84 61.39 85.59 82.22 NA [11,] 70.07 53.88 53.36 47.67 69.60 89.35 94.40 NA [12,] 68.55 52.95 46.23 49.16 79.23 89.42 102.18 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.00 -0.25 -1.62 -0.99 -0.99 -0.97 -1.26 -1.43 -1.74 9.62 -1.42 1.92 [2,] -0.19 -0.42 -9.39 -0.63 -0.82 3.36 -1.46 -0.94 1.43 14.68 -1.53 -1.70 [3,] -0.17 -0.38 -0.70 2.81 -0.72 -0.06 1.69 -0.03 -3.19 1.41 -1.45 0.82 [4,] 0.67 -0.36 1.29 -1.14 -0.78 -1.05 -0.35 -0.34 -2.79 1.00 -1.45 -0.66 [5,] 3.13 -0.34 3.89 -1.20 -0.96 -1.12 -1.39 -0.65 1.68 -1.04 -1.29 -1.66 [6,] -0.03 -0.32 1.85 -0.74 -0.95 -0.90 -1.36 -1.59 -2.02 -1.60 -1.31 -1.48 [7,] -0.39 -0.33 1.99 -0.88 -1.01 -1.11 -1.50 1.49 -1.68 -1.55 -1.16 -1.44 [8,] -0.75 -0.73 3.78 -0.90 -0.98 -0.93 2.60 0.54 1.82 -1.80 -1.04 -1.48 [9,] -0.61 -1.11 1.50 -1.11 3.79 5.96 0.52 -1.42 -1.87 -1.18 1.49 1.37 [10,] -0.53 -0.88 1.27 -1.13 2.43 -1.66 -0.37 -0.76 3.28 1.05 1.65 4.33 [11,] -0.59 -0.76 0.76 2.79 -0.79 -1.60 -1.46 0.77 -0.02 -1.79 -0.99 -1.43 [12,] -0.75 -0.71 0.15 -0.44 -1.16 -1.25 -0.33 -1.44 5.59 -1.59 0.12 0.07 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] -3.03 -1.80 -1.15 -3.76 -0.48 6.85 -0.90 2.60 18.99 7.02 -2.91 -5.09 [2,] -0.89 1.98 1.24 -1.86 -0.47 19.65 -2.08 -0.45 43.89 -15.06 1.17 1.09 [3,] 6.46 2.20 3.66 -0.66 2.61 -1.17 -0.78 -2.94 -9.95 1.68 -1.48 -4.17 [4,] -1.69 -1.69 -1.16 -0.65 0.33 -0.84 3.03 0.75 24.51 -11.49 -1.63 -1.15 [5,] 0.82 -1.61 -1.16 -0.66 -0.70 0.96 -3.38 -1.16 0.06 -1.57 0.06 -5.39 [6,] -0.53 -1.49 -1.11 -3.80 -0.50 -1.28 -1.96 2.70 -39.24 0.99 -5.97 -0.32 [7,] -1.50 -1.50 -1.17 -3.35 6.31 19.60 -3.66 15.79 -6.69 -1.23 -2.55 -3.99 [8,] -1.55 6.29 -1.50 -0.56 -0.55 14.20 -1.47 1.93 -4.84 -8.42 -4.20 -0.84 [9,] -0.37 -0.82 1.77 -0.56 -0.55 2.61 -3.25 2.53 -5.73 -1.64 -0.03 1.14 [10,] 0.51 -1.41 -0.61 -0.54 -0.28 -0.27 0.12 0.04 -4.97 5.43 8.84 1.31 [11,] -1.63 -1.49 -1.27 -0.56 -0.79 -4.08 -0.09 -7.76 7.37 1.08 12.23 -1.52 [12,] -1.75 -1.28 0.00 -0.50 -0.47 -1.97 -6.31 -4.81 -0.22 -1.22 -9.25 -3.25 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [1,] -6.34 -1.47 -1.38 0.28 13.68 -9.41 NA [2,] 0.21 5.17 -3.22 2.40 13.95 -6.05 NA [3,] 3.20 3.07 2.60 -2.03 -7.45 1.17 NA [4,] 3.91 -5.94 1.51 1.58 -2.88 -3.47 NA [5,] -10.66 -2.59 -3.43 5.42 -6.85 -6.99 NA [6,] -0.39 -6.91 3.98 -8.43 -9.81 1.24 NA [7,] 0.62 -1.84 -2.93 0.98 -0.06 6.13 NA [8,] 0.90 5.43 0.47 6.80 -4.76 -3.71 NA [9,] -5.86 0.00 -1.21 -1.15 9.77 -1.42 NA [10,] 2.99 3.36 0.83 8.21 3.76 12.18 NA [11,] -0.93 -7.13 1.49 9.63 0.07 7.78 NA [12,] 2.13 4.22 6.38 0.77 15.31 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/rcomp/tmp/18nsk1199713100.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/rcomp/tmp/2eim51199713100.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/rcomp/tmp/31z0n1199713100.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/rcomp/tmp/4mc5h1199713100.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] [1,] 47.240 46.760 45.85 48.450 49.230 46.530 46.41 44.57 48.050 [2,] 66.950 65.210 67.25 68.770 69.750 68.590 67.37 65.87 64.810 [3,] 78.085 77.985 79.04 78.315 77.405 76.875 75.97 75.25 74.435 [4,] 87.070 87.280 87.38 86.920 86.970 85.320 83.36 87.35 83.940 [5,] 105.960 113.440 109.59 111.270 105.230 101.850 99.89 102.52 109.670 [,10] [,11] [,12] [1,] 46.840 47.67 46.23 [2,] 66.580 67.62 68.39 [3,] 75.575 78.24 78.95 [4,] 85.590 88.41 89.42 [5,] 112.280 112.01 117.85 $n [1] 30 30 30 30 30 30 30 30 30 30 30 30 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 72.28104 71.61853 73.23315 73.07932 72.4376 72.04894 71.35741 69.05372 [2,] 83.88896 84.35147 84.84685 83.55068 82.3724 81.70106 80.58259 81.44628 [,9] [,10] [,11] [,12] [1,] 68.91662 70.09124 72.24277 72.88353 [2,] 79.95338 81.05876 84.23723 85.01647 $out [1] 117.63 124.65 157.33 147.38 171.89 171.95 132.71 126.02 121.18 115.45 $group [1] 1 2 3 4 5 6 7 8 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5uon11199713100.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] [,10] [,11] [1,] -3.76 -6.050 -4.170 -3.47 -3.43 -3.800 -3.990 -4.840 -3.250 -4.97 -4.080 [2,] -1.74 -1.530 -1.450 -1.63 -1.66 -1.960 -1.680 -1.500 -1.210 -0.54 -1.490 [3,] -1.07 -0.435 -0.115 -0.72 -1.14 -1.195 -1.165 -0.795 -0.555 0.67 -0.775 [4,] 0.28 1.980 2.200 0.75 0.06 -0.320 0.620 1.820 1.490 3.28 0.770 [5,] 2.60 5.170 6.460 3.91 1.68 1.850 1.990 6.800 3.790 8.84 2.790 [,12] [1,] -3.25 [2,] -1.44 [3,] -0.50 [4,] 0.12 [5,] 2.13 $n [1] 30 30 30 30 30 30 30 30 30 30 30 29 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -1.6527038 -1.4475199 -1.1679053 -1.40655197 -1.6361636 -1.6680862 [2,] -0.4872962 0.5775199 0.9379053 -0.03344803 -0.6438364 -0.7219138 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -1.8284746 -1.7527111 -1.3338615 -0.4319448 -1.4269359 -0.95770187 [2,] -0.5015254 0.1627111 0.2238615 1.7719448 -0.1230641 -0.04229813 $out [1] 9.62 6.85 18.99 7.02 -5.09 -6.34 13.68 -9.41 -9.39 14.68 [11] 19.65 43.89 -15.06 13.95 -9.95 -7.45 24.51 -11.49 -5.94 3.13 [21] 3.89 -5.39 -10.66 5.42 -6.85 -6.99 2.70 -39.24 -5.97 -6.91 [31] 3.98 -8.43 -9.81 6.31 19.60 15.79 -6.69 6.13 14.20 -8.42 [41] 5.96 -5.73 -5.86 9.77 12.18 -7.76 7.37 12.23 -7.13 9.63 [51] 7.78 5.59 -6.31 -4.81 -9.25 4.22 6.38 15.31 $group [1] 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 4 4 4 5 5 5 5 5 5 [26] 5 6 6 6 6 6 6 6 7 7 7 7 7 8 8 9 9 9 9 10 11 11 11 11 11 [51] 11 12 12 12 12 12 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-3.76000000000001, -1.74000000000001, -1.06999999999999, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/609st1199713100.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] [,10] [1,] 86.920 81.190 68.770 79.340 74.480 75.970 71.790 67.620 58.59 86.240 [2,] 87.280 83.385 72.005 81.855 76.945 77.895 73.840 68.385 60.39 86.610 [3,] 88.115 85.160 78.325 83.495 78.770 78.830 74.675 69.580 61.86 89.060 [4,] 89.925 86.210 82.320 84.755 80.305 80.675 76.100 70.365 64.33 92.425 [5,] 90.720 87.070 85.100 86.440 81.690 81.930 77.420 72.750 66.95 94.150 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [1,] 74.000 70.59 71.010 65.870 64.700 47.74 46.290 51.700 91.510 84.030 [2,] 75.820 72.79 72.275 68.210 65.990 49.12 47.635 76.030 93.195 85.210 [3,] 77.325 75.63 73.865 69.855 66.870 55.21 49.065 77.675 100.870 87.605 [4,] 80.975 76.57 75.925 70.905 68.035 58.75 53.100 108.800 104.020 103.485 [5,] 84.650 78.19 77.470 72.160 70.910 64.70 54.340 112.280 105.960 107.020 [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [1,] 94.450 87.910 75.660 67.620 50.890 44.570 45.85 54.760 75.820 79.98 [2,] 114.445 93.880 81.165 68.655 54.555 48.205 47.21 55.780 80.610 82.93 [3,] 123.600 98.705 88.380 72.610 56.300 53.340 48.25 57.995 89.385 88.31 [4,] 152.355 110.430 90.875 80.850 60.770 55.495 49.56 62.865 95.490 94.86 [5,] 171.950 124.650 96.730 87.480 66.280 61.850 50.51 69.600 107.630 104.73 [,31] [1,] NA [2,] NA [3,] NA [4,] NA [5,] 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 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 86.9086 83.8715 73.62026 82.17229 77.23748 77.56202 73.6442 68.67691 [2,] 89.3214 86.4485 83.02974 84.81771 80.30252 80.09798 75.7058 70.48309 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 60.06294 86.40774 74.97377 73.90592 72.20021 68.62579 65.93726 50.81769 [2,] 63.65706 91.71226 79.67623 77.35408 75.52979 71.08421 67.80274 59.60231 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 46.57238 62.72838 95.93264 79.26965 106.309 91.15643 83.9512 67.04778 [2,] 51.55762 92.62162 105.80736 95.94035 140.891 106.25357 92.8088 78.17222 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [1,] 53.4653 50.01498 47.17815 54.76348 82.59813 82.86865 NA [2,] 59.1347 56.66502 49.32185 61.22652 96.17187 93.75135 NA $out [1] 67.44 77.06 79.23 $group [1] 10 10 28 $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] NA Warning message: In bxp(list(stats = c(86.92, 87.28, 88.115, 89.925, 90.72, 81.19, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/73jmd1199713100.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,] 75.40833 74.4350 74.68875 [2,] 75.95567 75.7725 76.14125 [3,] 77.62100 77.6950 76.82250 [4,] 78.84733 78.2775 78.09875 [5,] 79.73967 79.0400 78.71125 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 76.30209 76.55245 75.92967 [2,] 78.93991 78.83755 77.71533 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(75.4083333333333, 75.9556666666667, 77.621, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/18nsk1199713100.ps tmp/18nsk1199713100.png") > system("convert tmp/2eim51199713100.ps tmp/2eim51199713100.png") > system("convert tmp/31z0n1199713100.ps tmp/31z0n1199713100.png") > system("convert tmp/4mc5h1199713100.ps tmp/4mc5h1199713100.png") > system("convert tmp/5uon11199713100.ps tmp/5uon11199713100.png") > system("convert tmp/609st1199713100.ps tmp/609st1199713100.png") > system("convert tmp/73jmd1199713100.ps tmp/73jmd1199713100.png") > > > proc.time() user system elapsed 1.595 1.049 2.790