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(294912,293488,290555,284736,281818,287854,316263,325412,326011,328282,317480,317539,313737,312276,309391,302950,300316,304035,333476,337698,335932,323931,313927,314485,313218,309664,302963,298989,298423,301631,329765,335083,327616,309119,295916,291413,291542,284678,276475,272566,264981,263290,296806,303598,286994,276427,266424,267153,268381,262522,255542,253158,243803,250741,280445,285257,270976,261076,255603,260376,263903,264291,263276,262572,256167,264221,293860) > 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] 67 > (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 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 294912 313737 313218 291542 268381 263903 [2,] 293488 312276 309664 284678 262522 264291 [3,] 290555 309391 302963 276475 255542 263276 [4,] 284736 302950 298989 272566 253158 262572 [5,] 281818 300316 298423 264981 243803 256167 [6,] 287854 304035 301631 263290 250741 264221 [7,] 316263 333476 329765 296806 280445 293860 [8,] 325412 337698 335083 303598 285257 NA [9,] 326011 335932 327616 286994 270976 NA [10,] 328282 323931 309119 276427 261076 NA [11,] 317480 313927 295916 266424 255603 NA [12,] 317539 314485 291413 267153 260376 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -1424 -1461 -3554 -6864 -5859 388 [2,] -2933 -2885 -6701 -8203 -6980 -1015 [3,] -5819 -6441 -3974 -3909 -2384 -704 [4,] -2918 -2634 -566 -7585 -9355 -6405 [5,] 6036 3719 3208 -1691 6938 8054 [6,] 28409 29441 28134 33516 29704 29639 [7,] 9149 4222 5318 6792 4812 NA [8,] 599 -1766 -7467 -16604 -14281 NA [9,] 2271 -12001 -18497 -10567 -9900 NA [10,] -10802 -10004 -13203 -10003 -5473 NA [11,] 59 558 -4503 729 4773 NA [12,] -3802 -1267 129 1228 3527 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/1m2um1258488443.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/2txu91258488443.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/3havn1258488443.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/4ou8q1258488443.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] [1,] 263903 262522 255542 253158 243803.0 250741.0 280445.0 285257 270976 [2,] 268381 264291 263276 262572 256167.0 263290.0 293860.0 303598 286994 [3,] 293227 289083 283515 278651 273399.5 276037.5 306534.5 325412 326011 [4,] 313218 309664 302963 298989 298423.0 301631.0 329765.0 335083 327616 [5,] 313737 312276 309391 302950 300316.0 304035.0 333476.0 337698 335932 [,10] [,11] [,12] [1,] 261076 255603 260376 [2,] 276427 266424 267153 [3,] 309119 295916 291413 [4,] 323931 313927 314485 [5,] 328282 317480 317539 $n [1] 6 6 6 6 6 6 6 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 264305.7 259815.9 257915.6 255160.9 246143 251306.3 283374.6 303164.8 [2,] 322148.3 318350.1 309114.4 302141.1 300656 300768.7 329694.4 347659.2 [,9] [,10] [,11] [,12] [1,] 297307.6 275552.8 262350.5 257968.3 [2,] 354714.4 342685.2 329481.5 324857.7 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(263903, 268381, 293227, 313218, 313737, 262522, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5n0os1258488443.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] [1,] -6864.0 -8203 -6441.0 -9355.0 -1691.0 28134 4222 -16604 -12001 -10802 [2,] -5859.0 -6980 -5819.0 -7585.0 3208.0 28409 4812 -14281 -12001 -10802 [3,] -2507.5 -4817 -3941.5 -4661.5 4877.5 29540 5318 -7467 -10567 -10004 [4,] -1424.0 -2885 -2384.0 -2634.0 6938.0 29704 6792 -1766 -9900 -10003 [5,] 388.0 -1015 -704.0 -566.0 8054.0 29704 9149 599 -9900 -10003 [,11] [,12] [1,] 59 -3802 [2,] 59 -1267 [3,] 558 129 [4,] 729 1228 [5,] 729 3527 $n [1] 6 6 6 6 6 6 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -5368.2182 -7458.407 -6157.186 -7855.055 2471.530 28704.68 3918.937 [2,] 353.2182 -2175.593 -1725.814 -1467.945 7283.470 30375.32 6717.063 [,8] [,9] [,10] [,11] [,12] [1,] -16310.067 -12051.561 -10568.571 84.57969 -1633.961 [2,] 1376.067 -9082.439 -9439.429 1031.42031 1891.961 $out [1] 33516 2271 -18497 -13203 -5473 -4503 4773 $group [1] 6 9 9 10 10 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-6864, -5859, -2507.5, -1424, 388, -8203, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6lnsz1258488443.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,] 281818.0 300316.0 291413 263290.0 243803.0 262572 [2,] 289204.5 306713.0 298706 266788.5 254350.0 262924 [3,] 305587.5 313832.0 306041 276451.0 260726.0 263903 [4,] 321475.5 328703.5 320417 289268.0 269678.5 264256 [5,] 328282.0 337698.0 335083 303598.0 285257.0 264291 $n [1] 12 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 290868.5 303802 296138.5 266197.9 253734.6 263107.6 [2,] 320306.5 323862 315943.5 286704.1 267717.4 264698.4 $out [1] 256167 293860 $group [1] 6 6 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(281818, 289204.5, 305587.5, 321475.5, 328282, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7r4km1258488443.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,] 274251.3 273399.5 276321.1 [2,] 281097.8 281083.0 282036.6 [3,] 290031.6 292320.0 290497.2 [4,] 304101.4 307826.8 303742.0 [5,] 317409.6 326011.0 319340.5 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 279539.5 280122 280597.3 [2,] 300523.7 304518 300397.2 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(274251.333333333, 281097.75, 290031.6, 304101.416666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > system("convert tmp/1m2um1258488443.ps tmp/1m2um1258488443.png") > system("convert tmp/2txu91258488443.ps tmp/2txu91258488443.png") > system("convert tmp/3havn1258488443.ps tmp/3havn1258488443.png") > system("convert tmp/4ou8q1258488443.ps tmp/4ou8q1258488443.png") > system("convert tmp/5n0os1258488443.ps tmp/5n0os1258488443.png") > system("convert tmp/6lnsz1258488443.ps tmp/6lnsz1258488443.png") > system("convert tmp/7r4km1258488443.ps tmp/7r4km1258488443.png") > > > proc.time() user system elapsed 1.227 0.964 1.625