R version 2.6.0 (2007-10-03) Copyright (C) 2007 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.0270241445914126,-4.15985203764392e-05,-4.71564807232549e-05,0.074608257466791,-8.78440892624426e-06,-0.672883349838305,0.0387110679828922,-9.76779782417108e-06,0.599095899504094,0.00467216910723649,-2.38649217857090e-05,-0.315798808127945,-0.00973174119077954,-0.0873269077991939,-0.268181140360116,0.0208091793953472,-0.0423688325156033,0.539373590787594,0.0987223883190398,-0.177238454207196,-0.352702516983866,-0.0414035274018648,-0.083111753344775,0.00536205734659646,-0.0505515109284127,-0.0312774511414949,0.172107145309883,0.0164817741582300,-0.105658742579198,0.292404389402516,0.0410125524051896,-0.0199766112752285,-0.686937519093817,-0.0135344606911849,-0.0760228103898655,-0.300918567991244,-0.0224037203718124,-0.252466202387083,0.165990420346424,-0.0779856129719014,-0.230732662402762,0.486912235946702,-0.194449943015247,-0.265898227594253,0.291780605444579,-0.101259838196459,-0.314743113441016,-0.194402363838198) > 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,] -2.702414e-02 -0.009731741 -0.05055151 -0.02240372 NA [2,] -4.159852e-05 -0.087326908 -0.03127745 -0.25246620 NA [3,] -4.715648e-05 -0.268181140 0.17210715 0.16599042 NA [4,] 7.460826e-02 0.020809179 0.01648177 -0.07798561 NA [5,] -8.784409e-06 -0.042368833 -0.10565874 -0.23073266 NA [6,] -6.728833e-01 0.539373591 0.29240439 0.48691224 NA [7,] 3.871107e-02 0.098722388 0.04101255 -0.19444994 NA [8,] -9.767798e-06 -0.177238454 -0.01997661 -0.26589823 NA [9,] 5.990959e-01 -0.352702517 -0.68693752 0.29178061 NA [10,] 4.672169e-03 -0.041403527 -0.01353446 -0.10125984 NA [11,] -2.386492e-05 -0.083111753 -0.07602281 -0.31474311 NA [12,] -3.157988e-01 0.005362057 -0.30091857 -0.19440236 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 2.698255e-02 -0.07759517 0.01927406 -0.23006248 NA [2,] -5.557960e-06 -0.18085423 0.20338460 0.41845662 NA [3,] 7.465541e-02 0.28899032 -0.15562537 -0.24397603 NA [4,] -7.461704e-02 -0.06317801 -0.12214052 -0.15274705 NA [5,] -6.728746e-01 0.58174242 0.39806313 0.71764490 NA [6,] 7.115944e-01 -0.44065120 -0.25139184 -0.68136218 NA [7,] -3.872084e-02 -0.27596084 -0.06098916 -0.07144828 NA [8,] 5.991057e-01 -0.17546406 -0.66696091 0.55767883 NA [9,] -5.944237e-01 0.31129899 0.67340306 -0.39304044 NA [10,] -4.696034e-03 -0.04170823 -0.06248835 -0.21348328 NA [11,] -3.157749e-01 0.08847381 -0.22489576 0.12034075 NA [12,] 3.060671e-01 -0.05591357 0.27851485 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/1uaqa1197556562.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/2gfav1197556562.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/3zdxx1197556562.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/4vo8z1197556562.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.050551511 -2.524662e-01 -0.26818114 -0.07798561 -2.307327e-01 [2,] -0.038787828 -1.698966e-01 -0.13411415 -0.03075192 -1.681957e-01 [3,] -0.024713932 -5.930218e-02 0.08297163 0.01864548 -7.401379e-02 [4,] -0.016067731 -1.565952e-02 0.16904878 0.04770872 -2.118881e-02 [5,] -0.009731741 -4.159852e-05 0.17210715 0.07460826 -8.784409e-06 [,6] [,7] [,8] [,9] [,10] [1,] -0.6728833 -0.19444994 -2.658982e-01 -0.68693752 -0.101259838 [2,] -0.1902395 -0.07786944 -2.215683e-01 -0.51982002 -0.071331683 [3,] 0.3896583 0.03986181 -9.860753e-02 -0.03046096 -0.027468994 [4,] 0.5131429 0.06986747 -9.993190e-03 0.44543825 -0.004431146 [5,] 0.5393736 0.09872239 -9.767798e-06 0.59909590 0.004672169 [,11] [,12] [1,] -3.147431e-01 -0.315798808 [2,] -1.989274e-01 -0.308358688 [3,] -7.956728e-02 -0.247660466 [4,] -3.802334e-02 -0.094520153 [5,] -2.386492e-05 0.005362057 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.042662809 -0.18114943 -0.1565271 -0.04333843 -0.19014923 -0.1660138 [2,] -0.006765056 0.06254507 0.3224703 0.08062938 0.04212166 0.9453304 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.07685035 -0.26575190 -0.7930150 -0.08032042 -0.20668152 -0.41659291 [2,] 0.15657397 0.06853684 0.7320931 0.02538243 0.04754695 -0.07872802 $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.0505515109284127, -0.0387878277599127, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5omlc1197556562.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.23006248 -0.1808542 -0.24397603 -0.15274705 -0.6728746 -0.6813622 [2,] -0.15382882 -0.0904299 -0.19980070 -0.13744378 -0.1374057 -0.5610067 [3,] -0.02916055 0.1016895 -0.04048498 -0.09837878 0.4899028 -0.3460215 [4,] 0.02312830 0.3109206 0.18182287 -0.06889753 0.6496937 0.2301013 [5,] 0.02698255 0.4184566 0.28899032 -0.06317801 0.7176449 0.7115944 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.27596084 -0.6669609 -0.59442373 -0.213483275 -0.31577494 -0.05591357 [2,] -0.17370456 -0.4212125 -0.49373209 -0.137985812 -0.27033535 0.11130064 [3,] -0.06621872 0.1911074 -0.04087073 -0.052098288 -0.06821097 0.27851485 [4,] -0.04985500 0.5783923 0.49235102 -0.023202130 0.10440728 0.29229096 [5,] -0.03872084 0.5991057 0.67340306 -0.004696034 0.12034075 0.30606707 $n [1] 4 4 4 4 4 4 4 4 4 4 4 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.1689567 -0.2153774 -0.3419676 -0.15253032 -0.1319057 -0.9709968 [2,] 0.1106356 0.4187564 0.2609976 -0.04422724 1.1117113 0.2789538 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.16405988 -0.5985804 -0.8198764 -0.14277740 -0.3642577 0.1134131 [2,] 0.03162243 0.9807951 0.7381349 0.03858082 0.2278357 0.4436166 $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.230062482015271, -0.153828824311842, -0.0291605534107483, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6cg9p1197556562.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,] -2.702414e-02 -0.26818114 -0.10565874 -0.31474311 NA [2,] -1.353565e-02 -0.13228268 -0.09084078 -0.24159943 NA [3,] -1.681636e-05 -0.04188618 -0.02562703 -0.14783110 NA [4,] 2.169162e-02 0.01308562 0.02874716 0.07179335 NA [5,] 3.871107e-02 0.09872239 0.17210715 0.48691224 NA $n [1] 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.01608421 -0.10818964 -0.08017189 -0.290771653 NA [2,] 0.01605058 0.02441728 0.02891783 -0.004890549 NA $out [1] 0.07460826 -0.67288335 0.59909590 -0.31579881 0.53937359 -0.35270252 [7] 0.29240439 -0.68693752 -0.30091857 $group [1] 1 1 1 1 2 2 3 3 3 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.0270241445914126, -0.0135356505360679, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7igwp1197556562.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.201439421 -0.09860753 -0.22454994 [2,] -0.105236510 -0.07679053 -0.09168718 [3,] -0.037536149 -0.02896497 -0.03325056 [4,] 0.002238708 0.02925364 0.01574618 [5,] 0.161451717 0.08297163 0.05021947 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.08655632 -0.07733244 -0.08225164 [2,] 0.01148402 0.01940249 0.01575051 $out [1] 0.3896583 -0.2476605 0.2755550 $group [1] 2 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.201439420652698, -0.105236510347624, -0.0375361485389103, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1uaqa1197556562.ps tmp/1uaqa1197556562.png") > system("convert tmp/2gfav1197556562.ps tmp/2gfav1197556562.png") > system("convert tmp/3zdxx1197556562.ps tmp/3zdxx1197556562.png") > system("convert tmp/4vo8z1197556562.ps tmp/4vo8z1197556562.png") > system("convert tmp/5omlc1197556562.ps tmp/5omlc1197556562.png") > system("convert tmp/6cg9p1197556562.ps tmp/6cg9p1197556562.png") > system("convert tmp/7igwp1197556562.ps tmp/7igwp1197556562.png") > > > proc.time() user system elapsed 1.448 0.944 1.768