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.0173048550615281,0.00970821286152053,-0.00079996540379326,-0.00892066473015527,-0.0141008132138289,0.00759339542659559,-0.00535785170826649,0.0122184275141973,0.0162622652747317,-0.00273973183483127,-0.0246245251960173,-0.00939313017474595,-0.00118110304103207,-0.0310066632805764,-0.0122299234249117,-0.0131151854993210,0.0547412915989603,-0.0262297232767024,-0.0177472894511915,0.00677361015670445,-0.0342326146132618,-0.0263917601072953,0.0345414780136525,0.00934818751456624,-0.04291732548152,0.0229871189215208,0.00583061899185308,0.033403678623598,0.00323736391037539,-0.00232346094165913,-0.0128854409654862,0.0132916562626398,-0.033685507227119,0.043812485956927,0.00677871591310374,-0.0269603234365068,0.00612710347230793,-0.00195884473559086,0.0366963052514831,0.0344551911103022,0.0386513683774856,0.0181809583929349,0.0375555506713487,0.00753090985514504,-0.0174691093854333,0.00463114466519787,-0.0278211608566847,0.0082690870152522) > 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.0173048551 -0.001181103 -0.042917325 0.006127103 NA [2,] 0.0097082129 -0.031006663 0.022987119 -0.001958845 NA [3,] -0.0007999654 -0.012229923 0.005830619 0.036696305 NA [4,] -0.0089206647 -0.013115185 0.033403679 0.034455191 NA [5,] -0.0141008132 0.054741292 0.003237364 0.038651368 NA [6,] 0.0075933954 -0.026229723 -0.002323461 0.018180958 NA [7,] -0.0053578517 -0.017747289 -0.012885441 0.037555551 NA [8,] 0.0122184275 0.006773610 0.013291656 0.007530910 NA [9,] 0.0162622653 -0.034232615 -0.033685507 -0.017469109 NA [10,] -0.0027397318 -0.026391760 0.043812486 0.004631145 NA [11,] -0.0246245252 0.034541478 0.006778716 -0.027821161 NA [12,] -0.0093931302 0.009348188 -0.026960323 0.008269087 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.027013068 -0.029825560 0.065904444 -0.008085948 NA [2,] -0.010508178 0.018776740 -0.017156500 0.038655150 NA [3,] -0.008120699 -0.000885262 0.027573060 -0.002241114 NA [4,] -0.005180148 0.067856477 -0.030166315 0.004196177 NA [5,] 0.021694209 -0.080971015 -0.005560825 -0.020470410 NA [6,] -0.012951247 0.008482434 -0.010561980 0.019374592 NA [7,] 0.017576279 0.024520900 0.026177097 -0.030024641 NA [8,] 0.004043838 -0.041006225 -0.046977163 -0.025000019 NA [9,] -0.019001997 0.007840855 0.077497993 0.022100254 NA [10,] -0.021884793 0.060933238 -0.037033770 -0.032452306 NA [11,] 0.015231395 -0.025193290 -0.033739039 0.036090248 NA [12,] 0.008212027 -0.052265513 0.033087427 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/1p9la1260466695.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/2x98d1260466696.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/31q8w1260466696.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/4ija31260466696.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.042917325 -0.031006663 -0.012229923 -0.01311519 -0.014100813 [2,] -0.030111090 -0.016482754 -0.006514944 -0.01101793 -0.005431725 [3,] -0.009242979 0.003874684 0.002515327 0.01224151 0.020944366 [4,] 0.002473000 0.016347666 0.021263462 0.03392943 0.046696330 [5,] 0.006127103 0.022987119 0.036696305 0.03445519 0.054741292 [,6] [,7] [,8] [,9] [,10] [1,] -0.026229723 -0.017747289 0.006773610 -0.034232615 -0.0263917601 [2,] -0.014276592 -0.015316365 0.007152260 -0.033959061 -0.0145657460 [3,] 0.002634967 -0.009121646 0.009874669 -0.025577308 0.0009457064 [4,] 0.012887177 0.016098849 0.012755042 -0.000603422 0.0242218153 [5,] 0.018180958 0.037555551 0.013291656 0.016262265 0.0438124860 [,11] [,12] [1,] -0.027821161 -0.0269603234 [2,] -0.026222843 -0.0181767268 [3,] -0.008922905 -0.0005620216 [4,] 0.020660097 0.0088086373 [5,] 0.034541478 0.0093481875 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.03498441 -0.02206135 -0.01942961 -0.02326691 -0.02023680 -0.01882441 [2,] 0.01649845 0.02981072 0.02446027 0.04774992 0.06212553 0.02409434 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.03393967 0.005448471 -0.0519282630 -0.02969647 -0.04596043 -0.02188046 [2,] 0.01569637 0.014300866 0.0007736464 0.03158788 0.02811462 0.02075642 $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.04291732548152, -0.0301110902715240, -0.00924297905128008, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5fuce1260466696.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.02982556 -0.017156500 -0.008120699 -0.0301663147 -0.080971015 [2,] -0.01895575 -0.013832339 -0.005180907 -0.0176732316 -0.050720712 [3,] 0.00946356 0.004134281 -0.001563188 -0.0004919856 -0.013015617 [4,] 0.04645876 0.028715945 0.013343899 0.0360263272 0.008066692 [5,] 0.06590444 0.038655150 0.027573060 0.0678564771 0.021694209 [,6] [,7] [,8] [,9] [,10] [1,] -0.012951247 -0.030024641 -0.046977163 -0.019001997 -0.03703377 [2,] -0.011756614 -0.006224181 -0.043991694 -0.005580571 -0.03474304 [3,] -0.001039773 0.021048589 -0.033003122 0.014970554 -0.02716855 [4,] 0.013928513 0.025348998 -0.010478091 0.049799124 0.01952422 [5,] 0.019374592 0.026177097 0.004043838 0.077497993 0.06093324 [,11] [,12] [1,] -0.033739039 -0.052265513 [2,] -0.029466165 -0.022026743 [3,] -0.004980948 0.008212027 [4,] 0.025660821 0.020649727 [5,] 0.036090248 0.033087427 $n [1] 4 4 4 4 4 4 4 4 4 4 4 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.04221390 -0.02947886 -0.01619778 -0.04291464 -0.05945767 -0.02133102 [2,] 0.06114102 0.03774743 0.01307141 0.04193067 0.03342643 0.01925148 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.003894222 -0.059478869 -0.02877940 -0.07003968 -0.04853127 -0.03071802 [2,] 0.045991401 -0.006527375 0.05872051 0.01570259 0.03856937 0.04714208 $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.0298255602395443, -0.0189557542237216, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6hrea1260466696.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.024624525 -0.034232615 -0.042917325 -0.027821161 NA [2,] -0.011746972 -0.026310742 -0.019922882 0.001336150 NA [3,] -0.004048792 -0.012672554 0.004533991 0.007899998 NA [4,] 0.008650804 0.008060899 0.018139388 0.035575748 NA [5,] 0.016262265 0.054741292 0.043812486 0.038651368 NA $n [1] 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.013352354 -0.028349690 -0.01282647 -0.007716912 NA [2,] 0.005254771 0.003004582 0.02189445 0.023516909 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.0246245251960173, -0.0117469716942874, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7r5yl1260466696.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.0172812415 -0.025577308 -0.021429275 [2,] -0.0037327089 -0.009022275 -0.005108670 [3,] 0.0001618490 0.001730517 0.001436850 [4,] 0.0086639549 0.006874676 0.007429476 [5,] 0.0206323027 0.020944366 0.020788334 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.005492353 -0.005520190 -0.004281883 [2,] 0.005816051 0.008981223 0.007155583 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.0172812414877706, -0.00373270890092251, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1p9la1260466695.ps tmp/1p9la1260466695.png") > system("convert tmp/2x98d1260466696.ps tmp/2x98d1260466696.png") > system("convert tmp/31q8w1260466696.ps tmp/31q8w1260466696.png") > system("convert tmp/4ija31260466696.ps tmp/4ija31260466696.png") > system("convert tmp/5fuce1260466696.ps tmp/5fuce1260466696.png") > system("convert tmp/6hrea1260466696.ps tmp/6hrea1260466696.png") > system("convert tmp/7r5yl1260466696.ps tmp/7r5yl1260466696.png") > > > proc.time() user system elapsed 1.212 0.960 1.486