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. Natural language support but running in an English locale 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(112.3,117.3,111.1,102.2,104.3,122.9,107.6,121.3,131.5,89,104.4,128.9,135.9,133.3,121.3,120.5,120.4,137.9,126.1,133.2,151.1,105,119,140.4,156.6,137.1,122.7,125.8,139.3,134.9,149.2,132.3,149,117.2,119.6,152,149.4,127.3,114.1,102.1,107.7,104.4,102.1,96,109.3,90,83.9,112,114.3,103.6,91.7,80.8,87.2,109.2,102.7,95.1,117.5,85.1,92.1,113.5) > 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,] 112.3 135.9 156.6 149.4 114.3 NA [2,] 117.3 133.3 137.1 127.3 103.6 NA [3,] 111.1 121.3 122.7 114.1 91.7 NA [4,] 102.2 120.5 125.8 102.1 80.8 NA [5,] 104.3 120.4 139.3 107.7 87.2 NA [6,] 122.9 137.9 134.9 104.4 109.2 NA [7,] 107.6 126.1 149.2 102.1 102.7 NA [8,] 121.3 133.2 132.3 96.0 95.1 NA [9,] 131.5 151.1 149.0 109.3 117.5 NA [10,] 89.0 105.0 117.2 90.0 85.1 NA [11,] 104.4 119.0 119.6 83.9 92.1 NA [12,] 128.9 140.4 152.0 112.0 113.5 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 5.0 -2.6 -19.5 -22.1 -10.7 NA [2,] -6.2 -12.0 -14.4 -13.2 -11.9 NA [3,] -8.9 -0.8 3.1 -12.0 -10.9 NA [4,] 2.1 -0.1 13.5 5.6 6.4 NA [5,] 18.6 17.5 -4.4 -3.3 22.0 NA [6,] -15.3 -11.8 14.3 -2.3 -6.5 NA [7,] 13.7 7.1 -16.9 -6.1 -7.6 NA [8,] 10.2 17.9 16.7 13.3 22.4 NA [9,] -42.5 -46.1 -31.8 -19.3 -32.4 NA [10,] 15.4 14.0 2.4 -6.1 7.0 NA [11,] 24.5 21.4 32.4 28.1 21.4 NA [12,] 7.0 16.2 -2.6 2.3 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/1yi6d1291147008.ps",horizontal=F,onefile=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/2yi6d1291147008.ps",horizontal=F,onefile=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/3q9ny1291147008.ps",horizontal=F,onefile=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/4q9ny1291147008.ps",horizontal=F,onefile=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,] 112.3 103.6 111.1 80.8 87.2 104.4 102.1 95.1 109.3 85.1 83.9 112.0 [2,] 114.3 117.3 111.1 102.1 104.3 109.2 102.7 96.0 117.5 89.0 92.1 113.5 [3,] 135.9 127.3 114.1 102.2 107.7 122.9 107.6 121.3 131.5 90.0 104.4 128.9 [4,] 149.4 133.3 121.3 120.5 120.4 134.9 126.1 132.3 149.0 105.0 119.0 140.4 [5,] 156.6 137.1 122.7 125.8 139.3 137.9 149.2 133.2 151.1 117.2 119.6 152.0 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 111.0984 115.9944 106.8927 89.1986 96.32378 104.7404 91.06562 95.65051 [2,] 160.7016 138.6056 121.3073 115.2014 119.07622 141.0596 124.13438 146.94949 [,9] [,10] [,11] [,12] [1,] 109.2422 78.69444 85.39253 109.8925 [2,] 153.7578 101.30556 123.40747 147.9075 $out [1] 91.7 $group [1] 3 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(112.3, 114.3, 135.9, 149.4, 156.6, 103.6, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5q9ny1291147008.ps",horizontal=F,onefile=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,] -22.1 -14.4 -12.0 -0.1 -4.4 -15.3 -16.9 10.2 -46.1 -6.1 21.4 -2.60 [2,] -19.5 -13.2 -10.9 2.1 -3.3 -11.8 -7.6 13.3 -42.5 2.4 21.4 -0.15 [3,] -10.7 -12.0 -8.9 5.6 17.5 -6.5 -6.1 16.7 -32.4 7.0 24.5 4.65 [4,] -2.6 -11.9 -0.8 6.4 18.6 -2.3 7.1 17.9 -31.8 14.0 28.1 11.60 [5,] 5.0 -11.9 3.1 6.4 22.0 -2.3 13.7 22.4 -19.3 15.4 32.4 16.20 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -22.641497 -12.91858 -16.036635 2.561631 2.025515 -13.2126761 -16.486983 [2,] 1.241497 -11.08142 -1.763365 8.638369 32.974485 0.2126761 4.286983 [,8] [,9] [,10] [,11] [,12] [1,] 13.44965 -39.96059 -1.196531 19.76580 -4.6325 [2,] 19.95035 -24.83941 15.196531 29.23420 13.9325 $out [1] -6.2 13.5 14.3 $group [1] 2 4 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-22.1, -19.5, -10.7, -2.59999999999999, 5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6j0mj1291147008.ps",horizontal=F,onefile=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,] 89.00 105.00 117.20 83.90 80.80 NA [2,] 104.35 120.45 124.25 99.05 89.45 NA [3,] 111.70 129.65 136.00 106.05 98.90 NA [4,] 122.10 136.90 149.10 113.05 111.35 NA [5,] 131.50 151.10 156.60 127.30 117.50 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 103.6041 122.1470 124.6657 99.6645 88.91126 NA [2,] 119.7959 137.1530 147.3343 112.4355 108.88874 NA $out [1] 149.4 $group [1] 4 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(89, 104.35, 111.7, 122.1, 131.5, 105, 120.45, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7j0mj1291147008.ps",horizontal=F,onefile=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,] 97.26 90.0 97.000 [2,] 109.03 106.0 111.825 [3,] 116.56 117.7 115.300 [4,] 126.54 128.1 126.125 [5,] 133.70 135.9 133.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 108.5736 107.6200 108.7777 [2,] 124.5464 127.7800 121.8223 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(97.26, 109.03, 116.56, 126.54, 133.7, 90, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1yi6d1291147008.ps tmp/1yi6d1291147008.png",intern=TRUE)) character(0) > try(system("convert tmp/2yi6d1291147008.ps tmp/2yi6d1291147008.png",intern=TRUE)) character(0) > try(system("convert tmp/3q9ny1291147008.ps tmp/3q9ny1291147008.png",intern=TRUE)) character(0) > try(system("convert tmp/4q9ny1291147008.ps tmp/4q9ny1291147008.png",intern=TRUE)) character(0) > try(system("convert tmp/5q9ny1291147008.ps tmp/5q9ny1291147008.png",intern=TRUE)) character(0) > try(system("convert tmp/6j0mj1291147008.ps tmp/6j0mj1291147008.png",intern=TRUE)) character(0) > try(system("convert tmp/7j0mj1291147008.ps tmp/7j0mj1291147008.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.944 1.467 2.189