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(428800 + ,424800 + ,403400 + ,398400 + ,393500 + ,380500 + ,398300 + ,387300 + ,370400 + ,372800 + ,444600 + ,449900 + ,458100 + ,424800 + ,420600 + ,400100 + ,393000 + ,387100 + ,377500 + ,400400 + ,391400 + ,363600 + ,431000 + ,441700 + ,448500 + ,415600 + ,408000 + ,416600 + ,409300 + ,387600 + ,394500 + ,407600 + ,378500 + ,359600 + ,435700 + ,433800 + ,427700 + ,413300 + ,379500 + ,379300 + ,353700 + ,378200 + ,380600 + ,394000 + ,374000 + ,375000 + ,437600 + ,443900 + ,488800 + ,463900 + ,440000 + ,453800 + ,451600 + ,453400 + ,461400 + ,509100 + ,540600 + ,555100 + ,677400 + ,694600 + ,750100 + ,733900 + ,709300 + ,720500 + ,693200 + ,687200 + ,686800 + ,720900 + ,653100 + ,624700 + ,690000 + ,717800 + ,736500 + ,699900 + ,675600 + ,635600 + ,632500 + ,594900 + ,604000 + ,620800 + ,578400 + ,571200 + ,627400 + ,657700 + ,674100 + ,672800 + ,615300 + ,609100 + ,607600 + ,566900 + ,572700 + ,589200 + ,534800 + ,543100 + ,591100 + ,624800 + ,665300 + ,642600 + ,608700 + ,594500 + ,563800 + ,596100 + ,597600 + ,633100 + ,591000 + ,584200 + ,655800 + ,670700 + ,699700 + ,712900 + ,652000 + ,635100 + ,603100 + ,610100 + ,602000 + ,597600 + ,585400 + ,567100 + ,620600 + ,646200 + ,644800 + ,645200 + ,644800 + ,593000 + ,569100 + ,518800 + ,538700 + ,554600 + ,507900 + ,488400 + ,563300 + ,592400 + ,598100 + ,546300 + ,516100 + ,518500 + ,477400 + ,483400 + ,469400 + ,501300 + ,457400 + ,446700 + ,501900 + ,550400 + ,593700 + ,548900 + ,534200 + ,550500 + ,541800 + ,569300 + ,587400 + ,627700 + ,607000 + ,629500 + ,704600 + ,767700 + ,812200 + ,824600 + ,856300 + ,812200 + ,764100 + ,801700 + ,806000 + ,867200 + ,801600 + ,817500 + ,920900 + ,959700 + ,997700 + ,949100 + ,910900 + ,920400 + ,914200 + ,926300 + ,906400 + ,926100 + ,902500 + ,895300 + ,979900 + ,1009700 + ,1043800 + ,979800 + ,921600 + ,923500 + ,914500 + ,891700 + ,916000 + ,931700 + ,902400 + ,893700 + ,941500 + ,980100 + ,1006900 + ,949200 + ,883200 + ,849900 + ,839200 + ,803900 + ,797900 + ,830800 + ,753300 + ,764100 + ,807600 + ,853700 + ,886200 + ,815700 + ,743000 + ,753600 + ,724800 + ,709600 + ,721900) > 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] 211 > (np <- floor(n / par1)) [1] 17 > 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] 18 18 18 18 18 18 18 17 17 17 17 17 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 428800 458100 448500 427700 488800 750100 736500 674100 665300 699700 [2,] 424800 424800 415600 413300 463900 733900 699900 672800 642600 712900 [3,] 403400 420600 408000 379500 440000 709300 675600 615300 608700 652000 [4,] 398400 400100 416600 379300 453800 720500 635600 609100 594500 635100 [5,] 393500 393000 409300 353700 451600 693200 632500 607600 563800 603100 [6,] 380500 387100 387600 378200 453400 687200 594900 566900 596100 610100 [7,] 398300 377500 394500 380600 461400 686800 604000 572700 597600 602000 [8,] 387300 400400 407600 394000 509100 720900 620800 589200 633100 597600 [9,] 370400 391400 378500 374000 540600 653100 578400 534800 591000 585400 [10,] 372800 363600 359600 375000 555100 624700 571200 543100 584200 567100 [11,] 444600 431000 435700 437600 677400 690000 627400 591100 655800 620600 [12,] 449900 441700 433800 443900 694600 717800 657700 624800 670700 646200 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] 644800 598100 593700 812200 997700 1043800 1006900 886200 [2,] 645200 546300 548900 824600 949100 979800 949200 815700 [3,] 644800 516100 534200 856300 910900 921600 883200 743000 [4,] 593000 518500 550500 812200 920400 923500 849900 753600 [5,] 569100 477400 541800 764100 914200 914500 839200 724800 [6,] 518800 483400 569300 801700 926300 891700 803900 709600 [7,] 538700 469400 587400 806000 906400 916000 797900 721900 [8,] 554600 501300 627700 867200 926100 931700 830800 NA [9,] 507900 457400 607000 801600 902500 902400 753300 NA [10,] 488400 446700 629500 817500 895300 893700 764100 NA [11,] 563300 501900 704600 920900 979900 941500 807600 NA [12,] 592400 550400 767700 959700 1009700 980100 853700 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] -4000 -33300 -32900 -14400 -24900 -16200 -36600 -1300 -22700 13200 [2,] -21400 -4200 -7600 -33800 -23900 -24600 -24300 -57500 -33900 -60900 [3,] -5000 -20500 8600 -200 13800 11200 -40000 -6200 -14200 -16900 [4,] -4900 -7100 -7300 -25600 -2200 -27300 -3100 -1500 -30700 -32000 [5,] -13000 -5900 -21700 24500 1800 -6000 -37600 -40700 32300 7000 [6,] 17800 -9600 6900 2400 8000 -400 9100 5800 1500 -8100 [7,] -11000 22900 13100 13400 47700 34100 16800 16500 35500 -4400 [8,] -16900 -9000 -29100 -20000 31500 -67800 -42400 -54400 -42100 -12200 [9,] 2400 -27800 -18900 1000 14500 -28400 -7200 8300 -6800 -18300 [10,] 71800 67400 76100 62600 122300 65300 56200 48000 71600 53500 [11,] 5300 10700 -1900 6300 17200 27800 30300 33700 14900 25600 [12,] 8200 6800 -6100 44900 55500 18700 16400 40500 29000 -1400 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] 400 -51800 -44800 12400 -48600 -64000 -57700 -70500 [2,] -400 -30200 -14700 31700 -38200 -58200 -66000 -72700 [3,] -51800 2400 16300 -44100 9500 1900 -33300 10600 [4,] -23900 -41100 -8700 -48100 -6200 -9000 -10700 -28800 [5,] -50300 6000 27500 37600 12100 -22800 -35300 -15200 [6,] 19900 -14000 18100 4300 -19900 24300 -6000 12300 [7,] 15900 31900 40300 61200 19700 15700 32900 NA [8,] -46700 -43900 -20700 -65600 -23600 -29300 -77500 NA [9,] -19500 -10700 22500 15900 -7200 -8700 10800 NA [10,] 74900 55200 75100 103400 84600 47800 43500 NA [11,] 29100 48500 63100 38800 29800 38600 46100 NA [12,] 5700 43300 44500 38000 34100 26800 32500 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/1i3k71241354940.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/2peea1241354940.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/3s8ke1241354940.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/4zf3f1241354940.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] [,10] [1,] 427700 413300 379500 379300 353700 378200 377500 387300 370400 359600 [2,] 488800 463900 440000 453800 451600 453400 461400 501300 457400 446700 [3,] 669700 659000 630050 601800 586100 582100 592500 597600 578400 567100 [4,] 812200 815700 743000 753600 724800 709600 721900 720900 653100 629500 [5,] 1043800 979800 921600 923500 914500 926300 916000 931700 902500 895300 [,11] [,12] [1,] 431000 433800 [2,] 501900 550400 [3,] 627400 657700 [4,] 704600 767700 [5,] 979900 1009700 $n [1] 18 18 18 18 18 18 18 17 17 17 17 17 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 549262.7 527986.3 517209.9 490151.6 484357.7 486688.7 495487.3 513447.9 [2,] 790137.3 790013.7 742890.1 713448.4 687842.3 677511.3 689512.7 681752.1 [,9] [,10] [,11] [,12] [1,] 503406.5 497049.9 549724.1 574429.3 [2,] 653393.5 637150.1 705075.9 740970.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(427700, 488800, 669700, 812200, 1043800, 413300, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/50zzb1241354940.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,] -70500 -72700 -51800 -48100 -50300 -19900 -11000 -77500 -28400 43500 [2,] -48600 -57500 -20500 -28800 -22800 -6000 15700 -46700 -18300 55200 [3,] -28900 -27400 -2600 -9850 -5950 5050 19700 -29300 -7200 67400 [4,] -4000 -14700 9500 -6200 12100 12300 34100 -20000 8300 75100 [5,] 13200 31700 16300 -1500 37600 24300 61200 -9000 22500 103400 [,11] [,12] [1,] -1900 -6100 [2,] 14900 8200 [3,] 29100 29000 [4,] 38600 40500 [5,] 63100 55500 $n [1] 18 18 18 18 18 18 17 17 17 17 17 17 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -45509.47 -43339.13 -13772.287 -18266.456 -18947.094 -1765.095 12649.00 [2,] -12290.53 -11460.87 8572.287 -1433.544 7047.094 11865.095 26751.00 [,8] [,9] [,10] [,11] [,12] [1,] -39531.61 -17393.287 59774.2 20018.01 16622.44 [2,] -19068.39 2993.287 75025.8 38181.99 41377.56 $out [1] 31500 122300 $group [1] 8 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-70500, -48600, -28900, -4000, 13200, -72700, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6a7081241354940.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] [,10] [1,] 370400 363600 359600 353700 440000 653100 571200 534800 563800 567100 [2,] 383900 389250 391050 376600 453600 687000 599450 569800 592750 599800 [3,] 398350 400250 408650 380050 476350 701250 629950 599350 603150 615350 [4,] 426800 427900 425200 420500 547850 720700 666650 620050 649200 649100 [5,] 449900 458100 448500 443900 677400 750100 736500 674100 670700 712900 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] 488400 446700 534200 764100 895300 891700 753300 709600 [2,] 528750 473400 549700 803850 908650 908450 800900 723350 [3,] 566200 501600 590550 814850 923250 922550 835000 743000 [4,] 618900 532400 628600 861750 964500 960650 868450 784650 [5,] 645200 598100 704600 920900 1009700 980100 949200 815700 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 378783 382621.5 393074.0 360026.9 433361.9 685879.2 599299.6 576430.6 [2,] 417917 417878.5 424226.0 400073.1 519338.1 716620.8 660600.4 622269.4 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 577402.8 592863.9 525082 474689.7 554563.2 788441.4 897776.4 898741.2 [2,] 628897.2 637836.1 607318 528510.3 626536.8 841258.6 948723.6 946358.8 [,17] [,18] [1,] 804190 706392.6 [2,] 865810 779607.4 $out [1] 694600 624700 767700 959700 1043800 1006900 886200 $group [1] 5 6 13 14 16 17 18 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" NA Warning message: In bxp(list(stats = c(370400, 383900, 398350, 426800, 449900, 363600, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/79e3w1241354940.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,] 579505.9 567100 538100.0 [2,] 599050.0 584100 584962.5 [3,] 618933.7 599700 600025.0 [4,] 653974.3 643875 625487.5 [5,] 686722.2 669700 659050.0 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 593882.3 572436.2 581541.3 [2,] 643985.0 626963.8 618508.7 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(579505.882352941, 599050, 618933.660130719, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1i3k71241354940.ps tmp/1i3k71241354940.png") > system("convert tmp/2peea1241354940.ps tmp/2peea1241354940.png") > system("convert tmp/3s8ke1241354940.ps tmp/3s8ke1241354940.png") > system("convert tmp/4zf3f1241354940.ps tmp/4zf3f1241354940.png") > system("convert tmp/50zzb1241354940.ps tmp/50zzb1241354940.png") > system("convert tmp/6a7081241354940.ps tmp/6a7081241354940.png") > system("convert tmp/79e3w1241354940.ps tmp/79e3w1241354940.png") > > > proc.time() user system elapsed 1.937 1.493 2.470