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(100 + ,211.3043478 + ,333.9130435 + ,187.826087 + ,254.7826087 + ,121.7391304 + ,113.9130435 + ,140.8695652 + ,142.6086957 + ,173.9130435 + ,121.7391304 + ,113.9130435 + ,115.6521739 + ,155.6521739 + ,256.5217391 + ,217.3913043 + ,230.4347826 + ,171.3043478 + ,155.6521739 + ,133.0434783 + ,140.8695652 + ,169.5652174 + ,183.4782609 + ,173.9130435 + ,176.5217391 + ,176.5217391 + ,306.9565217 + ,247.826087 + ,171.3043478 + ,140.8695652 + ,171.3043478 + ,153.0434783 + ,184.3478261 + ,236.5217391 + ,177.3913043 + ,164.3478261 + ,208.6956522 + ,211.3043478 + ,403.4782609 + ,328.6956522 + ,261.7391304 + ,254.7826087 + ,277.3913043 + ,282.6086957 + ,357.3913043 + ,284.3478261 + ,183.4782609 + ,186.0869565 + ,179.1304348 + ,175.6521739 + ,317.3913043 + ,375.6521739 + ,277.3913043 + ,313.9130435 + ,203.4782609 + ,199.1304348 + ,240.8695652 + ,238.2608696 + ,258.2608696 + ,335.6521739 + ,302.6086957 + ,373.9130435 + ,570.4347826 + ,382.6086957 + ,307.826087 + ,435.6521739 + ,306.0869565 + ,250.4347826 + ,350.4347826 + ,426.0869565 + ,477.3913043 + ,437.3913043 + ,334.7826087 + ,305.2173913 + ,571.3043478 + ,523.4782609 + ,426.9565217 + ,360.8695652 + ,267.826087 + ,254.7826087 + ,357.3913043 + ,348.6956522 + ,259.1304348 + ,290.4347826 + ,329.5652174 + ,304.3478261 + ,566.9565217 + ,409.5652174 + ,345.2173913 + ,445.2173913 + ,315.6521739 + ,270.4347826 + ,306.0869565 + ,326.9565217 + ,313.0434783 + ,300 + ,360 + ,334.7826087 + ,384.3478261 + ,531.3043478 + ,325.2173913 + ,399.1304348 + ,365.2173913 + ,216.5217391 + ,333.9130435 + ,348.6956522 + ,180.8695652 + ,292.173913) > 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] 108 > (np <- floor(n / par1)) [1] 9 > 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] 9 9 9 9 9 9 9 9 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 100.0000 115.6522 176.5217 208.6957 179.1304 302.6087 334.7826 329.5652 [2,] 211.3043 155.6522 176.5217 211.3043 175.6522 373.9130 305.2174 304.3478 [3,] 333.9130 256.5217 306.9565 403.4783 317.3913 570.4348 571.3043 566.9565 [4,] 187.8261 217.3913 247.8261 328.6957 375.6522 382.6087 523.4783 409.5652 [5,] 254.7826 230.4348 171.3043 261.7391 277.3913 307.8261 426.9565 345.2174 [6,] 121.7391 171.3043 140.8696 254.7826 313.9130 435.6522 360.8696 445.2174 [7,] 113.9130 155.6522 171.3043 277.3913 203.4783 306.0870 267.8261 315.6522 [8,] 140.8696 133.0435 153.0435 282.6087 199.1304 250.4348 254.7826 270.4348 [9,] 142.6087 140.8696 184.3478 357.3913 240.8696 350.4348 357.3913 306.0870 [10,] 173.9130 169.5652 236.5217 284.3478 238.2609 426.0870 348.6957 326.9565 [11,] 121.7391 183.4783 177.3913 183.4783 258.2609 477.3913 259.1304 313.0435 [12,] 113.9130 173.9130 164.3478 186.0870 335.6522 437.3913 290.4348 300.0000 [,9] [,10] [1,] 360.0000 NA [2,] 334.7826 NA [3,] 384.3478 NA [4,] 531.3043 NA [5,] 325.2174 NA [6,] 399.1304 NA [7,] 365.2174 NA [8,] 216.5217 NA [9,] 333.9130 NA [10,] 348.6957 NA [11,] 180.8696 NA [12,] 292.1739 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 111.304348 40.000000 0.00000 2.608696 -3.478261 71.30435 [2,] 122.608696 100.869565 130.43478 192.173913 141.739130 196.52174 [3,] -146.086957 -39.130435 -59.13043 -74.782609 58.260870 -187.82609 [4,] 66.956522 13.043478 -76.52174 -66.956522 -98.260870 -74.78261 [5,] -133.043478 -59.130435 -30.43478 -6.956522 36.521739 127.82609 [6,] -7.826087 -15.652174 30.43478 22.608696 -110.434783 -129.56522 [7,] 26.956522 -22.608696 -18.26087 5.217391 -4.347826 -55.65217 [8,] 1.739130 7.826087 31.30435 74.782609 41.739130 100.00000 [9,] 31.304348 28.695652 52.17391 -73.043478 -2.608696 75.65217 [10,] -52.173913 13.913044 -59.13043 -100.869565 20.000000 51.30435 [11,] -7.826087 -9.565217 -13.04348 2.608696 77.391304 -40.00000 [12,] 1.739130 2.608696 44.34783 -6.956522 -33.043478 -102.60870 [,7] [,8] [,9] [,10] [1,] -29.565217 -25.21739 -25.21739 NA [2,] 266.086956 262.60870 49.56522 NA [3,] -47.826087 -157.39130 146.95652 NA [4,] -96.521739 -64.34783 -206.08696 NA [5,] -66.086956 100.00000 73.91304 NA [6,] -93.043478 -129.56522 -33.91304 NA [7,] -13.043478 -45.21739 -148.69565 NA [8,] 102.608696 35.65217 117.39130 NA [9,] -8.695652 20.86957 14.78261 NA [10,] -89.565217 -13.91304 -167.82609 NA [11,] 31.304348 -13.04348 111.30435 NA [12,] 39.130435 60.00000 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/1w1pg1257854877.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/2jegr1257854877.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/3lzdg1257854877.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/4lk7z1257854877.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] [1,] 100.0000 155.6522 256.5217 187.8261 171.3043 121.7391 113.9130 133.0435 [2,] 176.5217 176.5217 317.3913 247.8261 254.7826 171.3043 171.3043 153.0435 [3,] 208.6957 211.3043 384.3478 375.6522 277.3913 313.9130 267.8261 216.5217 [4,] 329.5652 305.2174 566.9565 409.5652 325.2174 399.1304 306.0870 254.7826 [5,] 360.0000 373.9130 571.3043 531.3043 426.9565 445.2174 365.2174 282.6087 [,9] [,10] [,11] [,12] [1,] 140.8696 169.5652 121.7391 113.9130 [2,] 184.3478 236.5217 180.8696 173.9130 [3,] 306.0870 284.3478 183.4783 290.4348 [4,] 350.4348 348.6957 259.1304 300.0000 [5,] 357.3913 426.0870 313.0435 437.3913 $n [1] 9 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 128.0928 143.5246 252.9101 290.4696 240.2957 193.9246 196.8406 162.9391 [2,] 289.2986 279.0841 515.7855 460.8348 314.4870 433.9014 338.8116 270.1043 [,9] [,10] [,11] [,12] [1,] 218.6145 225.2696 142.2609 224.0290 [2,] 393.5594 343.4261 224.6957 356.8406 $out [1] 477.3913 $group [1] 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(100, 176.5217391, 208.6956522, 329.5652174, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5v3181257854877.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] [1,] -29.56522 49.56522 -187.82609 -98.26087 -133.043478 -129.565217 [2,] -25.21739 122.60870 -146.08696 -96.52174 -59.130435 -110.434783 [3,] 0.00000 141.73913 -59.13043 -74.78261 -6.956522 -33.913044 [4,] 40.00000 196.52174 -39.13043 -64.34783 73.913044 -7.826087 [5,] 111.30435 266.08696 58.26087 -64.34783 127.826087 30.434783 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -55.652174 1.739130 -8.695652 -167.82609 -40.000000 -102.608696 [2,] -45.217391 31.304348 -2.608696 -89.56522 -13.043478 -20.000000 [3,] -18.260869 41.739130 20.869565 -52.17391 -7.826087 2.173913 [4,] -4.347826 100.000000 31.304348 13.91304 31.304348 41.739130 [5,] 26.956522 117.391304 75.652174 51.30435 77.391304 60.000000 $n [1] 9 9 9 9 9 9 9 9 9 9 9 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -34.34783 102.8116 -115.4609 -91.72754 -77.02609 -87.95362 -39.785507 [2,] 34.34783 180.6667 -2.8000 -57.83768 63.11304 20.12754 3.263768 [,8] [,9] [,10] [,11] [,12] [1,] 5.55942 3.008696 -106.672464 -31.18261 -32.31445 [2,] 77.91884 38.730435 2.324638 15.53043 36.66228 $out [1] 146.95652 66.95652 13.04348 -206.08696 -148.69565 -73.04348 111.30435 $group [1] 3 4 4 4 7 9 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-29.5652174000001, -25.2173913, 0, 40, 111.3043478, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ze9b1257854877.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] [,7] [,8] [1,] 100.0000 115.6522 140.8696 183.4783 175.6522 250.4348 254.7826 270.4348 [2,] 117.8261 148.2609 167.8261 210.0000 201.3043 306.9565 279.1304 305.2174 [3,] 141.7391 170.4348 176.5217 269.5652 249.5652 378.2609 341.7391 321.3043 [4,] 199.5652 200.4348 210.4348 306.5217 315.6522 436.5217 393.9130 377.3913 [5,] 254.7826 256.5217 247.8261 403.4783 375.6522 570.4348 523.4783 445.2174 [,9] [,10] [1,] 216.5217 NA [2,] 308.6957 NA [3,] 341.7391 NA [4,] 374.7826 NA [5,] 399.1304 NA $n [1] 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 104.4574 146.6379 157.0876 225.5410 197.4104 319.1653 289.3860 288.3853 [2,] 179.0209 194.2317 195.9559 313.5894 301.7200 437.3564 394.0922 354.2234 [,9] [,10] [1,] 311.5964 NA [2,] 371.8818 NA $out [1] 333.9130 306.9565 571.3043 566.9565 531.3043 180.8696 $group [1] 1 3 7 8 9 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(100, 117.82608695, 141.73913045, 199.5652174, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7nqon1257854877.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,] 211.2077 183.4783 203.9130 [2,] 240.6280 213.9130 237.8261 [3,] 261.5459 280.8696 260.2174 [4,] 291.3527 310.0000 291.3043 [5,] 356.0386 384.3478 328.6957 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 238.4100 237.0437 235.8256 [2,] 284.6817 324.6955 284.6092 $out [1] 412.3671 442.1739 $group [1] 1 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(211.207729477778, 240.628019322222, 261.545893711111, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1w1pg1257854877.ps tmp/1w1pg1257854877.png") > system("convert tmp/2jegr1257854877.ps tmp/2jegr1257854877.png") > system("convert tmp/3lzdg1257854877.ps tmp/3lzdg1257854877.png") > system("convert tmp/4lk7z1257854877.ps tmp/4lk7z1257854877.png") > system("convert tmp/5v3181257854877.ps tmp/5v3181257854877.png") > system("convert tmp/6ze9b1257854877.ps tmp/6ze9b1257854877.png") > system("convert tmp/7nqon1257854877.ps tmp/7nqon1257854877.png") > > > proc.time() user system elapsed 1.216 0.956 1.707