R version 2.9.0 (2009-04-17) 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(8.82000 + ,8.80000 + ,8.82000 + ,8.58000 + ,8.54000 + ,8.42000 + ,8.43000 + ,8.44000 + ,8.09000 + ,7.69000 + ,7.56000 + ,7.54000 + ,7.40000 + ,7.39000 + ,7.37000 + ,7.31000 + ,7.35000 + ,7.26000 + ,7.37000 + ,7.35000 + ,7.33000 + ,7.32000 + ,7.31000 + ,7.33000 + ,7.32000 + ,7.27000 + ,7.48000 + ,7.70000 + ,7.77000 + ,7.80000 + ,7.84000 + ,7.81000 + ,7.78000 + ,7.82000 + ,7.80000 + ,7.81000 + ,7.80000 + ,7.66000 + ,7.41000 + ,7.35000 + ,7.39000 + ,7.32000 + ,7.32000 + ,7.30000 + ,7.29000 + ,7.26000 + ,7.22000 + ,7.21000 + ,7.21000 + ,7.21000 + ,7.20000 + ,7.19000 + ,7.18000 + ,7.12000 + ,7.12000 + ,7.07000 + ,7.08000 + ,7.05000 + ,7.06000 + ,7.07000 + ,7.08000 + ,7.08000 + ,7.09000 + ,7.07000 + ,7.06000 + ,6.99000 + ,6.99000 + ,6.99000 + ,6.98000 + ,6.96000 + ,6.95000 + ,6.91000 + ,6.91000 + ,6.87000 + ,6.91000 + ,6.89000 + ,6.88000 + ,6.90000 + ,6.91000 + ,6.85000 + ,6.86000 + ,6.82000 + ,6.80000 + ,6.83000 + ,6.84000 + ,6.89000 + ,7.14000 + ,7.21000 + ,7.25000 + ,7.31000 + ,7.30000 + ,7.48000 + ,7.49000 + ,7.40000 + ,7.44000 + ,7.42000 + ,7.14000 + ,7.24000 + ,7.33000 + ,7.61000 + ,7.66000 + ,7.69000 + ,7.70000 + ,7.68000 + ,7.71000 + ,7.71000 + ,7.72000 + ,7.68000 + ,7.72000 + ,7.74000 + ,7.76000 + ,7.90000 + ,7.97000 + ,7.96000 + ,7.95000 + ,7.97000 + ,7.93000 + ,7.99000 + ,7.96000 + ,7.92000 + ,7.97000 + ,7.98000 + ,8.00000 + ,8.04000 + ,8.17000 + ,8.29000 + ,8.26000 + ,8.30000 + ,8.32000 + ,8.28000 + ,8.27000 + ,8.32000) > 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] 132 > (np <- floor(n / par1)) [1] 11 > 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] 11 11 11 11 11 11 11 11 11 11 11 11 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 8.82 7.40 7.32 7.80 7.21 7.08 6.91 6.84 7.14 7.72 7.97 NA [2,] 8.80 7.39 7.27 7.66 7.21 7.08 6.87 6.89 7.24 7.74 7.98 NA [3,] 8.82 7.37 7.48 7.41 7.20 7.09 6.91 7.14 7.33 7.76 8.00 NA [4,] 8.58 7.31 7.70 7.35 7.19 7.07 6.89 7.21 7.61 7.90 8.04 NA [5,] 8.54 7.35 7.77 7.39 7.18 7.06 6.88 7.25 7.66 7.97 8.17 NA [6,] 8.42 7.26 7.80 7.32 7.12 6.99 6.90 7.31 7.69 7.96 8.29 NA [7,] 8.43 7.37 7.84 7.32 7.12 6.99 6.91 7.30 7.70 7.95 8.26 NA [8,] 8.44 7.35 7.81 7.30 7.07 6.99 6.85 7.48 7.68 7.97 8.30 NA [9,] 8.09 7.33 7.78 7.29 7.08 6.98 6.86 7.49 7.71 7.93 8.32 NA [10,] 7.69 7.32 7.82 7.26 7.05 6.96 6.82 7.40 7.71 7.99 8.28 NA [11,] 7.56 7.31 7.80 7.22 7.06 6.95 6.80 7.44 7.72 7.96 8.27 NA [12,] 7.54 7.33 7.81 7.21 7.07 6.91 6.83 7.42 7.68 7.92 8.32 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.02 -0.01 -0.05 -0.14 0.00 0.00 -0.04 0.05 0.10 0.02 0.01 NA [2,] 0.02 -0.02 0.21 -0.25 -0.01 0.01 0.04 0.25 0.09 0.02 0.02 NA [3,] -0.24 -0.06 0.22 -0.06 -0.01 -0.02 -0.02 0.07 0.28 0.14 0.04 NA [4,] -0.04 0.04 0.07 0.04 -0.01 -0.01 -0.01 0.04 0.05 0.07 0.13 NA [5,] -0.12 -0.09 0.03 -0.07 -0.06 -0.07 0.02 0.06 0.03 -0.01 0.12 NA [6,] 0.01 0.11 0.04 0.00 0.00 0.00 0.01 -0.01 0.01 -0.01 -0.03 NA [7,] 0.01 -0.02 -0.03 -0.02 -0.05 0.00 -0.06 0.18 -0.02 0.02 0.04 NA [8,] -0.35 -0.02 -0.03 -0.01 0.01 -0.01 0.01 0.01 0.03 -0.04 0.02 NA [9,] -0.40 -0.01 0.04 -0.03 -0.03 -0.02 -0.04 -0.09 0.00 0.06 -0.04 NA [10,] -0.13 -0.01 -0.02 -0.04 0.01 -0.01 -0.02 0.04 0.01 -0.03 -0.01 NA [11,] -0.02 0.02 0.01 -0.01 0.01 -0.04 0.03 -0.02 -0.04 -0.04 0.05 NA [12,] -0.14 -0.01 -0.01 0.00 0.01 0.00 0.01 -0.28 0.04 0.05 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/1ftwp1264519810.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/2k6qs1264519810.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/3j3ly1264519810.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/4zh8r1264519810.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] [,11] [,12] [1,] 6.84 6.870 6.91 6.89 6.880 6.90 6.910 6.850 6.860 6.820 6.80 6.830 [2,] 7.11 7.145 7.17 7.20 7.215 7.19 7.210 7.185 7.185 7.155 7.14 7.140 [3,] 7.32 7.270 7.37 7.35 7.390 7.32 7.370 7.480 7.490 7.400 7.44 7.420 [4,] 7.76 7.700 7.62 7.80 7.870 7.88 7.895 7.890 7.855 7.765 7.76 7.745 [5,] 7.97 7.980 8.00 8.58 8.540 8.42 8.430 8.440 8.320 8.280 8.27 8.320 $n [1] 11 11 11 11 11 11 11 11 11 11 11 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 7.010348 7.005605 7.155625 7.064167 7.077966 6.991292 7.043674 7.144147 [2,] 7.629652 7.534395 7.584375 7.635833 7.702034 7.648708 7.696326 7.815853 [,9] [,10] [,11] [,12] [1,] 7.17082 7.109403 7.14464 7.131785 [2,] 7.80918 7.690597 7.73536 7.708215 $out [1] 8.82 8.80 8.82 $group [1] 1 2 3 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(6.84, 7.11, 7.32, 7.76, 7.97, 6.87, 7.145, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5e49o1264519810.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,] -0.050 -0.020 -0.240 -0.04 -0.12 -0.010 -0.060 -0.040 -0.090 -0.040 -0.040 [2,] -0.030 0.000 -0.040 -0.01 -0.07 -0.005 -0.025 -0.025 -0.040 -0.025 -0.030 [3,] 0.000 0.020 -0.010 0.04 -0.01 0.000 -0.020 -0.010 -0.030 -0.010 -0.010 [4,] 0.015 0.065 0.105 0.06 0.03 0.010 0.015 0.010 -0.005 0.000 0.015 [5,] 0.050 0.090 0.280 0.13 0.12 0.010 0.040 0.030 0.040 0.010 0.050 [,12] [1,] -0.01 [2,] -0.01 [3,] 0.00 [4,] 0.01 [5,] 0.01 $n [1] 11 11 11 11 11 11 11 11 11 11 11 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.02143746 -0.01096522 -0.07907625 0.006652845 -0.05763879 -0.007145819 [2,] 0.02143746 0.05096522 0.05907625 0.073347155 0.03763879 0.007145819 [,7] [,8] [,9] [,10] [,11] [1,] -0.039055517 -0.026673577 -0.04667358 -0.021909698 -0.03143746 [2,] -0.000944483 0.006673577 -0.01332642 0.001909698 0.01143746 [,12] [1,] -0.009992797 [2,] 0.009992797 $out [1] -0.14 0.10 0.21 -0.25 0.25 0.11 0.04 -0.03 0.18 -0.35 -0.40 0.06 [13] -0.13 0.04 -0.14 -0.28 0.04 0.05 $group [1] 1 1 2 2 2 6 6 6 7 8 9 9 10 10 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(-0.0500000000000007, -0.0299999999999998, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/67svt1264519810.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] [,11] [,12] [1,] 7.540 7.260 7.27 7.210 7.050 6.910 6.800 6.840 7.140 7.720 7.970 NA [2,] 7.890 7.315 7.59 7.275 7.070 6.970 6.840 7.175 7.470 7.830 8.020 NA [3,] 8.435 7.340 7.79 7.320 7.120 6.990 6.875 7.305 7.680 7.940 8.265 NA [4,] 8.690 7.370 7.81 7.400 7.195 7.075 6.905 7.430 7.705 7.965 8.295 NA [5,] 8.820 7.400 7.84 7.410 7.210 7.090 6.910 7.490 7.720 7.990 8.320 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 8.070115 7.314914 7.689657 7.262987 7.062987 6.942109 6.845353 7.188693 [2,] 8.799885 7.365086 7.890343 7.377013 7.177013 7.037891 6.904647 7.421307 [,9] [,10] [,11] [,12] [1,] 7.572815 7.878426 8.13957 NA [2,] 7.787185 8.001574 8.39043 NA $out [1] 7.80 7.66 $group [1] 4 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" NA Warning message: In bxp(list(stats = c(7.54, 7.89, 8.435, 8.69, 8.82, 7.26, 7.315, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7znvr1264519810.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,] 7.458182 7.270 7.39500 [2,] 7.470000 7.335 7.43875 [3,] 7.516364 7.380 7.48000 [4,] 7.556818 7.430 7.53625 [5,] 7.567273 7.490 7.55250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 7.476765 7.33667 7.43553 [2,] 7.555962 7.42333 7.52447 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(7.45818181818182, 7.47, 7.51636363636364, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1ftwp1264519810.ps tmp/1ftwp1264519810.png",intern=TRUE)) character(0) > try(system("convert tmp/2k6qs1264519810.ps tmp/2k6qs1264519810.png",intern=TRUE)) character(0) > try(system("convert tmp/3j3ly1264519810.ps tmp/3j3ly1264519810.png",intern=TRUE)) character(0) > try(system("convert tmp/4zh8r1264519810.ps tmp/4zh8r1264519810.png",intern=TRUE)) character(0) > try(system("convert tmp/5e49o1264519810.ps tmp/5e49o1264519810.png",intern=TRUE)) character(0) > try(system("convert tmp/67svt1264519810.ps tmp/67svt1264519810.png",intern=TRUE)) character(0) > try(system("convert tmp/7znvr1264519810.ps tmp/7znvr1264519810.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.230 0.976 1.675