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(34671,34355,34035,33373,39924,39577,34671,31409,31724,31724,32075,32706,34671,34035,35017,36631,45813,45813,43853,41888,43502,45466,45813,46795,49742,47777,47777,50724,58893,59555,57911,53982,56928,56928,57244,58893,60191,60853,60853,62817,70355,72315,72630,67724,70355,69373,67408,71653,72630,70986,71333,73613,82133,86372,86372,84412,87355,84412,82764,89004,89981,87670,93559,95870,102741,107301,106670,106319,108950,108630,104706,110594,112559,110594,118763,122692,131839,135448,134470,132505,134150,136114,129559,134785,138079,136750,145265,148207,160652,162932,159990,161634,162616,163598,157358,163247,166509,163247,172749,175692,188451,190416,191047,194340,194340,195638,189749,192696,194656,191047,201527,203491,216567,218878,222140,225087,225402,225749,219860,225749) > 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] [,11] [1,] 34671 34671 49742 60191 72630 89981 112559 138079 166509 194656 NA [2,] 34355 34035 47777 60853 70986 87670 110594 136750 163247 191047 NA [3,] 34035 35017 47777 60853 71333 93559 118763 145265 172749 201527 NA [4,] 33373 36631 50724 62817 73613 95870 122692 148207 175692 203491 NA [5,] 39924 45813 58893 70355 82133 102741 131839 160652 188451 216567 NA [6,] 39577 45813 59555 72315 86372 107301 135448 162932 190416 218878 NA [7,] 34671 43853 57911 72630 86372 106670 134470 159990 191047 222140 NA [8,] 31409 41888 53982 67724 84412 106319 132505 161634 194340 225087 NA [9,] 31724 43502 56928 70355 87355 108950 134150 162616 194340 225402 NA [10,] 31724 45466 56928 69373 84412 108630 136114 163598 195638 225749 NA [11,] 32075 45813 57244 67408 82764 104706 129559 157358 189749 219860 NA [12,] 32706 46795 58893 71653 89004 110594 134785 163247 192696 225749 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -316 -636 -1965 662 -1644 -2311 -1965 -1329 -3262 -3609 NA [2,] -320 982 0 0 347 5889 8169 8515 9502 10480 NA [3,] -662 1614 2947 1964 2280 2311 3929 2942 2943 1964 NA [4,] 6551 9182 8169 7538 8520 6871 9147 12445 12759 13076 NA [5,] -347 0 662 1960 4239 4560 3609 2280 1965 2311 NA [6,] -4906 -1960 -1644 315 0 -631 -978 -2942 631 3262 NA [7,] -3262 -1965 -3929 -4906 -1960 -351 -1965 1644 3293 2947 NA [8,] 315 1614 2946 2631 2943 2631 1645 982 0 315 NA [9,] 0 1964 0 -982 -2943 -320 1964 982 1298 347 NA [10,] 351 347 316 -1965 -1648 -3924 -6555 -6240 -5889 -5889 NA [11,] 631 982 1649 4245 6240 5888 5226 5889 2947 5889 NA [12,] 1965 2947 1298 977 977 1965 3294 3262 1960 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/wessaorg/rcomp/tmp/18tul1313684018.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/2iqat1313684018.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/38vmk1313684018.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/4nozn1313684018.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] [1,] 34671.0 34035 34035 33373.0 39924 39577.0 34671 31409.0 31724.0 [2,] 49742.0 47777 47777 50724.0 58893 59555.0 57911 53982.0 56928.0 [3,] 81305.5 79328 82446 84741.5 92437 96836.5 96521 95365.5 98152.5 [4,] 138079.0 136750 145265 148207.0 160652 162932.0 159990 161634.0 162616.0 [5,] 194656.0 191047 201527 203491.0 216567 218878.0 222140 225087.0 225402.0 [,10] [,11] [,12] [1,] 31724 32075 32706 [2,] 56928 57244 58893 [3,] 96521 93735 99799 [4,] 163598 157358 163247 [5,] 225749 219860 225749 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 37168.81 34873.54 33737.11 36035.11 41594.15 45185.23 45518.26 [2,] 125442.19 123782.46 131154.89 133447.89 143279.85 148487.77 147523.74 [,8] [,9] [,10] [,11] [,12] [1,] 41578.27 45346.56 43224.42 43714.05 47659.58 [2,] 149152.73 150958.44 149817.58 143755.95 151938.42 $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(34671, 49742, 81305.5, 138079, 194656, 34035, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5u82h1313684018.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] [1,] -3609.0 -320.0 1614.0 6551.0 -347.0 -4906.0 -4906.0 0.0 -982.0 [2,] -2311.0 0.0 1964.0 7538.0 662.0 -1960.0 -3262.0 315.0 -320.0 [3,] -1804.5 3435.5 2295.5 8833.5 2122.5 -804.5 -1962.5 1629.5 173.5 [4,] -636.0 8515.0 2943.0 12445.0 3609.0 315.0 1644.0 2631.0 1298.0 [5,] 662.0 10480.0 3929.0 13076.0 4560.0 3262.0 3293.0 2946.0 1964.0 [,10] [,11] [,12] [1,] -6555.0 631.0 977 [2,] -5889.0 1649.0 1298 [3,] -2944.5 4735.5 1965 [4,] 316.0 5889.0 2947 [5,] 351.0 6240.0 3294 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2641.3968 -818.9335 1806.353 6381.767 650.0613 -1941.1807 -4413.7332 [2,] -967.6032 7689.9335 2784.647 11285.233 3594.9387 332.1807 488.7332 [,8] [,9] [,10] [,11] [,12] [1,] 472.3341 -634.9173 -6044.7654 2617.027 1096.527 [2,] 2786.6659 981.9173 155.7654 6853.973 2833.473 $out [1] -662 -2943 $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(-3609, -2311, -1804.5, -636, 662, -320, 0, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/638yp1313684018.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] [,9] [1,] 31409.0 34035.0 47777 60191.0 70986.0 87670.0 110594.0 136750 163247.0 [2,] 31899.5 35824.0 50233 61835.0 73121.5 94714.5 120727.5 146736 174220.5 [3,] 33704.0 43677.5 56928 68548.5 83588.0 105512.5 132172.0 160321 190082.5 [4,] 34671.0 45813.0 58402 71004.0 86372.0 107965.5 134627.5 162774 193518.0 [5,] 34671.0 46795.0 59555 72630.0 89004.0 110594.0 136114.0 163598 195638.0 [,10] [,11] [1,] 191047.0 NA [2,] 202509.0 NA [3,] 219369.0 NA [4,] 225244.5 NA [5,] 225749.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 32439.9 39121.45 53202.06 64366.46 77544.36 99468.63 125832.1 153006 [2,] 34968.1 48233.55 60653.94 72730.54 89631.64 111556.37 138511.9 167636 [,9] [,10] [,11] [1,] 181280.8 208999.2 NA [2,] 198884.2 229738.8 NA $out [1] 39924 39577 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(31409, 31899.5, 33704, 34671, 34671, 34035, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/73hz01313684018.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,] 93731.4 79328.00 90628.50 [2,] 99199.4 83593.75 96315.25 [3,] 109833.4 94550.25 106742.50 [4,] 111647.7 96678.75 108137.62 [5,] 112612.2 99799.00 109403.00 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 104155.6 88582.09 101350.2 [2,] 115511.2 100518.41 112134.8 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(93731.4, 99199.4, 109833.4, 111647.7, 112612.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/18tul1313684018.ps tmp/18tul1313684018.png",intern=TRUE)) character(0) > try(system("convert tmp/2iqat1313684018.ps tmp/2iqat1313684018.png",intern=TRUE)) character(0) > try(system("convert tmp/38vmk1313684018.ps tmp/38vmk1313684018.png",intern=TRUE)) character(0) > try(system("convert tmp/4nozn1313684018.ps tmp/4nozn1313684018.png",intern=TRUE)) character(0) > try(system("convert tmp/5u82h1313684018.ps tmp/5u82h1313684018.png",intern=TRUE)) character(0) > try(system("convert tmp/638yp1313684018.ps tmp/638yp1313684018.png",intern=TRUE)) character(0) > try(system("convert tmp/73hz01313684018.ps tmp/73hz01313684018.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.534 0.308 1.847