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(101.3,106.3,94,102.8,102,105.1,92.4,81.4,105.8,120.3,100.7,88.8,94.3,99.9,103.4,103.3,98.8,104.2,91.2,74.7,108.5,114.5,96.9,89.6,97.1,100.3,122.6,115.4,109,129.1,102.8,96.2,127.7,128.9,126.5,119.8,113.2,114.1,134.1,130,121.8,132.1,105.3,103,117.1,126.3,138.1,119.5,138,135.5,178.6,162.2,176.9,204.9,132.2,142.5,164.3,174.9,175.4,143) > 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,] 101.3 94.3 97.1 113.2 138.0 NA [2,] 106.3 99.9 100.3 114.1 135.5 NA [3,] 94.0 103.4 122.6 134.1 178.6 NA [4,] 102.8 103.3 115.4 130.0 162.2 NA [5,] 102.0 98.8 109.0 121.8 176.9 NA [6,] 105.1 104.2 129.1 132.1 204.9 NA [7,] 92.4 91.2 102.8 105.3 132.2 NA [8,] 81.4 74.7 96.2 103.0 142.5 NA [9,] 105.8 108.5 127.7 117.1 164.3 NA [10,] 120.3 114.5 128.9 126.3 174.9 NA [11,] 100.7 96.9 126.5 138.1 175.4 NA [12,] 88.8 89.6 119.8 119.5 143.0 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 5.0 5.6 3.2 0.9 -2.5 NA [2,] -12.3 3.5 22.3 20.0 43.1 NA [3,] 8.8 -0.1 -7.2 -4.1 -16.4 NA [4,] -0.8 -4.5 -6.4 -8.2 14.7 NA [5,] 3.1 5.4 20.1 10.3 28.0 NA [6,] -12.7 -13.0 -26.3 -26.8 -72.7 NA [7,] -11.0 -16.5 -6.6 -2.3 10.3 NA [8,] 24.4 33.8 31.5 14.1 21.8 NA [9,] 14.5 6.0 1.2 9.2 10.6 NA [10,] -19.6 -17.6 -2.4 11.8 0.5 NA [11,] -11.9 -7.3 -6.7 -18.6 -32.4 NA [12,] 5.5 7.5 -6.6 18.5 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/rcomp/tmp/1v6kn1258324714.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/27iax1258324714.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/3f41z1258324714.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/48zsi1258324714.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,] 94.3 99.9 94.0 102.8 98.8 104.2 91.2 74.7 105.8 114.5 96.9 88.8 [2,] 97.1 100.3 103.4 103.3 102.0 105.1 92.4 81.4 108.5 120.3 100.7 89.6 [3,] 101.3 106.3 122.6 115.4 109.0 129.1 102.8 96.2 117.1 126.3 126.5 119.5 [4,] 113.2 114.1 134.1 130.0 121.8 132.1 105.3 103.0 127.7 128.9 138.1 119.8 [5,] 113.2 114.1 178.6 162.2 121.8 132.1 105.3 103.0 127.7 128.9 175.4 143.0 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 89.92378 96.54895 100.9075 96.53385 95.00937 110.0219 93.6849 [2,] 112.67622 116.05105 144.2925 134.26615 122.99063 148.1781 111.9151 [,8] [,9] [,10] [,11] [,12] [1,] 80.9375 103.5333 120.2233 100.0733 98.16076 [2,] 111.4625 130.6667 132.3767 152.9267 140.83924 $out [1] 138.0 135.5 176.9 204.9 132.2 142.5 164.3 174.9 $group [1] 1 2 5 6 7 8 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(94.3, 97.1, 101.3, 113.2, 113.2, 99.9, 100.3, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5a3k21258324714.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.5 -12.3 -16.4 -8.2 3.1 -26.8 -16.5 14.1 1.2 -19.6 -32.4 -6.60 [2,] 0.9 3.5 -7.2 -6.4 5.4 -26.8 -11.0 21.8 6.0 -17.6 -18.6 -0.55 [3,] 3.2 20.0 -4.1 -4.5 10.3 -26.3 -6.6 24.4 9.2 -2.4 -11.9 6.50 [4,] 5.0 22.3 -0.1 -0.8 20.1 -13.0 -2.3 31.5 10.6 0.5 -7.3 13.00 [5,] 5.6 43.1 8.8 -0.8 28.0 -12.7 10.3 33.8 14.5 11.8 -6.7 18.50 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.3029503 6.715967 -9.1168421 -8.4569459 -0.08698297 -36.05105 [2,] 6.0970497 33.284033 0.9168421 -0.5430541 20.68698297 -16.54895 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -12.7473981 17.54600 5.949652 -15.18941 -19.884552 -4.2045 [2,] -0.4526019 31.25400 12.450348 10.38941 -3.915448 17.2045 $out [1] 14.7 -72.7 $group [1] 4 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-2.5, 0.899999999999991, 3.2, 5, 5.60000000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ezq61258324714.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,] 81.40 89.60 96.20 103.00 132.20 NA [2,] 93.20 92.75 101.55 113.65 140.25 NA [3,] 101.65 99.35 117.60 120.65 163.25 NA [4,] 105.45 103.80 127.10 131.05 176.15 NA [5,] 120.30 114.50 129.10 138.10 204.90 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 96.0627 94.31002 105.9465 112.7137 146.8758 NA [2,] 107.2373 104.38998 129.2535 128.5863 179.6242 NA $out [1] 74.7 $group [1] 2 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(81.4, 93.2, 101.65, 105.45, 120.3, 89.6, 92.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7p1ph1258324714.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,] 99.56 96.20 92.200 [2,] 110.00 104.55 104.925 [3,] 122.22 116.25 114.275 [4,] 127.03 124.45 118.675 [5,] 135.08 129.10 124.600 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 114.4525 107.1735 108.0035 [2,] 129.9875 125.3265 120.5465 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(99.56, 110, 122.22, 127.03, 135.08, 96.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > system("convert tmp/1v6kn1258324714.ps tmp/1v6kn1258324714.png") > system("convert tmp/27iax1258324714.ps tmp/27iax1258324714.png") > system("convert tmp/3f41z1258324714.ps tmp/3f41z1258324714.png") > system("convert tmp/48zsi1258324714.ps tmp/48zsi1258324714.png") > system("convert tmp/5a3k21258324714.ps tmp/5a3k21258324714.png") > system("convert tmp/6ezq61258324714.ps tmp/6ezq61258324714.png") > system("convert tmp/7p1ph1258324714.ps tmp/7p1ph1258324714.png") > > > proc.time() user system elapsed 1.194 0.931 1.645