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(100,96.21,96.31,107.18,114.91,92.56,115.00,107.12,117.78,107.37,106.30,114.51,98.00,103.06,100.29,104.61,111.15,104.99,109.93,111.54,132.50,100.34,123.10,114.24,104.57,109.08,106.98,133.68,124.85,122.51,116.80,116.01,129.76,125.20,143.79,127.95,130.30,108.44,129.37,143.68,131.88,117.62,118.96,104.82,134.62,140.40,143.80,153.43,153.29,127.31,153.55,136.93,131.77,144.34,107.42,113.62,124.22,102.06,96.37,111.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] 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,] 100.00 98.00 104.57 130.30 153.29 NA [2,] 96.21 103.06 109.08 108.44 127.31 NA [3,] 96.31 100.29 106.98 129.37 153.55 NA [4,] 107.18 104.61 133.68 143.68 136.93 NA [5,] 114.91 111.15 124.85 131.88 131.77 NA [6,] 92.56 104.99 122.51 117.62 144.34 NA [7,] 115.00 109.93 116.80 118.96 107.42 NA [8,] 107.12 111.54 116.01 104.82 113.62 NA [9,] 117.78 132.50 129.76 134.62 124.22 NA [10,] 107.37 100.34 125.20 140.40 102.06 NA [11,] 106.30 123.10 143.79 143.80 96.37 NA [12,] 114.51 114.24 127.95 153.43 111.68 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -3.79 5.06 4.51 -21.86 -25.98 NA [2,] 0.10 -2.77 -2.10 20.93 26.24 NA [3,] 10.87 4.32 26.70 14.31 -16.62 NA [4,] 7.73 6.54 -8.83 -11.80 -5.16 NA [5,] -22.35 -6.16 -2.34 -14.26 12.57 NA [6,] 22.44 4.94 -5.71 1.34 -36.92 NA [7,] -7.88 1.61 -0.79 -14.14 6.20 NA [8,] 10.66 20.96 13.75 29.80 10.60 NA [9,] -10.41 -32.16 -4.56 5.78 -22.16 NA [10,] -1.07 22.76 18.59 3.40 -5.69 NA [11,] 8.21 -8.86 -15.84 9.63 15.31 NA [12,] -16.51 -9.67 2.35 -0.14 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/1jlmi1257774990.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/2knd01257774990.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/3u2b01257774990.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/48iat1257774990.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,] 98.00 96.21 96.31 104.61 111.15 92.56 107.42 104.82 117.78 100.34 [2,] 100.00 103.06 100.29 107.18 114.91 104.99 109.93 107.12 124.22 102.06 [3,] 104.57 108.44 106.98 133.68 124.85 117.62 115.00 111.54 129.76 107.37 [4,] 130.30 109.08 129.37 136.93 131.77 122.51 116.80 113.62 132.50 125.20 [5,] 153.29 109.08 153.55 143.68 131.88 144.34 118.96 116.01 134.62 140.40 [,11] [,12] [1,] 96.37 111.68 [2,] 106.30 114.24 [3,] 123.10 114.51 [4,] 143.79 127.95 [5,] 143.80 127.95 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 83.1601 104.1863 86.43215 112.6587 112.9368 105.2404 110.1457 106.9471 [2,] 125.9799 112.6937 127.52785 154.7013 136.7632 129.9996 119.8543 116.1329 [,9] [,10] [,11] [,12] [1,] 123.9094 91.01933 96.60966 104.8225 [2,] 135.6106 123.72067 149.59034 124.1975 $out [1] 127.31 153.43 $group [1] 2 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(98, 100, 104.57, 130.3, 153.29, 96.21, 103.06, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5phrw1257774990.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,] -25.98 -2.77 4.32 -11.80 -22.35 -5.71 -14.14 10.60 -32.16 -5.69 -15.84 [2,] -21.86 -2.10 4.32 -8.83 -14.26 -5.71 -7.88 10.66 -22.16 -1.07 -8.86 [3,] -3.79 0.10 10.87 -5.16 -6.16 1.34 -0.79 13.75 -10.41 3.40 8.21 [4,] 4.51 20.93 14.31 6.54 -2.34 4.94 1.61 20.96 -4.56 18.59 9.63 [5,] 5.06 26.24 26.70 7.73 12.57 4.94 6.20 29.80 5.78 22.76 15.31 [,12] [1,] -16.510 [2,] -13.090 [3,] -4.905 [4,] 1.105 [5,] 2.350 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -22.42298 -16.17294 3.811091 -16.020403 -14.582642 -6.185263 -7.49561 [2,] 14.84298 16.37294 17.928909 5.700403 2.262642 8.865263 5.91561 [,8] [,9] [,10] [,11] [,12] [1,] 6.472046 -22.846116 -10.49171 -4.854987 -16.11905 [2,] 21.027954 2.026116 17.29171 21.274987 6.30905 $out [1] -16.62 22.44 -36.92 $group [1] 3 6 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-25.98, -21.86, -3.79000000000001, 4.51000000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6kxt01257774990.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,] 92.560 98.00 104.570 104.82 96.370 NA [2,] 98.155 101.70 112.545 118.29 109.550 NA [3,] 107.150 107.46 123.680 131.09 125.765 NA [4,] 114.710 112.89 128.855 142.04 140.635 NA [5,] 117.780 123.10 143.790 153.43 153.550 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 99.59915 102.3562 116.2409 120.2575 111.5869 NA [2,] 114.70085 112.5638 131.1191 141.9225 139.9431 NA $out [1] 132.5 $group [1] 2 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(92.56, 98.155, 107.15, 114.71, 117.78, 98, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7miyz1257774990.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,] 108.820 104.570 106.0700 [2,] 114.348 107.905 113.4975 [3,] 117.266 114.755 114.9900 [4,] 123.637 123.975 122.6975 [5,] 127.776 133.680 128.3600 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 113.0292 107.4254 110.7938 [2,] 121.5028 122.0846 119.1862 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(108.82, 114.348, 117.266, 123.637, 127.776, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1jlmi1257774990.ps tmp/1jlmi1257774990.png") > system("convert tmp/2knd01257774990.ps tmp/2knd01257774990.png") > system("convert tmp/3u2b01257774990.ps tmp/3u2b01257774990.png") > system("convert tmp/48iat1257774990.ps tmp/48iat1257774990.png") > system("convert tmp/5phrw1257774990.ps tmp/5phrw1257774990.png") > system("convert tmp/6kxt01257774990.ps tmp/6kxt01257774990.png") > system("convert tmp/7miyz1257774990.ps tmp/7miyz1257774990.png") > > > proc.time() user system elapsed 1.252 0.968 1.769