R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(76.3,71,85.9,109,83.4,101.3,81.6,84.8,88,88.4,85.2,139.5,85.3,79.4,131,91.4,88.2,113.2,81.9,89,95.7,93.2,98,140.5,93,82.1,86.3,84.8,85.7,110.8,87.4,100.4,100.4,99.4,108,161.7,109.3,99.9,110.3,97.5,102.4,122.3,104.4,122.7,127.3,128.4,125.3,187.1,131.4,125.6,142.9,116.6,129.7,155.4,138.7,167.7,155.8,157.1,159.9,244.6,154.3,143.2,147.3,155.1,152.2,166.5,161.2,180.1,169.2,164.5,166.2,267.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] 72 > (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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 76.3 85.3 93.0 109.3 131.4 154.3 NA [2,] 71.0 79.4 82.1 99.9 125.6 143.2 NA [3,] 85.9 131.0 86.3 110.3 142.9 147.3 NA [4,] 109.0 91.4 84.8 97.5 116.6 155.1 NA [5,] 83.4 88.2 85.7 102.4 129.7 152.2 NA [6,] 101.3 113.2 110.8 122.3 155.4 166.5 NA [7,] 81.6 81.9 87.4 104.4 138.7 161.2 NA [8,] 84.8 89.0 100.4 122.7 167.7 180.1 NA [9,] 88.0 95.7 100.4 127.3 155.8 169.2 NA [10,] 88.4 93.2 99.4 128.4 157.1 164.5 NA [11,] 85.2 98.0 108.0 125.3 159.9 166.2 NA [12,] 139.5 140.5 161.7 187.1 244.6 267.3 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -5.3 -5.9 -10.9 -9.4 -5.8 -11.1 NA [2,] 14.9 51.6 4.2 10.4 17.3 4.1 NA [3,] 23.1 -39.6 -1.5 -12.8 -26.3 7.8 NA [4,] -25.6 -3.2 0.9 4.9 13.1 -2.9 NA [5,] 17.9 25.0 25.1 19.9 25.7 14.3 NA [6,] -19.7 -31.3 -23.4 -17.9 -16.7 -5.3 NA [7,] 3.2 7.1 13.0 18.3 29.0 18.9 NA [8,] 3.2 6.7 0.0 4.6 -11.9 -10.9 NA [9,] 0.4 -2.5 -1.0 1.1 1.3 -4.7 NA [10,] -3.2 4.8 8.6 -3.1 2.8 1.7 NA [11,] 54.3 42.5 53.7 61.8 84.7 101.1 NA [12,] -54.2 -47.5 -52.4 -55.7 -90.3 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/1aims1259598714.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/2wacu1259598714.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/34d041259598714.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/4w59y1259598714.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,] 76.30 71.0 85.90 84.80 83.4 101.30 81.6 84.80 88.00 88.4 85.20 [2,] 85.30 79.4 86.30 91.40 85.7 110.80 81.9 89.00 95.70 93.2 98.00 [3,] 101.15 91.0 120.65 103.25 95.3 117.75 95.9 111.55 113.85 113.9 116.65 [4,] 131.40 125.6 142.90 116.60 129.7 155.40 138.7 167.70 155.80 157.1 159.90 [5,] 154.30 143.2 147.30 116.60 152.2 166.50 161.2 180.10 169.20 164.5 166.20 [,12] [1,] 139.5 [2,] 140.5 [3,] 174.4 [4,] 244.6 [5,] 267.3 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 71.41401 61.19951 84.14117 86.99519 66.91858 88.98156 59.26217 [2,] 130.88599 120.80049 157.15883 119.50481 123.68142 146.51844 132.53783 [,8] [,9] [,10] [,11] [,12] [1,] 60.78596 75.08356 72.68244 76.7225 107.2521 [2,] 162.31404 152.61644 155.11756 156.5775 241.5479 $out [1] 155.1 $group [1] 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(76.3, 85.3, 101.15, 131.4, 154.3, 71, 79.4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/52z311259598714.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,] -11.10 4.10 -39.60 -3.2 14.30 -31.3 3.20 -11.9 -4.7 -3.20 42.50 -55.7 [2,] -10.90 4.20 -26.30 -3.2 17.90 -23.4 7.10 -10.9 -2.5 -3.10 53.70 -55.7 [3,] -7.65 12.65 -7.15 -1.0 22.45 -18.8 15.65 1.6 -0.3 2.25 58.05 -54.2 [4,] -5.80 17.30 7.80 4.9 25.10 -16.7 18.90 4.6 1.1 4.80 84.70 -52.4 [5,] -5.30 17.30 23.10 13.1 25.70 -16.7 29.00 6.7 1.3 8.60 101.10 -47.5 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -10.939665 4.200077 -29.1456 -6.224762 17.80577 -23.12172 8.038619 [2,] -4.360335 21.099923 14.8456 4.224762 27.09423 -14.47828 23.261381 [,8] [,9] [,10] [,11] [,12] [1,] -8.398 -2.622116 -2.845755 38.054 -56.53177 [2,] 11.598 2.022116 7.345755 78.046 -51.86823 $out [1] 51.6 -25.6 -5.3 -90.3 $group [1] 2 4 6 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-11.1000000000000, -10.9, -7.64999999999999, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6n4in1259598714.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,] 71.00 79.40 82.1 97.5 116.60 143.20 NA [2,] 82.50 86.75 86.0 103.4 130.55 153.25 NA [3,] 85.55 92.30 96.2 116.3 149.15 162.85 NA [4,] 94.85 105.60 104.2 126.3 158.50 167.85 NA [5,] 109.00 131.00 110.8 128.4 167.70 180.10 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 79.91708 83.70239 87.89886 105.8552 136.4018 156.1908 NA [2,] 91.18292 100.89761 104.50114 126.7448 161.8982 169.5092 NA $out [1] 139.5 140.5 161.7 187.1 244.6 267.3 $group [1] 1 2 3 4 5 6 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(71, 82.5, 85.55, 94.85, 109, 79.4, 86.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/76pye1259598714.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,] 100.2000 91.000 99.6250 [2,] 108.6667 98.525 105.5750 [3,] 119.5583 112.700 119.2250 [4,] 123.9417 117.200 125.0125 [5,] 128.2500 120.650 129.2625 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 112.5913 104.1822 110.3594 [2,] 126.5254 121.2178 128.0906 $out [1] 190.1167 174.4000 188.0125 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(100.2, 108.666666666667, 119.558333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1aims1259598714.ps tmp/1aims1259598714.png") > system("convert tmp/2wacu1259598714.ps tmp/2wacu1259598714.png") > system("convert tmp/34d041259598714.ps tmp/34d041259598714.png") > system("convert tmp/4w59y1259598714.ps tmp/4w59y1259598714.png") > system("convert tmp/52z311259598714.ps tmp/52z311259598714.png") > system("convert tmp/6n4in1259598714.ps tmp/6n4in1259598714.png") > system("convert tmp/76pye1259598714.ps tmp/76pye1259598714.png") > > > proc.time() user system elapsed 1.194 0.940 1.606