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(0.24,0.23,0.23,0.24,0.23,0.23,0.25,0.21,0.26,0.25,0.24,0.24,0.27,0.25,0.26,0.29,0.24,0.26,0.24,0.26,0.25,0.26,0.24,0.21,0.20,0.22,0.20,0.21,0.20,0.19,0.20,0.20,0.21,0.24,0.22,0.19,0.23,0.23,0.23,0.22,0.23,0.25,0.25,0.22,0.25,0.25,0.24,0.19,0.24,0.26,0.24,0.24,0.25,0.23,0.27,0.24,0.26,0.27,0.29,0.28,0.32,0.29,0.27,0.26,0.28,0.31,0.29,0.31,0.31,0.32,0.32,0.26,0.31,0.31,0.31,0.31,0.29,0.27,0.30,0.27,0.27,0.30,0.28,0.24,0.28,0.28,0.33,0.28,0.29,0.25,0.31,0.29,0.37,0.31,0.29,0.28,0.30,0.32,0.31,0.28,0.29,0.29,0.28,0.26,0.28,0.30,0.33,0.31,0.37,0.36,0.37,0.37,0.36,0.33,0.33,0.40,0.32,0.39,0.39,0.37,0.37,0.30,0.33,0.33,0.34,0.35,0.34,0.37,0.37,0.37,0.36,0.32,0.33,0.35,0.36,0.35,0.37,0.35,0.32,0.33,0.28,0.32,0.35,0.30,0.32,0.32,0.32,0.32,0.36,0.31,0.26,0.33,0.31,0.34,0.33,0.38,0.32,0.30,0.32,0.33,0.34,0.29,0.33,0.36,0.32,0.32,0.32,0.31,0.30,0.34,0.34,0.30,0.28,0.25,0.27,0.33,0.28,0.33,0.32,0.27,0.27,0.28,0.27,0.27,0.25,0.25,0.22,0.27) > 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] 188 > (np <- floor(n / par1)) [1] 15 > 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] 16 16 16 16 16 16 16 16 15 15 15 15 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 0.24 0.27 0.20 0.23 0.24 0.32 0.31 0.28 0.30 0.37 0.37 0.33 0.32 [2,] 0.23 0.25 0.22 0.23 0.26 0.29 0.31 0.28 0.32 0.36 0.30 0.35 0.32 [3,] 0.23 0.26 0.20 0.23 0.24 0.27 0.31 0.33 0.31 0.37 0.33 0.36 0.32 [4,] 0.24 0.29 0.21 0.22 0.24 0.26 0.31 0.28 0.28 0.37 0.33 0.35 0.32 [5,] 0.23 0.24 0.20 0.23 0.25 0.28 0.29 0.29 0.29 0.36 0.34 0.37 0.36 [6,] 0.23 0.26 0.19 0.25 0.23 0.31 0.27 0.25 0.29 0.33 0.35 0.35 0.31 [7,] 0.25 0.24 0.20 0.25 0.27 0.29 0.30 0.31 0.28 0.33 0.34 0.32 0.26 [8,] 0.21 0.26 0.20 0.22 0.24 0.31 0.27 0.29 0.26 0.40 0.37 0.33 0.33 [9,] 0.26 0.25 0.21 0.25 0.26 0.31 0.27 0.37 0.28 0.32 0.37 0.28 0.31 [10,] 0.25 0.26 0.24 0.25 0.27 0.32 0.30 0.31 0.30 0.39 0.37 0.32 0.34 [11,] 0.24 0.24 0.22 0.24 0.29 0.32 0.28 0.29 0.33 0.39 0.36 0.35 0.33 [12,] 0.24 0.21 0.19 0.19 0.28 0.26 0.24 0.28 0.31 0.37 0.32 0.30 0.38 [,14] [,15] [,16] [1,] 0.32 0.30 0.27 [2,] 0.30 0.34 0.28 [3,] 0.32 0.34 0.27 [4,] 0.33 0.30 0.27 [5,] 0.34 0.28 0.25 [6,] 0.29 0.25 0.25 [7,] 0.33 0.27 0.22 [8,] 0.36 0.33 0.27 [9,] 0.32 0.28 NA [10,] 0.32 0.33 NA [11,] 0.32 0.32 NA [12,] 0.31 0.27 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.01 -0.02 0.02 0.00 0.02 -0.03 0.00 0.00 0.02 -0.01 -0.07 0.02 [2,] 0.00 0.01 -0.02 0.00 -0.02 -0.02 0.00 0.05 -0.01 0.01 0.03 0.01 [3,] 0.01 0.03 0.01 -0.01 0.00 -0.01 0.00 -0.05 -0.03 0.00 0.00 -0.01 [4,] -0.01 -0.05 -0.01 0.01 0.01 0.02 -0.02 0.01 0.01 -0.01 0.01 0.02 [5,] 0.00 0.02 -0.01 0.02 -0.02 0.03 -0.02 -0.04 0.00 -0.03 0.01 -0.02 [6,] 0.02 -0.02 0.01 0.00 0.04 -0.02 0.03 0.06 -0.01 0.00 -0.01 -0.03 [7,] -0.04 0.02 0.00 -0.03 -0.03 0.02 -0.03 -0.02 -0.02 0.07 0.03 0.01 [8,] 0.05 -0.01 0.01 0.03 0.02 0.00 0.00 0.08 0.02 -0.08 0.00 -0.05 [9,] -0.01 0.01 0.03 0.00 0.01 0.01 0.03 -0.06 0.02 0.07 0.00 0.04 [10,] -0.01 -0.02 -0.02 -0.01 0.02 0.00 -0.02 -0.02 0.03 0.00 -0.01 0.03 [11,] 0.00 -0.03 -0.03 -0.05 -0.01 -0.06 -0.04 -0.01 -0.02 -0.02 -0.04 -0.05 [12,] 0.03 -0.01 0.04 0.05 0.04 0.05 0.04 0.02 0.06 0.00 0.01 0.02 [,13] [,14] [,15] [,16] [1,] 0.00 -0.02 0.04 0.01 [2,] 0.00 0.02 0.00 -0.01 [3,] 0.00 0.01 -0.04 0.00 [4,] 0.04 0.01 -0.02 -0.02 [5,] -0.05 -0.05 -0.03 0.00 [6,] -0.05 0.04 0.02 -0.03 [7,] 0.07 0.03 0.06 0.05 [8,] -0.02 -0.04 -0.05 NA [9,] 0.03 0.00 0.05 NA [10,] -0.01 0.00 -0.01 NA [11,] 0.05 -0.01 -0.05 NA [12,] -0.06 -0.01 0.00 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/18me11228676790.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/2u63e1228676790.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/36w3h1228676790.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/40igh1228676790.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,] 0.200 0.220 0.20 0.210 0.200 0.190 0.200 0.20 0.210 0.240 0.22 0.19 [2,] 0.255 0.255 0.25 0.250 0.245 0.250 0.250 0.25 0.260 0.265 0.26 0.24 [3,] 0.300 0.295 0.31 0.285 0.285 0.265 0.275 0.28 0.280 0.310 0.32 0.28 [4,] 0.320 0.320 0.33 0.325 0.340 0.310 0.315 0.33 0.315 0.325 0.33 0.31 [5,] 0.370 0.360 0.37 0.370 0.370 0.350 0.340 0.40 0.370 0.390 0.39 0.38 $n [1] 16 16 16 16 16 16 16 16 15 15 15 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 0.274325 0.269325 0.2784 0.255375 0.247475 0.2413 0.249325 0.2484 [2,] 0.325675 0.320675 0.3416 0.314625 0.322525 0.2887 0.300675 0.3116 [,9] [,10] [,11] [,12] [1,] 0.2575625 0.2855227 0.2914432 0.2514432 [2,] 0.3024375 0.3344773 0.3485568 0.3085568 $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(0.2, 0.255, 0.3, 0.32, 0.37, 0.22, 0.255, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/537ty1228676790.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.030 -0.02 -0.030 -0.050 -0.050 -0.050 -0.040 -0.08 -0.01 -0.020 -0.060 [2,] -0.015 -0.01 -0.010 -0.015 -0.030 -0.020 -0.025 -0.03 0.00 -0.015 -0.045 [3,] 0.000 0.00 0.000 0.010 -0.015 0.000 0.015 0.00 0.01 -0.010 -0.030 [4,] 0.020 0.01 0.005 0.010 0.005 0.025 0.040 0.02 0.03 0.000 -0.010 [5,] 0.040 0.03 0.010 0.040 0.030 0.060 0.070 0.08 0.07 0.020 0.000 [,12] [1,] -0.06 [2,] 0.00 [3,] 0.02 [4,] 0.04 [5,] 0.06 $n [1] 16 16 16 16 16 16 16 15 15 15 15 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -0.013825 -0.0079 -0.005925 0.000125 -0.028825 -0.017775 -0.010675 [2,] 0.013825 0.0079 0.005925 0.019875 -0.001175 0.017775 0.040675 [,8] [,9] [,10] [,11] [,12] [1,] -0.02039771 -0.002238627 -0.016119314 -0.0442784 0.00368183 [2,] 0.02039771 0.022238627 -0.003880686 -0.0157216 0.03631817 $out [1] -0.07 0.05 0.03 -0.05 -0.04 -0.06 0.03 0.03 0.05 $group [1] 1 2 3 3 3 9 10 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.03, -0.015, 0, 0.0200000000000000, 0.04, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/68kya1228676790.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] [,10] [,11] [,12] [1,] 0.210 0.210 0.190 0.190 0.230 0.260 0.240 0.25 0.260 0.320 0.300 0.28 [2,] 0.230 0.240 0.200 0.225 0.240 0.275 0.270 0.28 0.280 0.345 0.330 0.32 [3,] 0.240 0.255 0.200 0.230 0.255 0.300 0.295 0.29 0.295 0.370 0.345 0.34 [4,] 0.245 0.260 0.215 0.250 0.270 0.315 0.310 0.31 0.310 0.380 0.370 0.35 [5,] 0.260 0.290 0.220 0.250 0.290 0.320 0.310 0.33 0.330 0.400 0.370 0.37 [,13] [,14] [,15] [,16] [1,] 0.310 0.300 0.250 0.22 [2,] 0.315 0.315 0.275 0.25 [3,] 0.320 0.320 0.300 0.27 [4,] 0.335 0.330 0.330 0.27 [5,] 0.360 0.340 0.340 0.28 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.2331584 0.2458779 0.1931584 0.2185973 0.2413168 0.2817557 0.2767557 [2,] 0.2468416 0.2641221 0.2068416 0.2414027 0.2686832 0.3182443 0.3132443 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 0.2763168 0.2813168 0.3540363 0.3267557 0.3263168 0.3108779 0.3131584 [2,] 0.3036832 0.3086832 0.3859637 0.3632443 0.3536832 0.3291221 0.3268416 [,15] [,16] [1,] 0.2749141 0.2588277 [2,] 0.3250859 0.2811723 $out [1] 0.24 0.37 0.26 0.38 0.29 0.36 $group [1] 3 8 13 13 14 14 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] NA Warning message: In bxp(list(stats = c(0.21, 0.23, 0.24, 0.245, 0.26, 0.21, 0.24, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7wati1228676790.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,] 0.2756250 0.265 0.275000 [2,] 0.2831250 0.280 0.284375 [3,] 0.2896667 0.285 0.290000 [4,] 0.2925000 0.305 0.293125 [5,] 0.3046667 0.320 0.295000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 0.2853907 0.2735973 0.2860091 [2,] 0.2939427 0.2964027 0.2939909 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(0.275625, 0.283125, 0.289666666666667, 0.2925, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/18me11228676790.ps tmp/18me11228676790.png") > system("convert tmp/2u63e1228676790.ps tmp/2u63e1228676790.png") > system("convert tmp/36w3h1228676790.ps tmp/36w3h1228676790.png") > system("convert tmp/40igh1228676790.ps tmp/40igh1228676790.png") > system("convert tmp/537ty1228676790.ps tmp/537ty1228676790.png") > system("convert tmp/68kya1228676790.ps tmp/68kya1228676790.png") > system("convert tmp/7wati1228676790.ps tmp/7wati1228676790.png") > > > proc.time() user system elapsed 1.893 1.488 2.110