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(493395.00,487190.00,519493.00,519453.00,538588.00,438224.00,542034.00,512027.00,619880.00,533737.00,573789.00,589213.00,532168.00,551102.00,593789.00,527106.00,547327.00,601305.00,610872.00,601325.00,642143.00,614216.00,657979.00,673098.00,602297.00,615381.00,703671.00,733852.00,716596.00,745798.00,742027.10,679181.20,739022.70,645410.60,729382.10,671052.70,744954.80,677639.30,778207.20,763316.20,658531.60,831700.10,664156.30,621402.10,683588.70,600023.80,643273.80,653615.90,620177.50,574128.80,599828.00,599369.40,596617.70,616114.60,510226.90,493960.10,634503.30,588556.20,603239.00,617458.20,646543.50,680125.60,731595.80,759600.30,785031.70,849573.30,762342.00,815346.60,929603.20,784057.50,944667.70,1007258.30,664292.70,873207.40,1146510.00,1417266.80,1089387.90,1373379.70,1009397.60,818175.10,1003458.10,961142.70,1121906.60,1141713.30,1042352.60,992223.60,920525.30,1076093.40,967880.40,1236416.10) > 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] 90 > (np <- floor(n / par1)) [1] 7 > 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] 8 8 8 8 8 8 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 493395 532168 602297.0 744954.8 620177.5 646543.5 664292.7 1042352.6 [2,] 487190 551102 615381.0 677639.3 574128.8 680125.6 873207.4 992223.6 [3,] 519493 593789 703671.0 778207.2 599828.0 731595.8 1146510.0 920525.3 [4,] 519453 527106 733852.0 763316.2 599369.4 759600.3 1417266.8 1076093.4 [5,] 538588 547327 716596.0 658531.6 596617.7 785031.7 1089387.9 967880.4 [6,] 438224 601305 745798.0 831700.1 616114.6 849573.3 1373379.7 1236416.1 [7,] 542034 610872 742027.1 664156.3 510226.9 762342.0 1009397.6 NA [8,] 512027 601325 679181.2 621402.1 493960.1 815346.6 818175.1 NA [9,] 619880 642143 739022.7 683588.7 634503.3 929603.2 1003458.1 NA [10,] 533737 614216 645410.6 600023.8 588556.2 784057.5 961142.7 NA [11,] 573789 657979 729382.1 643273.8 603239.0 944667.7 1121906.6 NA [12,] 589213 673098 671052.7 653615.9 617458.2 1007258.3 1141713.3 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -6205 18934 13084.0 -67315.5 -46048.7 33582.1 208914.7 -50129.0 [2,] 32303 42687 88290.0 100567.9 25699.2 51470.2 273302.6 -71698.3 [3,] -40 -66683 30181.0 -14891.0 -458.6 28004.5 270756.8 155568.1 [4,] 19135 20221 -17256.0 -104784.6 -2751.7 25431.4 -327878.9 -108213.0 [5,] -100364 53978 29202.0 173168.5 19496.9 64541.6 283991.8 268535.7 [6,] 103810 9567 -3770.9 -167543.8 -105887.7 -87231.3 -363982.1 NA [7,] -30007 -9547 -62845.9 -42754.2 -16266.8 53004.6 -191222.5 NA [8,] 107853 40818 59841.5 62186.6 140543.2 114256.6 185283.0 NA [9,] -86143 -27927 -93612.1 -83564.9 -45947.1 -145545.7 -42315.4 NA [10,] 40052 43763 83971.5 43250.0 14682.8 160610.2 160763.9 NA [11,] 15424 15119 -58329.4 10342.1 14219.2 62590.6 19806.7 NA [12,] -57045 -70801 73902.1 -33438.4 29085.3 -342965.6 -99360.7 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/132hm1244392911.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/2qrhh1244392911.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/3gepl1244392911.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/4jwry1244392911.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] [1,] 493395.0 487190.0 519493.0 519453.0 538588.0 438224.0 510226.9 [2,] 567232.5 562615.4 596808.5 563237.7 571972.3 608709.8 576453.0 [3,] 633360.5 646510.2 717633.4 746726.2 687563.8 788749.1 664156.3 [4,] 704623.8 776666.5 849366.2 919704.8 876456.1 1042994.7 752184.6 [5,] 744954.8 992223.6 1146510.0 1417266.8 1089387.9 1373379.7 1009397.6 [,8] [,9] [,10] [,11] [,12] [1,] 493960.1 619880.0 533737.0 573789.0 589213.0 [2,] 556676.0 638323.2 594290.0 623256.4 635537.1 [3,] 621402.1 683588.7 614216.0 657979.0 671052.7 [4,] 747263.9 834312.9 714734.1 837024.9 840178.2 [5,] 818175.1 1003458.1 784057.5 1121906.6 1141713.3 $n [1] 8 8 8 8 8 8 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 556611.8 526938.1 576551.0 547598.5 517474.8 546151.3 559212.3 507586.1 [2,] 710109.2 766082.2 858715.8 945853.8 857652.8 1031346.8 769100.3 735218.1 [,9] [,10] [,11] [,12] [1,] 566546.8 542288.8 530319.9 548844.3 [2,] 800630.6 686143.2 785638.1 793261.1 $out [1] 1042352.6 961142.7 $group [1] 1 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(493395, 567232.5, 633360.5, 704623.75, 744954.8, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/54fdo1244392911.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] [1,] -67315.50 25699.20 -66683.00 -108213.00 -100364.00 -167543.80 -62845.90 [2,] -48088.85 29001.10 -7674.80 -106498.80 24349.45 -136715.75 -52800.05 [3,] 3439.50 47078.60 13982.25 -10003.85 59259.80 -87231.30 -30007.00 [4,] 26258.05 94428.95 92874.55 19678.00 220852.10 2898.05 -12906.90 [5,] 33582.10 100567.90 155568.10 25431.40 283991.80 103810.00 -9547.00 [,8] [,9] [,10] [,11] [,12] [1,] 40818.00 -145545.70 14682.8 10342.10 -99360.70 [2,] 61014.05 -89877.55 41651.0 12280.65 -85080.85 [3,] 107853.00 -83564.90 43763.0 15119.00 -57045.00 [4,] 127399.90 -44131.25 122290.8 17615.35 -2176.55 [5,] 185283.00 -27927.00 160763.9 19806.70 73902.10 $n [1] 8 8 8 8 8 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -38091.75 10529.66 -42186.06 -80488.02 -50509.4 -170606.409 -53830.546 [2,] 44970.75 83627.54 70150.56 60480.32 169029.0 -3856.191 -6183.454 [,8] [,9] [,10] [,11] [,12] [1,] 68208.44 -110883.85 -4393.817 11933.20 -106554.11 [2,] 147497.56 -56245.95 91919.817 18304.80 -7535.89 $out [1] 208914.7 273302.6 -71698.3 270756.8 -327878.9 -363982.1 53004.6 [8] -191222.5 -58329.4 62590.6 -342965.6 $group [1] 1 2 2 3 4 6 7 7 11 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-67315.5, -48088.85, 3439.5, 26258.05, 33582.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6f2341244392911.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,] 438224.0 527106.0 602297.0 600023.8 574128.8 646543.5 664292.7 920525.3 [2,] 502711.0 549214.5 658231.6 648444.9 581342.5 745598.1 917175.1 967880.4 [3,] 526615.0 601315.0 710133.5 670897.8 599598.7 784544.6 1049392.8 1017288.1 [4,] 557911.5 628179.5 736437.3 754135.5 616786.4 889588.2 1144111.6 1076093.4 [5,] 619880.0 673098.0 745798.0 831700.1 634503.3 1007258.3 1417266.8 1236416.1 $n [1] 12 12 12 12 12 12 12 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 501437.7 565298.5 674463.4 622691.6 583432.5 718869.7 945885.4 947487.2 [2,] 551792.3 637331.5 745803.6 719104.0 615764.9 850219.5 1152900.1 1087089.0 $out [1] 510226.9 493960.1 $group [1] 5 5 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(438224, 502711, 526615, 557911.5, 619880, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/79mfr1244392911.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,] 648773.9 614216.0 634611.5 [2,] 678340.5 639935.3 653241.0 [3,] 743348.7 667604.5 706786.0 [4,] 759117.6 702598.6 733229.3 [5,] 836563.8 788749.1 779348.1 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 706505.7 639023.4 670302.7 [2,] 780191.7 696185.6 743269.2 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(648773.871428571, 678340.484821429, 743348.725, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/132hm1244392911.ps tmp/132hm1244392911.png") > system("convert tmp/2qrhh1244392911.ps tmp/2qrhh1244392911.png") > system("convert tmp/3gepl1244392911.ps tmp/3gepl1244392911.png") > system("convert tmp/4jwry1244392911.ps tmp/4jwry1244392911.png") > system("convert tmp/54fdo1244392911.ps tmp/54fdo1244392911.png") > system("convert tmp/6f2341244392911.ps tmp/6f2341244392911.png") > system("convert tmp/79mfr1244392911.ps tmp/79mfr1244392911.png") > > > proc.time() user system elapsed 1.213 0.944 2.068