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(80.9,113.8,71.2,55.8,61,60.8,74.3,66.5,71,61.6,66.4,74.7,72.4,76.4,86.4,73.9,133.1,83.3,85.8,76.8,92.7,68.8,69.4,97,84.6,84.2,201.8,83.1,69.4,78.8,77.1,90.7,107.5,81.4,83.3,105.7,85.8,87.2,87.5,77.5,78.5,89.7,94.3,104.4,115.9,90,91.1,119.8,101.9,99.2,115.7,94.9,99.8,109.6,109,119.6,119.2,94.6,102.5,105.4) > 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] 60 > (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] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 80.9 72.4 84.6 85.8 101.9 NA [2,] 113.8 76.4 84.2 87.2 99.2 NA [3,] 71.2 86.4 201.8 87.5 115.7 NA [4,] 55.8 73.9 83.1 77.5 94.9 NA [5,] 61.0 133.1 69.4 78.5 99.8 NA [6,] 60.8 83.3 78.8 89.7 109.6 NA [7,] 74.3 85.8 77.1 94.3 109.0 NA [8,] 66.5 76.8 90.7 104.4 119.6 NA [9,] 71.0 92.7 107.5 115.9 119.2 NA [10,] 61.6 68.8 81.4 90.0 94.6 NA [11,] 66.4 69.4 83.3 91.1 102.5 NA [12,] 74.7 97.0 105.7 119.8 105.4 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 32.9 4.0 -0.4 1.4 -2.7 NA [2,] -42.6 10.0 117.6 0.3 16.5 NA [3,] -15.4 -12.5 -118.7 -10.0 -20.8 NA [4,] 5.2 59.2 -13.7 1.0 4.9 NA [5,] -0.2 -49.8 9.4 11.2 9.8 NA [6,] 13.5 2.5 -1.7 4.6 -0.6 NA [7,] -7.8 -9.0 13.6 10.1 10.6 NA [8,] 4.5 15.9 16.8 11.5 -0.4 NA [9,] -9.4 -23.9 -26.1 -25.9 -24.6 NA [10,] 4.8 0.6 1.9 1.1 7.9 NA [11,] 8.3 27.6 22.4 28.7 2.9 NA [12,] -2.3 -12.4 -19.9 -17.9 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/1jfix1257968487.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/2d2ju1257968487.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/3zuez1257968487.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/4yc2k1257968487.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,] 80.9 76.4 71.2 73.9 61.0 78.8 74.3 66.5 71.0 61.6 66.4 97.0 [2,] 80.9 84.2 86.4 73.9 69.4 78.8 77.1 76.8 92.7 68.8 69.4 97.0 [3,] 84.6 87.2 87.5 77.5 78.5 83.3 85.8 90.7 107.5 81.4 83.3 105.4 [4,] 85.8 99.2 115.7 83.1 99.8 89.7 94.3 104.4 115.9 90.0 91.1 105.7 [5,] 85.8 113.8 115.7 94.9 133.1 89.7 109.0 119.6 119.2 94.6 102.5 105.7 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 81.13767 76.60104 66.7967 70.9993 57.01944 75.59809 73.64652 71.19791 [2,] 88.06233 97.79896 108.2033 84.0007 99.98056 91.00191 97.95348 110.20209 [,9] [,10] [,11] [,12] [1,] 91.10694 66.42013 67.96683 99.2526 [2,] 123.89306 96.37987 98.63317 111.5474 $out [1] 72.4 101.9 201.8 55.8 60.8 109.6 74.7 119.8 $group [1] 1 1 3 4 6 6 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(80.9, 80.9, 84.6, 85.8, 85.8, 76.4, 84.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/54k1w1257968487.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,] -2.7 0.3 -20.8 1.0 -0.2 -1.7 -9.0 -0.4 -26.1 0.6 2.9 -19.90 [2,] -0.4 0.3 -20.8 1.0 -0.2 -0.6 -7.8 4.5 -25.9 1.1 8.3 -18.90 [3,] 1.4 10.0 -15.4 4.9 9.4 2.5 10.1 11.5 -24.6 1.9 22.4 -15.15 [4,] 4.0 16.5 -12.5 5.2 9.8 4.6 10.6 15.9 -23.9 4.8 27.6 -7.35 [5,] 4.0 16.5 -10.0 5.2 11.2 4.6 13.6 16.8 -23.9 7.9 28.7 -2.30 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.709029 -1.446879 -21.26476 1.932291 2.334025 -1.174307 -2.901394 [2,] 4.509029 21.446879 -9.53524 7.867709 16.465975 6.174307 23.101394 [,8] [,9] [,10] [,11] [,12] [1,] 3.444789 -26.01319 -0.7144107 8.762669 -24.2745 [2,] 19.555211 -23.18681 4.5144107 36.037331 -6.0255 $out [1] 32.9 -42.6 117.6 -118.7 59.2 -13.7 -49.8 13.5 -9.4 $group [1] 1 2 2 3 4 4 5 6 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-2.70000000000000, -0.399999999999991, 1.40000000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6pxir1257968487.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,] 55.80 68.80 69.40 77.50 94.60 NA [2,] 61.30 73.15 80.10 86.50 99.50 NA [3,] 68.75 80.05 83.75 89.85 103.95 NA [4,] 74.50 89.55 98.20 99.35 112.65 NA [5,] 80.90 97.00 107.50 115.90 119.60 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 62.72939 72.56985 75.49447 83.98903 97.9522 NA [2,] 74.77061 87.53015 92.00553 95.71097 109.9478 NA $out [1] 113.8 133.1 201.8 119.8 $group [1] 1 2 3 4 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(55.8, 61.3, 68.75, 74.5, 80.9, 68.8, 73.15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7l2c41257968487.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,] 77.04 77.50 78.500 [2,] 83.49 82.35 81.800 [3,] 88.23 85.20 85.150 [4,] 96.34 89.10 96.375 [5,] 112.52 90.70 104.300 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 82.36903 82.12128 78.50224 [2,] 94.09097 88.27872 91.79776 $out [1] 107.5 105.4 $group [1] 2 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(77.04, 83.49, 88.23, 96.34, 112.52, 77.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1jfix1257968487.ps tmp/1jfix1257968487.png") > system("convert tmp/2d2ju1257968487.ps tmp/2d2ju1257968487.png") > system("convert tmp/3zuez1257968487.ps tmp/3zuez1257968487.png") > system("convert tmp/4yc2k1257968487.ps tmp/4yc2k1257968487.png") > system("convert tmp/54k1w1257968487.ps tmp/54k1w1257968487.png") > system("convert tmp/6pxir1257968487.ps tmp/6pxir1257968487.png") > system("convert tmp/7l2c41257968487.ps tmp/7l2c41257968487.png") > > > proc.time() user system elapsed 1.874 1.468 2.225