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(29.90,28.77,15.64,23.73,25.65,21.81,28.97,24.29,25.33,28.84,19.99,19.75,22.70,23.28,24.15,20.38,27.75,27.31,25.61,22.64,26.05,28.07,21.02,25.00,17.93,35.45,17.70,28.53,26.55,26.51,30.78,26.83,27.49,25.89,20.44,19.79,18.14,27.98,35.90,34.38,21.58,21.53,31.14,28.25,25.16,20.51,30.05,20.17,32.37,22.46,25.40,19.82,18.14,20.10,20.25,19.73,24.74,26.17,20.14,31.71,26.66,20.75,20.01,26.67,23.91,26.81,29.31,31.76,22.99,23.94,27.04,20.28,23.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] 73 > (np <- floor(n / par1)) [1] 6 > 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 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 29.90 22.70 17.93 18.14 32.37 26.66 23.32 [2,] 28.77 23.28 35.45 27.98 22.46 20.75 NA [3,] 15.64 24.15 17.70 35.90 25.40 20.01 NA [4,] 23.73 20.38 28.53 34.38 19.82 26.67 NA [5,] 25.65 27.75 26.55 21.58 18.14 23.91 NA [6,] 21.81 27.31 26.51 21.53 20.10 26.81 NA [7,] 28.97 25.61 30.78 31.14 20.25 29.31 NA [8,] 24.29 22.64 26.83 28.25 19.73 31.76 NA [9,] 25.33 26.05 27.49 25.16 24.74 22.99 NA [10,] 28.84 28.07 25.89 20.51 26.17 23.94 NA [11,] 19.99 21.02 20.44 30.05 20.14 27.04 NA [12,] 19.75 25.00 19.79 20.17 31.71 20.28 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.13 0.58 17.52 9.84 -9.91 -5.91 NA [2,] -13.13 0.87 -17.75 7.92 2.94 -0.74 NA [3,] 8.09 -3.77 10.83 -1.52 -5.58 6.66 NA [4,] 1.92 7.37 -1.98 -12.80 -1.68 -2.76 NA [5,] -3.84 -0.44 -0.04 -0.05 1.96 2.90 NA [6,] 7.16 -1.70 4.27 9.61 0.15 2.50 NA [7,] -4.68 -2.97 -3.95 -2.89 -0.52 2.45 NA [8,] 1.04 3.41 0.66 -3.09 5.01 -8.77 NA [9,] 3.51 2.02 -1.60 -4.65 1.43 0.95 NA [10,] -8.85 -7.05 -5.45 9.54 -6.03 3.10 NA [11,] -0.24 3.98 -0.65 -9.88 11.57 -6.76 NA [12,] 2.95 -7.07 -1.65 12.20 -5.05 3.04 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/1iht71209412532.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/2it1u1209412532.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/32hz31209412532.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/4tk3d1209412532.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,] 17.93 20.75 15.64 19.82 18.14 20.10 20.25 19.73 22.990 20.51 19.99 19.750 [2,] 20.42 22.46 17.70 20.38 21.58 21.53 25.61 22.64 24.740 23.94 20.14 19.790 [3,] 23.32 25.63 22.08 25.20 24.78 24.16 29.14 25.56 25.245 26.03 20.73 20.225 [4,] 28.28 28.77 25.40 28.53 26.55 26.81 30.78 28.25 26.050 28.07 27.04 25.000 [5,] 32.37 35.45 35.90 34.38 27.75 27.31 31.14 31.76 27.490 28.84 30.05 31.710 $n [1] 7 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 18.62613 21.55985 17.11325 19.94299 21.57419 20.75423 25.80518 21.94137 [2,] 28.01387 29.70015 27.04675 30.45701 27.98581 27.56577 32.47482 29.17863 [,9] [,10] [,11] [,12] [1,] 24.40001 23.36602 16.27928 16.86438 [2,] 26.08999 28.69398 25.18072 23.58562 $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.93, 20.42, 23.32, 28.28, 32.37, 20.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5x4ai1209412532.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,] -9.910 -17.750 -5.58 -2.76 -3.840 -1.700 -4.68 -8.77 -4.65 -8.85 -9.880 [2,] -5.910 -13.130 -3.77 -2.76 -0.440 0.150 -3.95 -3.09 -1.60 -7.05 -6.760 [3,] -0.275 0.065 2.57 -1.83 -0.045 3.385 -2.93 0.85 1.19 -5.74 -0.445 [4,] 9.840 2.940 8.09 1.92 1.960 7.160 -0.52 3.41 2.02 3.10 3.980 [5,] 17.520 7.920 10.83 7.37 2.900 9.610 2.45 5.01 3.51 9.54 11.570 [,12] [1,] -7.07 [2,] -5.05 [3,] 0.65 [4,] 3.04 [5,] 12.20 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -10.434259 -10.30067 -5.080083 -4.848751 -1.593078 -1.136676 -5.1424608 [2,] 9.884259 10.43067 10.220083 1.188751 1.503078 7.906676 -0.7175392 [,8] [,9] [,10] [,11] [,12] [1,] -3.34271 -1.145017 -12.2870778 -7.372647 -4.568311 [2,] 5.04271 3.525017 0.8070778 6.482647 5.868311 $out [1] -12.8 $group [1] 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-9.91, -5.91, -0.274999999999999, 9.84, 17.52, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6elx81209412532.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] [1,] 15.640 20.380 17.700 18.140 18.140 20.010 23.32 [2,] 20.900 22.670 20.115 21.020 19.960 21.870 23.32 [3,] 24.810 24.575 26.530 26.570 21.355 25.300 23.32 [4,] 28.805 26.680 28.010 30.595 25.785 26.925 23.32 [5,] 29.900 28.070 35.450 35.900 32.370 31.760 23.32 $n [1] 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 21.20448 22.74601 22.92904 22.20278 18.69818 22.99438 23.32 [2,] 28.41552 26.40399 30.13096 30.93722 24.01182 27.60562 23.32 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(15.64, 20.9, 24.81, 28.805, 29.9, 20.38, 22.67, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7dcow1209412532.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,] 22.78333 20.225 21.68250 [2,] 23.53167 22.700 23.52125 [3,] 24.86238 24.990 24.49563 [4,] 25.58417 25.595 25.54625 [5,] 27.67667 29.140 28.43125 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 23.92622 23.66957 23.57201 [2,] 25.79854 26.31043 25.41924 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(22.7833333333333, 23.5316666666667, 24.8623809523810, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1iht71209412532.ps tmp/1iht71209412532.png") > system("convert tmp/2it1u1209412532.ps tmp/2it1u1209412532.png") > system("convert tmp/32hz31209412532.ps tmp/32hz31209412532.png") > system("convert tmp/4tk3d1209412532.ps tmp/4tk3d1209412532.png") > system("convert tmp/5x4ai1209412532.ps tmp/5x4ai1209412532.png") > system("convert tmp/6elx81209412532.ps tmp/6elx81209412532.png") > system("convert tmp/7dcow1209412532.ps tmp/7dcow1209412532.png") > > > proc.time() user system elapsed 2.124 1.391 2.282