R version 2.8.1 (2008-12-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(115.47,103.34,102.6,100.69,105.67,123.61,113.08,106.46,123.38,109.87,95.74,123.06,123.39,120.28,115.33,110.4,114.49,132.03,123.16,118.82,128.32,112.24,104.53,132.57,122.52,131.8,124.55,120.96,122.6,145.52,118.57,134.25,136.7,121.37,111.63,134.42,137.65,137.86,119.77,130.69,128.28,147.45,128.42,136.9,143.95,135.64,122.48,136.83,153.04,142.71,123.46,144.37,146.15,147.61,158.51,147.4,165.05,154.64,126.2,157.36,154.15,123.21,113.07,110.45,113.57,122.44,114.93,111.85,126.04,121.34) > 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] 70 > (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 6 6 6 6 6 6 6 6 6 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 115.47 123.39 122.52 137.65 153.04 154.15 [2,] 103.34 120.28 131.80 137.86 142.71 123.21 [3,] 102.60 115.33 124.55 119.77 123.46 113.07 [4,] 100.69 110.40 120.96 130.69 144.37 110.45 [5,] 105.67 114.49 122.60 128.28 146.15 113.57 [6,] 123.61 132.03 145.52 147.45 147.61 122.44 [7,] 113.08 123.16 118.57 128.42 158.51 114.93 [8,] 106.46 118.82 134.25 136.90 147.40 111.85 [9,] 123.38 128.32 136.70 143.95 165.05 126.04 [10,] 109.87 112.24 121.37 135.64 154.64 121.34 [11,] 95.74 104.53 111.63 122.48 126.20 NA [12,] 123.06 132.57 134.42 136.83 157.36 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -12.13 -3.11 9.28 0.21 -10.33 -30.94 [2,] -0.74 -4.95 -7.25 -18.09 -19.25 -10.14 [3,] -1.91 -4.93 -3.59 10.92 20.91 -2.62 [4,] 4.98 4.09 1.64 -2.41 1.78 3.12 [5,] 17.94 17.54 22.92 19.17 1.46 8.87 [6,] -10.53 -8.87 -26.95 -19.03 10.90 -7.51 [7,] -6.62 -4.34 15.68 8.48 -11.11 -3.08 [8,] 16.92 9.50 2.45 7.05 17.65 14.19 [9,] -13.51 -16.08 -15.33 -8.31 -10.41 -4.70 [10,] -14.13 -7.71 -9.74 -13.16 -28.44 NA [11,] 27.32 28.04 22.79 14.35 31.16 NA [12,] 0.33 -10.05 3.23 16.21 -3.21 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/1y9o61258036858.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/2tdq71258036858.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/3aan61258036858.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/40cpi1258036858.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] [1,] 115.47 103.340 102.60 100.690 105.670 122.440 113.080 106.460 123.38 [2,] 122.52 120.280 113.07 110.400 113.570 123.610 114.930 111.850 126.04 [3,] 130.52 127.505 117.55 115.705 118.545 138.775 120.865 126.535 132.51 [4,] 153.04 137.860 123.46 130.690 128.280 147.450 128.420 136.900 143.95 [5,] 154.15 142.710 124.55 144.370 146.150 147.610 128.420 147.400 165.05 [,10] [,11] [,12] [1,] 109.870 95.74 132.57 [2,] 112.240 104.53 132.57 [3,] 121.355 111.63 134.42 [4,] 135.640 122.48 136.83 [5,] 154.640 126.20 136.83 $n [1] 6 6 6 6 6 6 6 6 6 6 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 110.8336 116.1653 110.8481 102.6173 109.0566 123.3974 112.1635 110.3769 [2,] 150.2064 138.8447 124.2519 128.7927 128.0334 154.1526 129.5665 142.6931 [,9] [,10] [,11] [,12] [1,] 120.9575 106.2612 98.94658 131.4099 [2,] 144.0625 136.4488 124.31342 137.4301 $out [1] 158.51 123.06 157.36 $group [1] 7 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(115.47, 122.52, 130.52, 153.04, 154.15, 103.34, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5bggf1258036858.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,] -12.13 -19.250 -4.930 1.64 1.46 -26.95 -11.11 2.450 -16.08 -14.13 22.79 [2,] -12.13 -18.090 -3.590 1.64 8.87 -19.03 -6.62 7.050 -15.33 -14.13 22.79 [3,] -6.72 -8.695 -2.265 2.45 17.74 -9.70 -3.71 11.845 -11.96 -13.16 27.32 [4,] 0.21 -4.950 10.920 4.09 19.17 -7.51 8.48 16.920 -8.31 -9.74 28.04 [5,] 9.28 -0.740 20.910 4.98 22.92 -7.51 15.68 17.650 -4.70 -7.71 31.16 [,12] [1,] -10.05 [2,] -3.21 [3,] 0.33 [4,] 3.23 [5,] 3.23 $n [1] 6 6 6 6 6 6 6 6 6 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -14.679699 -17.1707244 -11.624419 0.8696709 11.09617 -17.130772 -13.449988 [2,] 1.239699 -0.2192756 7.094419 4.0303291 24.38383 -2.269228 6.029988 [,8] [,9] [,10] [,11] [,12] [1,] 5.478531 -16.488127 -16.26196 23.61036 -4.220488 [2,] 18.211469 -7.431873 -10.05804 31.02964 4.880488 $out [1] -30.94 -2.41 10.90 -28.44 14.35 16.21 $group [1] 1 4 6 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-12.13, -12.13, -6.71999999999999, 0.210000000000008, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/645is1258036858.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,] 95.740 104.530 111.630 119.770 126.200 110.450 [2,] 102.970 113.365 121.165 128.350 143.540 113.070 [3,] 108.165 119.550 123.575 136.235 147.505 118.135 [4,] 119.265 125.855 134.335 137.755 156.000 123.210 [5,] 123.610 132.570 145.520 147.450 165.050 126.040 $n [1] 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 100.7327 113.8532 117.5681 131.9453 141.8219 113.0687 [2,] 115.5973 125.2468 129.5819 140.5247 153.1881 123.2013 $out [1] 123.46 154.15 $group [1] 5 6 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(95.74, 102.97, 108.165, 119.265, 123.61, 104.53, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7uq521258036858.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,] 112.1160 111.6300 113.5050 [2,] 120.6933 118.0475 119.8325 [3,] 126.0292 123.9450 124.1044 [4,] 135.4067 131.5150 134.5369 [5,] 137.2400 138.7750 136.3413 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 119.3183 117.8024 117.3976 [2,] 132.7400 130.0876 130.8111 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(112.116, 120.693333333333, 126.029166666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1y9o61258036858.ps tmp/1y9o61258036858.png") > system("convert tmp/2tdq71258036858.ps tmp/2tdq71258036858.png") > system("convert tmp/3aan61258036858.ps tmp/3aan61258036858.png") > system("convert tmp/40cpi1258036858.ps tmp/40cpi1258036858.png") > system("convert tmp/5bggf1258036858.ps tmp/5bggf1258036858.png") > system("convert tmp/645is1258036858.ps tmp/645is1258036858.png") > system("convert tmp/7uq521258036858.ps tmp/7uq521258036858.png") > > > proc.time() user system elapsed 1.430 1.440 2.897