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(665272,661735,621014,574889,677734,717075,653612,690697,665864,830701,789303,617808,805775,909449,599973,955874,799494,876097,823300,900079,860754,923882,1121084,741757,966066,901978,648659,852732,706036,835792,722489,714262,739459,816834,743082,683375,1006000,866000,644000,703000,699000,713000,688000,672000,600000,847000,697000,687000,973000,796000,658000,709000,798000,820000,776000,699000,828433,942131,792916,864942,982689,948143,874863,735794,854605,1284216,961585,818379,1079498,1095091,1008925,967118,1127715) > 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] 73 > (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] 7 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 665272 805775 966066 1006000 973000 982689 1127715 [2,] 661735 909449 901978 866000 796000 948143 NA [3,] 621014 599973 648659 644000 658000 874863 NA [4,] 574889 955874 852732 703000 709000 735794 NA [5,] 677734 799494 706036 699000 798000 854605 NA [6,] 717075 876097 835792 713000 820000 1284216 NA [7,] 653612 823300 722489 688000 776000 961585 NA [8,] 690697 900079 714262 672000 699000 818379 NA [9,] 665864 860754 739459 600000 828433 1079498 NA [10,] 830701 923882 816834 847000 942131 1095091 NA [11,] 789303 1121084 743082 697000 792916 1008925 NA [12,] 617808 741757 683375 687000 864942 967118 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3537 103674 -64088 -140000 -177000 -34546 NA [2,] -40721 -309476 -253319 -222000 -138000 -73280 NA [3,] -46125 355901 204073 59000 51000 -139069 NA [4,] 102845 -156380 -146696 -4000 89000 118811 NA [5,] 39341 76603 129756 14000 22000 429611 NA [6,] -63463 -52797 -113303 -25000 -44000 -322631 NA [7,] 37085 76779 -8227 -16000 -77000 -143206 NA [8,] -24833 -39325 25197 -72000 129433 261119 NA [9,] 164837 63128 77375 247000 113698 15593 NA [10,] -41398 197202 -73752 -150000 -149215 -86166 NA [11,] -171495 -379327 -59707 -10000 72026 -41807 NA [12,] 187967 224309 322625 286000 117747 160597 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/1k4w31244127511.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/2tbog1244127511.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/3yp2v1244127511.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/4l7uk1244127511.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,] 805775.0 661735 599973.0 574889 677734 713000 653612.0 672000 600000 [2,] 885920.5 796000 621014.0 703000 699000 717075 688000.0 690697 665864 [3,] 973000.0 883989 646329.5 722397 752018 827896 749244.5 706631 783946 [4,] 994344.5 909449 658000.0 852732 799494 876097 823300.0 818379 860754 [5,] 1127715.0 948143 658000.0 955874 854605 876097 961585.0 900079 1079498 [,10] [,11] [,12] [1,] 816834 697000.0 617808.0 [2,] 830701 743082.0 683375.0 [3,] 885441 791109.5 714378.5 [4,] 942131 1008925.0 864942.0 [5,] 1095091 1121084.0 967118.0 $n [1] 7 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 908251 810810.7 622472.3 625815 687196.1 725321.7 661971.6 624272 [2,] 1037749 957167.3 670186.7 818979 816839.9 930470.3 836517.4 788990 [,9] [,10] [,11] [,12] [1,] 658235.7 813565 619632.2 597261.9 [2,] 909656.3 957317 962586.8 831495.1 $out [1] 665272 874863 1284216 $group [1] 1 3 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(805775, 885920.5, 973000, 994344.5, 1127715, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5rrjg1244127511.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,] -177000 -309476 -139069 -156380 14000 -113303 -143206.0 -72000 15593.0 [2,] -140000 -253319 -46125 -146696 22000 -113303 -77000.0 -39325 63128.0 [3,] -49317 -180000 55000 42500 57972 -58130 -12113.5 182 95536.5 [4,] -3537 -73280 204073 102845 129756 -44000 37085.0 129433 164837.0 [5,] 103674 -40721 355901 118811 129756 -25000 76779.0 261119 247000.0 [,10] [,11] [,12] [1,] -150000 -379327 117747 [2,] -149215 -171495 160597 [3,] -79959 -50757 206138 [4,] -41398 -10000 286000 [5,] -41398 72026 322625 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -137340.04 -296130.97 -106385.8 -118462 -11534.1 -102832.67 -85702.01 [2,] 38706.04 -63869.03 216385.8 203462 127478.1 -13427.33 61475.01 [,8] [,9] [,10] [,11] [,12] [1,] -108672.4 29930.91 -149504.45 -154926.49 125249.0 [2,] 109036.4 161142.09 -10413.55 53412.49 287027.0 $out [1] 429611 -322631 197202 $group [1] 5 6 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-177000, -140000, -49317, -3537, 103674, -309476, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/62xhs1244127511.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] [1,] 574889 741757.0 648659.0 600000 658000.0 735794.0 1127715 [2,] 637313 802634.5 710149.0 679500 742500.0 864734.0 1127715 [3,] 665568 868425.5 741270.5 698000 797000.0 964351.5 1127715 [4,] 703886 916665.5 844262.0 780000 846687.5 1044211.5 1127715 [5,] 789303 955874.0 966066.0 866000 973000.0 1284216.0 1127715 $n [1] 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 635203.6 816415.2 680100.7 652161.3 749479.4 882490.6 1127715 [2,] 695932.4 920435.8 802440.3 743838.7 844520.6 1046212.4 1127715 $out [1] 830701 599973 1121084 1006000 $group [1] 1 2 2 4 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(574889, 637313, 665568, 703886, 789303, 741757, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7wym41244127511.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,] 674418.2 646329.5 641212.6 [2,] 755513.2 718387.8 751994.2 [3,] 783249.5 767982.0 766233.5 [4,] 866540.8 855942.5 857660.3 [5,] 932359.6 973000.0 940132.5 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 732609 705242.4 718038.5 [2,] 833890 830721.6 814428.5 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(674418.166666667, 755513.166666667, 783249.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1k4w31244127511.ps tmp/1k4w31244127511.png") > system("convert tmp/2tbog1244127511.ps tmp/2tbog1244127511.png") > system("convert tmp/3yp2v1244127511.ps tmp/3yp2v1244127511.png") > system("convert tmp/4l7uk1244127511.ps tmp/4l7uk1244127511.png") > system("convert tmp/5rrjg1244127511.ps tmp/5rrjg1244127511.png") > system("convert tmp/62xhs1244127511.ps tmp/62xhs1244127511.png") > system("convert tmp/7wym41244127511.ps tmp/7wym41244127511.png") > > > proc.time() user system elapsed 1.897 1.497 2.196