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(63.152,60.106,72.616,73.159,68.848,77.056,62.246,60.777,64.513,58.353,56.511,44.554,71.414,65.719,80.997,69.826,65.386,75.589,65.520,59.003,63.961,59.716,57.520,42.886,69.805,64.656,80.353,71.321,76.577,81.580,71.127,63.478,48.152,69.236,57.038,43.621,69.551,72.009,72.140,81.519,73.310,80.406,70.697,59.328,68.281,70.041,51.244,46.538,61.443,62.256,73.117,74.155,65.191,77.889,68.688,59.983,65.470,65.089,54.795,47.123) > 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] 60 > (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] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 63.152 71.414 69.805 69.551 61.443 NA [2,] 60.106 65.719 64.656 72.009 62.256 NA [3,] 72.616 80.997 80.353 72.140 73.117 NA [4,] 73.159 69.826 71.321 81.519 74.155 NA [5,] 68.848 65.386 76.577 73.310 65.191 NA [6,] 77.056 75.589 81.580 80.406 77.889 NA [7,] 62.246 65.520 71.127 70.697 68.688 NA [8,] 60.777 59.003 63.478 59.328 59.983 NA [9,] 64.513 63.961 48.152 68.281 65.470 NA [10,] 58.353 59.716 69.236 70.041 65.089 NA [11,] 56.511 57.520 57.038 51.244 54.795 NA [12,] 44.554 42.886 43.621 46.538 47.123 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -3.046 -5.695 -5.149 2.458 0.813 NA [2,] 12.510 15.278 15.697 0.131 10.861 NA [3,] 0.543 -11.171 -9.032 9.379 1.038 NA [4,] -4.311 -4.440 5.256 -8.209 -8.964 NA [5,] 8.208 10.203 5.003 7.096 12.698 NA [6,] -14.810 -10.069 -10.453 -9.709 -9.201 NA [7,] -1.469 -6.517 -7.649 -11.369 -8.705 NA [8,] 3.736 4.958 -15.326 8.953 5.487 NA [9,] -6.160 -4.245 21.084 1.760 -0.381 NA [10,] -1.842 -2.196 -12.198 -18.797 -10.294 NA [11,] -11.957 -14.634 -13.417 -4.706 -7.672 NA [12,] 26.860 26.919 25.930 14.905 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/1uinc1292673635.ps",horizontal=F,onefile=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/2uinc1292673635.ps",horizontal=F,onefile=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/3nr4f1292673635.ps",horizontal=F,onefile=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/4nr4f1292673635.ps",horizontal=F,onefile=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] [,7] [,8] [,9] [,10] [1,] 61.443 60.106 72.140 69.826 65.191 75.589 62.246 59.003 63.961 58.353 [2,] 63.152 62.256 72.616 71.321 65.386 77.056 65.520 59.328 63.961 59.716 [3,] 69.551 64.656 73.117 73.159 68.848 77.889 68.688 59.983 64.513 65.089 [4,] 69.805 65.719 80.353 74.155 73.310 80.406 70.697 60.777 65.470 69.236 [5,] 71.414 65.719 80.997 74.155 76.577 81.580 71.127 60.777 65.470 70.041 [,11] [,12] [1,] 54.795 42.886 [2,] 54.795 43.621 [3,] 56.511 44.554 [4,] 57.038 46.538 [5,] 57.520 47.123 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 64.85001 62.20905 67.65006 71.1565 63.24892 75.5219 65.02994 58.95914 [2,] 74.25199 67.10295 78.58394 75.1615 74.44708 80.2561 72.34606 61.00686 [,9] [,10] [,11] [,12] [1,] 63.44674 58.36219 54.9261 42.49286 [2,] 65.57926 71.81581 58.0959 46.61514 $out [1] 72.009 81.519 63.478 48.152 68.281 51.244 $group [1] 2 4 8 9 9 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(61.443, 63.152, 69.551, 69.805, 71.414, 60.106, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5nr4f1292673635.ps",horizontal=F,onefile=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] [,7] [,8] [,9] [,10] [1,] -5.695 10.861 -11.171 -8.964 5.003 -10.453 -11.369 3.736 -6.160 -18.797 [2,] -5.149 10.861 -9.032 -8.209 7.096 -10.453 -8.705 3.736 -4.245 -12.198 [3,] -3.046 12.510 0.543 -4.440 8.208 -10.069 -7.649 4.958 -0.381 -10.294 [4,] 0.813 15.278 1.038 -4.311 10.203 -9.709 -6.517 5.487 1.760 -2.196 [5,] 2.458 15.697 9.379 -4.311 12.698 -9.201 -6.517 5.487 1.760 -1.842 [,11] [,12] [1,] -14.634 14.9050 [2,] -13.417 20.4175 [3,] -11.957 26.3950 [4,] -7.672 26.8895 [5,] -4.706 26.9190 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -7.258734 9.388959 -6.572437 -7.194317 6.012602 -10.594709 -9.195035 [2,] 1.166734 15.631041 7.658437 -1.685683 10.403398 -9.543291 -6.102965 [,8] [,9] [,10] [,11] [,12] [1,] 3.720748 -4.624118 -17.361388 -16.016403 21.28212 [2,] 6.195252 3.862118 -3.226612 -7.897597 31.50788 $out [1] 0.131 5.256 -14.810 -1.469 -15.326 8.953 21.084 $group [1] 2 4 6 7 8 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-5.69500000000001, -5.149, -3.046, 0.813000000000002, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6gi401292673635.ps",horizontal=F,onefile=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,] 44.5540 42.8860 43.6210 51.2440 47.1230 NA [2,] 59.2295 59.3595 60.2580 63.8045 60.7130 NA [3,] 62.6990 65.4530 69.5205 70.3690 65.1400 NA [4,] 70.7320 70.6200 73.9490 72.7250 70.9025 NA [5,] 77.0560 80.9970 81.5800 81.5190 77.8890 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 57.45263 60.31701 63.27594 66.3003 60.4925 NA [2,] 67.94537 70.58899 75.76506 74.4377 69.7875 NA $out [1] 46.538 $group [1] 4 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(44.554, 59.2295, 62.699, 70.732, 77.056, 42.886, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7gi401292673635.ps",horizontal=F,onefile=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,] 55.4216 56.5110 55.9165 [2,] 61.2946 62.2480 62.0200 [3,] 66.0111 66.8885 65.5970 [4,] 71.9292 71.3340 71.0430 [5,] 78.5040 77.8890 78.7310 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 61.16059 62.74431 61.48155 [2,] 70.86161 71.03269 69.71245 $out [1] 44.9444 44.5540 45.0795 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(55.4216, 61.2946, 66.0111, 71.9292, 78.504, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1uinc1292673635.ps tmp/1uinc1292673635.png",intern=TRUE)) character(0) > try(system("convert tmp/2uinc1292673635.ps tmp/2uinc1292673635.png",intern=TRUE)) character(0) > try(system("convert tmp/3nr4f1292673635.ps tmp/3nr4f1292673635.png",intern=TRUE)) character(0) > try(system("convert tmp/4nr4f1292673635.ps tmp/4nr4f1292673635.png",intern=TRUE)) character(0) > try(system("convert tmp/5nr4f1292673635.ps tmp/5nr4f1292673635.png",intern=TRUE)) character(0) > try(system("convert tmp/6gi401292673635.ps tmp/6gi401292673635.png",intern=TRUE)) character(0) > try(system("convert tmp/7gi401292673635.ps tmp/7gi401292673635.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.221 0.971 2.889