R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(20.98,20.1,20.61,20.27,20.08,23.58,22.31,22.89,21.78,22.19,22.58,22.78,25.06,25.16,25.47,25.34,24.2,25.32,25.57,25.76,24.79,23.14,22.66,22.06,24.26,23.15,22.92,21.43,21.56,23.48,24.35,24.83,24.19,23.58,23.58,24.35,27.18,25.69,24.81,23.26,23.49,26.86,27.12,27.66,26.26,25.51,24.63,23.57,27.63,25.85,26.09,24.47,24.19,25.09,25.26,25.58,24.76,25.02,24.24,24.14,28.69,26.74,26.48,24.45,23.88,26.58,26.23,28.63,26.81,26.56,26.64,26.8,28.37,27.13,28.44,28.62,27.28,31.32,31.26,31.41,31.76,32.72,32.15,33.62,35.97,33.78,33.77,32.75,32.55,33.22,32.88,31.56,30.27,28.65,27.89,27.07,30.8,28.38,27.5,28,28.02,29.2,27.59,27.22,27.16,26.31,25.67,26.41,28.34,25.43,23.72,23.33,23.8,27.7,26.28,27.51,27.93,28.76,28.65,29.52,31.23,27.9,27.87,27.52,27.59,31.2,30.22,30.62,31.52,30.59,31.42,31.95) > 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,] 20.98 25.06 24.26 27.18 27.63 28.69 28.37 35.97 30.80 28.34 31.23 NA [2,] 20.10 25.16 23.15 25.69 25.85 26.74 27.13 33.78 28.38 25.43 27.90 NA [3,] 20.61 25.47 22.92 24.81 26.09 26.48 28.44 33.77 27.50 23.72 27.87 NA [4,] 20.27 25.34 21.43 23.26 24.47 24.45 28.62 32.75 28.00 23.33 27.52 NA [5,] 20.08 24.20 21.56 23.49 24.19 23.88 27.28 32.55 28.02 23.80 27.59 NA [6,] 23.58 25.32 23.48 26.86 25.09 26.58 31.32 33.22 29.20 27.70 31.20 NA [7,] 22.31 25.57 24.35 27.12 25.26 26.23 31.26 32.88 27.59 26.28 30.22 NA [8,] 22.89 25.76 24.83 27.66 25.58 28.63 31.41 31.56 27.22 27.51 30.62 NA [9,] 21.78 24.79 24.19 26.26 24.76 26.81 31.76 30.27 27.16 27.93 31.52 NA [10,] 22.19 23.14 23.58 25.51 25.02 26.56 32.72 28.65 26.31 28.76 30.59 NA [11,] 22.58 22.66 23.58 24.63 24.24 26.64 32.15 27.89 25.67 28.65 31.42 NA [12,] 22.78 22.06 24.35 23.57 24.14 26.80 33.62 27.07 26.41 29.52 31.95 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.88 0.10 -1.11 -1.49 -1.78 -1.95 -1.24 -2.19 -2.42 -2.91 -3.33 NA [2,] 0.51 0.31 -0.23 -0.88 0.24 -0.26 1.31 -0.01 -0.88 -1.71 -0.03 NA [3,] -0.34 -0.13 -1.49 -1.55 -1.62 -2.03 0.18 -1.02 0.50 -0.39 -0.35 NA [4,] -0.19 -1.14 0.13 0.23 -0.28 -0.57 -1.34 -0.20 0.02 0.47 0.07 NA [5,] 3.50 1.12 1.92 3.37 0.90 2.70 4.04 0.67 1.18 3.90 3.61 NA [6,] -1.27 0.25 0.87 0.26 0.17 -0.35 -0.06 -0.34 -1.61 -1.42 -0.98 NA [7,] 0.58 0.19 0.48 0.54 0.32 2.40 0.15 -1.32 -0.37 1.23 0.40 NA [8,] -1.11 -0.97 -0.64 -1.40 -0.82 -1.82 0.35 -1.29 -0.06 0.42 0.90 NA [9,] 0.41 -1.65 -0.61 -0.75 0.26 -0.25 0.96 -1.62 -0.85 0.83 -0.93 NA [10,] 0.39 -0.48 0.00 -0.88 -0.78 0.08 -0.57 -0.76 -0.64 -0.11 0.83 NA [11,] 0.20 -0.60 0.77 -1.06 -0.10 0.16 1.47 -0.82 0.74 0.87 0.53 NA [12,] 2.28 2.20 2.83 4.06 4.55 1.57 2.35 3.73 1.93 1.71 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/wessaorg/rcomp/tmp/1zjqv1319540337.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2r37n1319540337.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3ml2m1319540337.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4pu6n1319540337.ps",horizontal=F,onefile=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,] 20.980 23.150 20.610 20.270 20.080 23.480 22.310 22.890 21.780 22.190 [2,] 26.120 25.295 24.265 23.295 23.645 25.205 25.415 25.670 24.775 24.300 [3,] 28.340 25.850 26.090 24.470 24.190 26.860 26.280 27.510 26.810 26.310 [4,] 29.745 27.515 27.685 27.760 27.435 30.200 28.905 29.625 29.100 28.705 [5,] 31.230 28.380 28.440 32.750 32.550 33.220 32.880 31.560 31.760 32.720 [,11] [,12] [1,] 22.58 22.060 [2,] 23.91 23.855 [3,] 25.67 26.410 [4,] 28.27 28.295 [5,] 32.15 33.620 $n [1] 11 11 11 11 11 11 11 11 11 11 11 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 26.61309 24.79242 24.46075 22.34293 22.38449 24.48044 24.61741 25.62589 [2,] 30.06691 26.90758 27.71925 26.59707 25.99551 29.23956 27.94259 29.39411 [,9] [,10] [,11] [,12] [1,] 24.74962 24.21151 23.59295 24.29484 [2,] 28.87038 28.40849 27.74705 28.52516 $out [1] 35.97 20.10 33.78 33.77 $group [1] 1 2 2 3 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(20.98, 26.12, 28.34, 29.745, 31.23, 23.15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/54m6c1319540337.ps",horizontal=F,onefile=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.330 -1.710 -2.030 -1.140 0.670 -1.610 -0.37 -1.820 -1.650 -0.88 -1.060 [2,] -2.305 -0.570 -1.520 -0.425 1.150 -1.125 0.17 -1.200 -0.890 -0.70 -0.350 [3,] -1.780 -0.030 -0.390 -0.190 2.700 -0.340 0.40 -0.820 -0.610 -0.48 0.200 [4,] -1.175 0.275 -0.235 0.100 3.555 0.210 0.56 0.145 0.335 0.04 0.755 [5,] 0.100 1.310 0.500 0.470 4.040 0.870 0.58 0.900 0.960 0.83 1.470 [,12] [1,] 1.570 [2,] 1.930 [3,] 2.315 [4,] 3.730 [5,] 4.550 $n [1] 11 11 11 11 11 11 11 11 11 11 11 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2.318318 -0.4325478 -1.0021585 -0.44010366 1.554287 -0.9759779 0.2142087 [2,] -1.241682 0.3725478 0.2221585 0.06010366 3.845713 0.2959779 0.5857913 [,8] [,9] [,10] [,11] [,12] [1,] -1.4607418 -1.19357521 -0.8325271 -0.3264087 1.415648 [2,] -0.1792582 -0.02642479 -0.1274729 0.7264087 3.214352 $out [1] -1.34 2.40 -1.32 1.23 $group [1] 4 7 7 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-3.33, -2.305, -1.78, -1.175, 0.100000000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6bip51319540337.ps",horizontal=F,onefile=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,] 20.080 22.060 21.430 23.26 24.140 26.230 27.130 27.070 25.670 23.330 [2,] 20.440 23.670 23.035 24.10 24.355 26.355 28.405 29.460 26.785 24.615 [3,] 21.985 25.110 23.580 25.60 25.055 26.610 31.290 32.650 27.545 27.605 [4,] 22.680 25.405 24.305 26.99 25.715 26.805 31.955 33.495 28.200 28.495 [5,] 23.580 25.760 24.830 27.66 27.630 26.810 33.620 35.970 29.200 29.520 [,11] [,12] [1,] 27.520 NA [2,] 27.885 NA [3,] 30.605 NA [4,] 31.325 NA [5,] 31.950 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,] 20.96332 24.31865 23.00074 24.28185 24.43469 26.40475 29.67082 30.80961 [2,] 23.00668 25.90135 24.15926 26.91815 25.67531 26.81525 32.90918 34.49039 [,9] [,10] [,11] [,12] [1,] 26.89961 25.83531 29.03599 NA [2,] 28.19039 29.37469 32.17401 NA $out [1] 28.69 24.45 23.88 28.63 30.80 $group [1] 6 6 6 6 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" NA Warning message: In bxp(list(stats = c(20.08, 20.44, 21.985, 22.68, 23.58, 22.06, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7zn7j1319540337.ps",horizontal=F,onefile=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,] 25.14909 24.190 25.52750 [2,] 26.22682 25.760 26.02500 [3,] 26.60455 26.295 26.45375 [4,] 27.39182 26.835 27.40375 [5,] 28.04636 28.340 27.93250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 26.07318 25.80469 25.82489 [2,] 27.13591 26.78531 27.08261 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(25.1490909090909, 26.2268181818182, 26.6045454545455, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1zjqv1319540337.ps tmp/1zjqv1319540337.png",intern=TRUE)) character(0) > try(system("convert tmp/2r37n1319540337.ps tmp/2r37n1319540337.png",intern=TRUE)) character(0) > try(system("convert tmp/3ml2m1319540337.ps tmp/3ml2m1319540337.png",intern=TRUE)) character(0) > try(system("convert tmp/4pu6n1319540337.ps tmp/4pu6n1319540337.png",intern=TRUE)) character(0) > try(system("convert tmp/54m6c1319540337.ps tmp/54m6c1319540337.png",intern=TRUE)) character(0) > try(system("convert tmp/6bip51319540337.ps tmp/6bip51319540337.png",intern=TRUE)) character(0) > try(system("convert tmp/7zn7j1319540337.ps tmp/7zn7j1319540337.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.575 0.338 1.927