R version 2.6.2 (2008-02-08) 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. Natural language support but running in an English locale 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(120.05,120.05,120.08,120.12,120.18,120.2,120.25,120.25,120.24,120.29,120.25,120.26,120.32,120.31,120.36,120.4,120.4,120.39,120.44,120.5,120.53,120.64,120.78,120.94,121,121.05,121.15,121.07,121.18,121.46,121.71,121.71,121.74,121.76,121.76,121.82,121.82,121.82,121.94,121.99,122.18,122.41,122.48,122.52,122.62,122.63,122.74,122.58,122.59,122.61,122.63,122.37,122.36,122.47,122.46,122.45,122.49,122.5,122.37,122.37,122.51,122.51,122.55,122.56,122.72,122.97,123.03,123.05,123.08,123.08,123.12,123.07,123.04,123.06,123.39,124.02,124.05,123.99,124.46,124.46,124.6,124.84,124.84,124.99,125.02) > 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] 85 > (np <- floor(n / par1)) [1] 7 > 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 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 120.05 120.32 121.00 121.82 122.59 122.51 123.04 125.02 [2,] 120.05 120.31 121.05 121.82 122.61 122.51 123.06 NA [3,] 120.08 120.36 121.15 121.94 122.63 122.55 123.39 NA [4,] 120.12 120.40 121.07 121.99 122.37 122.56 124.02 NA [5,] 120.18 120.40 121.18 122.18 122.36 122.72 124.05 NA [6,] 120.20 120.39 121.46 122.41 122.47 122.97 123.99 NA [7,] 120.25 120.44 121.71 122.48 122.46 123.03 124.46 NA [8,] 120.25 120.50 121.71 122.52 122.45 123.05 124.46 NA [9,] 120.24 120.53 121.74 122.62 122.49 123.08 124.60 NA [10,] 120.29 120.64 121.76 122.63 122.50 123.08 124.84 NA [11,] 120.25 120.78 121.76 122.74 122.37 123.12 124.84 NA [12,] 120.26 120.94 121.82 122.58 122.37 123.07 124.99 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 0.00 -0.01 0.05 0.00 0.02 0.00 0.02 NA [2,] 0.03 0.05 0.10 0.12 0.02 0.04 0.33 NA [3,] 0.04 0.04 -0.08 0.05 -0.26 0.01 0.63 NA [4,] 0.06 0.00 0.11 0.19 -0.01 0.16 0.03 NA [5,] 0.02 -0.01 0.28 0.23 0.11 0.25 -0.06 NA [6,] 0.05 0.05 0.25 0.07 -0.01 0.06 0.47 NA [7,] 0.00 0.06 0.00 0.04 -0.01 0.02 0.00 NA [8,] -0.01 0.03 0.03 0.10 0.04 0.03 0.14 NA [9,] 0.05 0.11 0.02 0.01 0.01 0.00 0.24 NA [10,] -0.04 0.14 0.00 0.11 -0.13 0.04 0.00 NA [11,] 0.01 0.16 0.06 -0.16 0.00 -0.05 0.15 NA [12,] 0.06 0.06 0.00 0.01 0.14 -0.03 0.03 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/freestat/rcomp/tmp/1umwa1209302359.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/freestat/rcomp/tmp/22eek1209302359.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/freestat/rcomp/tmp/3onyn1209302359.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/freestat/rcomp/tmp/4dzbi1209302359.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,] 120.050 120.05 120.080 120.120 120.18 120.200 120.250 120.250 120.240 [2,] 120.660 120.68 120.755 120.735 120.79 120.925 121.075 121.105 121.135 [3,] 122.165 121.82 121.940 121.990 122.18 122.410 122.460 122.450 122.490 [4,] 122.815 122.56 122.590 122.465 122.54 122.720 122.755 122.785 122.850 [5,] 125.020 123.06 123.390 124.020 124.05 123.990 124.460 124.460 124.600 [,10] [,11] [,12] [1,] 120.290 120.25 120.260 [2,] 121.200 121.27 121.380 [3,] 122.500 122.37 122.370 [4,] 122.855 122.93 122.825 [5,] 124.840 124.84 124.990 $n [1] 8 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 120.9612 120.6973 120.8442 120.9569 121.1349 121.3381 121.4567 121.4467 [2,] 123.3688 122.9427 123.0358 123.0231 123.2251 123.4819 123.4633 123.4533 [,9] [,10] [,11] [,12] [1,] 121.4658 121.5117 121.3787 121.5071 [2,] 123.5142 123.4883 123.3613 123.2329 $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(120.05, 120.66, 122.165, 122.815, 125.02, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/58n0a1209302359.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.01 0.020 -0.080 -0.010 -0.060 -0.01 -0.01 -0.01 0.00 -0.130 -0.160 [2,] 0.00 0.035 -0.035 0.015 0.005 0.05 0.00 0.03 0.01 -0.020 -0.025 [3,] 0.00 0.050 0.040 0.060 0.110 0.06 0.00 0.03 0.02 0.000 0.010 [4,] 0.02 0.110 0.045 0.135 0.240 0.16 0.03 0.07 0.08 0.075 0.105 [5,] 0.02 0.120 0.050 0.190 0.280 0.25 0.06 0.10 0.11 0.140 0.160 [,12] [1,] -0.030 [2,] 0.005 [3,] 0.030 [4,] 0.060 [5,] 0.140 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.01194368 0.00521121 -0.00777471 -0.01166206 -0.03033821 -0.005690225 [2,] 0.01194368 0.09478879 0.08777471 0.13166206 0.25033821 0.125690225 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.01791552 0.006112645 -0.02180287 -0.05673247 -0.0676339 -0.002845113 [2,] 0.01791552 0.053887355 0.06180287 0.05673247 0.0876339 0.062845113 $out [1] 0.05 0.33 -0.26 0.63 0.47 0.14 0.24 $group [1] 1 2 3 3 6 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.0099999999999909, 0, 0, 0.0199999999999960, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6ecs11209302359.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] [1,] 120.05 120.310 121.000 121.820 122.360 122.510 123.040 125.02 [2,] 120.10 120.375 121.110 121.965 122.370 122.555 123.690 125.02 [3,] 120.22 120.420 121.585 122.445 122.465 123.000 124.255 125.02 [4,] 120.25 120.585 121.750 122.600 122.545 123.075 124.720 125.02 [5,] 120.29 120.780 121.820 122.740 122.630 123.120 124.990 125.02 $n [1] 12 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 120.1516 120.3242 121.2931 122.1554 122.3852 122.7628 123.7852 125.02 [2,] 120.2884 120.5158 121.8769 122.7346 122.5448 123.2372 124.7248 125.02 $out [1] 120.94 $group [1] 2 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(120.05, 120.1, 120.22, 120.25, 120.29, 120.31, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7jftq1209302359.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,] 121.6300 121.8200 121.6000 [2,] 121.8286 122.0775 121.6687 [3,] 122.0812 122.3700 121.8687 [4,] 122.2171 122.4550 122.0100 [5,] 122.2900 122.5000 122.1025 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 121.9039 122.1978 121.7131 [2,] 122.2584 122.5422 122.0244 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(121.63, 121.828571428571, 122.081160714286, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1umwa1209302359.ps tmp/1umwa1209302359.png") > system("convert tmp/22eek1209302359.ps tmp/22eek1209302359.png") > system("convert tmp/3onyn1209302359.ps tmp/3onyn1209302359.png") > system("convert tmp/4dzbi1209302359.ps tmp/4dzbi1209302359.png") > system("convert tmp/58n0a1209302359.ps tmp/58n0a1209302359.png") > system("convert tmp/6ecs11209302359.ps tmp/6ecs11209302359.png") > system("convert tmp/7jftq1209302359.ps tmp/7jftq1209302359.png") > > > proc.time() user system elapsed 2.221 1.497 2.390