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.048509187,0.004875268,0.040233239,-0.019762292,0.003750254,-0.030054106,-0.028300605,-0.082704171,-0.11302669,-0.109498238,-0.101817328,-0.118607274,-0.117509531,-0.14465196,-0.134191387,-0.112640144,-0.136506071,-0.106713217,-0.069545605,-0.035113962,-0.049813456,-0.046036192,-0.039874149,-0.052636033,-0.051174681,-0.005255431,-0.006168594,-0.026415186,8.8102E-11,0.011735286,0.041086939,0.020957834,0.022425278,0.047397174,0.044352605,0.094029705,0.121239498,0.151691,0.151877967,0.161512219,0.195347326,0.256849467,0.236134898,0.229416189,0.256698716,0.236711619,0.15261767,0.114467851,-0.043411072,-0.031464005,0.08604682,-0.033480636,-0.049140091,0.001907932,0.008216556,0.078480333,0.088515538,0.090288271,0.105740055) > 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] 59 > (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] 5 5 5 5 5 5 5 5 5 5 5 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] -0.048509187 -0.11750953 -5.117468e-02 0.1212395 -0.043411072 [2,] 0.004875268 -0.14465196 -5.255431e-03 0.1516910 -0.031464005 [3,] 0.040233239 -0.13419139 -6.168594e-03 0.1518780 0.086046820 [4,] -0.019762292 -0.11264014 -2.641519e-02 0.1615122 -0.033480636 [5,] 0.003750254 -0.13650607 8.810200e-11 0.1953473 -0.049140091 [6,] -0.030054106 -0.10671322 1.173529e-02 0.2568495 0.001907932 [7,] -0.028300605 -0.06954560 4.108694e-02 0.2361349 0.008216556 [8,] -0.082704171 -0.03511396 2.095783e-02 0.2294162 0.078480333 [9,] -0.113026690 -0.04981346 2.242528e-02 0.2566987 0.088515538 [10,] -0.109498238 -0.04603619 4.739717e-02 0.2367116 0.090288271 [11,] -0.101817328 -0.03987415 4.435261e-02 0.1526177 0.105740055 [12,] -0.118607274 -0.05263603 9.402971e-02 0.1144679 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.053384455 -0.027142429 0.045919250 0.030451502 0.011947067 [2,] 0.035357971 0.010460573 -0.000913163 0.000186967 0.117510825 [3,] -0.059995531 0.021551243 -0.020246592 0.009634252 -0.119527456 [4,] 0.023512546 -0.023865927 0.026415186 0.033835107 -0.015659455 [5,] -0.033804360 0.029792854 0.011735286 0.061502141 0.051048023 [6,] 0.001753501 0.037167612 0.029351653 -0.020714569 0.006308624 [7,] -0.054403566 0.034431643 -0.020129105 -0.006718709 0.070263777 [8,] -0.030322519 -0.014699494 0.001467444 0.027282527 0.010035205 [9,] 0.003528452 0.003777264 0.024971896 -0.019987097 0.001772733 [10,] 0.007680910 0.006162043 -0.003044569 -0.084093949 0.015451784 [11,] -0.016789946 -0.012761884 0.049677100 -0.038149819 NA [12,] 0.001097743 0.001461352 0.027209793 -0.157878923 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/16nj61257490540.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/22a501257490540.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/3c5a21257490540.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/4itfz1257490540.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.05117468 -0.031464005 -0.134191387 -0.03348064 -4.914009e-02 [2,] -0.05117468 -0.031464005 -0.006168594 -0.03348064 -4.914009e-02 [3,] -0.04850919 -0.005255431 0.040233239 -0.02641519 8.810200e-11 [4,] -0.04341107 0.004875268 0.086046820 -0.01976229 3.750254e-03 [5,] -0.04341107 0.004875268 0.151877967 -0.01976229 3.750254e-03 [,6] [,7] [,8] [,9] [,10] [,11] [1,] -0.030054106 -0.069545605 -0.08270417 -0.11302669 -0.10949824 -0.10181733 [2,] -0.030054106 -0.028300605 -0.03511396 -0.04981346 -0.04603619 -0.03987415 [3,] 0.001907932 0.008216556 0.02095783 0.02242528 0.04739717 0.04435261 [4,] 0.011735286 0.041086939 0.07848033 0.08851554 0.09028827 0.10574005 [5,] 0.011735286 0.041086939 0.22941619 0.25669872 0.23671162 0.15261767 [,12] [1,] -0.11860727 [2,] -0.08562165 [3,] 0.02069684 [4,] 0.10424878 [5,] 0.11446785 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.05399493 -0.03093267 -0.02492594 -0.03610853 -0.03737218 -0.02762035 [2,] -0.04302344 0.02042181 0.10539242 -0.01672184 0.03737218 0.03143621 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.04081251 -0.05930761 -0.07531764 -0.04892935 -0.05853802 -0.1293008 [2,] 0.05724562 0.10122328 0.12016820 0.14372370 0.14724323 0.1706945 $out [1] -0.1175095 0.1212395 -0.1446520 0.1516910 -0.1126401 0.1615122 [7] -0.1365061 0.1953473 -0.1067132 0.2568495 0.2361349 $group [1] 1 1 2 2 4 4 5 5 6 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.051174681, -0.051174681, -0.048509187, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5827w1257490540.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.02714243 -0.000913163 -0.119527456 -0.02386593 -0.03380436 -0.020714569 [2,] 0.01194707 0.000186967 -0.059995531 -0.01565945 0.01173529 0.001753501 [3,] 0.03045150 0.010460573 -0.020246592 0.02351255 0.02979285 0.006308624 [4,] 0.04591925 0.035357971 0.009634252 0.02641519 0.05104802 0.029351653 [5,] 0.05338446 0.035357971 0.021551243 0.03383511 0.06150214 0.037167612 [,7] [,8] [,9] [,10] [,11] [1,] -0.054403566 -0.030322519 0.001772733 -0.003044569 -0.03814982 [2,] -0.020129105 -0.014699494 0.001772733 -0.003044569 -0.02746988 [3,] -0.006718709 0.001467444 0.003528452 0.006162043 -0.01477592 [4,] 0.034431643 0.010035205 0.003777264 0.007680910 0.01845761 [5,] 0.070263777 0.027282527 0.003777264 0.015451784 0.04967710 [,12] [1,] -0.157878923 [2,] -0.078390590 [3,] 0.001279547 [4,] 0.014335572 [5,] 0.027209793 $n [1] 5 5 5 5 5 5 5 5 5 5 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.006446843 -0.01439117 -0.06944682 -0.00621729 0.002014573 -0.01319216 [2,] 0.054456161 0.03531232 0.02895364 0.05324238 0.057571135 0.02580941 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.04527120 -0.01601003 0.002112055 -0.001416553 -0.05105863 -0.07197412 [2,] 0.03183378 0.01894492 0.004944849 0.013740639 0.02150680 0.07453322 $out [1] 0.11751082 0.02497190 -0.01998710 -0.08409395 $group [1] 2 9 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.027142429, 0.011947067, 0.030451502, 0.04591925, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6qnhl1257490540.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.118607274 -0.14465196 -0.051174681 0.1144679 -0.049140091 [2,] -0.105657783 -0.12585046 -0.005712013 0.1517845 -0.032472320 [3,] -0.039281647 -0.08812941 0.016346560 0.1784298 0.008216556 [4,] -0.008006019 -0.04792482 0.042719772 0.2364233 0.087281179 [5,] 0.040233239 -0.03511396 0.094029705 0.2568495 0.105740055 $n [1] 12 12 12 12 11 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.083821272 -0.12367182 -0.005743502 0.1398255 -0.04883257 [2,] 0.005257979 -0.05258701 0.038436622 0.2170341 0.06526568 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.118607274, -0.105657783, -0.0392816465, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/75x1w1257490540.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.027872995 -0.048509187 -0.04729288 [2,] -0.001135371 -0.002627715 -0.01799464 [3,] 0.027152341 0.014456696 0.01069918 [4,] 0.039239157 0.031329258 0.02190461 [5,] 0.043772527 0.047397174 0.03993911 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 0.008737248 -0.001031308 -0.007499135 [2,] 0.045567434 0.029944700 0.028897502 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.0278729946, -0.00113537099118980, 0.0271523407, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/16nj61257490540.ps tmp/16nj61257490540.png") > system("convert tmp/22a501257490540.ps tmp/22a501257490540.png") > system("convert tmp/3c5a21257490540.ps tmp/3c5a21257490540.png") > system("convert tmp/4itfz1257490540.ps tmp/4itfz1257490540.png") > system("convert tmp/5827w1257490540.ps tmp/5827w1257490540.png") > system("convert tmp/6qnhl1257490540.ps tmp/6qnhl1257490540.png") > system("convert tmp/75x1w1257490540.ps tmp/75x1w1257490540.png") > > > proc.time() user system elapsed 1.187 0.951 1.436