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(-0.0175986496695149,4.88533001639715e-06,-0.0797142045361252,0.0380359724366437,0.00379990057404902,0.0592823368877745,-0.0401047028912839,-0.0903371585255128,0.0600435985327014,-0.0751519168673696,0.184740243449584,0.108923749072312,0.148749967789766,0.0299253269806808,-0.132723413231867,0.130073522884216,-0.0524172791064149,-0.0488662775462986,0.0210676927748520,-0.165146636794414,-0.0573375413155754,0.00334278216483919,-0.120892947016097,0.0140181973103131,-0.166333957140904,-0.204256136740815,-0.170162160910135,-0.40685370432818,0.089890889862059,0.121852449126500,0.0718684968241188,0.146330323762528,-0.254402885578116,-0.175254145082071,0.0883657966655966,0.0494861859369645,-0.179560126045166,0.0352805435009704,0.107242027471589,0.0159604480893762,0.229492599513338,-0.472622684816943,-0.0310496393734398,0.0823390463619987,-0.0459447059808648,-0.0671207054879601,0.062392997118939,0.0541581853860086,-0.000768654964170624) > 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] 49 > (np <- floor(n / par1)) [1] 4 > 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] 5 4 4 4 4 4 4 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] -1.759865e-02 0.148749968 -0.16633396 -0.17956013 -0.000768655 [2,] 4.885330e-06 0.029925327 -0.20425614 0.03528054 NA [3,] -7.971420e-02 -0.132723413 -0.17016216 0.10724203 NA [4,] 3.803597e-02 0.130073523 -0.40685370 0.01596045 NA [5,] 3.799901e-03 -0.052417279 0.08989089 0.22949260 NA [6,] 5.928234e-02 -0.048866278 0.12185245 -0.47262268 NA [7,] -4.010470e-02 0.021067693 0.07186850 -0.03104964 NA [8,] -9.033716e-02 -0.165146637 0.14633032 0.08233905 NA [9,] 6.004360e-02 -0.057337541 -0.25440289 -0.04594471 NA [10,] -7.515192e-02 0.003342782 -0.17525415 -0.06712071 NA [11,] 1.847402e-01 -0.120892947 0.08836580 0.06239300 NA [12,] 1.089237e-01 0.014018197 0.04948619 0.05415819 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.01760353 -0.118824641 -0.03792218 0.214840670 NA [2,] -0.07971909 -0.162648740 0.03409398 0.071961484 NA [3,] 0.11775018 0.262796936 -0.23669154 -0.091281579 NA [4,] -0.03423607 -0.182490802 0.49674459 0.213532151 NA [5,] 0.05548244 0.003551002 0.03196156 -0.702115284 NA [6,] -0.09938704 0.069933970 -0.04998395 0.441573045 NA [7,] -0.05023246 -0.186214330 0.07446183 0.113388686 NA [8,] 0.15038076 0.107809095 -0.40073321 -0.128283752 NA [9,] -0.13519552 0.060680323 0.07914874 -0.021176000 NA [10,] 0.25989216 -0.124235729 0.26361994 0.129513703 NA [11,] -0.07581649 0.134911144 -0.03887961 -0.008234812 NA [12,] 0.03982622 -0.180352154 -0.22904631 -0.054926840 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/1w7qw1260563588.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/2dqmw1260563588.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/39cr71260563588.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/4vjos1260563588.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] [1,] -0.179560126 -0.20425614 -0.17016216 -0.40685370 -0.05241728 -0.47262268 [2,] -0.166333957 -0.10212563 -0.15144279 -0.19544663 -0.02430869 -0.26074448 [3,] -0.017598650 0.01496511 -0.10621881 0.02699821 0.04684540 0.00520803 [4,] -0.000768655 0.03260294 0.01376391 0.08405475 0.15969174 0.09056739 [5,] 0.148749968 0.03528054 0.10724203 0.13007352 0.22949260 0.12185245 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.040104703 -0.165146637 -0.254402886 -0.175254145 -0.12089295 0.01401820 [2,] -0.035577171 -0.127741898 -0.155870213 -0.125203031 -0.02924997 0.03175219 [3,] -0.004990973 -0.003999056 -0.051641124 -0.071136311 0.07537940 0.05182219 [4,] 0.046468095 0.114334685 0.007049446 -0.031888962 0.13655302 0.08154097 [5,] 0.071868497 0.146330324 0.060043599 0.003342782 0.18474024 0.10892375 $n [1] 5 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.13458668 -0.09147046 -0.23673210 -0.1938079 -0.09851495 -0.2723284 [2,] 0.09938938 0.12140067 0.02429448 0.2478043 0.19220574 0.2827444 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.06980673 -0.1952396 -0.1803477 -0.144854426 -0.05560497 0.01248905 [2,] 0.05982479 0.1872414 0.0770654 0.002581804 0.20636376 0.09115532 $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(-0.179560126045166, -0.166333957140904, -0.0175986496695149, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5vudj1260563588.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,] -0.11882464 -0.16264874 -0.2366915 -0.18249080 -0.70211528 -0.09938704 [2,] -0.07837341 -0.12118392 -0.1639866 -0.10836344 -0.34928214 -0.07468550 [3,] -0.01015932 -0.02281256 0.0132343 0.08964804 0.01775628 0.00997501 [4,] 0.11622210 0.05302773 0.1902736 0.35513837 0.04372200 0.25575351 [5,] 0.21484067 0.07196148 0.2627969 0.49674459 0.05548244 0.44157305 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.18621433 -0.40073321 -0.13519552 -0.124235729 -0.07581649 -0.229046312 [2,] -0.11822339 -0.26450848 -0.07818576 0.002638987 -0.05734805 -0.204699233 [3,] 0.01211469 -0.01023733 0.01975216 0.194702931 -0.02355721 -0.117639497 [4,] 0.09392526 0.12909493 0.06991453 0.261756051 0.06333817 -0.007550311 [5,] 0.11338869 0.15038076 0.07914874 0.263619942 0.13491114 0.039826219 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.1638898 -0.1604398 -0.2666312 -0.2765184 -0.2927170 -0.2510718 [2,] 0.1435711 0.1148146 0.2930998 0.4558145 0.3282296 0.2710218 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.1554827 -0.3211840 -0.09724707 -0.00999955 -0.1188993 -0.27338715 [2,] 0.1797121 0.3007094 0.13675139 0.39940541 0.0717849 0.03810815 $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(-0.118824640809085, -0.0783734102044981, -0.0101593223001899, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ksqx1260563588.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] [1,] -0.090337159 -0.16514664 -0.40685370 -0.17956013 -0.000768655 [2,] -0.057628310 -0.08911524 -0.18975514 -0.05653271 -0.000768655 [3,] 0.001902393 -0.02276175 -0.05842389 0.02562050 -0.000768655 [4,] 0.059662968 0.02549651 0.08912834 0.07236602 -0.000768655 [5,] 0.184740243 0.14874997 0.14633032 0.22949260 -0.000768655 $n [1] 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.05159495 -0.07503694 -0.18562451 -0.03317108 -0.000768655 [2,] 0.05539973 0.02951344 0.06877674 0.08441207 -0.000768655 $out [1] -0.4726227 $group [1] 4 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.0903371585255128, -0.0576283098793268, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7vs6w1260563588.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,] -0.08508854 -0.1062188089 -0.08752912 [2,] -0.07162491 -0.0346198867 -0.06893345 [3,] -0.03893181 0.0006044868 -0.01212349 [4,] 0.02954849 0.0369218027 0.02723081 [5,] 0.06769153 0.0753793969 0.06451546 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.085077683 -0.03202616 -0.05598466 [2,] 0.007214054 0.03323513 0.03173768 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.0850885440872418, -0.0716249106935491, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > system("convert tmp/1w7qw1260563588.ps tmp/1w7qw1260563588.png") > system("convert tmp/2dqmw1260563588.ps tmp/2dqmw1260563588.png") > system("convert tmp/39cr71260563588.ps tmp/39cr71260563588.png") > system("convert tmp/4vjos1260563588.ps tmp/4vjos1260563588.png") > system("convert tmp/5vudj1260563588.ps tmp/5vudj1260563588.png") > system("convert tmp/6ksqx1260563588.ps tmp/6ksqx1260563588.png") > system("convert tmp/7vs6w1260563588.ps tmp/7vs6w1260563588.png") > > > proc.time() user system elapsed 1.198 0.925 1.585