R version 2.7.0 (2008-04-22) 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(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] 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,] 12300 13131 14243 13588 14993 NA [2,] 12093 12421 13824 15582 16044 NA [3,] 12381 14286 14963 15803 16930 NA [4,] 12197 12865 13203 14131 15921 NA [5,] 9455 11160 12199 12923 14417 NA [6,] 13168 14316 15509 15612 15961 NA [7,] 13428 14389 14200 16034 17852 NA [8,] 11981 14014 15170 16037 16484 NA [9,] 11885 13419 14058 14038 14216 NA [10,] 11692 12770 13786 15331 17430 NA [11,] 12234 13316 14148 15038 17840 NA [12,] 14341 15333 16542 17402 17629 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -207 -710 -419 1994 1051 NA [2,] 288 1865 1139 221 886 NA [3,] -184 -1421 -1760 -1672 -1009 NA [4,] -2742 -1705 -1004 -1208 -1504 NA [5,] 3713 3156 3310 2689 1544 NA [6,] 260 73 -1309 422 1891 NA [7,] -1447 -375 970 3 -1368 NA [8,] -96 -595 -1112 -1999 -2268 NA [9,] -193 -649 -272 1293 3214 NA [10,] 542 546 362 -293 410 NA [11,] 2107 2017 2394 2364 -211 NA [12,] -1210 -1090 -2954 -2409 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/14ww91225734981.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/28nwy1225734981.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/3wd751225734981.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/4q0up1225734981.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] [,11] [,12] [1,] 12300 12093 12381 12197 9455 13168 13428 11981 13419 11692 12234 14341 [2,] 13131 12421 14286 12865 11160 14316 14200 14014 13419 12770 13316 15333 [3,] 13588 13824 14963 13203 12199 15509 14389 15170 14038 13786 14148 16542 [4,] 14243 15582 15803 14131 12923 15612 16034 16037 14058 15331 15038 17402 [5,] 14993 16044 16930 15921 14417 15961 17852 16484 14216 17430 15038 17629 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 12802.26 11590.45 13891.09 12308.45 10953.27 14593.25 13093.1 13740.55 [2,] 14373.74 16057.55 16034.91 14097.55 13444.73 16424.75 15684.9 16599.45 [,9] [,10] [,11] [,12] [1,] 13586.48 11976.40 12931.24 15080.05 [2,] 14489.52 15595.60 15364.76 18003.95 $out [1] 11885 17840 $group [1] 9 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(12300, 13131, 13588, 14243, 14993, 12093, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/59g6a1225734981.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,] -710 221 -1760 -1705 2689 73 -1447 -2268 -649 362 2017 -2954.0 [2,] -419 288 -1672 -1705 2689 73 -1368 -1999 -272 362 2017 -2681.5 [3,] -207 886 -1421 -1504 3156 260 -375 -1112 -193 410 2107 -1809.5 [4,] 1051 1139 -1009 -1208 3310 422 3 -595 1293 542 2364 -1150.0 [5,] 1994 1865 -184 -1004 3713 422 970 -96 3214 546 2394 -1090.0 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1245.6983 284.6855 -1889.4741 -1855.179 2717.203 13.39748 -1343.7451 [2,] 831.6983 1487.3145 -952.5259 -1152.821 3594.797 506.60252 593.7451 [,8] [,9] [,10] [,11] [,12] [1,] -2104.0629 -1298.825 282.8125 1861.811 -3019.385 [2,] -119.9371 912.825 537.1875 2352.189 -599.615 $out [1] -2742 1544 -1309 1891 -293 -211 $group [1] 4 5 6 6 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-710, -419, -207, 1051, 1994, 221, 288, 886, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ysuj1225734981.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,] 11692.0 11160.0 12199.0 12923.0 14216.0 NA [2,] 11933.0 12817.5 13805.0 14084.5 15457.0 NA [3,] 12215.5 13367.5 14174.0 15456.5 16264.0 NA [4,] 12774.5 14301.0 15066.5 15918.5 17529.5 NA [5,] 13428.0 15333.0 16542.0 17402.0 17852.0 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 11831.69 12690.87 13598.62 14620 15318.72 NA [2,] 12599.31 14044.13 14749.38 16293 17209.28 NA $out [1] 9455 14341 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(11692, 11933, 12215.5, 12774.5, 13428, 11160, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7hh671225734981.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,] 12030.8 12199.0 12041.50 [2,] 13657.2 13687.0 13712.75 [3,] 14358.5 14093.0 14113.75 [4,] 14892.9 15066.5 15035.00 [5,] 16249.4 16542.0 16367.50 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 13794.89 13463.8 13510.66 [2,] 14922.11 14722.2 14716.84 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(12030.8, 13657.2, 14358.5, 14892.9, 16249.4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/14ww91225734981.ps tmp/14ww91225734981.png") > system("convert tmp/28nwy1225734981.ps tmp/28nwy1225734981.png") > system("convert tmp/3wd751225734981.ps tmp/3wd751225734981.png") > system("convert tmp/4q0up1225734981.ps tmp/4q0up1225734981.png") > system("convert tmp/59g6a1225734981.ps tmp/59g6a1225734981.png") > system("convert tmp/6ysuj1225734981.ps tmp/6ysuj1225734981.png") > system("convert tmp/7hh671225734981.ps tmp/7hh671225734981.png") > > > proc.time() user system elapsed 2.875 1.602 3.115