R version 2.8.0 (2008-10-20) 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(17.88,18.11,18.16,18.27,18.29,18.35,18.35,18.38,18.41,18.41,18.42,18.43,18.48,18.54,18.65,18.66,18.69,18.72,18.72,18.73,18.84,18.83,18.91,18.91,18.94,18.97,19,19.08,19.18,19.24,19.23,19.25,19.3,19.33,19.35,19.35,19.31,19.47,19.7,19.76,19.9,19.97,20.1,20.26,20.44,20.43,20.57,20.6,20.69,20.93,20.98,21.11,21.14,21.16,21.32,21.32,21.48,21.58,21.74,21.75,21.81,21.89,22.21,22.37,22.47,22.51,22.55,22.61,22.58,22.85,22.93,22.98,23.01,23.11,23.18,23.18,23.21,23.22,23.12,23.15,23.16,23.21,23.21,23.22) > 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] 84 > (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] 7 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 17.88 18.48 18.94 19.31 20.69 21.81 23.01 NA [2,] 18.11 18.54 18.97 19.47 20.93 21.89 23.11 NA [3,] 18.16 18.65 19.00 19.70 20.98 22.21 23.18 NA [4,] 18.27 18.66 19.08 19.76 21.11 22.37 23.18 NA [5,] 18.29 18.69 19.18 19.90 21.14 22.47 23.21 NA [6,] 18.35 18.72 19.24 19.97 21.16 22.51 23.22 NA [7,] 18.35 18.72 19.23 20.10 21.32 22.55 23.12 NA [8,] 18.38 18.73 19.25 20.26 21.32 22.61 23.15 NA [9,] 18.41 18.84 19.30 20.44 21.48 22.58 23.16 NA [10,] 18.41 18.83 19.33 20.43 21.58 22.85 23.21 NA [11,] 18.42 18.91 19.35 20.57 21.74 22.93 23.21 NA [12,] 18.43 18.91 19.35 20.60 21.75 22.98 23.22 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 0.23 0.06 0.03 0.16 0.24 0.08 0.10 NA [2,] 0.05 0.11 0.03 0.23 0.05 0.32 0.07 NA [3,] 0.11 0.01 0.08 0.06 0.13 0.16 0.00 NA [4,] 0.02 0.03 0.10 0.14 0.03 0.10 0.03 NA [5,] 0.06 0.03 0.06 0.07 0.02 0.04 0.01 NA [6,] 0.00 0.00 -0.01 0.13 0.16 0.04 -0.10 NA [7,] 0.03 0.01 0.02 0.16 0.00 0.06 0.03 NA [8,] 0.03 0.11 0.05 0.18 0.16 -0.03 0.01 NA [9,] 0.00 -0.01 0.03 -0.01 0.10 0.27 0.05 NA [10,] 0.01 0.08 0.02 0.14 0.16 0.08 0.00 NA [11,] 0.01 0.00 0.00 0.03 0.01 0.05 0.01 NA [12,] 0.05 0.03 -0.04 0.09 0.06 0.03 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/freestat/rcomp/tmp/1au6n1294088170.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/www/html/freestat/rcomp/tmp/2lu7c1294088170.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/www/html/freestat/rcomp/tmp/3g6zv1294088170.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/www/html/freestat/rcomp/tmp/476dg1294088170.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] [,11] [1,] 17.88 18.110 18.160 18.27 18.290 18.350 18.350 18.380 18.41 18.410 18.420 [2,] 18.71 18.755 18.825 18.87 18.935 18.980 18.975 18.990 19.07 19.080 19.130 [3,] 19.31 19.470 19.700 19.76 19.900 19.970 20.100 20.260 20.44 20.430 20.570 [4,] 21.25 21.410 21.595 21.74 21.805 21.835 21.935 21.965 22.03 22.215 22.335 [5,] 23.01 23.110 23.180 23.18 23.210 23.220 23.120 23.150 23.16 23.210 23.210 [,12] [1,] 18.430 [2,] 19.130 [3,] 20.600 [4,] 22.365 [5,] 23.220 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 17.79315 17.88448 18.0458 18.04608 18.18608 18.26504 18.33234 18.48338 [2,] 20.82685 21.05552 21.3542 21.47392 21.61392 21.67496 21.86766 22.03662 [,9] [,10] [,11] [,12] [1,] 18.67234 18.55783 18.65603 18.66811 [2,] 22.20766 22.30217 22.48397 22.53189 $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(17.88, 18.71, 19.31, 21.25, 23.01, 18.11, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5xrr31294088170.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] [,12] [1,] 0.030 0.03 0.000 0.02 0.010 -0.100 0.000 -0.030 -0.010 0.000 0.000 0.03 [2,] 0.070 0.05 0.035 0.03 0.025 -0.005 0.015 0.020 -0.005 0.015 0.005 0.03 [3,] 0.100 0.07 0.080 0.03 0.040 0.000 0.030 0.050 0.030 0.080 0.010 0.04 [4,] 0.195 0.17 0.120 0.10 0.060 0.085 0.045 0.135 0.075 0.110 0.020 0.06 [5,] 0.240 0.32 0.160 0.14 0.070 0.160 0.060 0.180 0.100 0.160 0.030 0.09 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.02535202 -0.001662064 0.02923937 -0.01180287 0.01909856 -0.05374655 [2,] 0.17464798 0.141662064 0.13076063 0.07180287 0.06090144 0.05374655 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.01208448 -0.01867614 -0.01777471 0.02326753 0.001042242 0.02064903 [2,] 0.04791552 0.11867614 0.07777471 0.13673247 0.018957758 0.05935097 $out [1] 0.16 0.27 0.05 -0.04 $group [1] 7 9 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(0.0299999999999976, 0.0700000000000003, 0.0999999999999979, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6y96b1294088170.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] [1,] 18.110 18.480 18.940 19.310 20.690 21.81 23.110 NA [2,] 18.215 18.655 19.040 19.730 21.045 22.29 23.135 NA [3,] 18.350 18.720 19.235 20.035 21.240 22.53 23.180 NA [4,] 18.410 18.835 19.315 20.435 21.530 22.73 23.210 NA [5,] 18.430 18.910 19.350 20.600 21.750 22.98 23.220 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 18.26106 18.6379 19.10957 19.71344 21.01879 22.32931 23.14579 NA [2,] 18.43894 18.8021 19.36043 20.35656 21.46121 22.73069 23.21421 NA $out [1] 17.88 23.01 $group [1] 1 7 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(18.11, 18.215, 18.35, 18.41, 18.43, 18.48, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7epkt1294088170.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,] 20.01714 19.310 19.98000 [2,] 20.30786 19.730 20.25750 [3,] 20.46857 20.035 20.43125 [4,] 20.63214 20.435 20.59875 [5,] 20.74857 20.600 20.74750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 20.32066 19.71344 20.27560 [2,] 20.61648 20.35656 20.58690 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(20.0171428571429, 20.3078571428571, 20.4685714285714, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1au6n1294088170.ps tmp/1au6n1294088170.png",intern=TRUE)) character(0) > try(system("convert tmp/2lu7c1294088170.ps tmp/2lu7c1294088170.png",intern=TRUE)) character(0) > try(system("convert tmp/3g6zv1294088170.ps tmp/3g6zv1294088170.png",intern=TRUE)) character(0) > try(system("convert tmp/476dg1294088170.ps tmp/476dg1294088170.png",intern=TRUE)) character(0) > try(system("convert tmp/5xrr31294088170.ps tmp/5xrr31294088170.png",intern=TRUE)) character(0) > try(system("convert tmp/6y96b1294088170.ps tmp/6y96b1294088170.png",intern=TRUE)) character(0) > try(system("convert tmp/7epkt1294088170.ps tmp/7epkt1294088170.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.946 1.463 2.359