R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(192528,190773,188996,185320,221698,219771,192528,174414,176163,176163,178112,181616,192528,188996,194449,203412,254400,254400,243516,232604,241567,252473,254400,259853,276217,265306,265306,281670,327033,330709,321580,299762,316121,316121,317876,327033,334241,337917,337917,348823,390681,401565,403314,376071,390681,385228,374317,397889,403314,394185,396112,408773,456085,479624,479624,468740,485082,468740,459588,494239,499664,486832,519533,532366,570521,595842,592339,590389,604999,603222,581432,614128,625040,614128,659491,681309,732102,752143,746712,735800,744935,755841,719441,748461,766753,759373,806657,822993,892101,904761,888424,897554,903007,908460,873809,906511,924624,906511,959275,975618,1046468,1057380,1060884,1079170,1079170,1086378,1053676,1070041,1080925,1060884,1119079,1129986,1202597,1215430,1233543,1249908,1251657,1253584,1220883,1253584) > 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] 120 > (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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 192528 192528 276217 334241 403314 499664 625040 766753 924624 1080925 [2,] 190773 188996 265306 337917 394185 486832 614128 759373 906511 1060884 [3,] 188996 194449 265306 337917 396112 519533 659491 806657 959275 1119079 [4,] 185320 203412 281670 348823 408773 532366 681309 822993 975618 1129986 [5,] 221698 254400 327033 390681 456085 570521 732102 892101 1046468 1202597 [6,] 219771 254400 330709 401565 479624 595842 752143 904761 1057380 1215430 [7,] 192528 243516 321580 403314 479624 592339 746712 888424 1060884 1233543 [8,] 174414 232604 299762 376071 468740 590389 735800 897554 1079170 1249908 [9,] 176163 241567 316121 390681 485082 604999 744935 903007 1079170 1251657 [10,] 176163 252473 316121 385228 468740 603222 755841 908460 1086378 1253584 [11,] 178112 254400 317876 374317 459588 581432 719441 873809 1053676 1220883 [12,] 181616 259853 327033 397889 494239 614128 748461 906511 1070041 1253584 [,11] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] -1755 -3532 -10911 3676 -9129 -12832 -10912 -7380 -18113 -20041 [2,] -1777 5453 0 0 1927 32701 45363 47284 52764 58195 [3,] -3676 8963 16364 10906 12661 12833 21818 16336 16343 10907 [4,] 36378 50988 45363 41858 47312 38155 50793 69108 70850 72611 [5,] -1927 0 3676 10884 23539 25321 20041 12660 10912 12833 [6,] -27243 -10884 -9129 1749 0 -3503 -5431 -16337 3504 18113 [7,] -18114 -10912 -21818 -27243 -10884 -1950 -10912 9130 18286 16365 [8,] 1749 8963 16359 14610 16342 14610 9135 5453 0 1749 [9,] 0 10906 0 -5453 -16342 -1777 10906 5453 7208 1927 [10,] 1949 1927 1755 -10911 -9152 -21790 -36400 -34651 -32702 -32701 [11,] 3504 5453 9157 23572 34651 32696 29020 32702 16365 32701 [12,] 10912 16364 7208 5425 5425 10912 18292 18113 10884 NA [,11] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [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/wessaorg/rcomp/tmp/13yen1312380831.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2qzic1312380831.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3psk71312380831.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4ku1k1312380831.ps",horizontal=F,onefile=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] [1,] 192528 188996.0 188996.0 185320.0 221698 219771 192528.0 174414.0 [2,] 276217 265306.0 265306.0 281670.0 327033 330709 321580.0 299762.0 [3,] 451489 440508.5 457822.5 470569.5 513303 537733 535981.5 529564.5 [4,] 766753 759373.0 806657.0 822993.0 892101 904761 888424.0 897554.0 [5,] 1080925 1060884.0 1119079.0 1129986.0 1202597 1215430 1233543.0 1249908.0 [,9] [,10] [,11] [,12] [1,] 176163.0 176163 178112 181616.0 [2,] 316121.0 316121 317876 327033.0 [3,] 545040.5 535981 520510 554183.5 [4,] 903007.0 908460 873809 906511.0 [5,] 1251657.0 1253584 1220883 1253584.0 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 206397.7 193652.9 187342 200102.9 230972.5 250913.7 252763.6 230883.8 [2,] 696580.3 687364.1 728303 741036.1 795633.5 824552.3 819199.4 828245.2 [,9] [,10] [,11] [,12] [1,] 251808.9 240024.8 242743.7 264653.2 [2,] 838272.1 831937.2 798276.3 843713.8 $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(192528, 276217, 451489, 766753, 1080925, 188996, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/580ch1312380831.ps",horizontal=F,onefile=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,] -20041 -1777 8963 36378.0 -1927 -27243 -27243 0 -5453.0 -36400.0 [2,] -12832 0 10906 41858.0 3676 -10884 -18114 1749 -1777.0 -32702.0 [3,] -10020 19077 12747 49052.5 11786 -4467 -10898 9049 963.5 -16350.5 [4,] -3532 47284 16343 69108.0 20041 1749 9130 14610 7208.0 1755.0 [5,] 3676 58195 21818 72611.0 25321 18113 18286 16359 10906.0 1949.0 [,11] [,12] [1,] 3504 5425 [2,] 9157 7208 [3,] 26296 10912 [4,] 32701 16364 [5,] 34651 18292 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -14666.651 -4547.972 10030.46 35437.31 3609.394 -10778.95 -24510.189 [2,] -5373.349 42701.972 15463.54 62667.69 19962.606 1844.95 2714.189 [,8] [,9] [,10] [,11] [,12] [1,] 2623.132 -3525.764 -33566.591 14532.48 6089.84 [2,] 15474.868 5452.764 865.591 38059.52 15734.16 $out [1] -3676 -16342 $group [1] 3 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-20041, -12832, -10020, -3532, 3676, -1777, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6lzpn1312380831.ps",horizontal=F,onefile=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,] 174414.0 188996.0 265306.0 334241.0 394185.0 486832.0 614128.0 759373.0 [2,] 177137.5 198930.5 278943.5 343370.0 406043.5 525949.5 670400.0 814825.0 [3,] 187158.0 242541.5 316121.0 380649.5 464164.0 585910.5 733951.0 890262.5 [4,] 192528.0 254400.0 324306.5 394285.0 479624.0 599532.0 747586.5 903884.0 [5,] 192528.0 259853.0 330709.0 403314.0 494239.0 614128.0 755841.0 908460.0 [,9] [,10] [,11] [1,] 906511.0 1060884 NA [2,] 967446.5 1124532 NA [3,] 1055528.0 1218156 NA [4,] 1074605.5 1250782 NA [5,] 1086378.0 1253584 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 180138.3 217241.5 295430.6 357426.8 430603.4 552349 698745.7 849642.1 [2,] 194177.7 267841.5 336811.4 403872.2 497724.6 619472 769156.3 930882.9 [,9] [,10] [,11] [1,] 1006652 1160573 NA [2,] 1104404 1275740 NA $out [1] 221698 219771 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(174414, 177137.5, 187158, 192528, 192528, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7ursy1312380831.ps",horizontal=F,onefile=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,] 520490.5 440508.5 503260.2 [2,] 550854.2 464196.0 534838.6 [3,] 609904.9 525037.2 592741.1 [4,] 619979.6 536857.2 600488.2 [5,] 625335.5 554183.5 607514.8 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 578376.4 491896.0 562797.8 [2,] 641433.4 558178.5 622684.3 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(520490.5, 550854.25, 609904.9, 619979.6, 625335.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/13yen1312380831.ps tmp/13yen1312380831.png",intern=TRUE)) character(0) > try(system("convert tmp/2qzic1312380831.ps tmp/2qzic1312380831.png",intern=TRUE)) character(0) > try(system("convert tmp/3psk71312380831.ps tmp/3psk71312380831.png",intern=TRUE)) character(0) > try(system("convert tmp/4ku1k1312380831.ps tmp/4ku1k1312380831.png",intern=TRUE)) character(0) > try(system("convert tmp/580ch1312380831.ps tmp/580ch1312380831.png",intern=TRUE)) character(0) > try(system("convert tmp/6lzpn1312380831.ps tmp/6lzpn1312380831.png",intern=TRUE)) character(0) > try(system("convert tmp/7ursy1312380831.ps tmp/7ursy1312380831.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.552 0.343 1.897