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(106.61,107.86,108.07,109.68,108.57,108.08,109.8,110.7,110.72,109.85,109.71,110,110.87,111.16,110.63,110.26,110.79,111.75,111.75,112.31,112.29,111.86,111.22,112.42,111.94,112.77,113.99,114.22,112.38,114.98,115.94,117.63,117.76,118.11,120.88,121.87,122.34,122.42,122.17,121.99,124.96,125.94,126.59,127.11,125.44,125.53,125.35,124.3,123.62,124,124.47,126.89,124.87,124.55,125.31,126.02,125.31,125.3,123.83,124.49,124.92,125.77,127.27,127.35,126.05,126.84,127.75,126.78,126.98,127.89,127.9,128.06,125.39,128.13,127.97,127.56,126.89,127.8,128.05,128.98,128.55,129.35,129.97,130.89,129.8,132.32,131.1,131.44,132.55,133.35,132.53,131.11,130.02,128.91,129.98,130.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] 96 > (np <- floor(n / par1)) [1] 8 > 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] 8 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 106.61 110.87 111.94 122.34 123.62 124.92 125.39 129.80 NA [2,] 107.86 111.16 112.77 122.42 124.00 125.77 128.13 132.32 NA [3,] 108.07 110.63 113.99 122.17 124.47 127.27 127.97 131.10 NA [4,] 109.68 110.26 114.22 121.99 126.89 127.35 127.56 131.44 NA [5,] 108.57 110.79 112.38 124.96 124.87 126.05 126.89 132.55 NA [6,] 108.08 111.75 114.98 125.94 124.55 126.84 127.80 133.35 NA [7,] 109.80 111.75 115.94 126.59 125.31 127.75 128.05 132.53 NA [8,] 110.70 112.31 117.63 127.11 126.02 126.78 128.98 131.11 NA [9,] 110.72 112.29 117.76 125.44 125.31 126.98 128.55 130.02 NA [10,] 109.85 111.86 118.11 125.53 125.30 127.89 129.35 128.91 NA [11,] 109.71 111.22 120.88 125.35 123.83 127.90 129.97 129.98 NA [12,] 110.00 112.42 121.87 124.30 124.49 128.06 130.89 130.91 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 1.25 0.29 0.83 0.08 0.38 0.85 2.74 2.52 NA [2,] 0.21 -0.53 1.22 -0.25 0.47 1.50 -0.16 -1.22 NA [3,] 1.61 -0.37 0.23 -0.18 2.42 0.08 -0.41 0.34 NA [4,] -1.11 0.53 -1.84 2.97 -2.02 -1.30 -0.67 1.11 NA [5,] -0.49 0.96 2.60 0.98 -0.32 0.79 0.91 0.80 NA [6,] 1.72 0.00 0.96 0.65 0.76 0.91 0.25 -0.82 NA [7,] 0.90 0.56 1.69 0.52 0.71 -0.97 0.93 -1.42 NA [8,] 0.02 -0.02 0.13 -1.67 -0.71 0.20 -0.43 -1.09 NA [9,] -0.87 -0.43 0.35 0.09 -0.01 0.91 0.80 -1.11 NA [10,] -0.14 -0.64 2.77 -0.18 -1.47 0.01 0.62 1.07 NA [11,] 0.29 1.20 0.99 -1.05 0.66 0.16 0.92 0.93 NA [12,] 0.87 -0.48 0.47 -0.68 0.43 -2.67 -1.09 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/161zk1249291886.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/2tkds1249291886.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/3r8441249291886.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/4rfdb1249291886.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,] 106.610 107.860 108.07 109.680 108.570 108.080 109.800 110.700 110.720 [2,] 111.405 111.965 112.31 112.240 111.585 113.365 113.845 114.970 115.025 [3,] 122.980 123.210 123.32 124.440 124.915 125.245 125.950 126.400 125.375 [4,] 125.155 126.950 127.62 127.455 126.470 127.320 127.900 128.045 127.765 [5,] 129.800 132.320 131.10 131.440 132.550 133.350 132.530 131.110 130.020 [,10] [,11] [,12] [1,] 109.850 109.710 110.000 [2,] 114.985 116.050 117.145 [3,] 125.415 124.590 124.395 [4,] 128.400 128.935 129.475 [5,] 129.350 129.980 130.910 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 115.2991 114.8392 114.7676 115.9407 116.6000 117.4495 118.0987 119.0961 [2,] 130.6609 131.5808 131.8724 132.9393 133.2300 133.0405 133.8013 133.7039 [,9] [,10] [,11] [,12] [1,] 118.2583 117.9212 117.3923 117.5073 [2,] 132.4917 132.9088 131.7877 131.2827 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(106.61, 111.405, 122.98, 125.155, 129.8, 107.86, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5fowx1249291886.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,] 0.080 -1.220 -0.410 -2.02 -0.490 -0.820 -1.420 -1.670 -1.110 -1.470 0.160 [2,] 0.335 -0.390 -0.275 -1.57 0.235 0.125 -0.225 -0.900 -0.650 -0.410 0.225 [3,] 0.840 0.025 0.155 -0.89 0.855 0.705 0.635 -0.225 0.040 -0.065 0.790 [4,] 1.885 0.845 0.975 0.82 0.970 0.935 0.915 0.075 0.575 0.845 0.960 [5,] 2.740 1.500 2.420 2.97 0.980 1.720 1.690 0.200 0.910 1.070 1.200 [,12] [1,] -2.670 [2,] -0.885 [3,] -0.480 [4,] 0.450 [5,] 0.870 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.02585225 -0.6648887 -0.543268 -2.2250883 0.4444184 0.2525224 [2,] 1.70585225 0.7148887 0.853268 0.4450883 1.2655816 1.1574776 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.001820367 -0.769649 -0.6443026 -0.766061 0.3794184 -1.2772405 [2,] 1.271820367 0.319649 0.7243026 0.636061 1.2005816 0.3172405 $out [1] 2.60 2.77 -1.05 $group [1] 5 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(0.0799999999999983, 0.334999999999994, 0.839999999999996, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6oetz1249291886.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] [,8] [,9] [1,] 106.610 110.260 111.940 121.990 123.620 124.920 126.890 128.910 NA [2,] 108.075 110.830 113.380 122.380 124.235 126.415 127.680 130.000 NA [3,] 109.695 111.485 115.460 125.155 124.710 127.125 128.090 131.105 NA [4,] 109.925 112.075 117.935 125.735 125.310 127.820 129.165 132.425 NA [5,] 110.720 112.420 121.870 127.110 126.890 128.060 130.890 133.350 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 108.8512 110.9171 113.3824 123.6248 124.2197 126.4842 127.4127 129.9989 [2,] 110.5388 112.0529 117.5376 126.6852 125.2003 127.7658 128.7673 132.2111 [,9] [1,] NA [2,] NA $out [1] 125.39 $group [1] 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(106.61, 108.075, 109.695, 109.925, 110.72, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7ocjo1249291886.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.4363 122.9800 118.3550 [2,] 120.7956 123.8575 119.8306 [3,] 121.8806 124.7525 120.9925 [4,] 122.2850 125.3950 122.1425 [5,] 122.8675 126.4000 124.1375 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 121.2013 124.0512 119.9380 [2,] 122.5599 125.4538 122.0470 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(119.43625, 120.795625, 121.880625, 122.285, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/161zk1249291886.ps tmp/161zk1249291886.png") > system("convert tmp/2tkds1249291886.ps tmp/2tkds1249291886.png") > system("convert tmp/3r8441249291886.ps tmp/3r8441249291886.png") > system("convert tmp/4rfdb1249291886.ps tmp/4rfdb1249291886.png") > system("convert tmp/5fowx1249291886.ps tmp/5fowx1249291886.png") > system("convert tmp/6oetz1249291886.ps tmp/6oetz1249291886.png") > system("convert tmp/7ocjo1249291886.ps tmp/7ocjo1249291886.png") > > > proc.time() user system elapsed 1.189 0.926 1.565