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(100 + ,108.1560276 + ,114.0150276 + ,102.1880309 + ,110.3672031 + ,96.8602511 + ,94.1944583 + ,99.51621961 + ,94.06333487 + ,97.5541476 + ,78.15062422 + ,81.2434643 + ,92.36262465 + ,96.06324371 + ,114.0523777 + ,110.6616666 + ,104.9171949 + ,90.00187193 + ,95.7008067 + ,86.02741157 + ,84.85287668 + ,100.04328 + ,80.91713823 + ,74.06539709 + ,77.30281369 + ,97.23043249 + ,90.75515676 + ,100.5614455 + ,92.01293267 + ,99.24012138 + ,105.8672755 + ,90.9920463 + ,93.30624423 + ,91.17419413 + ,77.33295039 + ,91.1277721 + ,85.01249943 + ,83.90390242 + ,104.8626302 + ,110.9039108 + ,95.43714373 + ,111.6238727 + ,108.8925403 + ,96.17511682 + ,101.9740205 + ,99.11953031 + ,86.78158147 + ,118.4195003 + ,118.7441447 + ,106.5296192 + ,134.7772694 + ,104.6778714 + ,105.2954304 + ,139.4139849 + ,103.6060491 + ,99.78182974 + ,103.4610301 + ,120.0594945 + ,96.71377168 + ,107.1308929 + ,105.3608372 + ,111.6942359 + ,132.0519998 + ,126.8037879 + ,154.4824253 + ,141.5570984 + ,109.9506882 + ,127.904198 + ,133.0888617 + ,120.0796299 + ,117.5557142 + ,143.0362309 + ,159.982927 + ,128.5991124 + ,149.7373327 + ,126.8169313 + ,140.9639674 + ,137.6691981 + ,117.9402337 + ,122.3095247 + ,127.7804207 + ,136.1677176 + ,116.2405856 + ,123.1576893 + ,116.3400234 + ,108.6119282 + ,125.8982264 + ,112.8003105 + ,107.5182447 + ,135.0955413 + ,115.5096488 + ,115.8640759 + ,104.5883906 + ,163.7213386 + ,113.4482275 + ,98.0428844 + ,116.7868521 + ,126.5330444 + ,113.0336597 + ,124.3392163 + ,109.8298759 + ,124.4434777 + ,111.5039454 + ,102.0350019 + ,116.8726598 + ,112.2073122 + ,101.1513902 + ,124.4255108) > 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] 108 > (np <- floor(n / par1)) [1] 9 > 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] 9 9 9 9 9 9 9 9 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 100.00000 92.36262 77.30281 85.01250 118.74414 105.3608 159.9829 [2,] 108.15603 96.06324 97.23043 83.90390 106.52962 111.6942 128.5991 [3,] 114.01503 114.05238 90.75516 104.86263 134.77727 132.0520 149.7373 [4,] 102.18803 110.66167 100.56145 110.90391 104.67787 126.8038 126.8169 [5,] 110.36720 104.91719 92.01293 95.43714 105.29543 154.4824 140.9640 [6,] 96.86025 90.00187 99.24012 111.62387 139.41398 141.5571 137.6692 [7,] 94.19446 95.70081 105.86728 108.89254 103.60605 109.9507 117.9402 [8,] 99.51622 86.02741 90.99205 96.17512 99.78183 127.9042 122.3095 [9,] 94.06333 84.85288 93.30624 101.97402 103.46103 133.0889 127.7804 [10,] 97.55415 100.04328 91.17419 99.11953 120.05949 120.0796 136.1677 [11,] 78.15062 80.91714 77.33295 86.78158 96.71377 117.5557 116.2406 [12,] 81.24346 74.06540 91.12777 118.41950 107.13089 143.0362 123.1577 [,8] [,9] [,10] [1,] 116.34002 116.7869 NA [2,] 108.61193 126.5330 NA [3,] 125.89823 113.0337 NA [4,] 112.80031 124.3392 NA [5,] 107.51824 109.8299 NA [6,] 135.09554 124.4435 NA [7,] 115.50965 111.5039 NA [8,] 115.86408 102.0350 NA [9,] 104.58839 116.8727 NA [10,] 163.72134 112.2073 NA [11,] 113.44823 101.1514 NA [12,] 98.04288 124.4255 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 8.156028 3.700619 19.927619 -1.1085970 -12.214526 6.333399 [2,] 5.859000 17.989134 -6.475276 20.9587278 28.247650 20.357764 [3,] -11.826997 -3.390711 9.806289 6.0412806 -30.099398 -5.248212 [4,] 8.179172 -5.744472 -8.548513 -15.4667671 0.617559 27.678637 [5,] -13.506952 -14.915323 7.227189 16.1867290 34.118555 -12.925327 [6,] -2.665793 5.698935 6.627154 -2.7313324 -35.807936 -31.606410 [7,] 5.321761 -9.673395 -14.875229 -12.7174235 -3.824219 17.953510 [8,] -5.452885 -1.174535 2.314198 5.7989037 3.679200 5.184664 [9,] 3.490813 15.190403 -2.132050 -2.8544902 16.598464 -13.009232 [10,] -19.403523 -19.126142 -13.841244 -12.3379488 -23.345723 -2.523916 [11,] 3.092840 -6.851741 13.794822 31.6379188 10.417121 25.480517 [12,] 11.119160 3.237417 -6.115273 0.3246444 -1.770056 16.946696 [,7] [,8] [,9] [,10] [1,] -31.383815 -7.7280952 9.746192 NA [2,] 21.138220 17.2862982 -13.499385 NA [3,] -22.920401 -13.0979159 11.305557 NA [4,] 14.147036 -5.2820658 -14.509340 NA [5,] -3.294769 27.5772966 14.613602 NA [6,] -19.728964 -19.5858925 -12.939532 NA [7,] 4.369291 0.3544271 -9.468944 NA [8,] 5.470896 -11.2756853 14.837658 NA [9,] 8.387297 59.1329480 -4.665348 NA [10,] -19.927132 -50.2731111 -11.055922 NA [11,] 6.917104 -15.4053431 23.274121 NA [12,] -6.817666 18.7439677 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/1dx2m1257854730.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/2tuu11257854730.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/36j7z1257854730.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/4v6f41257854730.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] [1,] 77.30281 83.90390 90.75516 100.5614 104.9172 90.00187 94.19446 [2,] 92.36262 97.23043 113.03366 104.6779 104.9172 99.24012 103.60605 [3,] 105.36084 108.15603 114.05238 110.9039 107.5182 124.44348 108.89254 [4,] 116.78685 111.69424 132.05200 124.3392 110.3672 137.66920 111.50395 [5,] 118.74414 128.59911 149.73733 126.8169 110.3672 141.55710 117.94023 [,8] [,9] [,10] [,11] [,12] [1,] 86.02741 84.85288 91.17419 77.33295 74.06540 [2,] 96.17512 94.06333 99.11953 80.91714 91.12777 [3,] 99.78183 103.46103 112.20731 96.71377 107.13089 [4,] 115.86408 116.87266 120.07963 113.44823 123.15769 [5,] 127.90420 133.08886 136.16772 117.55571 143.03623 $n [1] 9 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 92.49741 100.5384 104.0361 100.5489 104.6479 104.2042 104.7330 89.41231 [2,] 118.22426 115.7736 124.0687 121.2589 110.3886 144.6828 113.0521 110.15135 [,9] [,10] [,11] [,12] [1,] 91.44812 101.1683 79.58073 90.2618 [2,] 115.47394 123.2463 113.84681 124.0000 $out [1] 159.98293 92.01293 95.43714 154.48243 140.96397 163.72134 $group [1] 1 5 5 5 5 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(77.30281369, 92.36262465, 105.3608372, 116.7868521, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5xjmb1257854730.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] [1,] -31.383815 -13.49938 -30.099398 -15.466767 -14.915323 -35.807936 [2,] -7.728095 5.85900 -13.097916 -8.548513 -12.925327 -19.728964 [3,] 3.700619 17.98913 -5.248212 -5.282066 7.227189 -12.939532 [4,] 8.156028 20.95873 6.041281 8.179172 16.186729 -2.665793 [5,] 19.927619 28.24765 11.305557 27.678637 34.118555 6.627154 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -14.875229 -5.452885 -13.009232 -23.345723 -15.40534 -6.817666 [2,] -9.673395 -1.174535 -2.854490 -19.927132 3.09284 -3.942664 [3,] -3.824219 3.679200 3.490813 -19.126142 10.41712 1.781030 [4,] 4.369291 5.470896 15.190403 -12.337949 23.27412 14.032928 [5,] 17.953510 14.837658 16.598464 -2.523916 31.63792 18.743968 $n [1] 9 9 9 9 9 9 9 9 9 9 9 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -4.665019 10.03661 -15.328189 -14.091980 -8.10516 -21.926136 -11.220034 [2,] 12.066257 25.94166 4.831765 3.527848 22.55954 -3.952929 3.571595 [,8] [,9] [,10] [,11] [,12] [1,] 0.1792734 -6.012831 -23.12311 -0.2116865 -8.260393 [2,] 7.1791273 12.994457 -15.12917 21.0459290 11.822454 $out [1] -11.27569 59.13295 -50.27311 $group [1] 8 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-31.3838146, -7.72809520000001, 3.70061905999999, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6qdrc1257854730.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,] 78.15062 74.06540 90.75516 83.90390 96.71377 105.3608 116.2406 [2,] 94.12890 85.44014 90.87360 91.10936 103.53354 114.6250 122.7336 [3,] 98.53518 94.03172 91.59356 100.54678 105.91252 127.3540 128.1898 [4,] 105.17203 102.48024 98.23528 109.89823 119.40182 137.3230 139.3166 [5,] 114.01503 114.05238 105.86728 118.41950 139.41398 154.4824 159.9829 [,8] [,9] [,10] [1,] 98.04288 101.1514 NA [2,] 108.06509 110.6669 NA [3,] 114.47894 114.9103 NA [4,] 121.11912 124.3824 NA [5,] 135.09554 126.5330 NA $n [1] 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 93.49834 86.25961 88.23585 91.97705 98.6749 117.0013 120.6262 108.5249 [2,] 103.57203 101.80382 94.95127 109.11650 113.1502 137.7067 135.7534 120.4330 [,9] [,10] [1,] 108.6545 NA [2,] 121.1660 NA $out [1] 77.30281 77.33295 163.72134 $group [1] 3 3 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(78.15062422, 94.128896585, 98.535183605, 105.17202925, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/75wfz1257854730.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,] 96.47689 96.71377 97.18268 [2,] 106.70207 104.41093 105.02137 [3,] 107.73413 107.83714 107.34886 [4,] 114.49728 111.55561 112.05406 [5,] 119.90930 114.05238 122.54283 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 104.1787 104.5784 104.1412 [2,] 111.2896 111.0959 110.5565 $out [1] 124.4435 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(96.4768870544444, 106.702065626111, 107.734126027222, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1dx2m1257854730.ps tmp/1dx2m1257854730.png") > system("convert tmp/2tuu11257854730.ps tmp/2tuu11257854730.png") > system("convert tmp/36j7z1257854730.ps tmp/36j7z1257854730.png") > system("convert tmp/4v6f41257854730.ps tmp/4v6f41257854730.png") > system("convert tmp/5xjmb1257854730.ps tmp/5xjmb1257854730.png") > system("convert tmp/6qdrc1257854730.ps tmp/6qdrc1257854730.png") > system("convert tmp/75wfz1257854730.ps tmp/75wfz1257854730.png") > > > proc.time() user system elapsed 1.258 0.986 1.742