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(100,100.00,93.55,88.17,89.25,91.40,92.47,91.40,88.17,87.10,84.95,92.47,93.55,93.55,91.40,90.32,91.40,93.55,93.55,92.47,91.40,89.25,86.02,88.17,87.10,87.10,86.02,84.95,84.95,86.02,86.02,84.95,86.02,82.80,77.42,80.65,78.49,75.27,75.27,75.27,77.42,78.49,76.34,73.12,68.82,65.59,69.89,82.80,84.95,80.65,74.19,70.97,74.19,82.80,86.02,86.02,82.80,78.49,79.57,87.10,89.25) > 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] 61 > (np <- floor(n / par1)) [1] 5 > 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] 6 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 100.00 93.55 87.10 78.49 84.95 89.25 [2,] 100.00 93.55 87.10 75.27 80.65 NA [3,] 93.55 91.40 86.02 75.27 74.19 NA [4,] 88.17 90.32 84.95 75.27 70.97 NA [5,] 89.25 91.40 84.95 77.42 74.19 NA [6,] 91.40 93.55 86.02 78.49 82.80 NA [7,] 92.47 93.55 86.02 76.34 86.02 NA [8,] 91.40 92.47 84.95 73.12 86.02 NA [9,] 88.17 91.40 86.02 68.82 82.80 NA [10,] 87.10 89.25 82.80 65.59 78.49 NA [11,] 84.95 86.02 77.42 69.89 79.57 NA [12,] 92.47 88.17 80.65 82.80 87.10 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.00 0.00 0.00 -3.22 -4.30 NA [2,] -6.45 -2.15 -1.08 0.00 -6.46 NA [3,] -5.38 -1.08 -1.07 0.00 -3.22 NA [4,] 1.08 1.08 0.00 2.15 3.22 NA [5,] 2.15 2.15 1.07 1.07 8.61 NA [6,] 1.07 0.00 0.00 -2.15 3.22 NA [7,] -1.07 -1.08 -1.07 -3.22 0.00 NA [8,] -3.23 -1.07 1.07 -4.30 -3.22 NA [9,] -1.07 -2.15 -3.22 -3.23 -4.31 NA [10,] -2.15 -3.23 -5.38 4.30 1.08 NA [11,] 7.52 2.15 3.23 12.91 7.53 NA [12,] 1.08 -1.07 -2.16 2.15 2.15 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/17og51257888285.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/2gi4r1257888285.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/3hsbi1257888285.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/4b0ra1257888285.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,] 78.490 75.27 74.19 70.97 74.19 78.49 86.02 84.95 82.80 65.59 69.89 80.65 [2,] 84.950 80.65 75.27 75.27 77.42 82.80 86.02 84.95 82.80 78.49 77.42 82.80 [3,] 88.175 87.10 86.02 84.95 84.95 86.02 86.02 86.02 86.02 82.80 79.57 87.10 [4,] 93.550 93.55 91.40 88.17 89.25 91.40 92.47 91.40 88.17 87.10 84.95 88.17 [5,] 100.000 100.00 93.55 90.32 91.40 93.55 93.55 92.47 91.40 89.25 86.02 92.47 $n [1] 6 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 82.62772 77.9849 74.62258 75.8349 76.59095 79.94326 81.46245 81.46245 [2,] 93.72228 96.2151 97.41742 94.0651 93.30905 92.09674 90.57755 90.57755 [,9] [,10] [,11] [,12] [1,] 82.22557 76.7162 74.24932 83.30557 [2,] 89.81443 88.8838 84.89068 90.89443 $out [1] 76.34 73.12 68.82 $group [1] 7 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(78.49, 84.95, 88.175, 93.55, 100, 75.27, 80.65, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5tl1h1257888285.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] [,12] [1,] -4.30 -6.46 -5.38 0.00 1.07 0.00 -1.08 -4.30 -4.31 -5.38 2.15 -2.16 [2,] -3.22 -6.45 -3.22 1.08 1.07 0.00 -1.08 -3.23 -3.23 -3.23 3.23 -1.07 [3,] 0.00 -2.15 -1.08 1.08 2.15 0.00 -1.07 -3.22 -3.22 -2.15 7.52 1.08 [4,] 0.00 -1.08 -1.07 2.15 2.15 1.07 -1.07 -1.07 -2.15 1.08 7.53 2.15 [5,] 0.00 0.00 0.00 3.22 2.15 1.07 -1.07 1.07 -1.07 4.30 12.91 2.15 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2.275244 -5.944428 -2.5991846 0.3239407 1.386875 -0.7560593 -1.077066 [2,] 2.275244 1.644428 0.4391846 1.8360593 2.913125 0.7560593 -1.062934 [,8] [,9] [,10] [,11] [,12] [1,] -4.746251 -3.983125 -5.1954351 4.481631 -1.195244 [2,] -1.693749 -2.456875 0.8954351 10.558369 3.355244 $out [1] 8.61 -2.15 3.22 -3.22 0.00 $group [1] 5 6 6 7 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-4.3, -3.22, 0, 0, 0, -6.46000000000001, -6.45, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6ddw71257888285.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] [1,] 84.95 86.020 82.800 65.590 70.970 89.25 [2,] 88.17 89.785 83.875 71.505 76.340 89.25 [3,] 91.40 91.400 85.485 75.270 81.725 89.25 [4,] 93.01 93.550 86.020 77.955 85.485 89.25 [5,] 100.00 93.550 87.100 82.800 87.100 89.25 $n [1] 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 89.19244 89.68276 84.50665 72.32811 77.5539 89.25 [2,] 93.60756 93.11724 86.46335 78.21189 85.8961 89.25 $out [1] 77.42 80.65 $group [1] 3 3 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(84.95, 88.17, 91.4, 93.01, 100, 86.02, 89.785, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/75rw81257888285.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,] 79.570 82.800 81.1850 [2,] 82.689 84.950 83.0650 [3,] 84.839 86.020 85.4850 [4,] 86.666 86.560 87.6375 [5,] 88.890 88.175 89.2450 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 83.02506 85.28567 83.39945 [2,] 86.65294 86.75433 87.57055 $out [1] 79.57 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(79.57, 82.689, 84.839, 86.666, 88.89, 82.8, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/17og51257888285.ps tmp/17og51257888285.png") > system("convert tmp/2gi4r1257888285.ps tmp/2gi4r1257888285.png") > system("convert tmp/3hsbi1257888285.ps tmp/3hsbi1257888285.png") > system("convert tmp/4b0ra1257888285.ps tmp/4b0ra1257888285.png") > system("convert tmp/5tl1h1257888285.ps tmp/5tl1h1257888285.png") > system("convert tmp/6ddw71257888285.ps tmp/6ddw71257888285.png") > system("convert tmp/75rw81257888285.ps tmp/75rw81257888285.png") > > > proc.time() user system elapsed 1.919 1.507 2.247