R version 2.7.2 (2008-08-25) 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(12055,12113,9617,12646,13581,12162,10970,11880,11888,12927,12300,12093,12381,12197,9455,13168,13428,11981,11885,11692,12234,14341,13131,12421,14286,12865,11160,14316,14389,14014,13419,12770,13316,15333,14243,13824,14963,13203,12199,15509,14200,15170,14058,13786,14148,16542,13588,15582,15803,14131,12923,15612,16034,16037,14038,15331,15038,17402,14993,16044,16930,15921,14417,15961,17852,16484,14216,17430,17840,17629) > 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] 70 > (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 6 6 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 12055 12381 14286 14963 15803 16930 [2,] 12113 12197 12865 13203 14131 15921 [3,] 9617 9455 11160 12199 12923 14417 [4,] 12646 13168 14316 15509 15612 15961 [5,] 13581 13428 14389 14200 16034 17852 [6,] 12162 11981 14014 15170 16037 16484 [7,] 10970 11885 13419 14058 14038 14216 [8,] 11880 11692 12770 13786 15331 17430 [9,] 11888 12234 13316 14148 15038 17840 [10,] 12927 14341 15333 16542 17402 17629 [11,] 12300 13131 14243 13588 14993 NA [12,] 12093 12421 13824 15582 16044 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 58 -184 -1421 -1760 -1672 -1009 [2,] -2496 -2742 -1705 -1004 -1208 -1504 [3,] 3029 3713 3156 3310 2689 1544 [4,] 935 260 73 -1309 422 1891 [5,] -1419 -1447 -375 970 3 -1368 [6,] -1192 -96 -595 -1112 -1999 -2268 [7,] 910 -193 -649 -272 1293 3214 [8,] 8 542 546 362 -293 410 [9,] 1039 2107 2017 2394 2364 -211 [10,] -627 -1210 -1090 -2954 -2409 NA [11,] -207 -710 -419 1994 1051 NA [12,] 288 1865 1139 221 886 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/1sw571225641787.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/2tdzr1225641787.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/3pnke1225641787.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/4auc31225641787.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] [,10] [1,] 12055.0 12113 9455.0 12646.0 13428.0 11981 10970.0 11692 11888 12927.0 [2,] 12381.0 12197 9617.0 13168.0 13581.0 12162 11885.0 11880 12234 14341.0 [3,] 14624.5 13034 11679.5 14912.5 14294.5 14592 13728.5 13278 13732 15937.5 [4,] 15803.0 14131 12923.0 15612.0 16034.0 16037 14058.0 15331 15038 17402.0 [5,] 16930.0 15921 14417.0 15961.0 17852.0 16484 14216.0 17430 17840 17629.0 [,11] [,12] [1,] 12300 12093 [2,] 13131 12421 [3,] 13588 13824 [4,] 14243 15582 [5,] 14993 16044 $n [1] 6 6 6 6 6 6 6 6 6 6 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 12417.2 11786.51 9547.023 13336.04 12712.24 12092.5 12326.84 11051.99 [2,] 16831.8 14281.49 13811.977 16488.96 15876.76 17091.5 15130.16 15504.01 [,9] [,10] [,11] [,12] [1,] 11923.33 13963.06 12802.26 11590.45 [2,] 15540.67 17911.94 14373.74 16057.55 $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(12055, 12381, 14624.5, 15803, 16930, 12113, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5dhew1225641787.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] [,9] [,10] [,11] [,12] [1,] -1760 -2742.0 2689.0 73 -1447.0 -2268 -649.0 -293 -211 -2954 -710 221 [2,] -1672 -2496.0 2689.0 73 -1419.0 -1999 -272.0 8 1039 -2409 -419 288 [3,] -1215 -1604.5 3092.5 341 -871.5 -1152 358.5 386 2062 -1210 -207 886 [4,] -184 -1208.0 3310.0 935 3.0 -595 1293.0 542 2364 -1090 1051 1139 [5,] 58 -1004.0 3713.0 1891 970.0 -96 3214.0 546 2394 -627 1994 1865 $n [1] 6 6 6 6 6 6 6 6 6 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2174.8081 -2435.3016 2691.935 -215.0178 -1788.73593 -2057.6253 -650.9755 [2,] -255.1919 -773.6984 3493.065 897.0178 45.73593 -246.3747 1367.9755 [,8] [,9] [,10] [,11] [,12] [1,] 41.55275 1207.332 -2142.0021 -1245.6983 284.6855 [2,] 730.44725 2916.668 -277.9979 831.6983 1487.3145 $out [1] 1544 -1309 $group [1] 3 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1760, -1672, -1215, -184, 58, -2742, -2496, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6qr6a1225641787.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,] 11880 11692.0 12770.0 12199.0 12923.0 14216 [2,] 11884 11933.0 13090.5 13687.0 14562.0 15921 [3,] 12103 12307.5 13919.0 14174.0 15471.5 16707 [4,] 12473 13149.5 14301.0 15339.5 16035.5 17629 [5,] 12927 14341.0 15333.0 16542.0 17402.0 17852 $n [1] 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 11834.35 11752.65 13366.88 13420.28 14799.43 15853.62 [2,] 12371.65 12862.35 14471.12 14927.72 16143.57 17560.38 $out [1] 9617 13581 10970 9455 11160 $group [1] 1 1 1 2 3 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(11880, 11884, 12103, 12473, 12927, 11692, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7x2js1225641787.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,] 13097.67 11679.50 13131.50 [2,] 13528.00 13433.00 13342.19 [3,] 14035.07 13778.00 13844.25 [4,] 14469.17 14608.25 14372.88 [5,] 15695.67 15937.50 15888.00 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 13605.79 13241.96 13374.15 [2,] 14464.34 14314.04 14314.35 $out [1] 11628.50 11372.38 $group [1] 1 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(13097.6666666667, 13528, 14035.0666666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1sw571225641787.ps tmp/1sw571225641787.png") > system("convert tmp/2tdzr1225641787.ps tmp/2tdzr1225641787.png") > system("convert tmp/3pnke1225641787.ps tmp/3pnke1225641787.png") > system("convert tmp/4auc31225641787.ps tmp/4auc31225641787.png") > system("convert tmp/5dhew1225641787.ps tmp/5dhew1225641787.png") > system("convert tmp/6qr6a1225641787.ps tmp/6qr6a1225641787.png") > system("convert tmp/7x2js1225641787.ps tmp/7x2js1225641787.png") > > > proc.time() user system elapsed 1.449 0.952 1.704