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.0302547091228863,0.267612748598290,-0.05258221327681,-0.164815126360942,0.145251461046347,-0.159755000728184,0.0335569671830744,0.0880735117085241,-0.281966491596610,0.0474120142758863,-0.484987801933979,0.128597017186614,-0.0760862380652182,0.0847658982829797,-0.0514027132036804,-0.0566616141609282,0.069052111883504,0.00440676426158568,-0.0195453629146749,0.240643182454573,-0.33572474280097,0.0460275678865050,0.0696491220369025,-0.225076263667377,-0.161613617181671,0.306153117443247,-0.155979433289841,0.179413008162554,-0.0918257869952182,-0.0515872033510411,-0.0259692959341839,-0.341986993998746,0.248411293223152,0.610658133668028,0.207016437881814,0.0640914582469286,0.0207692245074071,-0.202335247943092,0.0784047280219788,0.0426875319465297,0.376752985339314,-0.0782742138517369,0.238382692897949,-0.0408923157917773,0.0426863373126538,-0.00863223279670068,-0.288046674979108,0.256083528841667) > 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,] 0.03025471 -0.076086238 -0.16161362 0.020769225 NA [2,] 0.26761275 0.084765898 0.30615312 -0.202335248 NA [3,] -0.05258221 -0.051402713 -0.15597943 0.078404728 NA [4,] -0.16481513 -0.056661614 0.17941301 0.042687532 NA [5,] 0.14525146 0.069052112 -0.09182579 0.376752985 NA [6,] -0.15975500 0.004406764 -0.05158720 -0.078274214 NA [7,] 0.03355697 -0.019545363 -0.02596930 0.238382693 NA [8,] 0.08807351 0.240643182 -0.34198699 -0.040892316 NA [9,] -0.28196649 -0.335724743 0.24841129 0.042686337 NA [10,] 0.04741201 0.046027568 0.61065813 -0.008632233 NA [11,] -0.48498780 0.069649122 0.20701644 -0.288046675 NA [12,] 0.12859702 -0.225076264 0.06409146 0.256083529 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.23735804 0.160852136 0.46776673 -0.22310447 NA [2,] -0.32019496 -0.136168611 -0.46213255 0.28073998 NA [3,] -0.11223291 -0.005258901 0.33539244 -0.03571720 NA [4,] 0.31006659 0.125713726 -0.27123880 0.33406545 NA [5,] -0.30500646 -0.064645348 0.04023858 -0.45502720 NA [6,] 0.19331197 -0.023952127 0.02561791 0.31665691 NA [7,] 0.05451654 0.260188545 -0.31601770 -0.27927501 NA [8,] -0.37004000 -0.576367925 0.59039829 0.08357865 NA [9,] 0.32937851 0.381752311 0.36224684 -0.05131857 NA [10,] -0.53239982 0.023621554 -0.40364170 -0.27941444 NA [11,] 0.61358482 -0.294725386 -0.14292498 0.54413020 NA [12,] -0.20468326 0.063462646 -0.04332223 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/189e21260227885.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/2opsl1260227885.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/3kly51260227886.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/44ebf1260227886.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.16161362 -0.20233525 -0.15597943 -0.164815126 -0.09182579 -0.159755001 [2,] -0.11884993 -0.05878467 -0.10428082 -0.110738370 -0.01138684 -0.119014607 [3,] -0.02765851 0.17618932 -0.05199246 -0.006987041 0.10715179 -0.064930709 [4,] 0.02551197 0.28688293 0.01350101 0.111050270 0.26100222 -0.023590220 [5,] 0.03025471 0.30615312 0.07840473 0.179413008 0.37675299 0.004406764 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.025969296 -0.3419870 -0.3357247 -0.008632233 -0.4849878 -0.22507626 [2,] -0.022757329 -0.1914397 -0.3088456 0.018697668 -0.3865172 -0.08049240 [3,] 0.007005802 0.0235906 -0.1196401 0.046719791 -0.1091988 0.09634424 [4,] 0.135969830 0.1643583 0.1455488 0.329035074 0.1383328 0.19234027 [5,] 0.238382693 0.2406432 0.2484113 0.610658134 0.2070164 0.25608353 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.14170440 -0.09688809 -0.14504011 -0.1822001 -0.1080356 -0.14031597 [2,] 0.08638739 0.44926673 0.04105518 0.1682260 0.3223391 0.01045456 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.1183887 -0.2574898 -0.4786117 -0.1984468 -0.5238303 -0.1191936 [2,] 0.1324003 0.3046710 0.2393315 0.2918863 0.3054327 0.3118821 $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.161613617181671, -0.118849927623445, -0.0276585067789056, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/52vbq1260227886.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.22310447 -0.46213255 -0.11223291 -0.27123880 -0.45502720 -0.0239521272 [2,] -0.03112617 -0.39116376 -0.07397505 -0.07276253 -0.38001683 0.0008328901 [3,] 0.19910509 -0.22818179 -0.02048805 0.21789016 -0.18482590 0.1094649377 [4,] 0.35256239 0.07228568 0.16506677 0.32206602 -0.01220338 0.2549844373 [5,] 0.46776673 0.28073998 0.33539244 0.33406545 0.04023858 0.3166569067 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.3160177 -0.5763679 -0.05131857 -0.53239982 -0.2947254 -0.20468326 [2,] -0.2976464 -0.4732040 0.13902997 -0.46802076 -0.2188252 -0.12400274 [3,] -0.1123792 -0.1432307 0.34581267 -0.34152807 0.2006026 -0.04332223 [4,] 0.1573525 0.3369885 0.37199958 -0.12789644 0.5788575 0.01007021 [5,] 0.2601885 0.5903983 0.38175231 0.02362155 0.6135848 0.06346265 $n [1] 4 4 4 4 4 4 4 4 4 4 4 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.1040089 -0.5943068 -0.2093311 -0.0940244 -0.4753985 -0.09131478 [2,] 0.5022190 0.1379433 0.1683550 0.5298047 0.1057467 0.31024466 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.4718284 -0.7832827 0.1617667 -0.61022628 -0.4295667 -0.16562538 [2,] 0.2470699 0.4968213 0.5298587 -0.07282986 0.8307719 0.07898091 $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.223104472450499, -0.0311261680511506, 0.199105087911801, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/68bon1260227886.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.48498780 -0.22507626 -0.34198699 -0.28804667 NA [2,] -0.16228506 -0.06637393 -0.12390261 -0.05958326 NA [3,] 0.03190584 -0.00756930 0.01906108 0.03172778 NA [4,] 0.10833526 0.06935062 0.22771387 0.15839371 NA [5,] 0.26761275 0.24064318 0.61065813 0.37675299 NA $n [1] 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.09152591 -0.06947417 -0.1413136 -0.06769298 NA [2,] 0.15533759 0.05433558 0.1794357 0.13114854 NA $out [1] -0.3357247 $group [1] 2 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.484987801933979, -0.162285063544563, 0.0319058381529803, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7zul91260227886.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.124092229 -1.196401e-01 -0.1166455029 [2,] -0.058985697 -5.846159e-02 -0.0584038733 [3,] -0.006692352 9.380514e-06 0.0008047132 [4,] 0.085327690 7.153201e-02 0.0932135837 [5,] 0.173866371 1.761893e-01 0.1451192263 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.07251466 -0.05928157 -0.06834903 [2,] 0.05912995 0.05930033 0.06995845 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.124092229248593, -0.0589856969107465, -0.00669235200502659, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/189e21260227885.ps tmp/189e21260227885.png") > system("convert tmp/2opsl1260227885.ps tmp/2opsl1260227885.png") > system("convert tmp/3kly51260227886.ps tmp/3kly51260227886.png") > system("convert tmp/44ebf1260227886.ps tmp/44ebf1260227886.png") > system("convert tmp/52vbq1260227886.ps tmp/52vbq1260227886.png") > system("convert tmp/68bon1260227886.ps tmp/68bon1260227886.png") > system("convert tmp/7zul91260227886.ps tmp/7zul91260227886.png") > > > proc.time() user system elapsed 1.163 0.899 1.453