R version 2.10.1 (2009-12-14) Copyright (C) 2009 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(70.5,70.65,70.71,15.80,170.99,15.26,17.21,16.97,16.80,18.37,18.53,18.708,18.73,18.63,19.15,19.76,20.02,20.26,20.32,20.52,20.81,21.26,22.15,23.05,270.01,270.69,70.38,26.33,27.36,28.55,30.10,32.20,36.28,700.91,7070.66,707.83,709.97,52.20,55.67,59.92,65.15,70.17,770.6,78.23,79.27,80.707,81.701,83.97,86.83,89.61,91.79,970.32,96.56,97.98,100.00,101.63,102.60,103.75,106.39,109.02,110.81,112.57,1170.93,118.13,120.7,122.707,127.97,127.87) > 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] 68 > (np <- floor(n / par1)) [1] 5 > 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] 6 6 6 6 6 6 6 6 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 70.500 18.73 270.01 709.970 86.83 110.810 [2,] 70.650 18.63 270.69 52.200 89.61 112.570 [3,] 70.710 19.15 70.38 55.670 91.79 1170.930 [4,] 15.800 19.76 26.33 59.920 970.32 118.130 [5,] 170.990 20.02 27.36 65.150 96.56 120.700 [6,] 15.260 20.26 28.55 70.170 97.98 122.707 [7,] 17.210 20.32 30.10 770.600 100.00 127.970 [8,] 16.970 20.52 32.20 78.230 101.63 127.870 [9,] 16.800 20.81 36.28 79.270 102.60 NA [10,] 18.370 21.26 700.91 80.707 103.75 NA [11,] 18.530 22.15 7070.66 81.701 106.39 NA [12,] 18.708 23.05 707.83 83.970 109.02 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.150 -0.10 0.68 -657.770 2.78 1.760 [2,] 0.060 0.52 -200.31 3.470 2.18 1058.360 [3,] -54.910 0.61 -44.05 4.250 878.53 -1052.800 [4,] 155.190 0.26 1.03 5.230 -873.76 2.570 [5,] -155.730 0.24 1.19 5.020 1.42 2.007 [6,] 1.950 0.06 1.55 700.430 2.02 5.263 [7,] -0.240 0.20 2.10 -692.370 1.63 -0.100 [8,] -0.170 0.29 4.08 1.040 0.97 NA [9,] 1.570 0.45 664.63 1.437 1.15 NA [10,] 0.160 0.89 6369.75 0.994 2.64 NA [11,] 0.178 0.90 -6362.83 2.269 2.63 NA [12,] 0.022 246.96 2.14 2.860 1.79 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/rcomp/tmp/1tc081272464741.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/rcomp/tmp/2tc081272464741.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/rcomp/tmp/3mlht1272464741.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/rcomp/tmp/4mlht1272464741.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,] 18.73 18.63 19.150 15.800 20.020 15.260 17.21 16.970 16.80 18.370 [2,] 70.50 52.20 55.670 19.760 27.360 20.260 20.32 20.520 20.81 21.260 [3,] 98.82 80.13 70.545 43.125 80.855 49.360 65.05 55.215 36.28 80.707 [4,] 270.01 112.57 91.790 118.130 120.700 97.980 127.97 101.630 79.27 103.750 [5,] 270.01 112.57 91.790 118.130 170.990 122.707 127.97 127.870 102.60 103.750 [,11] [,12] [1,] 18.530 18.708 [2,] 22.150 23.050 [3,] 81.701 83.970 [4,] 106.390 109.020 [5,] 106.390 109.020 $n [1] 6 6 6 6 6 6 6 6 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -29.87039 41.1894 47.24643 -20.32683 20.64769 -0.7719103 -4.387727 [2,] 227.51039 119.0706 93.84357 106.57683 141.06231 99.4919103 134.487727 [,8] [,9] [,10] [,11] [,12] [1,] 2.89643 -5.027689 22.41977 22.17723 23.22381 [2,] 107.53357 77.587689 138.99423 141.22477 144.71619 $out [1] 709.97 270.69 1170.93 970.32 770.60 700.91 7070.66 707.83 $group [1] 1 2 3 4 7 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(18.73, 70.5, 98.82, 270.01, 270.01, 18.63, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5wuye1272464741.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] [,12] [1,] -0.100 0.06 -54.91 0.26 0.240 0.060 -0.24 -0.17 1.150 0.160 0.178 1.79 [2,] -0.100 0.06 -54.91 0.26 0.240 1.550 -0.24 0.29 1.150 0.890 0.178 1.79 [3,] 0.415 1.35 -21.72 1.80 1.305 1.985 0.05 0.97 1.437 0.994 0.900 2.14 [4,] 1.760 3.47 4.25 5.23 2.007 5.263 1.63 1.04 1.570 2.640 2.269 2.86 [5,] 2.780 3.47 4.25 5.23 2.007 5.263 2.10 1.04 1.570 2.640 2.630 2.86 $n [1] 6 6 6 6 6 6 6 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -0.78476 -0.8495601 -59.88011 -1.405811 0.1652279 -0.4100049 -1.156210 [2,] 1.61476 3.5495601 16.44011 5.005811 2.4447721 4.3800049 1.256210 [,8] [,9] [,10] [,11] [,12] [1,] 0.4400519 1.140229 -0.2425456 -0.5774953 1.383941 [2,] 1.4999481 1.733771 2.2305456 2.3774953 2.896059 $out [1] -657.770 -200.310 1058.360 878.530 -1052.800 155.190 -873.760 [8] -155.730 5.020 700.430 -692.370 4.080 0.450 664.630 [15] 6369.750 -6362.830 0.022 246.960 $group [1] 1 2 2 3 3 4 4 5 5 6 7 8 9 9 10 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.100000000000001, -0.100000000000001, 0.415000000000006, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6wuye1272464741.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] [1,] 15.260 18.630 26.330 52.2000 86.830 110.8100 [2,] 16.885 19.455 29.325 62.5350 94.175 115.3500 [3,] 18.450 20.290 53.330 78.7500 100.815 121.7035 [4,] 70.575 21.035 485.800 82.8355 105.070 127.9200 [5,] 70.710 23.050 707.830 83.9700 109.020 127.9700 $n [1] 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -6.03837 19.56935 -154.8713 69.4908 95.84572 114.6817 [2,] 42.93837 21.01065 261.5313 88.0092 105.78428 128.7253 $out [1] 170.99 7070.66 709.97 770.60 970.32 1170.93 $group [1] 1 3 4 4 5 6 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(15.26, 16.885, 18.45, 70.575, 70.71, 18.63, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7wuye1272464741.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,] 51.15200 36.2800 50.04000 [2,] 73.18333 52.2875 61.05000 [3,] 181.34970 75.3375 65.15250 [4,] 206.42583 81.2780 74.33500 [5,] 246.43833 98.8200 81.82125 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 120.5769 62.11474 59.09312 [2,] 242.1225 88.56026 71.21188 $out [1] 1459.8862 152.3962 $group [1] 1 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(51.152, 73.1833333333333, 181.3497, 206.425833333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1tc081272464741.ps tmp/1tc081272464741.png",intern=TRUE)) character(0) > try(system("convert tmp/2tc081272464741.ps tmp/2tc081272464741.png",intern=TRUE)) character(0) > try(system("convert tmp/3mlht1272464741.ps tmp/3mlht1272464741.png",intern=TRUE)) character(0) > try(system("convert tmp/4mlht1272464741.ps tmp/4mlht1272464741.png",intern=TRUE)) character(0) > try(system("convert tmp/5wuye1272464741.ps tmp/5wuye1272464741.png",intern=TRUE)) character(0) > try(system("convert tmp/6wuye1272464741.ps tmp/6wuye1272464741.png",intern=TRUE)) character(0) > try(system("convert tmp/7wuye1272464741.ps tmp/7wuye1272464741.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.570 1.500 2.719