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(100,96.21064363,96.31280765,107.1793443,114.9066592,92.56060184,114.9995356,107.1236185,117.7765394,107.3650971,106.2970187,114.5072908,98.0031578,103.0649206,100.2879168,104.6066685,111.1544534,104.9874617,109.9284852,111.5352466,132.4974459,100.3436426,123.0983561,114.2379493,104.569518,109.0833101,106.9843039,133.6769759,124.8537197,122.5132349,116.8013374,116.0118882,129.7575926,125.1973623,143.7912139,127.9465032,130.2962757,108.4424631,129.3675118,143.6797622,131.8844618,117.6186496,118.9560695,104.8202842,134.624315,140.401226,143.8005015,153.4317823,153.2924677,127.3149438,153.5525216,136.9276493,131.7730101,144.3391845,107.4208229,113.6249652,124.2221603,102.0618557,96.36853348,111.6838488) > 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,] 100.00000 98.00316 104.5695 130.2963 153.29247 NA [2,] 96.21064 103.06492 109.0833 108.4425 127.31494 NA [3,] 96.31281 100.28792 106.9843 129.3675 153.55252 NA [4,] 107.17934 104.60667 133.6770 143.6798 136.92765 NA [5,] 114.90666 111.15445 124.8537 131.8845 131.77301 NA [6,] 92.56060 104.98746 122.5132 117.6186 144.33918 NA [7,] 114.99954 109.92849 116.8013 118.9561 107.42082 NA [8,] 107.12362 111.53525 116.0119 104.8203 113.62497 NA [9,] 117.77654 132.49745 129.7576 134.6243 124.22216 NA [10,] 107.36510 100.34364 125.1974 140.4012 102.06186 NA [11,] 106.29702 123.09836 143.7912 143.8005 96.36853 NA [12,] 114.50729 114.23795 127.9465 153.4318 111.68385 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -3.7893564 5.061763 4.5137921 -21.8538126 -25.977524 NA [2,] 0.1021640 -2.777004 -2.0990062 20.9250487 26.237578 NA [3,] 10.8665367 4.318752 26.6926720 14.3122504 -16.624872 NA [4,] 7.7273149 6.547785 -8.8232562 -11.7953004 -5.154639 NA [5,] -22.3460574 -6.166992 -2.3404848 -14.2658122 12.566174 NA [6,] 22.4389338 4.941024 -5.7118975 1.3374199 -36.918362 NA [7,] -7.8759171 1.606761 -0.7894492 -14.1357853 6.204142 NA [8,] 10.6529209 20.962199 13.7457044 29.8040308 10.597195 NA [9,] -10.4114423 -32.153803 -4.5602303 5.7769110 -22.160305 NA [10,] -1.0680784 22.754714 18.5938516 3.3992755 -5.693322 NA [11,] 8.2102721 -8.860407 -15.8447107 9.6312808 15.315315 NA [12,] -16.5041330 -9.668431 2.3497725 -0.1393146 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/162281257774792.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/2pmy41257774792.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/3vdrw1257774792.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/413241257774792.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] [1,] 98.00316 96.21064 96.3128 104.6067 111.1545 92.5606 107.4208 104.8203 [2,] 100.00000 103.06492 100.2879 107.1793 114.9067 104.9875 109.9285 107.1236 [3,] 104.56952 108.44246 106.9843 133.6770 124.8537 117.6186 114.9995 111.5352 [4,] 130.29628 109.08331 129.3675 136.9276 131.7730 122.5132 116.8013 113.6250 [5,] 153.29247 109.08331 153.5525 143.6798 131.8845 144.3392 118.9561 116.0119 [,9] [,10] [,11] [,12] [1,] 117.7765 100.3436 96.36853 111.6838 [2,] 124.2222 102.0619 106.29702 114.2379 [3,] 129.7576 107.3651 123.09836 114.5073 [4,] 132.4974 125.1974 143.79121 127.9465 [5,] 134.6243 140.4012 143.80050 127.9465 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 83.16225 104.1899 86.43674 112.6569 112.9360 105.2350 110.1432 106.9414 [2,] 125.97679 112.6950 127.53187 154.6971 136.7714 130.0023 119.8559 116.1291 [,9] [,10] [,11] [,12] [1,] 123.9103 91.0176 96.60505 104.8209 [2,] 135.6049 123.7126 149.59166 124.1937 $out [1] 127.3149 153.4318 $group [1] 2 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(98.0031578, 100, 104.569518, 130.2962757, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5xpi91257774792.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,] -25.977524 -2.7770038 4.318752 -11.795300 -22.346057 -5.711898 [2,] -21.853813 -2.0990062 4.318752 -8.823256 -14.265812 -5.711898 [3,] -3.789356 0.1021640 10.866537 -5.154639 -6.166992 1.337420 [4,] 4.513792 20.9250487 14.312250 6.547785 -2.340485 4.941024 [5,] 5.061763 26.2375778 26.692672 7.727315 12.566174 4.941024 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -14.1357853 10.59720 -32.153803 -5.693322 -15.844711 -16.504133 [2,] -7.8759171 10.65292 -22.160305 -1.068078 -8.860407 -13.086282 [3,] -0.7894492 13.74570 -10.411442 3.399275 8.210272 -4.903873 [4,] 1.6067614 20.96220 -4.560230 18.593852 9.631281 1.105229 [5,] 6.2041423 29.80403 5.776911 22.754714 15.315315 2.349773 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -22.42064 -16.16658 3.805156 -16.01578 -14.593398 -6.189907 -7.489886 [2,] 14.84193 16.37090 17.927918 5.70650 2.259415 8.864747 5.910988 [,8] [,9] [,10] [,11] [,12] [1,] 6.461194 -22.847610 -10.49379 -4.855908 -16.115167 [2,] 21.030215 2.024726 17.29235 21.276452 6.307421 $out [1] -16.62487 22.43893 -36.91836 $group [1] 3 6 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-25.9775239, -21.8538126, -3.78935636999999, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6868t1257774792.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,] 92.5606 98.00316 104.5695 104.8203 96.36853 NA [2,] 98.1564 101.70428 112.5476 118.2874 109.55234 NA [3,] 107.1515 107.45797 123.6835 131.0904 125.76855 NA [4,] 114.7070 112.88660 128.8520 142.0405 140.63342 NA [5,] 117.7765 123.09836 143.7912 153.4318 153.55252 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 99.60265 102.3576 116.2469 120.2564 111.5923 NA [2,] 114.70031 112.5583 131.1200 141.9243 139.9448 NA $out [1] 132.4974 $group [1] 2 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(92.56060184, 98.156403825, 107.1514814, 114.706975, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/704jr1257774792.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,] 108.8233 104.5695 106.0741 [2,] 114.3475 107.9038 113.4973 [3,] 117.2666 114.7534 114.9879 [4,] 123.6380 123.9760 122.6967 [5,] 127.7756 133.6770 128.3598 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 113.0292 107.4227 110.7920 [2,] 121.5041 122.0841 119.1838 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(108.823256246, 114.34754343, 117.266648095, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/162281257774792.ps tmp/162281257774792.png") > system("convert tmp/2pmy41257774792.ps tmp/2pmy41257774792.png") > system("convert tmp/3vdrw1257774792.ps tmp/3vdrw1257774792.png") > system("convert tmp/413241257774792.ps tmp/413241257774792.png") > system("convert tmp/5xpi91257774792.ps tmp/5xpi91257774792.png") > system("convert tmp/6868t1257774792.ps tmp/6868t1257774792.png") > system("convert tmp/704jr1257774792.ps tmp/704jr1257774792.png") > > > proc.time() user system elapsed 1.233 0.977 3.707