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. 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(510401 + ,508642 + ,501805 + ,495476 + ,490336 + ,490042 + ,553155 + ,569999 + ,573170 + ,571687 + ,575453 + ,580177 + ,579849 + ,574346 + ,563325 + ,555604 + ,545544 + ,545109 + ,605181 + ,627856 + ,631421 + ,625671 + ,613577 + ,606463 + ,601676 + ,589121 + ,573559 + ,558487 + ,552148 + ,545720 + ,606569 + ,636067 + ,630704 + ,623275 + ,617771 + ,605401 + ,619393 + ,596019 + ,569977 + ,546213 + ,528492 + ,505944 + ,554910 + ,567831 + ,564021 + ,552800 + ,541102 + ,542378 + ,540380 + ,521219 + ,504652 + ,490626 + ,481686 + ,477930 + ,522605 + ,531432 + ,532355 + ,539954 + ,524987 + ,533307 + ,530541 + ,508392 + ,495208 + ,482223 + ,470495 + ,466106 + ,515037 + ,517752 + ,515565 + ,510727 + ,499725 + ,498369 + ,493756 + ,476141 + ,458458 + ,443182 + ,429597 + ,424476 + ,476257 + ,480555 + ,469762 + ,459820 + ,451028 + ,450065 + ,444385 + ,428846 + ,421020 + ,399778 + ,389005 + ,384018 + ,431933 + ,445844 + ,431464 + ,423263 + ,415881 + ,416208 + ,413491 + ,399153 + ,385939 + ,373917 + ,364635 + ,364696 + ,418358 + ,428212 + ,423730 + ,420677 + ,417428 + ,423245 + ,423113 + ,418873 + ,405733 + ,397812 + ,389918 + ,391116 + ,443814 + ,460373 + ,455422 + ,456288 + ,452233 + ,459256 + ,461146 + ,451391 + ,443101 + ,438810 + ,430457 + ,435721 + ,488280 + ,505814 + ,502338 + ,500910 + ,501434 + ,515476 + ,520862 + ,519517 + ,511805 + ,508607 + ,505327 + ,511435 + ,570158 + ,591665 + ,593572 + ,586346 + ,586063 + ,591504 + ,594033 + ,585597 + ,572450 + ,562917 + ,554675 + ,553997 + ,601310 + ,622255 + ,616735 + ,606480 + ,595079 + ,598588 + ,599917 + ,591573 + ,575489 + ,567223 + ,555338 + ,555252 + ,608249 + ,630859 + ,628632 + ,624435 + ,609670 + ,615830 + ,621170 + ,604212 + ,584348 + ,573717 + ,555234 + ,544897 + ,598866 + ,620081 + ,607699 + ,589960 + ,578665 + ,580166 + ,579457 + ,571560 + ,560460 + ,551397 + ,536763 + ,540562 + ,588184 + ,607049 + ,598968 + ,577644 + ,562640 + ,565867 + ,561274 + ,554144 + ,539900 + ,526271 + ,511841 + ,505282 + ,554083 + ,584225 + ,568858 + ,539516 + ,521612 + ,525562 + ,526519 + ,515713 + ,503454 + ,489301 + ,479020 + ,475102 + ,523682 + ,551528 + ,531626 + ,511037 + ,492417 + ,492188 + ,492865 + ,480961 + ,461935 + ,456608 + ,441977 + ,439148 + ,488180 + ,520564 + ,501492 + ,485025 + ,464196 + ,460170 + ,467037 + ,460070 + ,447988 + ,442867 + ,436087 + ,431328 + ,484015 + ,509673 + ,512927 + ,502831 + ,470984 + ,471067 + ,476049 + ,474605 + ,470439 + ,461251 + ,454724 + ,455626 + ,516847 + ,525192 + ,522975 + ,518585 + ,509239 + ,512238 + ,519164 + ,517009 + ,509933 + ,509127 + ,500857 + ,506971 + ,569323 + ,579714 + ,577992 + ,565464 + ,547344 + ,554788 + ,562325 + ,560854 + ,555332 + ,543599 + ,536662 + ,542722 + ,593530 + ,610763 + ,612613 + ,611324 + ,594167 + ,595454 + ,590865 + ,589379 + ,584428 + ,573100 + ,567456 + ,569028 + ,620735 + ,628884 + ,628232 + ,612117 + ,595404 + ,597141 + ,593408 + ,590072 + ,579799 + ,574205 + ,572775 + ,572942 + ,619567 + ,625809 + ,619916 + ,587625 + ,565742 + ,557274 + ,560576 + ,548854 + ,531673 + ,525919 + ,511038 + ,498662 + ,555362 + ,564591 + ,541657 + ,527070 + ,509846 + ,514258 + ,516922 + ,507561 + ,492622 + ,490243 + ,469357 + ,477580 + ,528379 + ,533590 + ,517945 + ,506174) > 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] 322 > (np <- floor(n / par1)) [1] 26 > 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] 27 27 27 27 27 27 27 27 27 27 26 26 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 510401 579849 601676 619393 540380 530541 493756 444385 413491 423113 [2,] 508642 574346 589121 596019 521219 508392 476141 428846 399153 418873 [3,] 501805 563325 573559 569977 504652 495208 458458 421020 385939 405733 [4,] 495476 555604 558487 546213 490626 482223 443182 399778 373917 397812 [5,] 490336 545544 552148 528492 481686 470495 429597 389005 364635 389918 [6,] 490042 545109 545720 505944 477930 466106 424476 384018 364696 391116 [7,] 553155 605181 606569 554910 522605 515037 476257 431933 418358 443814 [8,] 569999 627856 636067 567831 531432 517752 480555 445844 428212 460373 [9,] 573170 631421 630704 564021 532355 515565 469762 431464 423730 455422 [10,] 571687 625671 623275 552800 539954 510727 459820 423263 420677 456288 [11,] 575453 613577 617771 541102 524987 499725 451028 415881 417428 452233 [12,] 580177 606463 605401 542378 533307 498369 450065 416208 423245 459256 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [1,] 461146 520862 594033 599917 621170 579457 561274 526519 492865 467037 [2,] 451391 519517 585597 591573 604212 571560 554144 515713 480961 460070 [3,] 443101 511805 572450 575489 584348 560460 539900 503454 461935 447988 [4,] 438810 508607 562917 567223 573717 551397 526271 489301 456608 442867 [5,] 430457 505327 554675 555338 555234 536763 511841 479020 441977 436087 [6,] 435721 511435 553997 555252 544897 540562 505282 475102 439148 431328 [7,] 488280 570158 601310 608249 598866 588184 554083 523682 488180 484015 [8,] 505814 591665 622255 630859 620081 607049 584225 551528 520564 509673 [9,] 502338 593572 616735 628632 607699 598968 568858 531626 501492 512927 [10,] 500910 586346 606480 624435 589960 577644 539516 511037 485025 502831 [11,] 501434 586063 595079 609670 578665 562640 521612 492417 464196 470984 [12,] 515476 591504 598588 615830 580166 565867 525562 492188 460170 471067 [,21] [,22] [,23] [,24] [,25] [,26] [,27] [1,] 476049 519164 562325 590865 593408 560576 516922 [2,] 474605 517009 560854 589379 590072 548854 507561 [3,] 470439 509933 555332 584428 579799 531673 492622 [4,] 461251 509127 543599 573100 574205 525919 490243 [5,] 454724 500857 536662 567456 572775 511038 469357 [6,] 455626 506971 542722 569028 572942 498662 477580 [7,] 516847 569323 593530 620735 619567 555362 528379 [8,] 525192 579714 610763 628884 625809 564591 533590 [9,] 522975 577992 612613 628232 619916 541657 517945 [10,] 518585 565464 611324 612117 587625 527070 506174 [11,] 509239 547344 594167 595404 565742 509846 NA [12,] 512238 554788 595454 597141 557274 514258 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] -1759 -5503 -12555 -23374 -19161 -22149 -17615 -15539 -14338 -4240 [2,] -6837 -11021 -15562 -26042 -16567 -13184 -17683 -7826 -13214 -13140 [3,] -6329 -7721 -15072 -23764 -14026 -12985 -15276 -21242 -12022 -7921 [4,] -5140 -10060 -6339 -17721 -8940 -11728 -13585 -10773 -9282 -7894 [5,] -294 -435 -6428 -22548 -3756 -4389 -5121 -4987 61 1198 [6,] 63113 60072 60849 48966 44675 48931 51781 47915 53662 52698 [7,] 16844 22675 29498 12921 8827 2715 4298 13911 9854 16559 [8,] 3171 3565 -5363 -3810 923 -2187 -10793 -14380 -4482 -4951 [9,] -1483 -5750 -7429 -11221 7599 -4838 -9942 -8201 -3053 866 [10,] 3766 -12094 -5504 -11698 -14967 -11002 -8792 -7382 -3249 -4055 [11,] 4724 -7114 -12370 1276 8320 -1356 -963 327 5817 7023 [12,] -328 -4787 13992 -1998 -2766 -4613 -5680 -2717 -132 1890 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] -9755 -1345 -8436 -8344 -16958 -7897 -7130 -10806 -11904 -6967 -1444 [2,] -8290 -7712 -13147 -16084 -19864 -11100 -14244 -12259 -19026 -12082 -4166 [3,] -4291 -3198 -9533 -8266 -10631 -9063 -13629 -14153 -5327 -5121 -9188 [4,] -8353 -3280 -8242 -11885 -18483 -14634 -14430 -10281 -14631 -6780 -6527 [5,] 5264 6108 -678 -86 -10337 3799 -6559 -3918 -2829 -4759 902 [6,] 52559 58723 47313 52997 53969 47622 48801 48580 49032 52687 61221 [7,] 17534 21507 20945 22610 21215 18865 30142 27846 32384 25658 8345 [8,] -3476 1907 -5520 -2227 -12382 -8081 -15367 -19902 -19072 3254 -2217 [9,] -1428 -7226 -10255 -4197 -17739 -21324 -29342 -20589 -16467 -10096 -4390 [10,] 524 -283 -11401 -14765 -11295 -15004 -17904 -18620 -20829 -31847 -9346 [11,] 14042 5441 3509 6160 1501 3227 3950 -229 -4026 83 2999 [12,] 5386 2529 1329 5340 -709 -4593 957 677 6867 4982 6926 [,22] [,23] [,24] [,25] [,26] [,27] [1,] -2155 -1471 -1486 -3336 -11722 -9361 [2,] -7076 -5522 -4951 -10273 -17181 -14939 [3,] -806 -11733 -11328 -5594 -5754 -2379 [4,] -8270 -6937 -5644 -1430 -14881 -20886 [5,] 6114 6060 1572 167 -12376 8223 [6,] 62352 50808 51707 46625 56700 50799 [7,] 10391 17233 8149 6242 9229 5211 [8,] -1722 1850 -652 -5893 -22934 -15645 [9,] -12528 -1289 -16115 -32291 -14587 -11771 [10,] -18120 -17157 -16713 -21883 -17224 NA [11,] 7444 1287 1737 -8468 4412 NA [12,] 7537 -4589 -3733 3302 2664 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/109p51228935362.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/2v0st1228935362.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/3grz91228935362.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/44n6c1228935362.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,] 413491.0 399153.0 385939 373917.0 364635.0 364696.0 418358.0 428212 423730 [2,] 493310.5 478551.0 466187 458929.5 448350.5 447387.0 501658.5 519158 514246 [3,] 530541.0 519517.0 509933 508607.0 500857.0 498662.0 554083.0 567831 564021 [4,] 585357.0 579971.5 566651 553500.5 541153.5 543809.5 596198.0 615422 610156 [5,] 621170.0 604212.0 584428 574205.0 572775.0 572942.0 620735.0 636067 631421 [,10] [,11] [,12] [1,] 420677.0 415881.0 416208.0 [2,] 504502.5 492417.0 492188.0 [3,] 539954.0 533044.5 537842.5 [4,] 588792.5 586063.0 591504.0 [5,] 625671.0 617771.0 615830.0 $n [1] 27 27 27 27 27 27 27 27 27 27 26 26 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 502552.3 488678 479384.8 479850.7 472638.3 469342.7 525336.3 538559.9 [2,] 558529.7 550356 540481.2 537363.3 529075.7 527981.3 582829.7 597102.1 [,9] [,10] [,11] [,12] [1,] 534857.5 514323.8 504027 507068.1 [2,] 593184.5 565584.2 562062 568616.9 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5nuie1228935362.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,] -23374.0 -26042 -23764 -20886.0 -12376 44675.0 2715.0 -22934.0 -29342.0 [2,] -13446.5 -15823 -13307 -14007.5 -4873 48866.0 9028.0 -11587.5 -15351.0 [3,] -8436.0 -13140 -9188 -9282.0 -435 51781.0 16844.0 -4482.0 -9942.0 [4,] -3788.0 -8058 -5674 -6858.5 1385 55334.5 22058.5 -1187.0 -4293.5 [5,] -1345.0 -4166 -806 -1430.0 8223 63113.0 32384.0 3565.0 7599.0 [,10] [,11] [,12] [1,] -31847 -8468 -5680 [2,] -17224 -229 -2766 [3,] -11896 2368 817 [4,] -7382 5441 4982 [5,] 3766 8320 13992 $n [1] 27 27 27 27 27 27 27 27 27 26 26 26 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -11372.871 -15501.11 -11508.975 -11455.805 -2337.877 49814.12 12881.8 [2,] -5499.129 -10778.89 -6867.025 -7108.195 1467.877 53747.88 20806.2 [,8] [,9] [,10] [,11] [,12] [1,] -7644.492 -13304.267 -14945.677 611.074 -1583.822 [2,] -1319.508 -6579.733 -8846.323 4124.926 3217.822 $out [1] -22548 -32291 -12370 14042 $group [1] 5 9 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-23374, -13446.5, -8436, -3788, -1345, -26042, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ds5b1228935362.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,] 490042.0 545109.0 545720.0 505944 477930 466106.0 424476.0 384018.0 [2,] 498640.5 559464.5 566023.0 541740 497639 488715.5 446623.5 407829.5 [3,] 531778.0 592515.0 603538.5 553855 523796 504058.5 459139.0 422141.5 [4,] 572428.5 619624.0 620523.0 568904 532831 515301.0 476199.0 431698.5 [5,] 580177.0 631421.0 636067.0 596019 540380 530541.0 493756.0 445844.0 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 364635.0 389918.0 430457.0 505327 553997.0 555252.0 544897.0 536763.0 [2,] 379928.0 401772.5 440955.5 511620 567683.5 571356.0 576191.0 555928.5 [3,] 415459.5 433463.5 474713.0 545510 594556.0 604083.0 587154.0 568713.5 [4,] 421961.0 455855.0 501886.0 588925 603895.0 620132.5 605955.5 583820.5 [5,] 428212.0 460373.0 515476.0 593572 622255.0 630859.0 621170.0 607049.0 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 505282 475102.0 439148.0 431328.0 454724 500857.0 536662.0 567456.0 [2,] 523587 490744.5 458389.0 445427.5 465845 509530.0 549465.5 578764.0 [3,] 539708 507245.5 472578.5 469010.5 492644 533254.0 577927.5 593134.5 [4,] 557709 525100.5 490522.5 493423.0 517716 567393.5 603108.5 616426.0 [5,] 584225 551528.0 520564.0 512927.0 525192 579714.0 612613.0 628884.0 [,25] [,26] [,27] [1,] 557274.0 498662.0 469357.0 [2,] 572858.5 512648.0 490243.0 [3,] 583712.0 529371.5 506867.5 [4,] 606487.5 552108.0 517945.0 [5,] 625809.0 564591.0 533590.0 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 [26] 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 498122.8 565075.8 578680.7 541465.3 507744.7 491932.7 445649.4 411254.7 [2,] 565433.2 619954.2 628396.3 566244.7 539847.3 516184.3 472628.6 433028.3 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 396288 408796.1 446922.2 510250.7 578039.7 581835.7 573578.2 555991.8 [2,] 434631 458130.9 502503.8 580769.3 611072.3 626330.3 600729.8 581435.2 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 524144.7 491575.5 457922.2 447119.4 468985.3 506862.1 553460.6 575956.6 [2,] 555271.3 522915.5 487234.8 490901.6 516302.7 559645.9 602394.4 610312.4 [,25] [,26] [,27] [1,] 568373.6 511373.5 493026.5 [2,] 599050.4 547369.5 520708.5 $out [1] 619393 $group [1] 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" NA Warning message: In bxp(list(stats = c(490042, 498640.5, 531778, 572428.5, 580177, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7t6ih1228935362.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,] 489311.6 498662.0 494752.0 [2,] 506357.6 509270.0 511317.0 [3,] 532233.2 531792.8 539281.2 [4,] 542097.7 547018.5 547787.9 [5,] 558451.0 567831.0 567290.0 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 515931.9 514575.4 522646.6 [2,] 548534.5 549010.1 555915.9 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(489311.555555556, 506357.62962963, 532233.211538462, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/109p51228935362.ps tmp/109p51228935362.png") > system("convert tmp/2v0st1228935362.ps tmp/2v0st1228935362.png") > system("convert tmp/3grz91228935362.ps tmp/3grz91228935362.png") > system("convert tmp/44n6c1228935362.ps tmp/44n6c1228935362.png") > system("convert tmp/5nuie1228935362.ps tmp/5nuie1228935362.png") > system("convert tmp/6ds5b1228935362.ps tmp/6ds5b1228935362.png") > system("convert tmp/7t6ih1228935362.ps tmp/7t6ih1228935362.png") > > > proc.time() user system elapsed 1.257 0.990 1.592