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.486292969646479 + ,-0.362917336110185 + ,-0.536348990777285 + ,-0.555623800824641 + ,-0.495834828403081 + ,-0.241674178268230 + ,0.0662024989408287 + ,0.0624686744880833 + ,0.131634556956347 + ,-0.081239572203563 + ,-0.542159033439704 + ,-0.917158420735452 + ,-0.640742361222955 + ,-0.247524779302203 + ,0.249546137870057 + ,0.469687002407094 + ,0.683505850858074 + ,0.599266846562666 + ,0.325917110911554 + ,0.128360909067523 + ,0.0428308525653205 + ,0.192221522252290 + ,0.317269218228926 + ,0.0135272792073247 + ,0.00925906096951977 + ,-0.160948944588365 + ,-0.000173407413776173 + ,0.120811485508215 + ,0.140348479268687 + ,0.30644371235379 + ,0.225399015556098 + ,0.152876286808116 + ,0.131795430931023 + ,0.369191902291354 + ,0.588814505036284 + ,0.777802395266439 + ,0.783221145955551 + ,0.586451320545133 + ,0.216295832728025 + ,0.0217444288285005 + ,0.144111863643870 + ,0.0582537457860351 + ,0.111620653582971 + ,0.194803324471820 + ,0.0839897045834115 + ,0.178865938694431 + ,0.221819257268543 + ,0.714491599311593 + ,0.823884191371536 + ,0.309683273192494 + ,-0.197099071242002 + ,-0.573850728645969 + ,-0.513098671100029 + ,-0.409419085955675 + ,-0.386611433873316 + ,-0.186734091800279 + ,-0.390250545036102 + ,-0.659039791034511 + ,-0.585743947094049 + ,-0.588662853049904 + ,-0.489329067427174 + ,-0.124743533736874 + ,0.26777949883498 + ,0.517231612726801 + ,0.0409673057324789 + ,-0.312871040478587 + ,-0.342527845118136 + ,-0.351775103035264) > 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] 68 > (np <- floor(n / par1)) [1] 5 > 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] 6 6 6 6 6 6 6 6 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.48629297 -0.64074236 0.0092590610 0.78322115 0.8238842 -0.48932907 [2,] -0.36291734 -0.24752478 -0.1609489446 0.58645132 0.3096833 -0.12474353 [3,] -0.53634899 0.24954614 -0.0001734074 0.21629583 -0.1970991 0.26777950 [4,] -0.55562380 0.46968700 0.1208114855 0.02174443 -0.5738507 0.51723161 [5,] -0.49583483 0.68350585 0.1403484793 0.14411186 -0.5130987 0.04096731 [6,] -0.24167418 0.59926685 0.3064437124 0.05825375 -0.4094191 -0.31287104 [7,] 0.06620250 0.32591711 0.2253990156 0.11162065 -0.3866114 -0.34252785 [8,] 0.06246867 0.12836091 0.1528762868 0.19480332 -0.1867341 -0.35177510 [9,] 0.13163456 0.04283085 0.1317954309 0.08398970 -0.3902505 NA [10,] -0.08123957 0.19222152 0.3691919023 0.17886594 -0.6590398 NA [11,] -0.54215903 0.31726922 0.5888145050 0.22181926 -0.5857439 NA [12,] -0.91715842 0.01352728 0.7778023953 0.71449160 -0.5886629 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.123375634 0.393217582 -0.170208006 -0.19676983 -0.514200918 [2,] -0.173431655 0.497070917 0.160775537 -0.37015549 -0.506782344 [3,] -0.019274810 0.220140865 0.120984893 -0.19455140 -0.376751657 [4,] 0.059788972 0.213818848 0.019536994 0.12236743 0.060752058 [5,] 0.254160650 -0.084239004 0.166095233 -0.08585812 0.103679585 [6,] 0.307876677 -0.273349736 -0.081044697 0.05336691 0.022807652 [7,] -0.003733824 -0.197556202 -0.072522729 0.08318267 0.199877342 [8,] 0.069165882 -0.085530057 -0.021080856 -0.11081362 -0.203516453 [9,] -0.212874129 0.149390670 0.237396471 0.09487623 -0.268789246 [10,] -0.460919461 0.125047696 0.219622603 0.04295332 0.073295844 [11,] -0.374999387 -0.303741939 0.188987890 0.49267234 -0.002918906 [12,] 0.276416060 -0.004268218 0.005418751 0.10939259 0.099333786 [,6] [1,] 0.364585534 [2,] 0.392523033 [3,] 0.249452114 [4,] -0.476264307 [5,] -0.353838346 [6,] -0.029656805 [7,] -0.009247258 [8,] NA [9,] NA [10,] NA [11,] NA [12,] 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/1dot81260716102.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/2x86m1260716102.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/39kzo1260716102.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/4q3vn1260716102.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.6407424 -0.3629173 -0.5363490 -0.57385073 -0.5130987 -0.40941909 [2,] -0.4893291 -0.2475248 -0.1970991 -0.55562380 -0.4958348 -0.31287104 [3,] -0.2385170 -0.1428462 0.1080612 0.07127796 0.0906579 -0.09171022 [4,] 0.7832211 0.3096833 0.2495461 0.46968700 0.1441119 0.30644371 [5,] 0.8238842 0.5864513 0.2677795 0.51723161 0.6835059 0.59926685 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.38661143 -0.3517751 0.04283085 -0.08123957 -0.5857439 -0.91715842 [2,] -0.34252785 -0.1867341 0.04283085 -0.08123957 -0.5421590 -0.58866285 [3,] 0.08891158 0.0954148 0.08398970 0.17886594 0.2218193 0.01352728 [4,] 0.22539902 0.1528763 0.13163456 0.19222152 0.3172692 0.71449160 [5,] 0.32591711 0.1948033 0.13179543 0.36919190 0.5888145 0.77780240 $n [1] 6 6 6 6 6 6 6 6 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -1.059353 -0.5022634 -0.1800394 -0.5900806 -0.3221284 -0.4911882 [2,] 0.582319 0.2165710 0.3961618 0.7326365 0.5034442 0.3077678 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.2774196 -0.1236449 0.02124123 -0.01436098 -0.3854506 -0.9072784 [2,] 0.4552427 0.3144745 0.14673818 0.37209286 0.8290891 0.9343329 $out [1] -0.3902505 -0.6590398 $group [1] 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.640742361222955, -0.489329067427174, -0.238516954338480, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5sl9h1260716102.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.51420092 -0.506782344 -0.37675166 0.01953699 -0.35383835 -0.273349736 [2,] -0.19676983 -0.370155488 -0.19455140 0.01953699 -0.08585812 -0.081044697 [3,] -0.02341619 -0.006328059 0.05085504 0.06027051 0.00972029 -0.003424576 [4,] 0.36458553 0.392523033 0.22014086 0.12236743 0.16609523 0.053366908 [5,] 0.39321758 0.497070917 0.24945211 0.21381885 0.25416065 0.053366908 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.197556202 -0.20351645 -0.26878925 0.04295332 -0.374999387 -0.004268218 [2,] -0.072522729 -0.11081362 -0.21287413 0.04295332 -0.303741939 0.005418751 [3,] -0.006490541 -0.08553006 0.09487623 0.07329584 -0.002918906 0.099333786 [4,] 0.083182671 -0.02108086 0.14939067 0.12504770 0.188987890 0.109392592 [5,] 0.199877342 0.06916588 0.23739647 0.21962260 0.492672342 0.109392592 $n [1] 6 6 6 6 6 6 6 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.3855085 -0.4982803 -0.2166349 -0.006058441 -0.1527978 -0.09012440 [2,] 0.3386762 0.4856242 0.3183449 0.126599471 0.1722383 0.08327525 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.10692555 -0.14893500 -0.1610992 0.01528816 -0.3510806 0.02586613 [2,] 0.09394447 -0.02212511 0.3508516 0.13130352 0.3452428 0.17280144 $out [1] -0.4762643 0.3078767 -0.4609195 0.2764161 $group [1] 4 6 10 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.514200918179042, -0.196769825410418, -0.0234161860107954, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6q15w1260716102.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] [,6] [1,] -0.917158421 -0.24752478 -0.16094894 0.02174443 -0.6590398 -0.4893291 [2,] -0.539254012 0.02817907 0.06503527 0.09780518 -0.5797973 -0.3471515 [3,] -0.424605153 0.22088383 0.14661238 0.18683463 -0.3998348 -0.2188073 [4,] -0.009385449 0.39780206 0.33781781 0.40413529 -0.1919166 0.1543734 [5,] 0.131634557 0.68350585 0.58881451 0.78322115 0.3096833 0.5172316 $n [1] 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.6662818 0.0522963 0.02219444 0.04711541 -0.5767498 -0.49896628 [2,] -0.1829285 0.3894714 0.27103033 0.32655385 -0.2229198 0.06135171 $out [1] -0.6407424 0.7778024 0.8238842 $group [1] 2 3 5 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(-0.917158420735452, -0.539254012108495, -0.424605152878332, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7na4d1260716102.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,] -2.058584e-16 -0.23851695 -0.1124449076 [2,] -1.312630e-16 -0.03909147 -0.0235166096 [3,] -1.606929e-17 0.08645064 -0.0006225518 [4,] 4.255765e-17 0.10173800 0.0530356329 [5,] 1.915135e-16 0.22181926 0.0872327048 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -9.535005e-17 0.02221737 -0.03553854 [2,] 6.321148e-17 0.15068391 0.03429344 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-2.05858369991633e-16, -1.31263001770104e-16, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1dot81260716102.ps tmp/1dot81260716102.png",intern=TRUE)) character(0) > try(system("convert tmp/2x86m1260716102.ps tmp/2x86m1260716102.png",intern=TRUE)) character(0) > try(system("convert tmp/39kzo1260716102.ps tmp/39kzo1260716102.png",intern=TRUE)) character(0) > try(system("convert tmp/4q3vn1260716102.ps tmp/4q3vn1260716102.png",intern=TRUE)) character(0) > try(system("convert tmp/5sl9h1260716102.ps tmp/5sl9h1260716102.png",intern=TRUE)) character(0) > try(system("convert tmp/6q15w1260716102.ps tmp/6q15w1260716102.png",intern=TRUE)) character(0) > try(system("convert tmp/7na4d1260716102.ps tmp/7na4d1260716102.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.219 0.938 1.460