R version 2.8.0 (2008-10-20) 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(2.13,1.87,2.23,3,2.12,1.6,1.17,1.02,1.22,1.8,2.13,2.21,2.38,1.99,1.82,2.47,1.94,1.39,1.11,0.97,1.38,2.39,1.88,2.11,2.11,2.17,2.54,3.13,2.25,1.39,1.36,1.33,1.6,1.95,2.23,2.53,2.36,1.95,2.16,2.76,2.09,1.49,1.17,1.3,1.26,2.17,2.03,2.18,2.61,2.58,3.86,3.81,2.41,1.47,1.33,1.38,1.57,2.6,2.18,2.36,2.24,2.41,2.51,2.98,1.87,1.9,1.47,1.45,2.71,2.9,2.11,2.18,2.24,2.05,2.42,2.77,1.99,1.47,1.09,0.93,1.32,2.03,2.04,2.78,2.8,3.03,3.11,2.75,2.78,1.76,1.29,1.28,1.43,1.71,1.89,1.84) > 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,] 2.13 2.38 2.11 2.36 2.61 2.24 2.24 2.80 NA [2,] 1.87 1.99 2.17 1.95 2.58 2.41 2.05 3.03 NA [3,] 2.23 1.82 2.54 2.16 3.86 2.51 2.42 3.11 NA [4,] 3.00 2.47 3.13 2.76 3.81 2.98 2.77 2.75 NA [5,] 2.12 1.94 2.25 2.09 2.41 1.87 1.99 2.78 NA [6,] 1.60 1.39 1.39 1.49 1.47 1.90 1.47 1.76 NA [7,] 1.17 1.11 1.36 1.17 1.33 1.47 1.09 1.29 NA [8,] 1.02 0.97 1.33 1.30 1.38 1.45 0.93 1.28 NA [9,] 1.22 1.38 1.60 1.26 1.57 2.71 1.32 1.43 NA [10,] 1.80 2.39 1.95 2.17 2.60 2.90 2.03 1.71 NA [11,] 2.13 1.88 2.23 2.03 2.18 2.11 2.04 1.89 NA [12,] 2.21 2.11 2.53 2.18 2.36 2.18 2.78 1.84 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -0.26 -0.39 0.06 -0.41 -0.03 0.17 -0.19 0.23 NA [2,] 0.36 -0.17 0.37 0.21 1.28 0.10 0.37 0.08 NA [3,] 0.77 0.65 0.59 0.60 -0.05 0.47 0.35 -0.36 NA [4,] -0.88 -0.53 -0.88 -0.67 -1.40 -1.11 -0.78 0.03 NA [5,] -0.52 -0.55 -0.86 -0.60 -0.94 0.03 -0.52 -1.02 NA [6,] -0.43 -0.28 -0.03 -0.32 -0.14 -0.43 -0.38 -0.47 NA [7,] -0.15 -0.14 -0.03 0.13 0.05 -0.02 -0.16 -0.01 NA [8,] 0.20 0.41 0.27 -0.04 0.19 1.26 0.39 0.15 NA [9,] 0.58 1.01 0.35 0.91 1.03 0.19 0.71 0.28 NA [10,] 0.33 -0.51 0.28 -0.14 -0.42 -0.79 0.01 0.18 NA [11,] 0.08 0.23 0.30 0.15 0.18 0.07 0.74 -0.05 NA [12,] 0.17 0.00 -0.17 0.43 -0.12 0.06 0.02 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/freestat/rcomp/tmp/13ey91228049992.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/2jb3u1228049992.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/3m4641228049992.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/4ogo61228049992.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] [,12] [1,] 2.110 1.870 1.820 2.470 1.870 1.39 1.090 0.930 1.220 1.710 1.880 1.840 [2,] 2.185 1.970 2.195 2.755 1.965 1.43 1.140 0.995 1.290 1.875 1.960 2.145 [3,] 2.300 2.110 2.465 2.875 2.105 1.48 1.230 1.290 1.405 2.100 2.075 2.195 [4,] 2.495 2.495 2.825 3.065 2.330 1.68 1.345 1.355 1.585 2.495 2.155 2.445 [5,] 2.800 3.030 3.110 3.130 2.780 1.90 1.470 1.450 1.600 2.900 2.230 2.780 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2.126830 1.816727 2.113073 2.701830 1.901106 1.340346 1.115484 1.088899 [2,] 2.473170 2.403273 2.816927 3.048170 2.308894 1.619654 1.344516 1.491101 [,9] [,10] [,11] [,12] [1,] 1.240209 1.753659 1.966070 2.027416 [2,] 1.569791 2.446341 2.183930 2.362584 $out [1] 3.86 3.81 2.71 $group [1] 3 4 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(2.11, 2.185, 2.3, 2.495, 2.8, 1.87, 1.97, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5acpg1228049992.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.410 -0.170 -0.360 -1.400 -1.020 -0.47 -0.160 -0.040 0.190 -0.790 -0.050 [2,] -0.325 0.090 0.150 -0.995 -0.900 -0.43 -0.145 0.170 0.315 -0.465 0.075 [3,] -0.110 0.285 0.530 -0.830 -0.575 -0.35 -0.025 0.235 0.645 -0.065 0.165 [4,] 0.115 0.370 0.625 -0.600 -0.520 -0.21 0.020 0.400 0.960 0.230 0.265 [5,] 0.230 0.370 0.770 -0.530 0.030 -0.03 0.130 0.410 1.030 0.330 0.300 [,12] [1,] -0.170 [2,] -0.060 [3,] 0.020 [4,] 0.115 [5,] 0.170 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.3557903 0.1285880 0.2646582 -1.0506527 -0.7872735 -0.4728952 [2,] 0.1357903 0.4414120 0.7953418 -0.6093473 -0.3627265 -0.2271048 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.11717137 0.1065187 0.2846937 -0.453237 0.05886327 -0.08450718 [2,] 0.06717137 0.3634813 1.0053063 0.323237 0.27113673 0.12450718 $out [1] 1.28 0.03 1.26 0.74 0.43 $group [1] 2 4 8 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.41, -0.325, -0.11, 0.115, 0.23, -0.17, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6j8ie1228049992.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,] 1.020 0.970 1.330 1.170 1.330 1.450 0.930 1.280 NA [2,] 1.410 1.385 1.495 1.395 1.520 1.885 1.395 1.570 NA [3,] 1.995 1.910 2.140 2.060 2.385 2.210 2.035 1.865 NA [4,] 2.170 2.245 2.390 2.175 2.605 2.610 2.330 2.790 NA [5,] 3.000 2.470 3.130 2.760 3.860 2.980 2.780 3.110 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1.648359 1.517748 1.731784 1.704237 1.890124 1.879323 1.608540 1.308550 [2,] 2.341641 2.302252 2.548216 2.415763 2.879876 2.540677 2.461460 2.421450 [,9] [1,] NA [2,] NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(1.02, 1.41, 1.995, 2.17, 3, 0.97, 1.385, 1.91, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7mmv51228049992.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,] 1.20750 1.2300 1.175000 [2,] 1.56000 1.4425 1.493125 [3,] 2.18750 2.1025 2.155625 [4,] 2.31625 2.2475 2.303750 [5,] 2.95875 2.8750 2.895000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1.842569 1.735334 1.785893 [2,] 2.532431 2.469666 2.525357 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1.2075, 1.56, 2.1875, 2.31625, 2.95875, 1.23, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/13ey91228049992.ps tmp/13ey91228049992.png") > system("convert tmp/2jb3u1228049992.ps tmp/2jb3u1228049992.png") > system("convert tmp/3m4641228049992.ps tmp/3m4641228049992.png") > system("convert tmp/4ogo61228049992.ps tmp/4ogo61228049992.png") > system("convert tmp/5acpg1228049992.ps tmp/5acpg1228049992.png") > system("convert tmp/6j8ie1228049992.ps tmp/6j8ie1228049992.png") > system("convert tmp/7mmv51228049992.ps tmp/7mmv51228049992.png") > > > proc.time() user system elapsed 1.867 1.463 2.224