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(36.80,35.40,33.00,28.73,26.70,26.46,24.60,28.00,31.60,33.50,34.50,35.00,34.76,33.50,32.74,34.40,31.93,29.24,25.75,26.03,26.08,23.80,20.61,19.70,18.18,19.60,20.60,20.03,23.00,23.60,22.56,22.55,23.75,24.92,24.50,30.58,28.07,27.70,27.00,25.23,26.86,25.60,24.55,23.96,23.50,23.64,21.55,21.05,21.89,21.98,21.45,22.15,22.58,23.80,23.30,22.38,23.00,21.96,22.40,20.80,20.40,16.00,12.78,9.75,7.50,11.24,12.24,12.75,12.52,14.49,14.21,14.32,22.15,22.58,23.80,23.30,22.38,23.00,21.96,22.40,20.80,20.40,16.00,12.78,9.75,7.50,11.24,12.24,12.75,12.52,14.49,14.21,14.32) > 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] 93 > (np <- floor(n / par1)) [1] 7 > 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 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 36.80 34.76 18.18 28.07 21.89 20.40 22.15 9.75 [2,] 35.40 33.50 19.60 27.70 21.98 16.00 22.58 7.50 [3,] 33.00 32.74 20.60 27.00 21.45 12.78 23.80 11.24 [4,] 28.73 34.40 20.03 25.23 22.15 9.75 23.30 12.24 [5,] 26.70 31.93 23.00 26.86 22.58 7.50 22.38 12.75 [6,] 26.46 29.24 23.60 25.60 23.80 11.24 23.00 12.52 [7,] 24.60 25.75 22.56 24.55 23.30 12.24 21.96 14.49 [8,] 28.00 26.03 22.55 23.96 22.38 12.75 22.40 14.21 [9,] 31.60 26.08 23.75 23.50 23.00 12.52 20.80 14.32 [10,] 33.50 23.80 24.92 23.64 21.96 14.49 20.40 NA [11,] 34.50 20.61 24.50 21.55 22.40 14.21 16.00 NA [12,] 35.00 19.70 30.58 21.05 20.80 14.32 12.78 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1.40 -1.26 1.42 -0.37 0.09 -4.40 0.43 -2.25 [2,] -2.40 -0.76 1.00 -0.70 -0.53 -3.22 1.22 3.74 [3,] -4.27 1.66 -0.57 -1.77 0.70 -3.03 -0.50 1.00 [4,] -2.03 -2.47 2.97 1.63 0.43 -2.25 -0.92 0.51 [5,] -0.24 -2.69 0.60 -1.26 1.22 3.74 0.62 -0.23 [6,] -1.86 -3.49 -1.04 -1.05 -0.50 1.00 -1.04 1.97 [7,] 3.40 0.28 -0.01 -0.59 -0.92 0.51 0.44 -0.28 [8,] 3.60 0.05 1.20 -0.46 0.62 -0.23 -1.60 0.11 [9,] 1.90 -2.28 1.17 0.14 -1.04 1.97 -0.40 NA [10,] 1.00 -3.19 -0.42 -2.09 0.44 -0.28 -4.40 NA [11,] 0.50 -0.91 6.08 -0.50 -1.60 0.11 -3.22 NA [12,] -0.24 -1.52 -2.51 0.84 -0.40 7.83 -3.03 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/1yon01244393654.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/2dhyv1244393654.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/3g4im1244393654.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/48yft1244393654.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,] 9.750 7.50 11.240 9.750 7.500 11.24 12.240 12.750 12.520 20.40 14.210 [2,] 19.290 17.80 16.690 16.135 17.565 17.76 18.225 18.295 17.560 21.18 18.305 [3,] 22.020 22.28 22.625 22.725 22.790 23.70 22.930 22.475 23.250 23.64 21.550 [4,] 31.415 30.60 29.870 26.980 26.780 26.03 24.575 24.995 24.915 24.36 23.450 [5,] 36.800 35.40 33.000 34.400 31.930 29.24 25.750 28.000 31.600 24.92 24.500 [,12] [1,] 12.780 [2,] 17.010 [3,] 20.800 [4,] 25.815 [5,] 35.000 $n [1] 8 8 8 8 8 8 8 8 8 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 15.2468 15.12974 15.26246 16.66683 17.64237 19.08026 19.3828 18.73228 [2,] 28.7932 29.43026 29.98754 28.78317 27.93763 28.31974 26.4772 26.21772 [,9] [,10] [,11] [,12] [1,] 19.14139 21.74096 18.47749 15.54180 [2,] 27.35861 25.53904 24.62251 26.05820 $out [1] 33.50 14.49 34.50 $group [1] 10 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(9.75, 19.29, 22.02, 31.415, 36.8, 7.5, 17.8, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/53bix1244393654.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,] -4.400 -3.220 -4.270 -2.470 -2.690 -3.490 -0.920 -1.600 -2.280 -4.40 [2,] -1.825 -1.580 -2.400 -2.140 -0.750 -1.455 -0.435 -0.345 -0.720 -2.64 [3,] -0.815 -0.615 -0.535 -0.245 0.185 -1.040 0.135 0.080 0.140 -0.42 [4,] 0.260 1.110 0.850 1.070 0.920 0.250 0.475 0.910 1.535 0.08 [5,] 1.420 3.740 1.660 2.970 1.220 1.970 0.510 1.200 1.970 1.00 [,11] [,12] [1,] -3.220 -3.030 [2,] -1.255 -2.015 [3,] -0.500 -0.400 [4,] 0.305 0.300 [5,] 0.500 0.840 $n [1] 8 8 8 8 8 8 8 8 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.9797109 -2.1176726 -2.350497 -2.038152 -0.747886 -1.99243748 -0.3733391 [2,] 0.3497109 0.8876726 1.280497 1.548152 1.117886 -0.08756252 0.6433391 [,8] [,9] [,10] [,11] [,12] [1,] -0.621061 -1.206650 -2.04434 -1.4316068 -1.7824807 [2,] 0.781061 1.486650 1.20434 0.4316068 0.9824807 $out [1] 3.74 3.40 3.60 6.08 7.83 $group [1] 5 7 8 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-4.4, -1.825, -0.815, 0.26, 1.42000000000000, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6bm001244393654.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] [1,] 24.60 19.700 18.180 21.05 20.800 9.750 20.400 7.50 [2,] 27.35 24.775 20.315 23.57 21.925 11.740 20.600 11.24 [3,] 32.30 27.660 22.780 24.89 22.265 12.765 22.265 12.52 [4,] 34.75 33.120 24.125 26.93 22.790 14.405 22.790 14.21 [5,] 36.80 34.760 24.920 28.07 23.800 16.000 23.800 14.49 $n [1] 12 12 12 12 12 12 12 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 28.92481 23.85379 21.04223 23.35748 21.87047 11.54948 21.26613 10.9558 [2,] 35.67519 31.46621 24.51777 26.42252 22.65953 13.98052 23.26387 14.0842 $out [1] 30.58 20.40 7.50 16.00 12.78 $group [1] 3 6 6 7 7 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(24.6, 27.35, 32.3, 34.75, 36.8, 19.7, 24.775, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/70r841244393654.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,] 21.18125 20.800 20.87750 [2,] 21.82250 22.150 21.92500 [3,] 21.97295 22.675 22.58875 [4,] 22.92938 23.090 23.44812 [5,] 24.00000 23.700 24.79375 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 21.46809 22.24626 21.89404 [2,] 22.47780 23.10374 23.28346 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(21.18125, 21.8225, 21.9729464285714, 22.929375, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1yon01244393654.ps tmp/1yon01244393654.png") > system("convert tmp/2dhyv1244393654.ps tmp/2dhyv1244393654.png") > system("convert tmp/3g4im1244393654.ps tmp/3g4im1244393654.png") > system("convert tmp/48yft1244393654.ps tmp/48yft1244393654.png") > system("convert tmp/53bix1244393654.ps tmp/53bix1244393654.png") > system("convert tmp/6bm001244393654.ps tmp/6bm001244393654.png") > system("convert tmp/70r841244393654.ps tmp/70r841244393654.png") > > > proc.time() user system elapsed 1.860 1.451 2.199