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.0361910293346112,-0.130573884966384,0.0508918469276857,0.285401362098854,-0.144383808058777,-0.476406414178597,0.278108962203585,-0.0138995480760208,-0.222087591687949,0.209175033072651,0.0192739441821673,0.125857350082273,0.168126817416534,0.274716043721393,-0.115012818239837,-0.0528152166781983,0.0480099460498225,-0.271042927666295,-0.0342951777172496,0.0627605675580158,0.307448877628319,-0.0499740104559924,0.217631447221166,0.588245078469366,-0.104976668923445,0.0161271326299674,0.360593922797499,0.34885308223225,0.549125366865315,-0.149732507116426,-9.69800498562644e-05,-0.110243693744648,0.519793940008804,-0.222103155113966,-0.532213370012485,-0.373612664474175,-1.5811601864973,-0.226230409909392,-0.351379591715873,0.415088956092972,-0.551569595958837,0.00975519186112273,0.440089664424697,0.484637284027762,-0.101777761975125,-0.635031201521745,0.562248596196737,-0.52856053825146,0.56712813530769) > 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,] -0.03619103 0.16812682 -1.049767e-01 -1.581160186 0.5671281 [2,] -0.13057388 0.27471604 1.612713e-02 -0.226230410 NA [3,] 0.05089185 -0.11501282 3.605939e-01 -0.351379592 NA [4,] 0.28540136 -0.05281522 3.488531e-01 0.415088956 NA [5,] -0.14438381 0.04800995 5.491254e-01 -0.551569596 NA [6,] -0.47640641 -0.27104293 -1.497325e-01 0.009755192 NA [7,] 0.27810896 -0.03429518 -9.698005e-05 0.440089664 NA [8,] -0.01389955 0.06276057 -1.102437e-01 0.484637284 NA [9,] -0.22208759 0.30744888 5.197939e-01 -0.101777762 NA [10,] 0.20917503 -0.04997401 -2.221032e-01 -0.635031202 NA [11,] 0.01927394 0.21763145 -5.322134e-01 0.562248596 NA [12,] 0.12585735 0.58824508 -3.736127e-01 -0.528560538 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] -0.09438286 0.10658923 0.12110380 1.35492978 NA [2,] 0.18146573 -0.38972886 0.34446679 -0.12514918 NA [3,] 0.23450952 0.06219760 -0.01174084 0.76646855 NA [4,] -0.42978517 0.10082516 0.20027228 -0.96665855 NA [5,] -0.33202261 -0.31905287 -0.69885787 0.56132479 NA [6,] 0.75451538 0.23674775 0.14963553 0.43033447 NA [7,] -0.29200851 0.09705575 -0.11014671 0.04454762 NA [8,] -0.20818804 0.24468831 0.63003763 -0.58641505 NA [9,] 0.43126262 -0.35742289 -0.74189710 -0.53325344 NA [10,] -0.18990109 0.26760546 -0.31011021 1.19727980 NA [11,] 0.10658341 0.37061363 0.15860071 -1.09080913 NA [12,] 0.04226947 -0.69322175 -1.20754752 1.09568867 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/1ybnm1259831679.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/2z5531259831679.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/33xu61259831679.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/4or4l1259831679.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.10497667 -0.22623041 -0.35137959 -0.05281522 -0.55156960 -0.476406414 [2,] -0.10497667 -0.17840215 -0.23319620 0.11629307 -0.34797670 -0.373724671 [3,] -0.03619103 -0.05722338 -0.03206049 0.31712722 -0.04818693 -0.210387717 [4,] 0.16812682 0.14542159 0.20574288 0.38197102 0.29856766 -0.069988658 [5,] 0.56712814 0.27471604 0.36059392 0.41508896 0.54912537 0.009755192 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.03429518 -0.11024369 -0.2220876 -0.63503120 -0.5322134 -0.5285605 [2,] -0.01719608 -0.06207162 -0.1619327 -0.42856718 -0.2564697 -0.4510866 [3,] 0.13900599 0.02443051 0.1028356 -0.13603858 0.1184527 -0.1238777 [4,] 0.35909931 0.27369893 0.4136214 0.07960051 0.3899400 0.3570512 [5,] 0.44008966 0.48463728 0.5197939 0.20917503 0.5622486 0.5882451 $n [1] 5 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.2291653 -0.3130441 -0.3788224 0.1072416 -0.5589570 -0.45033917 [2,] 0.1567832 0.1985974 0.3147014 0.5270128 0.4625831 0.02956373 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.1582674 -0.2408282 -0.3518522 -0.5374911 -0.3922110 -0.7623065 [2,] 0.4362794 0.2896892 0.5575233 0.2654139 0.6291164 0.5145512 $out [1] -1.581160 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.104976668923445, -0.104976668923445, -0.0361910293346112, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5kqv31259831679.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.094382856 -0.38972886 -0.01174084 -0.9666586 -0.6988579 0.1496355 [2,] 0.006103185 -0.25743902 0.02522838 -0.6982219 -0.5154402 0.1931916 [3,] 0.113846514 0.02815828 0.14835356 -0.1644800 -0.3255377 0.3335411 [4,] 0.738016789 0.26296626 0.50048903 0.1505487 0.1211360 0.5924249 [5,] 1.354929777 0.34446679 0.76646855 0.2002723 0.5613248 0.7545154 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.29200851 -0.58641505 -0.74189710 -0.31011021 -1.0908091 -1.2075475 [2,] -0.20107761 -0.39730154 -0.63757527 -0.25000565 -0.4921129 -0.9503846 [3,] -0.03279955 0.01825013 -0.44533816 0.03885218 0.1325921 -0.3254761 [4,] 0.07080168 0.43736297 0.03691987 0.73244263 0.2646072 0.5689791 [5,] 0.09705575 0.63003763 0.43126262 1.19727980 0.3706136 1.0956887 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.4643652 -0.3829619 -0.2271024 -0.8350088 -0.8284329 0.01814682 [2,] 0.6920583 0.4392784 0.5238095 0.5060488 0.1773575 0.64893541 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.2475842 -0.6411348 -0.9781893 -0.7372820 -0.4652168 -1.5257735 [2,] 0.1819851 0.6776351 0.0875130 0.8149863 0.7304009 0.8748212 $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.0943828556317728, 0.0061031853365431, 0.113846513929136, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6tf361259831679.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.476406414 -0.27104293 -0.53221337 -1.5811602 0.5671281 [2,] -0.137478847 -0.05139461 -0.18591783 -0.5400651 0.5671281 [3,] 0.002687198 0.05538526 -0.05253682 -0.1640041 0.5671281 [4,] 0.167516192 0.24617375 0.35472350 0.4275893 0.5671281 [5,] 0.285401362 0.58824508 0.54912537 0.5622486 0.5671281 $n [1] 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.1364231 -0.08033767 -0.2991270 -0.6053577 0.5671281 [2,] 0.1417975 0.19110818 0.1940533 0.2773496 0.5671281 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.476406414178597, -0.137478846512581, 0.00268719805307325, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7733n1259831679.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.22185666 -0.21038772 -0.216122191 [2,] -0.11075051 -0.09055052 -0.061152252 [3,] -0.01510847 -0.03412576 0.004340751 [4,] 0.11582901 0.11064413 0.103466943 [5,] 0.24913205 0.31712722 0.283129634 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.11845291 -0.12589198 -0.07074317 [2,] 0.08823597 0.05764047 0.07942467 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.221856664275049, -0.110750513524131, -0.0151084698443676, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1ybnm1259831679.ps tmp/1ybnm1259831679.png") > system("convert tmp/2z5531259831679.ps tmp/2z5531259831679.png") > system("convert tmp/33xu61259831679.ps tmp/33xu61259831679.png") > system("convert tmp/4or4l1259831679.ps tmp/4or4l1259831679.png") > system("convert tmp/5kqv31259831679.ps tmp/5kqv31259831679.png") > system("convert tmp/6tf361259831679.ps tmp/6tf361259831679.png") > system("convert tmp/7733n1259831679.ps tmp/7733n1259831679.png") > > > proc.time() user system elapsed 1.164 0.909 1.728