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(9487,8700,9627,8947,9283,8829,9947,9628,9318,9605,8640,9214,9567,8547,9185,9470,9123,9278,10170,9434,9655,9429,8739,9552,9784,9089,9763,9330,9144,9895,10404,10195,9987,9789,9437,10096,9776,9106,10258,9766,9826,9957,10036,10508,10146,10166,9365,9968,10123,9144,10447,9699,10451,10192,10404,10597,10633,10727,9784,9667,10297,9426,10274,9598,10400,9985,10761,11081,10297,10751,9760,10133,10806,9734,10083,10691,10446,10517,11353,10436,10721,10701,9793,10142) > 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] 84 > (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] 7 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 9487 9567 9784 9776 10123 10297 10806 NA [2,] 8700 8547 9089 9106 9144 9426 9734 NA [3,] 9627 9185 9763 10258 10447 10274 10083 NA [4,] 8947 9470 9330 9766 9699 9598 10691 NA [5,] 9283 9123 9144 9826 10451 10400 10446 NA [6,] 8829 9278 9895 9957 10192 9985 10517 NA [7,] 9947 10170 10404 10036 10404 10761 11353 NA [8,] 9628 9434 10195 10508 10597 11081 10436 NA [9,] 9318 9655 9987 10146 10633 10297 10721 NA [10,] 9605 9429 9789 10166 10727 10751 10701 NA [11,] 8640 8739 9437 9365 9784 9760 9793 NA [12,] 9214 9552 10096 9968 9667 10133 10142 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -787 -1020 -695 -670 -979 -871 -1072 NA [2,] 927 638 674 1152 1303 848 349 NA [3,] -680 285 -433 -492 -748 -676 608 NA [4,] 336 -347 -186 60 752 802 -245 NA [5,] -454 155 751 131 -259 -415 71 NA [6,] 1118 892 509 79 212 776 836 NA [7,] -319 -736 -209 472 193 320 -917 NA [8,] -310 221 -208 -362 36 -784 285 NA [9,] 287 -226 -198 20 94 454 -20 NA [10,] -965 -690 -352 -801 -943 -991 -908 NA [11,] 574 813 659 603 -117 373 349 NA [12,] 353 232 -320 155 630 673 NA 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/1bxaf1261223781.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/2pif61261223781.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/312hr1261223781.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/444c11261223781.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,] 9487.0 8547.0 9185 8947.0 9123.0 9278.0 9947.0 9434.0 9318 9429 [2,] 9671.5 8894.5 9695 9400.0 9213.5 9586.5 10103.0 9911.5 9821 9697 [3,] 9784.0 9106.0 10083 9598.0 9826.0 9957.0 10404.0 10436.0 10146 10166 [4,] 10210.0 9285.0 10266 9732.5 10423.0 10088.5 10582.5 10552.5 10465 10714 [5,] 10806.0 9734.0 10447 9766.0 10451.0 10517.0 10761.0 11081.0 10721 10751 [,11] [,12] [1,] 8640 9214.0 [2,] 9052 9609.5 [3,] 9437 9968.0 [4,] 9772 10114.5 [5,] 9793 10142.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 9462.416 8872.8 9742.008 9399.436 9103.706 9657.214 10117.65 10053.21 [2,] 10105.584 9339.2 10423.992 9796.564 10548.294 10256.786 10690.35 10818.79 [,9] [,10] [,11] [,12] [1,] 9761.414 9558.664 9007.028 9666.422 [2,] 10530.586 10773.336 9866.972 10269.578 $out [1] 10691 8829 11353 $group [1] 4 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(9487, 9671.5, 9784, 10210, 10806, 8547, 8894.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/56ar81261223781.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] [,11] [1,] -1072.0 349.0 -748 -347.0 -454 79.0 -917.0 -784.0 -226.0 -991.0 349 [2,] -999.5 656.0 -678 -215.5 -337 360.5 -527.5 -336.0 -109.0 -954.0 361 [3,] -871.0 848.0 -492 60.0 71 776.0 -209.0 -208.0 20.0 -908.0 574 [4,] -741.0 1039.5 -74 544.0 143 864.0 256.5 128.5 190.5 -745.5 631 [5,] -670.0 1303.0 608 802.0 751 1118.0 472.0 285.0 454.0 -690.0 813 [,12] [1,] -320.0 [2,] 155.0 [3,] 292.5 [4,] 630.0 [5,] 673.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1025.372 618.98 -852.6991 -393.5611 -215.6483 475.3179 -677.1922 [2,] -716.628 1077.02 -131.3009 513.5611 357.6483 1076.6821 259.1922 [,8] [,9] [,10] [,11] [,12] [1,] -485.3919 -158.8566 -1032.5128 412.7604 -13.89034 [2,] 69.3919 198.8566 -783.4872 735.2396 598.89034 $out [1] -352 -117 $group [1] 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1072, -999.5, -871, -741, -670, 349, 656, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6yclh1261223781.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,] 8640.0 8547.0 9089.0 9365.0 9144.0 9426.0 9734.0 NA [2,] 8888.0 9154.0 9383.5 9771.0 9741.5 9872.5 10112.5 NA [3,] 9300.5 9431.5 9786.5 9962.5 10298.0 10285.5 10481.5 NA [4,] 9616.0 9559.5 10041.5 10156.0 10524.0 10575.5 10711.0 NA [5,] 9947.0 9655.0 10404.0 10508.0 10727.0 11081.0 11353.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 8968.454 9246.549 9486.382 9786.899 9941.096 9964.857 10208.52 NA [2,] 9632.546 9616.451 10086.618 10138.101 10654.904 10606.143 10754.48 NA $out [1] 10170 9106 $group [1] 2 4 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(8640, 8888, 9300.5, 9616, 9947, 8547, 9154, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7i3wx1261223781.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,] 9359.714 9106.0 9089.750 [2,] 9725.286 9691.0 9692.250 [3,] 9886.357 9962.5 9901.375 [4,] 10137.500 10156.0 10174.250 [5,] 10439.286 10436.0 10342.750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 9698.343 9750.41 9681.532 [2,] 10074.371 10174.59 10121.218 $out [1] 9106.571 $group [1] 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(9359.7142857143, 9725.28571428571, 9886.35714285714, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1bxaf1261223781.ps tmp/1bxaf1261223781.png",intern=TRUE)) character(0) > try(system("convert tmp/2pif61261223781.ps tmp/2pif61261223781.png",intern=TRUE)) character(0) > try(system("convert tmp/312hr1261223781.ps tmp/312hr1261223781.png",intern=TRUE)) character(0) > try(system("convert tmp/444c11261223781.ps tmp/444c11261223781.png",intern=TRUE)) character(0) > try(system("convert tmp/56ar81261223781.ps tmp/56ar81261223781.png",intern=TRUE)) character(0) > try(system("convert tmp/6yclh1261223781.ps tmp/6yclh1261223781.png",intern=TRUE)) character(0) > try(system("convert tmp/7i3wx1261223781.ps tmp/7i3wx1261223781.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.217 0.950 2.329