R version 2.12.1 (2010-12-16) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(293896,295705,339828,336278,346017,351623,352478,356391,333962,336828,344530,406516,319235,314750,362781,352440,374399,367418,362980,376600,346981,349571,357797,419221,329877,324252,375221,359533,392530,377686,373303,388904,354829,369553,378740,427251,343705,345062,374186,370241,399458,379886,385254,384375,352107,351566,337330,386331,311953,301261,330481,331632,349725,346615,350251,355782,326844,341207,342127,403845,318619,315067,365498,362038,371518,364774,368462,369199,351696,361750,372533,434288) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 293896 319235 329877 343705 311953 318619 NA [2,] 295705 314750 324252 345062 301261 315067 NA [3,] 339828 362781 375221 374186 330481 365498 NA [4,] 336278 352440 359533 370241 331632 362038 NA [5,] 346017 374399 392530 399458 349725 371518 NA [6,] 351623 367418 377686 379886 346615 364774 NA [7,] 352478 362980 373303 385254 350251 368462 NA [8,] 356391 376600 388904 384375 355782 369199 NA [9,] 333962 346981 354829 352107 326844 351696 NA [10,] 336828 349571 369553 351566 341207 361750 NA [11,] 344530 357797 378740 337330 342127 372533 NA [12,] 406516 419221 427251 386331 403845 434288 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1809 -4485 -5625 1357 -10692 -3552 NA [2,] 44123 48031 50969 29124 29220 50431 NA [3,] -3550 -10341 -15688 -3945 1151 -3460 NA [4,] 9739 21959 32997 29217 18093 9480 NA [5,] 5606 -6981 -14844 -19572 -3110 -6744 NA [6,] 855 -4438 -4383 5368 3636 3688 NA [7,] 3913 13620 15601 -879 5531 737 NA [8,] -22429 -29619 -34075 -32268 -28938 -17503 NA [9,] 2866 2590 14724 -541 14363 10054 NA [10,] 7702 8226 9187 -14236 920 10783 NA [11,] 61986 61424 48511 49001 61718 61755 NA [12,] -87281 -89344 -83546 -74378 -85226 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/rcomp/tmp/106pw1319547345.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/rcomp/tmp/2rmk11319547345.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/rcomp/tmp/3lkrr1319547345.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/rcomp/tmp/4x7b61319547345.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] [1,] 293896 295705.0 330481.0 331632.0 346017.0 346615 350251 355782.0 326844.0 [2,] 311953 301261.0 339828.0 336278.0 349725.0 351623 352478 356391.0 333962.0 [3,] 318927 314908.5 364139.5 355986.5 372958.5 366096 365721 372899.5 349338.5 [4,] 329877 324252.0 374186.0 362038.0 392530.0 377686 373303 384375.0 352107.0 [5,] 343705 345062.0 375221.0 370241.0 399458.0 379886 385254 388904.0 354829.0 [,10] [,11] [,12] [1,] 336828.0 337330.0 386331.0 [2,] 341207.0 342127.0 403845.0 [3,] 350568.5 351163.5 412868.5 [4,] 361750.0 372533.0 427251.0 [5,] 369553.0 378740.0 434288.0 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 307365.4 300078.6 341977.5 339370.5 345347.9 349284.5 352288.2 354848.9 [2,] 330488.6 329738.4 386301.5 372602.5 400569.1 382907.5 379153.8 390950.1 [,9] [,10] [,11] [,12] [1,] 337634.4 337317.6 331550.6 397770.9 [2,] 361042.6 363819.4 370776.4 427966.1 $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(293896, 311953, 318927, 329877, 343705, 295705, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5utrw1319547345.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] [1,] -10692.0 29124 -15688.0 9480 -19572.0 -4438.0 -879 -34075.0 -541 920 [2,] -5625.0 29220 -10341.0 9739 -14844.0 -4383.0 737 -32268.0 2590 920 [3,] -4018.5 46077 -3747.5 20026 -6862.5 2245.5 4722 -29278.5 6460 7964 [4,] 1357.0 50431 -3460.0 29217 -3110.0 3688.0 13620 -22429.0 14363 9187 [5,] 1809.0 50969 1151.0 32997 5606.0 5368.0 15601 -17503.0 14724 10783 [,11] [,12] [1,] 48511 -89344 [2,] 49001 -87281 [3,] 61571 -85226 [4,] 61755 -83546 [5,] 61986 -83546 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -8522.1155 32395.22 -8185.9672 7462.061 -14431.309 -2960.556 -3587.951 [2,] 485.1155 59758.78 690.9672 32589.939 706.309 7451.556 13031.951 [,8] [,9] [,10] [,11] [,12] [1,] -35624.97 -1133.965 2631.518 53344.26 -87865.14 [2,] -22932.03 14053.965 13296.482 69797.74 -82586.86 $out [1] -14236 -74378 $group [1] 10 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-10692, -5625, -4018.5, 1357, 1809, 29124, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6tz1t1319547345.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] [,7] [1,] 333962.0 314750.0 324252 337330.0 301261.0 351696.0 NA [2,] 335120.0 348276.0 357181 348314.0 328662.5 356723.0 NA [3,] 342179.0 360289.0 374262 372213.5 341667.0 365136.0 NA [4,] 352050.5 370908.5 383822 384814.5 349988.0 370358.5 NA [5,] 356391.0 376600.0 392530 399458.0 355782.0 372533.0 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 334456.9 349966.2 362110.9 355565.4 331940.3 358916.8 NA [2,] 349901.1 370611.8 386413.1 388861.6 351393.7 371355.2 NA $out [1] 293896 295705 406516 419221 427251 403845 318619 315067 434288 $group [1] 1 1 1 2 3 5 6 6 6 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(333962, 335120, 342179, 352050.5, 356391, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7eg4t1319547345.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,] 319547.5 349338.5 320418.0 [2,] 348074.5 349953.5 347737.8 [3,] 356754.3 360063.0 357289.2 [4,] 368664.9 369497.8 368013.5 [5,] 372274.5 372958.5 371585.2 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 347362.9 351148.7 348041.4 [2,] 366145.8 368977.3 366537.1 $out [1] 316016.2 412908.7 318927.0 314908.5 412868.5 313294.5 414878.1 $group [1] 1 1 2 2 2 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(319547.5, 348074.5, 356754.333333333, 368664.916666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/106pw1319547345.ps tmp/106pw1319547345.png",intern=TRUE)) character(0) > try(system("convert tmp/2rmk11319547345.ps tmp/2rmk11319547345.png",intern=TRUE)) character(0) > try(system("convert tmp/3lkrr1319547345.ps tmp/3lkrr1319547345.png",intern=TRUE)) character(0) > try(system("convert tmp/4x7b61319547345.ps tmp/4x7b61319547345.png",intern=TRUE)) character(0) > try(system("convert tmp/5utrw1319547345.ps tmp/5utrw1319547345.png",intern=TRUE)) character(0) > try(system("convert tmp/6tz1t1319547345.ps tmp/6tz1t1319547345.png",intern=TRUE)) character(0) > try(system("convert tmp/7eg4t1319547345.ps tmp/7eg4t1319547345.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.848 0.448 2.384