R version 2.7.2 (2008-08-25) 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(118.63,121.83,119.97,124.98,129.99,126.60,121.71,119.28,122.63,116.74,114.23,113.23,112.75,113.54,115.30,121.05,119.51,116.78,117.17,117.50,119.65,120.97,117.18,116.87,119.46,122.52,124.10,118.39,113.10,113.94,114.58,118.79,120.44,118.37,118.44,117.93,117.76,118.29,121.11,124.86,131.17,130.16,131.76,134.70,135.32,140.23,136.31,131.62,128.90,133.89,138.21,146.12,144.69,149.18,156.60,158.87,164.85,162.89,153.31,150.91) > 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] 60 > (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] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 118.63 112.75 119.46 117.76 128.90 NA [2,] 121.83 113.54 122.52 118.29 133.89 NA [3,] 119.97 115.30 124.10 121.11 138.21 NA [4,] 124.98 121.05 118.39 124.86 146.12 NA [5,] 129.99 119.51 113.10 131.17 144.69 NA [6,] 126.60 116.78 113.94 130.16 149.18 NA [7,] 121.71 117.17 114.58 131.76 156.60 NA [8,] 119.28 117.50 118.79 134.70 158.87 NA [9,] 122.63 119.65 120.44 135.32 164.85 NA [10,] 116.74 120.97 118.37 140.23 162.89 NA [11,] 114.23 117.18 118.44 136.31 153.31 NA [12,] 113.23 116.87 117.93 131.62 150.91 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 3.20 0.79 3.06 0.53 4.99 NA [2,] -1.86 1.76 1.58 2.82 4.32 NA [3,] 5.01 5.75 -5.71 3.75 7.91 NA [4,] 5.01 -1.54 -5.29 6.31 -1.43 NA [5,] -3.39 -2.73 0.84 -1.01 4.49 NA [6,] -4.89 0.39 0.64 1.60 7.42 NA [7,] -2.43 0.33 4.21 2.94 2.27 NA [8,] 3.35 2.15 1.65 0.62 5.98 NA [9,] -5.89 1.32 -2.07 4.91 -1.96 NA [10,] -2.51 -3.79 0.07 -3.92 -9.58 NA [11,] -1.00 -0.31 -0.51 -4.69 -2.40 NA [12,] -0.48 2.59 -0.17 -2.72 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/16fb91225746503.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/2u0191225746503.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/3co8h1225746503.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/48x461225746503.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,] 117.76 113.54 115.30 118.39 113.10 113.94 114.58 117.50 119.65 116.74 [2,] 117.76 118.29 119.97 121.05 119.51 116.78 117.17 118.79 120.44 118.37 [3,] 118.63 121.83 121.11 124.86 129.99 126.60 121.71 119.28 122.63 120.97 [4,] 119.46 122.52 124.10 124.98 131.17 130.16 131.76 134.70 135.32 140.23 [5,] 119.46 122.52 124.10 124.98 144.69 149.18 131.76 134.70 135.32 162.89 [,11] [,12] [1,] 114.23 113.23 [2,] 117.18 116.87 [3,] 118.44 117.93 [4,] 136.31 131.62 [5,] 153.31 150.91 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 117.4288 118.8411 118.1918 122.0831 121.7511 117.1457 111.4007 108.0380 [2,] 119.8312 124.8189 124.0282 127.6369 138.2289 136.0543 132.0193 130.5220 [,9] [,10] [,11] [,12] [1,] 112.1158 105.5238 104.9228 107.5077 [2,] 133.1442 136.4162 131.9572 128.3523 $out [1] 112.75 128.90 133.89 138.21 146.12 156.60 158.87 164.85 $group [1] 1 1 2 3 4 7 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(117.76, 117.76, 118.63, 119.46, 119.46, 113.54, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5qfwh1225746503.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.53 1.58 3.75 -5.29 -3.39 0.39 -2.43 0.62 -5.89 -3.92 -4.69 -2.720 [2,] 0.79 1.58 3.75 -1.54 -2.73 0.39 0.33 1.65 -2.07 -3.92 -2.40 -1.600 [3,] 3.06 1.76 5.01 -1.43 -1.01 0.64 2.27 2.15 -1.96 -3.79 -1.00 -0.325 [4,] 3.20 2.82 5.75 5.01 0.84 1.60 2.94 3.35 1.32 -2.51 -0.51 1.210 [5,] 4.99 4.32 7.91 6.31 4.49 1.60 4.21 3.35 4.91 -2.51 -0.31 2.590 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.3571 0.8838191 3.596805 -6.058213 -3.532553 -0.2149830 0.4257806 [2,] 4.7629 2.6361809 6.423195 3.198213 1.512553 1.4949830 4.1142194 [,8] [,9] [,10] [,11] [,12] [1,] 0.9487843 -4.3553655 -4.786302 -2.3354692 -2.5449 [2,] 3.3512157 0.4353655 -2.793698 0.3354692 1.8949 $out [1] -1.86 -5.71 -4.89 7.42 5.98 0.07 -9.58 $group [1] 2 3 6 6 8 10 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(0.530000000000001, 0.790000000000006, 3.06, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6hrg31225746503.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,] 113.230 112.750 113.100 117.760 128.900 NA [2,] 117.685 116.040 116.255 122.985 141.450 NA [3,] 120.840 117.175 118.415 131.395 150.045 NA [4,] 123.805 119.580 119.950 135.010 157.735 NA [5,] 129.990 121.050 124.100 140.230 164.850 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 118.0486 115.5604 116.7297 125.9103 142.6173 NA [2,] 123.6314 118.7896 120.1003 136.8797 157.4727 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(113.23, 117.685, 120.84, 123.805, 129.99, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7v6p21225746503.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,] 119.500 117.930 118.610 [2,] 124.925 118.955 122.525 [3,] 127.512 121.410 124.355 [4,] 129.096 123.745 126.745 [5,] 132.578 129.990 129.300 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 125.6096 119.2252 122.4302 [2,] 129.4144 123.5948 126.2798 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(119.5, 124.925, 127.512, 129.096, 132.578, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/16fb91225746503.ps tmp/16fb91225746503.png") > system("convert tmp/2u0191225746503.ps tmp/2u0191225746503.png") > system("convert tmp/3co8h1225746503.ps tmp/3co8h1225746503.png") > system("convert tmp/48x461225746503.ps tmp/48x461225746503.png") > system("convert tmp/5qfwh1225746503.ps tmp/5qfwh1225746503.png") > system("convert tmp/6hrg31225746503.ps tmp/6hrg31225746503.png") > system("convert tmp/7v6p21225746503.ps tmp/7v6p21225746503.png") > > > proc.time() user system elapsed 1.420 0.933 2.985