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.000746099568952258 + ,0.0292111461738571 + ,-0.00280267596084322 + ,-0.00184368563002967 + ,0.0107457257770334 + ,-0.0175566443253344 + ,0.0131165820445929 + ,-0.0324353168918367 + ,0.0224382020922924 + ,-0.0196217959970839 + ,0.00495492412026389 + ,0.00168574453846993 + ,-0.00506586405101019 + ,-0.0116484658020629 + ,-0.00266112877646428 + ,0.0179412217929573 + ,-0.0168120805195169 + ,-0.0147697069296257 + ,0.035665172596409 + ,0.00293868331101443 + ,-0.0145835439671676 + ,0.0184871974956284 + ,-0.00672423964360273 + ,-0.00834464807562108 + ,0.0169769773776811 + ,0.0141585136611232 + ,0.00520840412118636 + ,-0.0225065086691335 + ,0.0221256937859677 + ,0.008638135296458 + ,0.0259853418092942 + ,-0.014962664383958 + ,0.0336773377780083 + ,-0.00849413990510539 + ,-0.0229882773676451 + ,0.0569000711376357 + ,0.025309869239244 + ,-0.029512999542613 + ,-0.00750220761300227 + ,0.00584857420117411 + ,0.061296189365771 + ,-0.0463306298649301 + ,0.0380278085379684 + ,-0.00825029078426785 + ,0.0306190401585776 + ,-0.04367173164859 + ,-0.0556117529677963 + ,-0.0645860417354424 + ,-0.111874591444632 + ,0.0165546744538422 + ,0.0221480269493989 + ,-0.0314854358112557 + ,0.000731004014562505 + ,0.038079969509555 + ,0.0370058066179545 + ,0.0559163673005735 + ,0.00562659064346704 + ,0.0194175429404750 + ,-0.00318505854805297) > 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] 59 > (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 5 5 5 5 5 5 5 5 5 5 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] 0.0007460996 -0.005065864 0.016976977 0.025309869 -0.111874591 [2,] 0.0292111462 -0.011648466 0.014158514 -0.029513000 0.016554674 [3,] -0.0028026760 -0.002661129 0.005208404 -0.007502208 0.022148027 [4,] -0.0018436856 0.017941222 -0.022506509 0.005848574 -0.031485436 [5,] 0.0107457258 -0.016812081 0.022125694 0.061296189 0.000731004 [6,] -0.0175566443 -0.014769707 0.008638135 -0.046330630 0.038079970 [7,] 0.0131165820 0.035665173 0.025985342 0.038027809 0.037005807 [8,] -0.0324353169 0.002938683 -0.014962664 -0.008250291 0.055916367 [9,] 0.0224382021 -0.014583544 0.033677338 0.030619040 0.005626591 [10,] -0.0196217960 0.018487197 -0.008494140 -0.043671732 0.019417543 [11,] 0.0049549241 -0.006724240 -0.022988277 -0.055611753 -0.003185059 [12,] 0.0016857445 -0.008344648 0.056900071 -0.064586042 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.0284650466 -0.006582602 -0.002818464 -0.054822869 0.128429266 [2,] -0.0320138221 0.008987337 -0.008950110 0.022010792 0.005593352 [3,] 0.0009589903 0.020602351 -0.027714913 0.013350782 -0.053633463 [4,] 0.0125894114 -0.034753302 0.044632202 0.055447615 0.032216440 [5,] -0.0283023701 0.002042374 -0.013487558 -0.107626819 0.037348965 [6,] 0.0306732264 0.050434880 0.017347207 0.084358438 -0.001074163 [7,] -0.0455518989 -0.032726489 -0.040948006 -0.046278099 0.018910561 [8,] 0.0548735190 -0.017522227 0.048640002 0.038869331 -0.050289777 [9,] -0.0420599981 0.033070741 -0.042171478 -0.074290772 0.013790952 [10,] 0.0245767201 -0.025211437 -0.014494137 -0.011940021 -0.022602601 [11,] -0.0032691796 -0.001620408 0.079888349 -0.008974289 NA [12,] -0.0067516086 0.025321625 -0.031590202 -0.047288550 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/17e2s1259953031.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/2sts81259953031.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/3s4zi1259953031.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/4egc61259953031.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] [1,] -0.0050658641 -0.02951300 -0.007502208 -0.031485436 -0.016812081 [2,] -0.0050658641 -0.01164847 -0.002802676 -0.022506509 0.000731004 [3,] 0.0007460996 0.01415851 -0.002661129 -0.001843686 0.010745726 [4,] 0.0169769774 0.01655467 0.005208404 0.005848574 0.022125694 [5,] 0.0253098692 0.02921115 0.005208404 0.017941222 0.022125694 [,6] [,7] [,8] [,9] [,10] [,11] [1,] -0.046330630 0.01311658 -0.032435317 -0.014583544 -0.04367173 -0.022988277 [2,] -0.017556644 0.02598534 -0.014962664 0.005626591 -0.01962180 -0.022988277 [3,] -0.014769707 0.03566517 -0.008250291 0.022438202 -0.00849414 -0.006724240 [4,] 0.008638135 0.03700581 0.002938683 0.030619040 0.01848720 -0.003185059 [5,] 0.038079970 0.03802781 0.002938683 0.033677338 0.01941754 0.004954924 [,12] [1,] -0.064586042 [2,] -0.036465345 [3,] -0.003329452 [4,] 0.029292908 [5,] 0.056900071 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.01482932 -0.005769754 -0.008321738 -0.02187932 -0.004371708 [2,] 0.01632152 0.034086782 0.002999480 0.01819194 0.025863160 [,6] [,7] [,8] [,9] [,10] [,11] [1,] -0.033278872 0.02787814 -0.020899338 0.0047786 -0.03542186 -0.020717144 [2,] 0.003739458 0.04345221 0.004398756 0.0400978 0.01843358 0.007268665 [,12] [1,] -0.05527847 [2,] 0.04861957 $out [1] -0.11187459 0.02214803 0.06129619 0.05591637 -0.05561175 $group [1] 1 3 5 8 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.00506586405101019, -0.00506586405101019, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5lbcc1259953031.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] [1,] -0.054822869 -0.032013822 -0.0536334628 -0.03475330 -0.028302370 [2,] -0.006582602 -0.008950110 -0.0277149128 0.01258941 -0.028302370 [3,] -0.002818464 0.005593352 0.0009589903 0.03221644 -0.013487558 [4,] 0.028465047 0.008987337 0.0133507818 0.04463220 0.002042374 [5,] 0.028465047 0.022010792 0.0206023506 0.05544762 0.037348965 [,6] [,7] [,8] [,9] [,10] [,11] [1,] -0.001074163 -0.04627810 -0.05028978 -0.07429077 -0.02521144 -0.008974289 [2,] 0.017347207 -0.04555190 -0.01752223 -0.04217148 -0.02260260 -0.006121734 [3,] 0.030673226 -0.04094801 0.03886933 -0.04206000 -0.01449414 -0.002444794 [4,] 0.050434880 -0.03272649 0.04864000 0.01379095 -0.01194002 0.039133970 [5,] 0.084358438 -0.03272649 0.05487352 0.03307074 -0.01194002 0.079888349 [,12] [1,] -0.047288550 [2,] -0.039439376 [3,] -0.019170905 [4,] 0.009285008 [5,] 0.025321625 $n [1] 5 5 5 5 5 5 5 5 5 5 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.02758304 -0.007081202 -0.02805793 0.009575084 -0.034929078 0.00729356 [2,] 0.02194612 0.018267907 0.02997591 0.054857795 0.007953961 0.05405289 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.05001041 -0.007880734 -0.081602910 -0.022028290 -0.03819680 -0.05766317 [2,] -0.03188560 0.085619396 -0.002517086 -0.006959985 0.03330721 0.01932136 $out [1] 0.12842927 -0.10762682 0.01891056 0.02457672 $group [1] 1 5 7 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.054822868781857, -0.00658260175105271, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6vo261259953031.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.032435317 -0.016812081 -0.02298828 -0.064586042 -0.031485436 [2,] -0.010179660 -0.013116005 -0.01172840 -0.045001181 -0.001227027 [3,] 0.001215922 -0.005895052 0.01139832 -0.007876249 0.016554674 [4,] 0.011931154 0.010439953 0.02405552 0.027964455 0.029576917 [5,] 0.029211146 0.035665173 0.05690007 0.061296189 0.055916367 $n [1] 12 12 12 12 11 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.008868969 -0.016639082 -0.004922962 -0.04115637 0.001880047 [2,] 0.011300813 0.004848978 0.027719611 0.02540387 0.031229301 $out [1] -0.1118746 $group [1] 5 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.0324353168918367, -0.0101796601430888, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7z6ob1259953031.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.016710881 -0.014769707 -0.0130866680 [2,] -0.006592876 -0.007487265 -0.0052356225 [3,] -0.001472431 -0.002252407 0.0003177824 [4,] 0.009654050 0.012452120 0.0086919528 [5,] 0.029960142 0.035665173 0.0181228154 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.008882763 -0.01134689 -0.006034678 [2,] 0.005937900 0.00684208 0.006670243 $out [1] 0.03149557 $group [1] 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.0167108808813666, -0.00659287612309633, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/17e2s1259953031.ps tmp/17e2s1259953031.png") > system("convert tmp/2sts81259953031.ps tmp/2sts81259953031.png") > system("convert tmp/3s4zi1259953031.ps tmp/3s4zi1259953031.png") > system("convert tmp/4egc61259953031.ps tmp/4egc61259953031.png") > system("convert tmp/5lbcc1259953031.ps tmp/5lbcc1259953031.png") > system("convert tmp/6vo261259953031.ps tmp/6vo261259953031.png") > system("convert tmp/7z6ob1259953031.ps tmp/7z6ob1259953031.png") > > > proc.time() user system elapsed 1.213 0.966 1.483