R version 2.6.2 (2008-02-08) 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. Natural language support but running in an English locale 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(23.11,18.64,14.94,16.90,15.46,11.15,13.13,12.48,12.95,12.59,10.58,10.58,12.39,15.53,13.06,10.22,16.33,19.72,21.31,18.84,24.84,15.67,15.57,12.73,13.56,15.54,17.22,12.14,11.07,12.02,11.55,6.92,10.33,8.38,12.11,11.46,12.75,13.32,13.00,11.90,11.79,12.55,11.84,11.25,11.15,10.99,11.70,14.01,17.51,17.27,16.90,15.79,15.45,16.24,16.71,16.77,16.64,17.80,16.87,16.13,15.76,15.66,15.54,15.30,15.05,14.69,14.39,14.18,13.70,13.66,13.27,13.56,13.14,14.19,22.57,23.09,23.31,22.91,22.36,43.06,64.67,64.68,56.90,48.79,45.21,41.40,22.17,25.52,20.28,22.87,27.63,22.95,21.35,18.38,17.15,18.27,19.40,20.52) > 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] 98 > (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] 9 9 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 23.11 12.39 13.56 12.75 17.51 15.76 13.14 45.21 19.40 [2,] 18.64 15.53 15.54 13.32 17.27 15.66 14.19 41.40 20.52 [3,] 14.94 13.06 17.22 13.00 16.90 15.54 22.57 22.17 NA [4,] 16.90 10.22 12.14 11.90 15.79 15.30 23.09 25.52 NA [5,] 15.46 16.33 11.07 11.79 15.45 15.05 23.31 20.28 NA [6,] 11.15 19.72 12.02 12.55 16.24 14.69 22.91 22.87 NA [7,] 13.13 21.31 11.55 11.84 16.71 14.39 22.36 27.63 NA [8,] 12.48 18.84 6.92 11.25 16.77 14.18 43.06 22.95 NA [9,] 12.95 24.84 10.33 11.15 16.64 13.70 64.67 21.35 NA [10,] 12.59 15.67 8.38 10.99 17.80 13.66 64.68 18.38 NA [11,] 10.58 15.57 12.11 11.70 16.87 13.27 56.90 17.15 NA [12,] 10.58 12.73 11.46 14.01 16.13 13.56 48.79 18.27 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -4.47 3.14 1.98 0.57 -0.24 -0.10 1.05 -3.81 1.12 [2,] -3.70 -2.47 1.68 -0.32 -0.37 -0.12 8.38 -19.23 NA [3,] 1.96 -2.84 -5.08 -1.10 -1.11 -0.24 0.52 3.35 NA [4,] -1.44 6.11 -1.07 -0.11 -0.34 -0.25 0.22 -5.24 NA [5,] -4.31 3.39 0.95 0.76 0.79 -0.36 -0.40 2.59 NA [6,] 1.98 1.59 -0.47 -0.71 0.47 -0.30 -0.55 4.76 NA [7,] -0.65 -2.47 -4.63 -0.59 0.06 -0.21 20.70 -4.68 NA [8,] 0.47 6.00 3.41 -0.10 -0.13 -0.48 21.61 -1.60 NA [9,] -0.36 -9.17 -1.95 -0.16 1.16 -0.04 0.01 -2.97 NA [10,] -2.01 -0.10 3.73 0.71 -0.93 -0.39 -7.78 -1.23 NA [11,] 0.00 -2.84 -0.65 2.31 -0.74 0.29 -8.11 1.12 NA [12,] 1.81 0.83 1.29 3.50 -0.37 -0.42 -3.58 1.13 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/freestat/rcomp/tmp/1s2z91209200371.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/freestat/rcomp/tmp/26mk31209200371.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/freestat/rcomp/tmp/3m21u1209200372.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/freestat/rcomp/tmp/4ox791209200372.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] [,11] [1,] 12.39 13.32 13.000 10.220 11.070 11.150 11.550 6.920 10.330 8.380 10.580 [2,] 13.14 15.53 14.000 12.020 13.420 12.285 12.485 11.865 12.050 11.790 11.905 [3,] 15.76 15.66 16.220 15.545 15.455 15.465 15.550 15.475 15.170 14.665 14.420 [4,] 19.40 18.64 19.695 19.995 18.305 21.295 21.835 20.895 23.095 18.090 17.010 [5,] 23.11 20.52 22.570 25.520 23.310 22.910 27.630 22.950 24.840 18.380 17.150 [,12] [1,] 10.580 [2,] 12.095 [3,] 13.785 [4,] 17.200 [5,] 18.270 $n [1] 9 9 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 12.46307 14.02207 13.03869 11.09005 12.72617 10.43188 10.32696 10.43071 [2,] 19.05693 17.29793 19.40131 19.99995 18.18383 20.49812 20.77304 20.51929 [,9] [,10] [,11] [,12] [1,] 9.000104 11.14573 11.56827 10.93327 [2,] 21.339896 18.18427 17.27173 16.63673 $out [1] 45.21 41.40 43.06 64.67 64.68 56.90 48.79 $group [1] 1 2 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(12.39, 13.14, 15.76, 19.4, 23.11, 13.32, 15.53, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5c0mp1209200372.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] [1,] -0.24 -3.700 -5.080 -1.440 -0.400 -0.710 -4.680 -1.600 -2.970 -2.010 [2,] -0.24 -3.085 -1.975 -1.255 -0.380 -0.510 -3.550 -0.305 -2.460 -1.620 [3,] 0.57 -0.345 -0.670 -0.295 0.775 0.085 -0.620 0.185 -0.260 -0.660 [4,] 1.12 0.780 1.240 0.055 1.770 1.785 -0.075 4.705 -0.015 0.305 [5,] 3.14 1.680 3.350 0.220 3.390 4.760 0.060 6.000 1.160 0.710 [,11] [,12] [1,] -2.840 -0.420 [2,] -1.790 -0.395 [3,] -0.325 0.980 [4,] 0.705 1.550 [5,] 2.310 3.500 $n [1] 9 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -0.1462667 -2.504044 -2.465945 -1.0267848 -0.4260209 -1.19702 -2.561185 [2,] 1.2862667 1.814044 1.125945 0.4367848 1.9760209 1.36702 1.321185 [,8] [,9] [,10] [,11] [,12] [1,] -2.613658 -1.625812 -1.7353326 -1.718743 -0.1065049 [2,] 2.983658 1.105812 0.4153326 1.068743 2.0665049 $out [1] -4.47 -3.81 8.38 -19.23 6.11 -5.24 -4.31 20.70 21.61 -9.17 [11] 3.73 -7.78 -8.11 -3.58 $group [1] 1 1 2 2 4 4 5 7 8 9 10 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.240000000000002, -0.240000000000002, 0.57, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/690f81209200372.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,] 10.580 10.220 8.380 10.990 15.450 13.27 13.140 17.150 19.40 [2,] 11.815 12.895 10.700 11.475 16.185 13.68 22.465 19.330 19.40 [3,] 13.040 15.620 11.785 11.870 16.740 14.54 23.200 22.520 19.96 [4,] 16.180 19.280 12.850 12.875 17.085 15.42 52.845 26.575 20.52 [5,] 18.640 24.840 15.540 14.010 17.800 15.76 64.680 27.630 20.52 $n [1] 12 12 12 12 12 12 12 12 2 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 11.04909 12.70776 10.80437 11.23145 16.32950 13.74637 9.343478 19.21551 [2,] 15.03091 18.53224 12.76563 12.50855 17.15050 15.33363 37.056522 25.82449 [,9] [1,] 18.70870 [2,] 21.21130 $out [1] 23.11 17.22 6.92 45.21 41.40 $group [1] 1 3 3 8 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(10.58, 11.815, 13.04, 16.18, 18.64, 10.22, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7cqr81209200372.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,] 16.09250 14.4200 14.47375 [2,] 16.72188 14.9175 15.16562 [3,] 18.24875 15.4700 16.14500 [4,] 19.23604 15.6050 16.77437 [5,] 21.95375 16.2200 17.36125 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 17.10202 15.15643 15.41124 [2,] 19.39548 15.78357 16.87876 $out [1] 13.785 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(16.0925, 16.721875, 18.24875, 19.2360416666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1s2z91209200371.ps tmp/1s2z91209200371.png") > system("convert tmp/26mk31209200371.ps tmp/26mk31209200371.png") > system("convert tmp/3m21u1209200372.ps tmp/3m21u1209200372.png") > system("convert tmp/4ox791209200372.ps tmp/4ox791209200372.png") > system("convert tmp/5c0mp1209200372.ps tmp/5c0mp1209200372.png") > system("convert tmp/690f81209200372.ps tmp/690f81209200372.png") > system("convert tmp/7cqr81209200372.ps tmp/7cqr81209200372.png") > > > proc.time() user system elapsed 2.115 1.368 2.275