R version 2.6.2 (2008-02-08) Copyright (C) 2008 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.70291,0.6885,0.67127,0.66502,0.65825,0.65025,0.65779,0.66014,0.64683,0.64587,0.63702,0.62651,0.61834,0.61466,0.61063,0.59802,0.60151,0.62927,0.62304,0.6071,0.60773,0.58933,0.60039,0.61342,0.6348,0.634,0.62915,0.62168,0.61328,0.6089,0.60857,0.62672,0.62291,0.62393,0.61838,0.62012,0.61659,0.6116,0.61573,0.61407,0.62823,0.64405,0.6387,0.63633,0.63059,0.62994,0.63709,0.64217,0.65711,0.66977,0.68255,0.68902,0.71322,0.70224,0.70045,0.69919,0.69693,0.69763,0.69278,0.70196,0.69215,0.6769,0.67124,0.66532,0.67157,0.66428,0.66576,0.66942,0.6813,0.69144,0.69862,0.695,0.69867,0.68968,0.69233,0.68293,0.68399,0.66895,0.68756,0.68527,0.6776,0.68137,0.67933,0.67922,0.68598,0.68297,0.68935,0.69463,0.6833,0.68666,0.68782,0.67669,0.67511,0.67254,0.67397,0.67286,0.66341,0.668,0.68021,0.67934,0.68136,0.67562,0.6744,0.67766,0.68887,0.69614,0.70896,0.72064,0.74725) > 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] 109 > (np <- floor(n / par1)) [1] 9 > 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] 10 9 9 9 9 9 9 9 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 0.70291 0.61834 0.63480 0.61659 0.65711 0.69215 0.69867 0.68598 0.66341 [2,] 0.68850 0.61466 0.63400 0.61160 0.66977 0.67690 0.68968 0.68297 0.66800 [3,] 0.67127 0.61063 0.62915 0.61573 0.68255 0.67124 0.69233 0.68935 0.68021 [4,] 0.66502 0.59802 0.62168 0.61407 0.68902 0.66532 0.68293 0.69463 0.67934 [5,] 0.65825 0.60151 0.61328 0.62823 0.71322 0.67157 0.68399 0.68330 0.68136 [6,] 0.65025 0.62927 0.60890 0.64405 0.70224 0.66428 0.66895 0.68666 0.67562 [7,] 0.65779 0.62304 0.60857 0.63870 0.70045 0.66576 0.68756 0.68782 0.67440 [8,] 0.66014 0.60710 0.62672 0.63633 0.69919 0.66942 0.68527 0.67669 0.67766 [9,] 0.64683 0.60773 0.62291 0.63059 0.69693 0.68130 0.67760 0.67511 0.68887 [10,] 0.64587 0.58933 0.62393 0.62994 0.69763 0.69144 0.68137 0.67254 0.69614 [11,] 0.63702 0.60039 0.61838 0.63709 0.69278 0.69862 0.67933 0.67397 0.70896 [12,] 0.62651 0.61342 0.62012 0.64217 0.70196 0.69500 0.67922 0.67286 0.72064 [,10] [1,] 0.74725 [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -0.01441 -0.00368 -0.00080 -0.00499 0.01266 -0.01525 -0.00899 -0.00301 [2,] -0.01723 -0.00403 -0.00485 0.00413 0.01278 -0.00566 0.00265 0.00638 [3,] -0.00625 -0.01261 -0.00747 -0.00166 0.00647 -0.00592 -0.00940 0.00528 [4,] -0.00677 0.00349 -0.00840 0.01416 0.02420 0.00625 0.00106 -0.01133 [5,] -0.00800 0.02776 -0.00438 0.01582 -0.01098 -0.00729 -0.01504 0.00336 [6,] 0.00754 -0.00623 -0.00033 -0.00535 -0.00179 0.00148 0.01861 0.00116 [7,] 0.00235 -0.01594 0.01815 -0.00237 -0.00126 0.00366 -0.00229 -0.01113 [8,] -0.01331 0.00063 -0.00381 -0.00574 -0.00226 0.01188 -0.00767 -0.00158 [9,] -0.00096 -0.01840 0.00102 -0.00065 0.00070 0.01014 0.00377 -0.00257 [10,] -0.00885 0.01106 -0.00555 0.00715 -0.00485 0.00718 -0.00204 0.00143 [11,] -0.01051 0.01303 0.00174 0.00508 0.00918 -0.00362 -0.00011 -0.00111 [12,] -0.00817 0.02138 -0.00353 0.01494 -0.00981 0.00367 0.00676 -0.00945 [,9] [,10] [1,] 0.00459 NA [2,] 0.01221 NA [3,] -0.00087 NA [4,] 0.00202 NA [5,] -0.00574 NA [6,] -0.00122 NA [7,] 0.00326 NA [8,] 0.01121 NA [9,] 0.00727 NA [10,] 0.01282 NA [11,] 0.01168 NA [12,] 0.02661 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/1yi5f1209329241.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/2ntin1209329241.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/3mlop1209329241.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/4dwh51209329241.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] [,7] [,8] [,9] [1,] 0.616590 0.61160 0.61063 0.59802 0.60151 0.60890 0.60857 0.60710 0.60773 [2,] 0.634800 0.63400 0.62915 0.62168 0.62823 0.64405 0.63870 0.63633 0.63059 [3,] 0.674695 0.66977 0.67127 0.66532 0.67157 0.66428 0.66576 0.66942 0.67511 [4,] 0.698670 0.68297 0.68255 0.68293 0.68330 0.67562 0.68756 0.67766 0.68130 [5,] 0.747250 0.68968 0.69233 0.69463 0.71322 0.70224 0.70045 0.69919 0.69693 [,10] [,11] [,12] [1,] 0.58933 0.60039 0.61342 [2,] 0.62994 0.63702 0.62651 [3,] 0.67254 0.67397 0.67286 [4,] 0.69144 0.69278 0.69500 [5,] 0.69763 0.70896 0.72064 $n [1] 10 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.642783 0.6439791 0.643146 0.6330617 0.6425665 0.6476531 0.6400271 [2,] 0.706607 0.6955609 0.699394 0.6975783 0.7005735 0.6809069 0.6914929 [,8] [,9] [,10] [,11] [,12] [1,] 0.6476529 0.6484027 0.64015 0.6446031 0.6367886 [2,] 0.6911871 0.7018173 0.70493 0.7033369 0.7089314 $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.61659, 0.6348, 0.674695, 0.69867, 0.74725, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5z5171209329241.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] [,7] [,8] [1,] -0.01525 -0.01723 -0.01261 -0.01133 -0.01504 -0.00623 -0.00237 -0.01331 [2,] -0.00899 -0.00485 -0.00747 -0.00677 -0.00800 -0.00179 -0.00237 -0.00574 [3,] -0.00368 0.00265 -0.00592 0.00202 -0.00574 -0.00033 -0.00126 -0.00226 [4,] -0.00080 0.00638 -0.00087 0.00625 0.00336 0.00148 0.00326 0.00063 [5,] 0.00459 0.01278 0.00647 0.02420 0.01582 0.00148 0.00366 0.00063 [,9] [,10] [,11] [,12] [1,] -0.00257 -0.00885 -0.01051 -0.00981 [2,] -0.00096 -0.00485 -0.00111 -0.00817 [3,] 0.00070 0.00143 0.00174 0.00367 [4,] 0.00377 0.00718 0.00918 0.01494 [5,] 0.01014 0.01282 0.01303 0.02661 $n [1] 9 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.0079934 -0.003264467 -0.009396 -0.0048372 -0.0117229333 -0.0020522 [2,] 0.0006334 0.008564467 -0.002444 0.0088772 0.0002429333 0.0013922 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.004225133 -0.005614867 -0.001791133 -0.0049058 -0.0036794 -0.008501267 [2,] 0.001705133 0.001094867 0.003191133 0.0077658 0.0071594 0.015841267 $out [1] 0.01266 0.02776 0.00754 0.01861 -0.01594 0.01815 -0.01113 0.01188 [9] 0.01121 -0.01840 $group [1] 1 5 6 6 7 7 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(-0.0152500000000001, -0.00899000000000005, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6xjng1209329241.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] [,7] [,8] [1,] 0.626510 0.58933 0.608570 0.611600 0.669770 0.664280 0.668950 0.672540 [2,] 0.646350 0.60095 0.615830 0.616160 0.685785 0.667590 0.679275 0.674540 [3,] 0.658020 0.60918 0.622295 0.630265 0.697280 0.674235 0.683460 0.683135 [4,] 0.668145 0.61650 0.627935 0.637895 0.701205 0.691795 0.688620 0.687240 [5,] 0.688500 0.62927 0.634800 0.644050 0.713220 0.698620 0.698670 0.694630 [,9] [,10] [1,] 0.663410 0.74725 [2,] 0.675010 0.74725 [3,] 0.679775 0.74725 [4,] 0.692505 0.74725 [5,] 0.708960 0.74725 $n [1] 12 12 12 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.6480792 0.6020875 0.6167738 0.6203515 0.6902468 0.663195 0.6791977 [2,] 0.6679608 0.6162725 0.6278162 0.6401785 0.7043132 0.685275 0.6877223 [,8] [,9] [,10] [1,] 0.6773424 0.6717954 0.74725 [2,] 0.6889276 0.6877546 0.74725 $out [1] 0.70291 0.65711 0.72064 $group [1] 1 5 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(0.62651, 0.64635, 0.65802, 0.668145, 0.6885, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7nlks1209329241.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.6566700 0.664280 0.6523050 [2,] 0.6588006 0.667590 0.6558975 [3,] 0.6597000 0.671420 0.6591600 [4,] 0.6605906 0.673415 0.6619425 [5,] 0.6607267 0.675110 0.6687088 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 0.6588836 0.6687632 0.6564028 [2,] 0.6605164 0.6740768 0.6619172 $out [1] 0.6717210 0.6635444 $group [1] 1 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(0.65667, 0.658800555555556, 0.6597, 0.660590555555556, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1yi5f1209329241.ps tmp/1yi5f1209329241.png") > system("convert tmp/2ntin1209329241.ps tmp/2ntin1209329241.png") > system("convert tmp/3mlop1209329241.ps tmp/3mlop1209329241.png") > system("convert tmp/4dwh51209329241.ps tmp/4dwh51209329241.png") > system("convert tmp/5z5171209329241.ps tmp/5z5171209329241.png") > system("convert tmp/6xjng1209329241.ps tmp/6xjng1209329241.png") > system("convert tmp/7nlks1209329241.ps tmp/7nlks1209329241.png") > > > proc.time() user system elapsed 2.561 1.585 2.818