R version 2.11.1 (2010-05-31) Copyright (C) 2010 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(132,131.4,132.7,130.9,126,109.7,68.3,70.6,75.3,74.1,74.9,74,74.2,76,76.2,74.9,74.1,76.5,57.8,59.2,57.3,57.5,60.4,59.9,59.9,60,60.2,65.4,62.4,78.8,65.6,64.4,67.4,65.3,66.7,66.8,69.4,71.7,77.1,81.1,82.1,92.1,77.1,78.2,77.7,77.3,78.5,78.8,78.7,79.8,82.2,84,81.7,77.6,64.3,72.6,73.8,73.8,70.1,70,72.3,72.1,73.3,79.1,77,76.1,66.4,72.7,73.2,70.7,73.6,74.2,72.6,73.6,79.1,79.6,78,85.4,82,91.9,89.4,92.1,93.8,93.6,95.6,99.9,103.7,99.2,93.7,93.5,80.7,91.8,105.8,111.3,110.3,109.4,111.4,111.6,111.8,106.6,104.3,105.5,98.5,108.5,106,101.8,101.3,92.4,88.9,84.9,86.4,90.7,86.8,90.6,88.3,95.4,93.6,91.3,91.3,89.3) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 132.0 74.2 59.9 69.4 78.7 72.3 72.6 95.6 111.4 88.9 NA [2,] 131.4 76.0 60.0 71.7 79.8 72.1 73.6 99.9 111.6 84.9 NA [3,] 132.7 76.2 60.2 77.1 82.2 73.3 79.1 103.7 111.8 86.4 NA [4,] 130.9 74.9 65.4 81.1 84.0 79.1 79.6 99.2 106.6 90.7 NA [5,] 126.0 74.1 62.4 82.1 81.7 77.0 78.0 93.7 104.3 86.8 NA [6,] 109.7 76.5 78.8 92.1 77.6 76.1 85.4 93.5 105.5 90.6 NA [7,] 68.3 57.8 65.6 77.1 64.3 66.4 82.0 80.7 98.5 88.3 NA [8,] 70.6 59.2 64.4 78.2 72.6 72.7 91.9 91.8 108.5 95.4 NA [9,] 75.3 57.3 67.4 77.7 73.8 73.2 89.4 105.8 106.0 93.6 NA [10,] 74.1 57.5 65.3 77.3 73.8 70.7 92.1 111.3 101.8 91.3 NA [11,] 74.9 60.4 66.7 78.5 70.1 73.6 93.8 110.3 101.3 91.3 NA [12,] 74.0 59.9 66.8 78.8 70.0 74.2 93.6 109.4 92.4 89.3 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -0.6 1.8 0.1 2.3 1.1 -0.2 1.0 4.3 0.2 -4.0 NA [2,] 1.3 0.2 0.2 5.4 2.4 1.2 5.5 3.8 0.2 1.5 NA [3,] -1.8 -1.3 5.2 4.0 1.8 5.8 0.5 -4.5 -5.2 4.3 NA [4,] -4.9 -0.8 -3.0 1.0 -2.3 -2.1 -1.6 -5.5 -2.3 -3.9 NA [5,] -16.3 2.4 16.4 10.0 -4.1 -0.9 7.4 -0.2 1.2 3.8 NA [6,] -41.4 -18.7 -13.2 -15.0 -13.3 -9.7 -3.4 -12.8 -7.0 -2.3 NA [7,] 2.3 1.4 -1.2 1.1 8.3 6.3 9.9 11.1 10.0 7.1 NA [8,] 4.7 -1.9 3.0 -0.5 1.2 0.5 -2.5 14.0 -2.5 -1.8 NA [9,] -1.2 0.2 -2.1 -0.4 0.0 -2.5 2.7 5.5 -4.2 -2.3 NA [10,] 0.8 2.9 1.4 1.2 -3.7 2.9 1.7 -1.0 -0.5 0.0 NA [11,] -0.9 -0.5 0.1 0.3 -0.1 0.6 -0.2 -0.9 -8.9 -2.0 NA [12,] 0.2 0.0 2.6 -0.1 2.3 -1.6 2.0 2.0 -3.5 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/rcomp/tmp/1hlyk1290791942.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/rcomp/tmp/2hlyk1290791942.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/rcomp/tmp/3hlyk1290791942.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/rcomp/tmp/49dgn1290791942.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,] 59.90 60.0 60.20 65.40 62.4 76.1 57.8 59.20 57.3 57.5 60.4 59.9 [2,] 72.30 72.1 76.20 79.10 77.0 77.6 65.6 70.60 73.2 70.7 70.1 70.0 [3,] 76.45 77.9 80.65 82.55 81.9 88.0 72.7 75.45 76.5 75.7 76.7 76.5 [4,] 95.60 99.9 103.70 99.20 93.7 93.5 82.0 91.90 93.6 92.1 93.8 92.4 [5,] 111.40 131.4 132.70 106.60 104.3 109.7 98.5 108.50 106.0 111.3 110.3 109.4 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 64.80839 64.01001 66.9099 72.50724 73.55601 80.05573 64.5059 64.80767 [2,] 88.09161 91.78999 94.3901 92.59276 90.24399 95.94427 80.8941 86.09233 [,9] [,10] [,11] [,12] [1,] 66.30735 65.00771 64.85854 65.30807 [2,] 86.69265 86.39229 88.54146 87.69193 $out [1] 132.0 130.9 126.0 $group [1] 1 4 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(59.9, 72.3, 76.45, 95.6, 111.4, 60, 72.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/59dgn1290791942.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] [,12] [1,] -0.6 0.2 -5.20 -5.5 -4.1 -18.7 -1.2 -2.5 -4.2 -3.7 -2.00 -1.6 [2,] -0.2 0.2 -1.80 -3.9 -0.9 -15.0 1.4 -1.9 -2.3 -0.5 -0.90 -0.1 [3,] 0.6 1.4 1.15 -2.3 1.8 -13.0 6.7 0.0 -0.8 1.0 -0.35 0.2 [4,] 1.8 3.8 4.30 -1.6 7.4 -7.0 9.9 3.0 0.2 1.7 0.10 2.0 [5,] 4.3 5.5 5.80 1.0 16.4 -2.3 11.1 4.7 2.7 2.9 0.60 2.6 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -0.3992797 -0.3987035 -1.897803 -3.449172 -2.347011 -16.997119 2.453061 [2,] 1.5992797 3.1987035 4.197803 -1.150828 5.947011 -9.002881 10.946939 [,8] [,9] [,10] [,11] [,12] [1,] -2.448235 -2.0490997 -0.09920771 -0.8496399 -0.906 [2,] 2.448235 0.4490997 2.09920771 0.1496399 1.306 $out [1] -4.0 -16.3 -41.4 14.0 5.5 -8.9 -3.5 $group [1] 1 5 6 8 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.599999999999994, -0.200000000000003, 0.599999999999994, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6kmxp1290791942.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] [,10] [,11] [1,] 68.30 57.30 59.90 77.10 64.30 70.70 72.60 80.70 92.40 84.90 NA [2,] 74.05 58.50 61.30 77.10 71.35 72.20 78.55 93.60 101.55 87.55 NA [3,] 92.50 67.25 65.35 77.95 75.70 73.25 83.70 99.55 105.75 89.95 NA [4,] 131.15 75.45 66.75 79.95 80.75 75.15 92.00 107.60 109.95 91.30 NA [5,] 132.70 76.50 67.40 82.10 84.00 79.10 93.80 111.30 111.80 95.40 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 66.45631 59.51899 62.86422 76.6501 71.4126 71.90449 77.56536 93.1645 [2,] 118.54369 74.98101 67.83578 79.2499 79.9874 74.59551 89.83464 105.9355 [,9] [,10] [,11] [1,] 101.9187 88.2396 NA [2,] 109.5813 91.6604 NA $out [1] 78.8 69.4 71.7 92.1 66.4 $group [1] 3 4 4 4 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(68.3, 74.05, 92.5, 131.15, 132.7, 57.3, 58.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7kmxp1290791942.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,] 74.900 72.700 81.31250 [2,] 81.180 76.075 81.58750 [3,] 83.795 76.600 83.05000 [4,] 87.440 81.275 85.06875 [5,] 89.150 88.000 88.15000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 80.93977 74.22825 81.46218 [2,] 86.65023 78.97175 84.63782 $out [1] 73.7375 $group [1] 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(74.9, 81.18, 83.795, 87.44, 89.15, 72.7, 76.075, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1hlyk1290791942.ps tmp/1hlyk1290791942.png",intern=TRUE)) character(0) > try(system("convert tmp/2hlyk1290791942.ps tmp/2hlyk1290791942.png",intern=TRUE)) character(0) > try(system("convert tmp/3hlyk1290791942.ps tmp/3hlyk1290791942.png",intern=TRUE)) character(0) > try(system("convert tmp/49dgn1290791942.ps tmp/49dgn1290791942.png",intern=TRUE)) character(0) > try(system("convert tmp/59dgn1290791942.ps tmp/59dgn1290791942.png",intern=TRUE)) character(0) > try(system("convert tmp/6kmxp1290791942.ps tmp/6kmxp1290791942.png",intern=TRUE)) character(0) > try(system("convert tmp/7kmxp1290791942.ps tmp/7kmxp1290791942.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.00 0.62 2.58