R version 2.7.2 (2008-08-25) 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(109.1,113.8,97.4,72.5,82.7,88.9,105.9,100.8,94,105,58.5,87.6,113.1,112.5,89.6,74.5,82.7,90.1,109.4,96,89.2,109.1,49.1,92.9,107.7,103.5,91.1,79.8,71.9,82.9,90.1,100.7,90.7,108.8,44.1,93.6,107.4,96.5,93.6,76.5,76.7,84,103.3,88.5,99,105.9,44.7,94,107.1,104.8,102.5,77.7,85.2,91.3,106.5,92.4,97.5,107,51.1,98.6,102.2,114.3,99.4,72.5,92.3,99.4,85.9,109.4,97.6,104.7) > 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] 70 > (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 6 6 6 6 6 6 6 6 6 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 109.1 113.1 107.7 107.4 107.1 102.2 [2,] 113.8 112.5 103.5 96.5 104.8 114.3 [3,] 97.4 89.6 91.1 93.6 102.5 99.4 [4,] 72.5 74.5 79.8 76.5 77.7 72.5 [5,] 82.7 82.7 71.9 76.7 85.2 92.3 [6,] 88.9 90.1 82.9 84.0 91.3 99.4 [7,] 105.9 109.4 90.1 103.3 106.5 85.9 [8,] 100.8 96.0 100.7 88.5 92.4 109.4 [9,] 94.0 89.2 90.7 99.0 97.5 97.6 [10,] 105.0 109.1 108.8 105.9 107.0 104.7 [11,] 58.5 49.1 44.1 44.7 51.1 NA [12,] 87.6 92.9 93.6 94.0 98.6 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 4.7 -0.6 -4.2 -10.9 -2.3 12.1 [2,] -16.4 -22.9 -12.4 -2.9 -2.3 -14.9 [3,] -24.9 -15.1 -11.3 -17.1 -24.8 -26.9 [4,] 10.2 8.2 -7.9 0.2 7.5 19.8 [5,] 6.2 7.4 11.0 7.3 6.1 7.1 [6,] 17.0 19.3 7.2 19.3 15.2 -13.5 [7,] -5.1 -13.4 10.6 -14.8 -14.1 23.5 [8,] -6.8 -6.8 -10.0 10.5 5.1 -11.8 [9,] 11.0 19.9 18.1 6.9 9.5 7.1 [10,] -46.5 -60.0 -64.7 -61.2 -55.9 NA [11,] 29.1 43.8 49.5 49.3 47.5 NA [12,] 25.5 14.8 13.8 13.1 3.6 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/rcomp/tmp/12mbb1225740070.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/rcomp/tmp/2ijud1225740070.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/rcomp/tmp/3qlv31225740070.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/rcomp/tmp/41un21225740070.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,] 107.10 96.50 89.6 72.5 71.9 82.9 85.9 88.50 89.20 104.70 44.1 92.9 [2,] 107.10 103.50 91.1 72.5 76.7 84.0 90.1 92.40 90.70 105.00 44.7 92.9 [3,] 107.55 108.65 95.5 75.5 82.7 89.5 104.6 98.35 95.75 106.45 49.1 93.6 [4,] 109.10 113.80 99.4 77.7 85.2 91.3 106.5 100.80 97.60 108.80 51.1 94.0 [5,] 109.10 114.30 102.5 79.8 92.3 99.4 109.4 109.40 99.00 109.10 58.5 94.0 $n [1] 6 6 6 6 6 6 6 6 6 6 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 106.2599 102.0062 90.14623 72.14583 77.21723 84.79126 94.02147 92.93173 [2,] 108.8401 115.2938 100.85377 78.85417 88.18277 94.20874 115.17853 103.76827 [,9] [,10] [,11] [,12] [1,] 91.29928 103.9989 44.57778 92.82274 [2,] 100.20072 108.9011 53.62222 94.37726 $out [1] 113.1 102.2 87.6 98.6 $group [1] 1 1 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(107.1, 107.1, 107.55, 109.1, 109.1, 96.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5wrks1225740070.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,] -10.90 -22.90 -26.90 -7.90 6.1 7.2 -14.80 -11.8 6.90 -64.7 43.8 13.1 [2,] -4.20 -16.40 -24.90 0.20 6.2 7.2 -14.10 -10.0 7.10 -61.2 43.8 13.1 [3,] -1.45 -13.65 -20.95 7.85 7.2 16.1 -9.25 -6.8 10.25 -60.0 47.5 13.8 [4,] 4.70 -2.90 -15.10 10.20 7.4 19.3 10.60 5.1 18.10 -55.9 49.3 14.8 [5,] 12.10 -2.30 -11.30 19.80 7.4 19.3 23.50 10.5 19.90 -55.9 49.5 14.8 $n [1] 6 6 6 6 6 6 6 6 6 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -7.190787 -22.357936 -27.27132 1.399677 6.425961 8.29511 -25.182298 [2,] 4.290787 -4.942064 -14.62868 14.300323 7.974039 23.90489 6.682298 [,8] [,9] [,10] [,11] [,12] [1,] -16.539988 3.154645 -63.74497 43.61371 12.59878 [2,] 2.939988 17.345355 -56.25503 51.38629 15.00122 $out [1] 11.0 -13.5 -46.5 29.1 25.5 3.6 $group [1] 5 6 10 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-10.9, -4.2, -1.45000000000000, 4.7, 12.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6fj521225740070.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,] 58.50 74.50 71.90 76.50 77.70 85.9 [2,] 85.15 85.95 81.35 80.35 88.25 92.3 [3,] 95.70 91.50 90.90 93.80 98.05 99.4 [4,] 105.45 109.25 102.10 101.15 105.65 104.7 [5,] 113.80 113.10 108.80 107.40 107.10 114.3 $n [1] 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 86.44103 80.87271 81.43579 84.31298 90.11374 93.20447 [2,] 104.95897 102.12729 100.36421 103.28702 105.98626 105.59553 $out [1] 49.1 44.1 44.7 51.1 72.5 $group [1] 2 3 4 5 6 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(58.5, 85.15, 95.7, 105.45, 113.8, 74.5, 85.95, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7tx7u1225740070.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,] 75.58333 75.500 75.20000 [2,] 85.67500 86.100 84.75000 [3,] 95.13333 95.625 94.93125 [4,] 103.46667 105.525 103.33125 [5,] 107.76667 108.650 108.65000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 87.01843 86.76513 86.45622 [2,] 103.24823 104.48487 103.40628 $out [1] 49.5 49.1 47.9 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(75.5833333333333, 85.675, 95.1333333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/12mbb1225740070.ps tmp/12mbb1225740070.png") > system("convert tmp/2ijud1225740070.ps tmp/2ijud1225740070.png") > system("convert tmp/3qlv31225740070.ps tmp/3qlv31225740070.png") > system("convert tmp/41un21225740070.ps tmp/41un21225740070.png") > system("convert tmp/5wrks1225740070.ps tmp/5wrks1225740070.png") > system("convert tmp/6fj521225740070.ps tmp/6fj521225740070.png") > system("convert tmp/7tx7u1225740070.ps tmp/7tx7u1225740070.png") > > > proc.time() user system elapsed 1.451 0.982 1.767