R version 2.7.0 (2008-04-22) 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. 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(127.96,127.47,126.47,125.75,125.42,125.14,125.15,125.51,125.63,126.22,126.88,127.96,128.74,129.6,131.2,132.72,134.67,135.94,136.39,136.74,137.2,137.36,138.63,141.07,143.32,147.91,152.56,151.61,156.56,157.45,158.13,159.18,159.47,159.79,161.65,162.77,163.48,166.16,163.86,162.12,149.08,145.32,141.21,134.68,133.65,139.17,138.61,144.96,157.99,167.18,174.48,182.77,190.00,189.70,188.90,198.28,201.18,204.14,221.02,221.12,220.68) > 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] 61 > (np <- floor(n / par1)) [1] 5 > 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 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 127.96 128.74 143.32 163.48 157.99 220.68 [2,] 127.47 129.60 147.91 166.16 167.18 NA [3,] 126.47 131.20 152.56 163.86 174.48 NA [4,] 125.75 132.72 151.61 162.12 182.77 NA [5,] 125.42 134.67 156.56 149.08 190.00 NA [6,] 125.14 135.94 157.45 145.32 189.70 NA [7,] 125.15 136.39 158.13 141.21 188.90 NA [8,] 125.51 136.74 159.18 134.68 198.28 NA [9,] 125.63 137.20 159.47 133.65 201.18 NA [10,] 126.22 137.36 159.79 139.17 204.14 NA [11,] 126.88 138.63 161.65 138.61 221.02 NA [12,] 127.96 141.07 162.77 144.96 221.12 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.49 0.86 4.59 2.68 9.19 NA [2,] -1.00 1.60 4.65 -2.30 7.30 NA [3,] -0.72 1.52 -0.95 -1.74 8.29 NA [4,] -0.33 1.95 4.95 -13.04 7.23 NA [5,] -0.28 1.27 0.89 -3.76 -0.30 NA [6,] 0.01 0.45 0.68 -4.11 -0.80 NA [7,] 0.36 0.35 1.05 -6.53 9.38 NA [8,] 0.12 0.46 0.29 -1.03 2.90 NA [9,] 0.59 0.16 0.32 5.52 2.96 NA [10,] 0.66 1.27 1.86 -0.56 16.88 NA [11,] 1.08 2.44 1.12 6.35 0.10 NA [12,] 0.78 2.25 0.71 13.03 -0.44 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/103bi1229601445.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/2h8hg1229601445.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/3jwu41229601445.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/44hfj1229601445.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] [1,] 127.960 127.47 126.47 125.75 125.42 125.14 125.15 125.51 125.63 126.22 [2,] 128.740 129.60 131.20 132.72 134.67 135.94 136.39 134.68 133.65 137.36 [3,] 150.655 147.91 152.56 151.61 149.08 145.32 141.21 136.74 137.20 139.17 [4,] 163.480 166.16 163.86 162.12 156.56 157.45 158.13 159.18 159.47 159.79 [5,] 163.480 167.18 174.48 182.77 156.56 189.70 188.90 159.18 159.47 159.79 [,11] [,12] [1,] 126.88 127.96 [2,] 138.61 141.07 [3,] 138.63 144.96 [4,] 161.65 162.77 [5,] 161.65 162.77 $n [1] 6 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 128.2466 122.0768 129.4825 130.8360 133.6126 130.1211 125.8486 119.4284 [2,] 173.0634 173.7432 175.6375 172.3840 164.5474 160.5189 156.5714 154.0516 [,9] [,10] [,11] [,12] [1,] 118.9557 123.3210 122.35 129.6268 [2,] 155.4443 155.0190 154.91 160.2932 $out [1] 220.68 190.00 198.28 201.18 204.14 221.02 221.12 $group [1] 1 5 8 9 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(127.96, 128.74, 150.655, 163.48, 163.48, 127.47, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5idlq1229601445.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.49 -2.30 -1.74 -0.33 -0.30 -0.80 0.35 0.12 0.16 -0.56 0.10 -0.44 [2,] 0.86 -1.00 -0.95 -0.33 -0.30 -0.80 0.35 0.12 0.32 0.66 1.08 0.71 [3,] 2.68 1.60 -0.72 1.95 -0.28 0.01 0.36 0.29 0.59 1.27 1.12 0.78 [4,] 4.59 4.65 1.52 4.95 0.89 0.45 1.05 0.46 2.96 1.86 2.44 2.25 [5,] 9.19 7.30 1.52 7.23 1.27 0.68 1.05 0.46 5.52 1.86 2.44 2.25 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.0443914 -2.392276 -2.465296 -1.780835 -1.120851 -0.8732469 -0.1346182 [2,] 5.3156086 5.592276 1.025296 5.680835 0.560851 0.8932469 0.8546182 [,8] [,9] [,10] [,11] [,12] [1,] 0.04975686 -1.275417 0.422083 0.1590274 -0.3081601 [2,] 0.53024314 2.455417 2.117917 2.0809726 1.8681601 $out [1] 8.29 -13.04 -3.76 -4.11 -6.53 9.38 -1.03 2.90 16.88 6.35 [11] 13.03 $group [1] 3 4 5 6 7 7 8 8 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.489999999999995, 0.859999999999985, 2.68000000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/68raw1229601445.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] [1,] 125.140 128.740 143.320 133.65 157.990 220.68 [2,] 125.465 131.960 152.085 138.89 178.625 220.68 [3,] 125.985 136.165 157.790 145.14 189.850 220.68 [4,] 127.175 137.280 159.630 162.80 202.660 220.68 [5,] 127.960 141.070 162.770 166.16 221.120 220.68 $n [1] 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 125.2051 133.7385 154.3487 134.2345 178.8875 220.68 [2,] 126.7649 138.5915 161.2313 156.0455 200.8125 220.68 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(125.14, 125.465, 125.985, 127.175, 127.96, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7kovw1229601445.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,] 147.6640 136.7400 145.6150 [2,] 150.3330 138.9000 146.8125 [3,] 151.0700 145.1400 147.3400 [4,] 155.1822 149.8675 148.2275 [5,] 159.5760 152.5600 150.1300 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 148.8583 140.1376 146.6946 [2,] 153.2817 150.1424 147.9854 $out [1] 151.92 $group [1] 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(147.664, 150.333, 151.07, 155.182166666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/103bi1229601445.ps tmp/103bi1229601445.png") > system("convert tmp/2h8hg1229601445.ps tmp/2h8hg1229601445.png") > system("convert tmp/3jwu41229601445.ps tmp/3jwu41229601445.png") > system("convert tmp/44hfj1229601445.ps tmp/44hfj1229601445.png") > system("convert tmp/5idlq1229601445.ps tmp/5idlq1229601445.png") > system("convert tmp/68raw1229601445.ps tmp/68raw1229601445.png") > system("convert tmp/7kovw1229601445.ps tmp/7kovw1229601445.png") > > > proc.time() user system elapsed 3.508 2.523 3.828