R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(96.8,114.1,110.3,103.9,101.6,94.6,95.9,104.7,102.8,98.1,113.9,80.9,95.7,113.2,105.9,108.8,102.3,99,100.7,115.5,100.7,109.9,114.6,85.4,100.5,114.8,116.5,112.9,102,106,105.3,118.8,106.1,109.3,117.2,92.5,104.2,112.5,122.4,113.3,100,110.7,112.8,109.8,117.3,109.1,115.9,96,99.8,116.8,115.7,99.4,94.3,91,93.2,103.1,94.1,91.8,102.7,82.6,89.1) > 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,] 96.8 95.7 100.5 104.2 99.8 89.1 [2,] 114.1 113.2 114.8 112.5 116.8 NA [3,] 110.3 105.9 116.5 122.4 115.7 NA [4,] 103.9 108.8 112.9 113.3 99.4 NA [5,] 101.6 102.3 102.0 100.0 94.3 NA [6,] 94.6 99.0 106.0 110.7 91.0 NA [7,] 95.9 100.7 105.3 112.8 93.2 NA [8,] 104.7 115.5 118.8 109.8 103.1 NA [9,] 102.8 100.7 106.1 117.3 94.1 NA [10,] 98.1 109.9 109.3 109.1 91.8 NA [11,] 113.9 114.6 117.2 115.9 102.7 NA [12,] 80.9 85.4 92.5 96.0 82.6 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 17.3 17.5 14.3 8.3 17.0 NA [2,] -3.8 -7.3 1.7 9.9 -1.1 NA [3,] -6.4 2.9 -3.6 -9.1 -16.3 NA [4,] -2.3 -6.5 -10.9 -13.3 -5.1 NA [5,] -7.0 -3.3 4.0 10.7 -3.3 NA [6,] 1.3 1.7 -0.7 2.1 2.2 NA [7,] 8.8 14.8 13.5 -3.0 9.9 NA [8,] -1.9 -14.8 -12.7 7.5 -9.0 NA [9,] -4.7 9.2 3.2 -8.2 -2.3 NA [10,] 15.8 4.7 7.9 6.8 10.9 NA [11,] -33.0 -29.2 -24.7 -19.9 -20.1 NA [12,] 14.8 15.1 11.7 3.8 6.5 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/1atqa1258124341.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/2ztq71258124341.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/3tjn01258124341.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/4am9n1258124341.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,] 89.1 112.5 105.9 99.4 100.0 91.0 93.2 103.1 94.1 91.8 113.9 80.9 [2,] 95.7 113.2 110.3 103.9 100.0 94.6 95.9 104.7 100.7 98.1 113.9 82.6 [3,] 98.3 114.1 115.7 108.8 101.6 99.0 100.7 109.8 102.8 109.1 114.6 85.4 [4,] 100.5 114.8 116.5 112.9 102.0 106.0 105.3 115.5 106.1 109.3 115.9 92.5 [5,] 104.2 116.8 122.4 113.3 102.3 110.7 112.8 118.8 106.1 109.9 117.2 96.0 $n [1] 6 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 95.20384 112.9694 111.3191 102.4406 100.1868 90.94479 94.05798 102.1687 [2,] 101.39616 115.2306 120.0809 115.1594 103.0132 107.05521 107.34202 117.4313 [,9] [,10] [,11] [,12] [1,] 98.98437 101.1861 113.1868 78.40468 [2,] 106.61563 117.0139 116.0132 92.39532 $out [1] 94.3 117.3 102.7 $group [1] 5 9 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(89.1, 95.7, 98.3, 100.5, 104.2, 112.5, 113.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/532rs1258124341.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,] 14.3 -7.3 -16.3 -13.3 -7.0 1.3 8.8 -14.8 -8.2 4.7 -33.0 3.8 [2,] 14.3 -3.8 -9.1 -10.9 -3.3 1.3 8.8 -12.7 -4.7 6.8 -29.2 6.5 [3,] 17.0 -1.1 -6.4 -6.5 -3.3 1.7 9.9 -9.0 -2.3 7.9 -24.7 11.7 [4,] 17.3 1.7 -3.6 -5.1 4.0 2.1 13.5 -1.9 3.2 10.9 -20.1 14.8 [5,] 17.5 9.9 2.9 -2.3 10.7 2.2 14.8 7.5 9.2 15.8 -19.9 15.1 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 14.88021 -4.986286 -10.286286 -10.598265 -8.458162 1.134722 6.578992 [2,] 19.11979 2.786286 -2.513714 -2.401735 1.858162 2.265278 13.221008 [,8] [,9] [,10] [,11] [,12] [1,] -16.631253 -7.88212 5.00295 -31.13004 5.835241 [2,] -1.368747 3.28212 10.79705 -18.26996 17.564759 $out [1] 8.3 -0.7 -3.0 $group [1] 1 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(14.3, 14.3, 17, 17.3, 17.5, -7.3, -3.8, -1.09999999999999, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/67rsn1258124341.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,] 80.90 85.40 92.50 96.00 82.60 89.1 [2,] 96.35 99.85 103.65 106.65 92.50 89.1 [3,] 102.20 104.10 107.70 111.60 96.85 89.1 [4,] 107.50 111.55 115.65 114.60 102.90 89.1 [5,] 114.10 115.50 118.80 122.40 116.80 89.1 $n [1] 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 97.11441 98.76355 102.2267 107.9740 92.10649 89.1 [2,] 107.28559 109.43645 113.1733 115.2260 101.59351 89.1 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(80.9, 96.35, 102.2, 107.5, 114.1, 85.4, 99.85, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/77p1x1258124341.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,] 87.48 85.40 87.55 [2,] 100.15 99.85 100.45 [3,] 103.92 105.80 103.55 [4,] 111.62 111.95 111.75 [5,] 114.28 115.70 114.90 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 98.68846 100.2811 98.396 [2,] 109.15154 111.3189 108.704 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(87.48, 100.15, 103.92, 111.62, 114.28, 85.4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1atqa1258124341.ps tmp/1atqa1258124341.png") > system("convert tmp/2ztq71258124341.ps tmp/2ztq71258124341.png") > system("convert tmp/3tjn01258124341.ps tmp/3tjn01258124341.png") > system("convert tmp/4am9n1258124341.ps tmp/4am9n1258124341.png") > system("convert tmp/532rs1258124341.ps tmp/532rs1258124341.png") > system("convert tmp/67rsn1258124341.ps tmp/67rsn1258124341.png") > system("convert tmp/77p1x1258124341.ps tmp/77p1x1258124341.png") > > > proc.time() user system elapsed 1.219 0.954 1.708