R version 2.7.0 (2008-04-22) 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(15.14,15.09,15.17,15.18,15.21,15.27,15.28,15.3,15.34,15.33,15.27,15.35,15.38,15.38,15.39,15.4,15.39,15.43,15.43,15.47,15.48,15.47,15.58,15.56,15.61,15.56,15.62,15.63,15.58,15.65,15.79,15.76,15.77,15.79,15.87,15.79,15.9,15.96,16.05,16.18,16.29,16.43,16.38,16.39,16.35,16.48,16.52,16.44,16.46,16.52,16.47,16.59,16.59,16.59,16.54,16.48,16.47,16.56,16.61,16.57,16.72,16.69,16.72,16.81,16.75,16.85,16.84,16.92,17.02,17.11,17.2,17.3,17.37,17.42,17.51,17.56,17.62,17.59,17.78,17.73,17.79,17.85,17.86,17.79,17.97,17.96,18.03,18.02,18.03,18.14,18.16,18.24,18.28,18.18,18.19,18.32) > 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] 96 > (np <- floor(n / par1)) [1] 8 > 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] 8 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 15.14 15.38 15.61 15.90 16.46 16.72 17.37 17.97 NA [2,] 15.09 15.38 15.56 15.96 16.52 16.69 17.42 17.96 NA [3,] 15.17 15.39 15.62 16.05 16.47 16.72 17.51 18.03 NA [4,] 15.18 15.40 15.63 16.18 16.59 16.81 17.56 18.02 NA [5,] 15.21 15.39 15.58 16.29 16.59 16.75 17.62 18.03 NA [6,] 15.27 15.43 15.65 16.43 16.59 16.85 17.59 18.14 NA [7,] 15.28 15.43 15.79 16.38 16.54 16.84 17.78 18.16 NA [8,] 15.30 15.47 15.76 16.39 16.48 16.92 17.73 18.24 NA [9,] 15.34 15.48 15.77 16.35 16.47 17.02 17.79 18.28 NA [10,] 15.33 15.47 15.79 16.48 16.56 17.11 17.85 18.18 NA [11,] 15.27 15.58 15.87 16.52 16.61 17.20 17.86 18.19 NA [12,] 15.35 15.56 15.79 16.44 16.57 17.30 17.79 18.32 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -0.05 0.00 -0.05 0.06 0.06 -0.03 0.05 -0.01 NA [2,] 0.08 0.01 0.06 0.09 -0.05 0.03 0.09 0.07 NA [3,] 0.01 0.01 0.01 0.13 0.12 0.09 0.05 -0.01 NA [4,] 0.03 -0.01 -0.05 0.11 0.00 -0.06 0.06 0.01 NA [5,] 0.06 0.04 0.07 0.14 0.00 0.10 -0.03 0.11 NA [6,] 0.01 0.00 0.14 -0.05 -0.05 -0.01 0.19 0.02 NA [7,] 0.02 0.04 -0.03 0.01 -0.06 0.08 -0.05 0.08 NA [8,] 0.04 0.01 0.01 -0.04 -0.01 0.10 0.06 0.04 NA [9,] -0.01 -0.01 0.02 0.13 0.09 0.09 0.06 -0.10 NA [10,] -0.06 0.11 0.08 0.04 0.05 0.09 0.01 0.01 NA [11,] 0.08 -0.02 -0.08 -0.08 -0.04 0.10 -0.07 0.13 NA [12,] 0.03 0.05 0.11 0.02 0.15 0.07 0.18 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/1xmvj1218463902.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/2x76n1218463902.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/3zi9s1218463902.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/4k5ix1218463902.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] [1,] 15.140 15.090 15.170 15.180 15.210 15.27 15.28 15.300 15.340 15.33 15.270 [2,] 15.495 15.470 15.505 15.515 15.485 15.54 15.61 15.615 15.625 15.63 15.725 [3,] 16.180 16.240 16.260 16.385 16.440 16.51 16.46 16.435 16.410 16.52 16.565 [4,] 17.045 17.055 17.115 17.185 17.185 17.22 17.31 17.325 17.405 17.48 17.530 [5,] 17.970 17.960 18.030 18.020 18.030 18.14 18.16 18.240 18.280 18.18 18.190 [,12] [1,] 15.350 [2,] 15.675 [3,] 16.505 [4,] 17.545 [5,] 18.320 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 15.31415 15.35460 15.36063 15.45211 15.49036 15.57153 15.51036 15.47977 [2,] 17.04585 17.12540 17.15937 17.31789 17.38964 17.44847 17.40964 17.39023 [,9] [,10] [,11] [,12] [1,] 15.41567 15.48656 15.5567 15.46039 [2,] 17.40433 17.55344 17.5733 17.54961 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(15.14, 15.495, 16.18, 17.045, 17.97, 15.09, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5ot9u1218463902.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] [1,] -0.050 -0.050 -0.010 -0.060 -0.030 -0.050 -0.060 -4.000000e-02 -0.10 [2,] -0.040 0.020 0.010 -0.030 0.020 -0.030 -0.040 -8.881784e-16 -0.01 [3,] -0.005 0.065 0.030 0.005 0.065 0.005 0.015 2.500000e-02 0.04 [4,] 0.055 0.085 0.105 0.045 0.105 0.080 0.060 5.000000e-02 0.09 [5,] 0.060 0.090 0.130 0.110 0.140 0.190 0.080 1.000000e-01 0.13 [,10] [,11] [,12] [1,] -0.060 -0.080 0.02 [2,] 0.010 -0.075 0.04 [3,] 0.045 -0.030 0.07 [4,] 0.085 0.090 0.13 [5,] 0.110 0.130 0.18 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.05806836 0.02869007 -0.02306836 -0.03689608 0.01751778 -0.05644758 [2,] 0.04806836 0.10130993 0.08306836 0.04689608 0.11248222 0.06644758 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.04086144 -0.002930718 -0.01586144 0.003103923 -0.12217137 0.01625345 [2,] 0.07086144 0.052930718 0.09586144 0.086896077 0.06217137 0.12374655 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.0500000000000007, -0.0399999999999983, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6styy1218463902.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,] 15.090 15.380 15.560 15.900 16.460 16.690 17.370 17.960 NA [2,] 15.175 15.390 15.615 16.115 16.475 16.735 17.535 18.025 NA [3,] 15.270 15.430 15.705 16.365 16.550 16.845 17.675 18.150 NA [4,] 15.315 15.475 15.790 16.435 16.590 17.065 17.790 18.215 NA [5,] 15.350 15.580 15.870 16.520 16.610 17.300 17.860 18.320 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 15.20615 15.39123 15.62518 16.21905 16.49755 16.69448 17.55869 18.06334 [2,] 15.33385 15.46877 15.78482 16.51095 16.60245 16.99552 17.79131 18.23666 [,9] [1,] NA [2,] NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(15.09, 15.175, 15.27, 15.315, 15.35, 15.38, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/734zy1218463902.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,] 16.31875 16.1800 16.19375 [2,] 16.39563 16.3225 16.24500 [3,] 16.50937 16.4375 16.35125 [4,] 16.57937 16.5075 16.47875 [5,] 16.64000 16.5650 16.58125 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 16.42557 16.35312 16.24464 [2,] 16.59318 16.52188 16.45786 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(16.31875, 16.395625, 16.509375, 16.579375, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1xmvj1218463902.ps tmp/1xmvj1218463902.png") > system("convert tmp/2x76n1218463902.ps tmp/2x76n1218463902.png") > system("convert tmp/3zi9s1218463902.ps tmp/3zi9s1218463902.png") > system("convert tmp/4k5ix1218463902.ps tmp/4k5ix1218463902.png") > system("convert tmp/5ot9u1218463902.ps tmp/5ot9u1218463902.png") > system("convert tmp/6styy1218463902.ps tmp/6styy1218463902.png") > system("convert tmp/734zy1218463902.ps tmp/734zy1218463902.png") > > > proc.time() user system elapsed 2.569 1.582 2.824