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.182516799296172,-0.94944747492426,0.325776406812306,1.34843784634604,1.75851178494723,-4.42943737703812,-2.31767034001953,-2.79684843911972,1.12860275120424,0.217683629453875,-0.216346021448786,0.255354934382147,2.08016887397565,1.98907623974419,-3.52512220833211,1.92157710024745,-1.23357233233317,-0.578487095948195,5.06860933422247,-4.07661330236988,-1.08022426396436,-2.92628938204488,1.85793956152257,-0.724378886250608,-1.18764122154349,-4.95476521135668,-2.32628727835623,-0.861616376891362,0.240169843746047,-1.29423362584277,-0.405011745027942,0.234630784545567,2.02245409411938,-1.53833721759871,-1.22520941203817,-1.34782503211567,-1.09231372313534,-2.23500733525864,2.01390498049689,-0.25916025754455,-0.251975811363012,3.19306194411413,-0.520346345270884,0.263552410490847,3.6273589881661,0.757824778633043,-2.20481542891670,-0.350907792837528) > 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.1825168 2.0801689 -1.1876412 -1.0923137 NA [2,] -0.9494475 1.9890762 -4.9547652 -2.2350073 NA [3,] 0.3257764 -3.5251222 -2.3262873 2.0139050 NA [4,] 1.3484378 1.9215771 -0.8616164 -0.2591603 NA [5,] 1.7585118 -1.2335723 0.2401698 -0.2519758 NA [6,] -4.4294374 -0.5784871 -1.2942336 3.1930619 NA [7,] -2.3176703 5.0686093 -0.4050117 -0.5203463 NA [8,] -2.7968484 -4.0766133 0.2346308 0.2635524 NA [9,] 1.1286028 -1.0802243 2.0224541 3.6273590 NA [10,] 0.2176836 -2.9262894 -1.5383372 0.7578248 NA [11,] -0.2163460 1.8579396 -1.2252094 -2.2048154 NA [12,] 0.2553549 -0.7243789 -1.3478250 -0.3509078 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] -0.7669307 -0.09109263 -3.7671240 -1.142693612 NA [2,] 1.2752239 -5.51419845 2.6284779 4.248912316 NA [3,] 1.0226614 5.44669931 1.4646709 -2.273065238 NA [4,] 0.4100739 -3.15514943 1.1017862 0.007184446 NA [5,] -6.1879492 0.65508524 -1.5344035 3.445037755 NA [6,] 2.1117670 5.64709643 0.8892219 -3.713408289 NA [7,] -0.4791781 -9.14522264 0.6396425 0.783898756 NA [8,] 3.9254512 2.99638904 1.7878233 3.363806578 NA [9,] -0.9109191 -1.84606512 -3.5607913 -2.869534210 NA [10,] -0.4340297 4.78422894 0.3131278 -2.962640208 NA [11,] 0.4717010 -2.58231845 -0.1226156 1.853907636 NA [12,] 1.8248139 -0.46326234 0.2555113 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/1m4u81259933628.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/21koz1259933628.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/31dbg1259933628.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/4nmiy1259933628.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,] -1.1876412 -4.9547652 -3.525122 -0.8616164 -1.233572332 -4.4294374 [2,] -1.1399775 -3.5948863 -2.925705 -0.5603883 -0.742774072 -2.8618355 [3,] -0.6374153 -1.5922274 -1.000255 0.5446388 -0.005902984 -0.9363604 [4,] 0.9488260 0.5198144 1.169841 1.6350075 0.999340814 1.3072874 [5,] 2.0801689 1.9890762 2.013905 1.9215771 1.758511785 3.1930619 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -2.317670 -4.0766133 -1.08022426 -2.9262894 -2.2048154 -1.34782503 [2,] -1.419008 -3.4367309 0.02418924 -2.2323133 -1.7150124 -1.03610196 [3,] -0.462679 -1.2811088 1.57552842 -0.6603268 -0.7207777 -0.53764334 [4,] 2.331799 0.2490916 2.82490654 0.4877542 0.8207968 -0.04777643 [5,] 5.068609 0.2635524 3.62735899 0.7578248 1.8579396 0.25535493 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2.287570 -4.842841 -4.235736 -1.189724 -1.382174 -4.229967 -3.425817 [2,] 1.012740 1.658386 2.235225 2.279001 1.370368 2.357247 2.500459 [,8] [,9] [,10] [,11] [,12] [1,] -4.192909 -0.6370382 -2.809180 -2.724067 -1.3184205 [2,] 1.630691 3.7880951 1.488527 1.282512 0.2431338 $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(-1.18764122154349, -1.13997747233941, -0.637415261215756, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5z68q1259933628.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,] -3.76712399 -5.514198 -2.2730652 -3.1551494 -6.1879492 -3.713408 [2,] -2.45490880 -2.119487 -0.6252019 -1.5739825 -3.8611763 -1.412093 [3,] -0.95481214 1.951851 1.2436662 0.2086292 -0.4396591 1.500494 [4,] -0.42901165 3.438695 3.4556851 0.7559301 2.0500615 3.879432 [5,] -0.09109263 4.248912 5.4466993 1.1017862 3.4450378 5.647096 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -9.14522264 1.787823 -3.5607913 -2.96264021 -2.5823184 -0.4632623 [2,] -4.81220037 2.392106 -3.2151628 -1.69833493 -1.3524670 -0.1038755 [3,] 0.08023222 3.180098 -2.3577997 -0.06045092 0.1745427 0.2555113 [4,] 0.71177064 3.644629 -1.3784921 2.54867837 1.1628043 1.0401626 [5,] 0.78389876 3.925451 -0.9109191 4.78422894 1.8539076 1.8248139 $n [1] 4 4 4 4 4 4 4 4 4 4 4 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2.5552709 -2.439113 -1.980235 -1.632002 -5.109537 -2.679810 -4.283705 [2,] 0.6456466 6.342815 4.467567 2.049260 4.230219 5.680799 4.444169 [,8] [,9] [,10] [,11] [,12] [1,] 2.190605 -3.8087695 -3.415591 -1.812522 -0.7880956 [2,] 4.169591 -0.9068299 3.294690 2.161607 1.2991183 $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(-3.76712398981319, -2.45490880096825, -0.954812143875694, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6708a1259933628.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,] -4.42943738 -4.076613 -2.32628728 -2.2350073 NA [2,] -1.63355891 -2.079931 -1.44308112 -0.8063300 NA [3,] 0.01758342 -0.651433 -1.20642532 -0.2555680 NA [4,] 0.72718958 1.955327 -0.08519048 1.3858649 NA [5,] 1.75851178 5.068609 0.24016984 3.6273590 NA $n [1] 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [1,] -1.059170 -2.491941 -1.8257684 -1.2554429 NA [2,] 1.094337 1.189075 -0.5870823 0.7443068 NA $out [1] -4.954765 2.022454 $group [1] 3 3 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-4.42943737703812, -1.63355890747190, 0.0175834150788515, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/712m71259933628.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,] -1.5938196 -1.5922274 -1.56488168 [2,] -0.8751058 -0.9683079 -0.89795547 [3,] -0.4945235 -0.6488710 -0.56186702 [4,] 0.2923393 -0.2342910 0.02902414 [5,] 1.4245479 0.5446388 0.54097419 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -1.02700305 -0.983661 -0.9846686 [2,] 0.03795603 -0.314081 -0.1390654 $out [1] 1.575528 1.500038 $group [1] 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-1.59381963661330, -0.875105786366977, -0.494523509712843, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > system("convert tmp/1m4u81259933628.ps tmp/1m4u81259933628.png") > system("convert tmp/21koz1259933628.ps tmp/21koz1259933628.png") > system("convert tmp/31dbg1259933628.ps tmp/31dbg1259933628.png") > system("convert tmp/4nmiy1259933628.ps tmp/4nmiy1259933628.png") > system("convert tmp/5z68q1259933628.ps tmp/5z68q1259933628.png") > system("convert tmp/6708a1259933628.ps tmp/6708a1259933628.png") > system("convert tmp/712m71259933628.ps tmp/712m71259933628.png") > > > proc.time() user system elapsed 1.172 0.910 1.438