R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. Natural language support but running in an English locale 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(586111,604378,600991,544686,537034,551531,563250,574761,580112,575093,557560,564478,580523,596594,586570,536214,523597,536535,536322,532638,528222,516141,501866,506174,517945,533590,528379,477580,469357,490243,492622,507561,516922,514258,509846,527070,541657,564591,555362,498662,511038,525919,531673,548854,560576,557274,565742,587625,619916,625809,619567,572942,572775,574205,579799,590072,593408,597141,595404,612117,628232,628884,620735,569028,567456,573100,584428,589379,590865,595454,594167,611324,612613,610763,593530,542722,536662,543599,555332,560854,562325,554788,547344,565464,577992,579714,569323,506971,500857,509127,509933,517009,519164,512238,509239,518585,522975,525192,516847,455626,454724,461251,470439,474605,476049,471067,470984,502831,512927,509673,484015,431328,436087,442867,447988,460070,467037,460170,464196,485025,501492,520564,488180,439148,441977,456608,461935,480961,492865) > 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] 129 > (np <- floor(n / par1)) [1] 10 > 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] 11 11 11 11 11 11 11 11 11 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 586111 580523 517945 541657 619916 628232 612613 577992 522975 512927 [2,] 604378 596594 533590 564591 625809 628884 610763 579714 525192 509673 [3,] 600991 586570 528379 555362 619567 620735 593530 569323 516847 484015 [4,] 544686 536214 477580 498662 572942 569028 542722 506971 455626 431328 [5,] 537034 523597 469357 511038 572775 567456 536662 500857 454724 436087 [6,] 551531 536535 490243 525919 574205 573100 543599 509127 461251 442867 [7,] 563250 536322 492622 531673 579799 584428 555332 509933 470439 447988 [8,] 574761 532638 507561 548854 590072 589379 560854 517009 474605 460070 [9,] 580112 528222 516922 560576 593408 590865 562325 519164 476049 467037 [10,] 575093 516141 514258 557274 597141 595454 554788 512238 471067 460170 [11,] 557560 501866 509846 565742 595404 594167 547344 509239 470984 464196 [12,] 564478 506174 527070 587625 612117 611324 565464 518585 502831 485025 [,11] [1,] 501492 [2,] 520564 [3,] 488180 [4,] 439148 [5,] 441977 [6,] 456608 [7,] 461935 [8,] 480961 [9,] 492865 [10,] NA [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 18267 16071 15645 22934 5893 652 -1850 1722 2217 -3254 [2,] -3387 -10024 -5211 -9229 -6242 -8149 -17233 -10391 -8345 -25658 [3,] -56305 -50356 -50799 -56700 -46625 -51707 -50808 -62352 -61221 -52687 [4,] -7652 -12617 -8223 12376 -167 -1572 -6060 -6114 -902 4759 [5,] 14497 12938 20886 14881 1430 5644 6937 8270 6527 6780 [6,] 11719 -213 2379 5754 5594 11328 11733 806 9188 5121 [7,] 11511 -3684 14939 17181 10273 4951 5522 7076 4166 12082 [8,] 5351 -4416 9361 11722 3336 1486 1471 2155 1444 6967 [9,] -5019 -12081 -2664 -3302 3733 4589 -7537 -6926 -4982 -6867 [10,] -17533 -14275 -4412 8468 -1737 -1287 -7444 -2999 -83 4026 [11,] 6918 4308 17224 21883 16713 17157 18120 9346 31847 20829 [12,] 16045 11771 14587 32291 16115 1289 12528 4390 10096 16467 [,11] [1,] 19072 [2,] -32384 [3,] -49032 [4,] 2829 [5,] 14631 [6,] 5327 [7,] 19026 [8,] 11904 [9,] NA [10,] NA [11,] NA [12,] 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/freestat/rcomp/tmp/186mi1289562180.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/freestat/rcomp/tmp/286mi1289562180.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/freestat/rcomp/tmp/386mi1289562180.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/freestat/rcomp/tmp/4if4l1289562180.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,] 501492 509673.0 484015.0 431328 436087.0 442867 447988.0 460070.0 467037.0 [2,] 520460 529391.0 522613.0 466603 462040.5 475747 481530.5 494261.0 504893.5 [3,] 577992 579714.0 569323.0 506971 511038.0 525919 531673.0 532638.0 528222.0 [4,] 599362 607570.5 597260.5 543704 536848.0 547565 559291.0 567807.5 571218.5 [5,] 628232 628884.0 620735.0 572942 572775.0 574205 584428.0 590072.0 593408.0 [,10] [,11] [,12] [1,] 460170.0 464196 485025 [2,] 512238.0 501866 506174 [3,] 535464.5 528595 545774 [4,] 575093.0 565742 587625 [5,] 597141.0 595404 612117 $n [1] 11 11 11 11 11 11 11 11 11 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 540404 542470.2 533761.8 470241 475400.6 491705.8 494628.8 497601.3 [2,] 615580 616957.8 604884.2 543701 546675.4 560132.2 568717.2 567674.7 [,9] [,10] [,11] [,12] [1,] 496625.6 504059.6 496680 505077.8 [2,] 559818.4 566869.4 560510 586470.2 $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(501492, 520460, 577992, 599362, 628232, 509673, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5if4l1289562180.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] [1,] -3254 -17233.0 -62352.0 -12617 1430.0 -213 -3684.0 -4416.0 -12081.0 [2,] 1187 -13812.0 -56502.5 -6883 6653.5 3750 5236.5 1478.5 -6926.0 [3,] 5893 -9229.0 -51707.0 -1572 8270.0 5594 10273.0 3336.0 -5000.5 [4,] 17169 -7195.5 -50577.5 1331 14564.0 10258 13510.5 8164.0 -2664.0 [5,] 22934 -3387.0 -46625.0 12376 20886.0 11733 19026.0 11904.0 -2664.0 [,10] [,11] [,12] [1,] -17533 4308.0 1289.0 [2,] -7444 9346.0 10096.0 [3,] -2368 17190.5 13557.5 [4,] -83 20829.0 16115.0 [5,] 8468 31847.0 16467.0 $n [1] 11 11 11 11 11 11 11 11 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1720.632 -12381.021 -54529.6 -5485.050 4501.533 2493.667 6331.366 [2,] 13506.632 -6076.979 -48884.4 2341.050 12038.467 8694.333 14214.634 [,8] [,9] [,10] [,11] [,12] [1,] 151.1085 -7129.965 -6045.849 11453.14 10550.17 [2,] 6520.8915 -2871.035 1309.849 22927.86 16564.83 $out [1] -25658 -32384 3733 4589 32291 $group [1] 2 2 9 9 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-3254, 1187, 5893, 17169, 22934, -17233, -13812, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6t63o1289562180.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] [,9] [1,] 537034.0 501866 469357.0 498662.0 572775 567456.0 536662.0 500857.0 454724 [2,] 554545.5 519869 491432.5 528796.0 577002 578764.0 545471.5 509183.0 465845 [3,] 569619.5 534426 512052.0 552108.0 594406 592516.0 558093.0 514623.5 472836 [4,] 583111.5 558529 522507.5 562583.5 615842 616029.5 579497.0 544243.5 509839 [5,] 604378.0 596594 533590.0 587625.0 625809 628884.0 612613.0 579714.0 525192 [,10] [,11] [1,] 431328.0 439148 [2,] 445427.5 456608 [3,] 462183.0 480961 [4,] 484520.0 492865 [5,] 512927.0 520564 $n [1] 12 12 12 12 12 12 12 12 12 12 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 556590.4 516792.9 497878.5 536697.3 576690.8 575519 542573.7 498632.2 [2,] 582648.6 552059.1 526225.5 567518.7 612121.2 609513 573612.3 530614.8 [,9] [,10] [,11] [1,] 452770 444352.6 461865.6 [2,] 492902 480013.4 500056.4 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(537034, 554545.5, 569619.5, 583111.5, 604378, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7t63o1289562180.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,] 504687.6 506971.0 499444.2 [2,] 518123.0 527070.5 516033.4 [3,] 533433.1 532155.5 535879.4 [4,] 554193.7 557548.5 552795.9 [5,] 572704.7 579714.0 568480.8 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 516981.0 518254.3 519111.8 [2,] 549885.2 546056.7 552647.1 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(504687.636363636, 518123, 533433.081818182, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/186mi1289562180.ps tmp/186mi1289562180.png",intern=TRUE)) character(0) > try(system("convert tmp/286mi1289562180.ps tmp/286mi1289562180.png",intern=TRUE)) character(0) > try(system("convert tmp/386mi1289562180.ps tmp/386mi1289562180.png",intern=TRUE)) character(0) > try(system("convert tmp/4if4l1289562180.ps tmp/4if4l1289562180.png",intern=TRUE)) character(0) > try(system("convert tmp/5if4l1289562180.ps tmp/5if4l1289562180.png",intern=TRUE)) character(0) > try(system("convert tmp/6t63o1289562180.ps tmp/6t63o1289562180.png",intern=TRUE)) character(0) > try(system("convert tmp/7t63o1289562180.ps tmp/7t63o1289562180.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.923 1.473 2.126