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.0322477563269516,-3.10259255310256e-05,0.301112186970538,0.0548018747461919,-0.133474622792253,0.225291776937562,0.000852895716111257,0.084341638102194,0.175638628134376,-0.245187524575567,0.0183476060840533,-0.440793105929656,0.05602794751164,-0.0510397796031188,-0.00363943386617446,-0.0951080672001434,-0.103351045622180,0.0742974455944183,0.0310716492059746,-0.039678164485403,0.227962566341790,-0.324561783178065,-0.039986746601846,0.115065195584166,-0.241826610194115,-0.265010543596649,0.264312619508534,-0.145748811773052,0.074984846215623,-0.0792037460618129,-0.0897138675350125,-0.0318891439749827,-0.379520044986777,0.147081259023469,0.644573158872214,0.231406353542911,0.0257968317972824,0.102825006570345,-0.0644471569186438,0.143283137022815,0.088575110656249,0.345773687682914,-0.0556153301369342,0.236619255674046,0.109662493707481,0.137836003726974,0.0941980971503112,-0.255678058936345) > 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] 48 > (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] 4 4 4 4 4 4 4 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] -3.224776e-02 0.056027948 -0.24182661 0.02579683 NA [2,] -3.102593e-05 -0.051039780 -0.26501054 0.10282501 NA [3,] 3.011122e-01 -0.003639434 0.26431262 -0.06444716 NA [4,] 5.480187e-02 -0.095108067 -0.14574881 0.14328314 NA [5,] -1.334746e-01 -0.103351046 0.07498485 0.08857511 NA [6,] 2.252918e-01 0.074297446 -0.07920375 0.34577369 NA [7,] 8.528957e-04 0.031071649 -0.08971387 -0.05561533 NA [8,] 8.434164e-02 -0.039678164 -0.03188914 0.23661926 NA [9,] 1.756386e-01 0.227962566 -0.37952004 0.10966249 NA [10,] -2.451875e-01 -0.324561783 0.14708126 0.13783600 NA [11,] 1.834761e-02 -0.039986747 0.64457316 0.09419810 NA [12,] -4.407931e-01 0.115065196 0.23140635 -0.25567806 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.03221673 -0.107067727 -0.02318393 0.07702817 NA [2,] 0.30114321 0.047400346 0.52932316 -0.16727216 NA [3,] -0.24631031 -0.091468633 -0.41006143 0.20773029 NA [4,] -0.18827650 -0.008242978 0.22073366 -0.05470803 NA [5,] 0.35876640 0.177648491 -0.15418859 0.25719858 NA [6,] -0.22443888 -0.043225796 -0.01051012 -0.40138902 NA [7,] 0.08348874 -0.070749814 0.05782472 0.29223459 NA [8,] 0.09129699 0.267640731 -0.34763090 -0.12695676 NA [9,] -0.42082615 -0.552524350 0.52660130 0.02817351 NA [10,] 0.26353513 0.284575037 0.49749190 -0.04363791 NA [11,] -0.45914071 0.155051942 -0.41316681 -0.34987616 NA [12,] 0.49682105 -0.356891806 -0.20560952 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/19gth1260392963.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/25rtr1260392963.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/31no21260392963.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/4g8hu1260392963.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.241826610 -0.26501054 -0.06444716 -0.14574881 -0.13347462 -0.07920375 [2,] -0.137037183 -0.15802516 -0.03404330 -0.12042844 -0.11841283 -0.00245315 [3,] -0.003225462 -0.02553540 0.13033659 -0.02015310 -0.01418310 0.14979461 [4,] 0.040912390 0.05139699 0.28271240 0.09904251 0.08177998 0.28553273 [5,] 0.056027948 0.10282501 0.30111219 0.14328314 0.08857511 0.34577369 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.08971387 -0.03967816 -0.3795200 -0.32456178 -0.03998675 -0.44079311 [2,] -0.07266460 -0.03578365 -0.1349288 -0.28487465 -0.01081957 -0.34823558 [3,] -0.02738122 0.02622625 0.1426506 -0.05367576 0.05627285 -0.07030643 [4,] 0.01596227 0.16048045 0.2018006 0.14245863 0.36938563 0.17323577 [5,] 0.03107165 0.23661926 0.2279626 0.14708126 0.64457316 0.23140635 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.1438056 -0.1909789 -0.1199004 -0.1935351 -0.1723354 -0.07771424 [2,] 0.1373547 0.1399081 0.3805736 0.1532290 0.1439692 0.37730346 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.09739645 -0.1288224 -0.1233656 -0.3912691 -0.2440893 -0.4822688 [2,] 0.04263401 0.1812749 0.4086668 0.2839175 0.3566350 0.3416559 $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.241826610194115, -0.137037183260533, -0.0032254622648346, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5jr3k1260392963.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.107067727 -0.16727216 -0.41006143 -0.1882765 -0.15418859 -0.40138902 [2,] -0.065125830 -0.05993591 -0.32818587 -0.1214923 0.01172995 -0.31291395 [3,] 0.004516398 0.17427178 -0.16888947 -0.0314755 0.21742353 -0.13383234 [4,] 0.054622453 0.41523319 0.05813083 0.1062453 0.30798249 -0.02686796 [5,] 0.077028175 0.52932316 0.20773029 0.2207337 0.35876640 -0.01051012 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.070749814 -0.34763090 -0.5525243 -0.04363791 -0.4591407 -0.3568918 [2,] -0.006462545 -0.23729383 -0.4866753 0.10994861 -0.4361538 -0.2812507 [3,] 0.070656733 -0.01782989 -0.1963263 0.27405508 -0.3815215 -0.2056095 [4,] 0.187861664 0.17946886 0.2773874 0.39103347 -0.0974121 0.1456058 [5,] 0.292234586 0.26764073 0.5266013 0.49749190 0.1550519 0.4968211 $n [1] 4 4 4 4 4 4 4 4 4 4 4 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.09008474 -0.2011118 -0.4740797 -0.2113882 -0.01661597 -0.3598087 [2,] 0.09911754 0.5496554 0.1363007 0.1484372 0.45146304 0.0921440 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.0828594 -0.3470724 -0.7999358 0.05199805 -0.6491274 -0.5949937 [2,] 0.2241729 0.3114126 0.4072832 0.49611212 -0.1139156 0.1837746 $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.107067727114759, -0.0651258302586464, 0.00451639849944329, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/675k91260392963.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.24518752 -0.10335105 -0.37952004 -0.06444716 NA [2,] -0.08286119 -0.07307392 -0.19378771 -0.01490925 NA [3,] 0.00960025 -0.02165880 -0.05554645 0.09851155 NA [4,] 0.12999013 0.06516270 0.18924381 0.14055957 NA [5,] 0.30111219 0.22796257 0.64457316 0.34577369 NA $n [1] 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.08748267 -0.08470945 -0.2302497 0.02760118 NA [2,] 0.10668317 0.04139185 0.1191568 0.16942192 NA $out [1] -0.4407931 -0.3245618 -0.2556781 $group [1] 1 2 4 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.245187524575567, -0.0828611895596023, 0.00960025090008228, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7fjwr1260392963.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.08749990 -0.070306432 -0.07890317 [2,] -0.05068824 -0.026458310 -0.03364547 [3,] -0.01450470 -0.008704281 -0.01583640 [4,] 0.09334148 0.093304722 0.10291059 [5,] 0.17928303 0.149794611 0.14566720 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.08019762 -0.06332900 -0.07812053 [2,] 0.05118822 0.04592044 0.04644774 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.087499903934731, -0.0506882412208872, -0.0145046973433437, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/19gth1260392963.ps tmp/19gth1260392963.png") > system("convert tmp/25rtr1260392963.ps tmp/25rtr1260392963.png") > system("convert tmp/31no21260392963.ps tmp/31no21260392963.png") > system("convert tmp/4g8hu1260392963.ps tmp/4g8hu1260392963.png") > system("convert tmp/5jr3k1260392963.ps tmp/5jr3k1260392963.png") > system("convert tmp/675k91260392963.ps tmp/675k91260392963.png") > system("convert tmp/7fjwr1260392963.ps tmp/7fjwr1260392963.png") > > > proc.time() user system elapsed 1.232 0.973 2.925