R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-redhat-linux-gnu (64-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. Natural language support but running in an English locale 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(10407,10463,10556,10646,10702,11353,11346,11451,11964,12574,13031,13812,14544,14931,14886,16005,17064,15168,16050,15839,15137,14954,15648,15305,15579,16348,15928,16171,15937,15713,15594,15683,16438,17032,17696,17745,19394,20148,20108,18584,18441,18391,19178,18079,18483,19644,19195,19650,20830,23595,22937,21814,21928,21777,21383,21467,22052,22680,24320,24977,25204,25739,26434,27525,30695,32436,30160,30236,31293,31077,32226,33865,32810) > 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,] 10407 14544 15579 19394 20830 25204 32810 [2,] 10463 14931 16348 20148 23595 25739 NA [3,] 10556 14886 15928 20108 22937 26434 NA [4,] 10646 16005 16171 18584 21814 27525 NA [5,] 10702 17064 15937 18441 21928 30695 NA [6,] 11353 15168 15713 18391 21777 32436 NA [7,] 11346 16050 15594 19178 21383 30160 NA [8,] 11451 15839 15683 18079 21467 30236 NA [9,] 11964 15137 16438 18483 22052 31293 NA [10,] 12574 14954 17032 19644 22680 31077 NA [11,] 13031 15648 17696 19195 24320 32226 NA [12,] 13812 15305 17745 19650 24977 33865 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 56 387 769 754 2765 535 NA [2,] 93 -45 -420 -40 -658 695 NA [3,] 90 1119 243 -1524 -1123 1091 NA [4,] 56 1059 -234 -143 114 3170 NA [5,] 651 -1896 -224 -50 -151 1741 NA [6,] -7 882 -119 787 -394 -2276 NA [7,] 105 -211 89 -1099 84 76 NA [8,] 513 -702 755 404 585 1057 NA [9,] 610 -183 594 1161 628 -216 NA [10,] 457 694 664 -449 1640 1149 NA [11,] 781 -343 49 455 657 1639 NA [12,] 732 274 1649 1180 227 -1055 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/wessaorg/rcomp/tmp/13sj11305801133.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/www/wessaorg/rcomp/tmp/2w6s71305801133.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/www/wessaorg/rcomp/tmp/3r9jm1305801133.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/www/wessaorg/rcomp/tmp/41qb91305801133.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] [,9] [,10] [1,] 10407.0 10463 10556 10646.0 10702.0 11353 11346 11451 11964.0 12574 [2,] 15061.5 14931 14886 16005.0 15937.0 15168 15594 15683 15137.0 14954 [3,] 19394.0 18248 18018 17377.5 17752.5 17052 17614 16959 17460.5 18338 [4,] 23017.0 23595 22937 21814.0 21928.0 21777 21383 21467 22052.0 22680 [5,] 32810.0 25739 26434 27525.0 30695.0 21777 21383 21467 31293.0 31077 [,11] [,12] [1,] 13031.0 13812.0 [2,] 15648.0 15305.0 [3,] 18445.5 18697.5 [4,] 24320.0 24977.0 [5,] 32226.0 33865.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,] 14643.10 12659.44 12824.84 13630.51 13888.11 12788.98 13879.91 13228.13 [2,] 24144.90 23836.56 23211.16 21124.49 21616.89 21315.02 21348.09 20689.87 [,9] [,10] [,11] [,12] [1,] 13000.10 13354.48 12851.78 12458.75 [2,] 21920.90 23321.52 24039.22 24936.25 $out [1] 32436 30160 30236 $group [1] 6 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(10407, 15061.5, 19394, 23017, 32810, 10463, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/5sysm1305801133.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] [,11] [,12] [1,] 56.0 -658.0 -1524.0 -234 -224.0 -394 -211 404 -216 -449 -343 -1055 [2,] 387.0 -420.0 -1123.0 -143 -224.0 -394 -211 404 -183 457 49 227 [3,] 644.5 -42.5 166.5 85 -100.5 -63 80 549 602 679 556 503 [4,] 769.0 93.0 1091.0 1059 651.0 787 89 755 628 1149 781 1180 [5,] 769.0 695.0 1119.0 1059 1741.0 882 105 1057 1161 1640 1639 1649 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 398.0977 -373.4016 -1261.602 -690.3288 -664.9033 -824.7831 -113.5097 [2,] 890.9023 288.4016 1594.602 860.3288 463.9033 698.7831 273.5097 [,8] [,9] [,10] [,11] [,12] [1,] 322.5937 78.8788 232.6376 83.83636 -111.7158 [2,] 775.4063 1125.1212 1125.3624 1028.16364 1117.7158 $out [1] 2765 3170 -1896 -2276 -1099 -702 $group [1] 1 4 5 6 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(56, 387, 644.5, 769, 769, -658, -420, -42.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/6riq21305801133.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] [1,] 10407.0 14544.0 15579 18079.0 20830 25204.0 32810 [2,] 10601.0 14942.5 15698 18462.0 21622 26979.5 32810 [3,] 11349.5 15236.5 16054 19186.5 21990 30465.5 32810 [4,] 12269.0 15922.0 16735 19647.0 23266 31759.5 32810 [5,] 13812.0 17064.0 17745 20148.0 24977 33865.0 32810 $n [1] 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 10588.71 14789.74 15581.02 18646.01 21240.16 28285.31 32810 [2,] 12110.29 15683.26 16526.98 19726.99 22739.84 32645.69 32810 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(10407, 10601, 11349.5, 12269, 13812, 14544, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/786gc1305801134.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,] 18457.50 16959.00 18117.38 [2,] 18664.92 17419.00 18290.44 [3,] 19133.75 17885.25 18662.81 [4,] 19742.08 18391.75 19024.25 [5,] 20892.33 19394.00 19780.12 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 18642.45 17441.57 18328.12 [2,] 19625.05 18328.93 18997.51 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(18457.5, 18664.9166666667, 19133.75, 19742.0833333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/13sj11305801133.ps tmp/13sj11305801133.png",intern=TRUE)) character(0) > try(system("convert tmp/2w6s71305801133.ps tmp/2w6s71305801133.png",intern=TRUE)) character(0) > try(system("convert tmp/3r9jm1305801133.ps tmp/3r9jm1305801133.png",intern=TRUE)) character(0) > try(system("convert tmp/41qb91305801133.ps tmp/41qb91305801133.png",intern=TRUE)) character(0) > try(system("convert tmp/5sysm1305801133.ps tmp/5sysm1305801133.png",intern=TRUE)) character(0) > try(system("convert tmp/6riq21305801133.ps tmp/6riq21305801133.png",intern=TRUE)) character(0) > try(system("convert tmp/786gc1305801134.ps tmp/786gc1305801134.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.500 0.530 5.152