R version 2.8.0 (2008-10-20) 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. Natural language support but running in an English locale 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(41086,39690,43129,37863,35953,29133,24693,22205,21725,27192,21790,13253,37702,30364,32609,30212,29965,28352,25814,22414,20506,28806,22228,13971,36845,35338,35022,34777,26887,23970,22780,17351,21382,24561,17409,11514,31514,27071,29462,26105,22397,23843,21705,18089,20764,25316,17704,15548,28029,29383,36438,32034,22679,24319,18004,17537,20366,22782,19169,13807,29743,25591,29096,26482,22405,27044,17970,18730,19684,19785,18479,10698) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 41086 37702 36845 31514 28029 29743 NA [2,] 39690 30364 35338 27071 29383 25591 NA [3,] 43129 32609 35022 29462 36438 29096 NA [4,] 37863 30212 34777 26105 32034 26482 NA [5,] 35953 29965 26887 22397 22679 22405 NA [6,] 29133 28352 23970 23843 24319 27044 NA [7,] 24693 25814 22780 21705 18004 17970 NA [8,] 22205 22414 17351 18089 17537 18730 NA [9,] 21725 20506 21382 20764 20366 19684 NA [10,] 27192 28806 24561 25316 22782 19785 NA [11,] 21790 22228 17409 17704 19169 18479 NA [12,] 13253 13971 11514 15548 13807 10698 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1396 -7338 -1507 -4443 1354 -4152 NA [2,] 3439 2245 -316 2391 7055 3505 NA [3,] -5266 -2397 -245 -3357 -4404 -2614 NA [4,] -1910 -247 -7890 -3708 -9355 -4077 NA [5,] -6820 -1613 -2917 1446 1640 4639 NA [6,] -4440 -2538 -1190 -2138 -6315 -9074 NA [7,] -2488 -3400 -5429 -3616 -467 760 NA [8,] -480 -1908 4031 2675 2829 954 NA [9,] 5467 8300 3179 4552 2416 101 NA [10,] -5402 -6578 -7152 -7612 -3613 -1306 NA [11,] -8537 -8257 -5895 -2156 -5362 -7781 NA [12,] 24449 22874 20000 12481 15936 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/freestat/rcomp/tmp/1wofm1289748013.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/freestat/rcomp/tmp/2oyfp1289748013.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/freestat/rcomp/tmp/3oyfp1289748013.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/freestat/rcomp/tmp/4zpws1289748013.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,] 28029.0 25591.0 29096.0 26105 22397 23843.0 17970.0 17351.0 19684 19785.0 [2,] 29743.0 27071.0 29462.0 26482 22405 23970.0 18004.0 17537.0 20366 22782.0 [3,] 34179.5 29873.5 33815.5 31123 24783 25681.5 22242.5 18409.5 20635 24938.5 [4,] 37702.0 35338.0 36438.0 34777 29965 28352.0 24693.0 22205.0 21382 27192.0 [5,] 41086.0 39690.0 43129.0 37863 35953 29133.0 25814.0 22414.0 21725 28806.0 [,11] [,12] [1,] 17409 10698 [2,] 17704 11514 [3,] 18824 13530 [4,] 21790 13971 [5,] 22228 15548 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 29045.69 24541.02 29315.75 25772.46 19906.56 22854.97 17927.88 15398.49 [2,] 39313.31 35205.98 38315.25 36473.54 29659.44 28508.03 26557.12 21420.51 [,9] [,10] [,11] [,12] [1,] 19979.65 22093.91 16188.40 11945.16 [2,] 21290.35 27783.09 21459.60 15114.84 $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(28029, 29743, 34179.5, 37702, 41086, 25591, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5zpws1289748013.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] [1,] -7338.0 2245 -5266.0 -9355.0 -6820.0 -9074 -5429 -1908.0 101.0 -7612 [2,] -4443.0 2245 -4404.0 -7890.0 -2917.0 -6315 -3616 -480.0 2416.0 -7152 [3,] -2829.5 2915 -2985.5 -3892.5 -83.5 -3489 -2944 1814.5 3865.5 -5990 [4,] -1396.0 3505 -2397.0 -1910.0 1640.0 -2138 -467 2829.0 5467.0 -3613 [5,] 1354.0 3505 -245.0 -247.0 4639.0 -1190 760 4031.0 8300.0 -1306 [,11] [,12] [1,] -8537 12481 [2,] -8257 15936 [3,] -6838 20000 [4,] -5362 22874 [5,] -2156 24449 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -4794.9134 2102.259 -4280.080 -7749.79315 -3022.912 -6183.2999 -4975.2067 [2,] -864.0866 3727.741 -1690.920 -35.20685 2855.912 -794.7001 -912.7933 [,8] [,9] [,10] [,11] [,12] [1,] -319.9119 1897.506 -8272.769 -8705.369 15097.63 [2,] 3948.9119 5833.494 -3707.231 -4970.631 24902.37 $out [1] -316 7055 $group [1] 2 2 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-7338, -4443, -2829.5, -1396, 1354, 2245, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6zpws1289748013.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] [1,] 13253.0 13971 11514.0 15548.0 13807.0 10698.0 NA [2,] 21997.5 22321 19395.5 19426.5 18586.5 18604.5 NA [3,] 28162.5 28579 24265.5 23120.0 22730.5 21095.0 NA [4,] 38776.5 30288 34899.5 26588.0 28706.0 26763.0 NA [5,] 43129.0 37702 36845.0 31514.0 36438.0 29743.0 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 20509.49 24945.2 17194.02 19853.59 18114.93 17373.85 NA [2,] 35815.51 32212.8 31336.98 26386.41 27346.07 24816.15 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(13253, 21997.5, 28162.5, 38776.5, 43129, 13971, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7sgvd1289748013.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,] 13131.83 13530.00 12939.38 [2,] 20100.50 19729.50 20165.25 [3,] 25425.25 24860.75 25404.69 [4,] 31242.50 30498.25 30812.31 [5,] 34292.67 34179.50 33836.75 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 20343.31 19949.05 20548.49 [2,] 30507.19 29772.45 30260.88 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > try(system("convert tmp/1wofm1289748013.ps tmp/1wofm1289748013.png",intern=TRUE)) character(0) > try(system("convert tmp/2oyfp1289748013.ps tmp/2oyfp1289748013.png",intern=TRUE)) character(0) > try(system("convert tmp/3oyfp1289748013.ps tmp/3oyfp1289748013.png",intern=TRUE)) character(0) > try(system("convert tmp/4zpws1289748013.ps tmp/4zpws1289748013.png",intern=TRUE)) character(0) > try(system("convert tmp/5zpws1289748013.ps tmp/5zpws1289748013.png",intern=TRUE)) character(0) > try(system("convert tmp/6zpws1289748013.ps tmp/6zpws1289748013.png",intern=TRUE)) character(0) > try(system("convert tmp/7sgvd1289748013.ps tmp/7sgvd1289748013.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.894 1.432 2.117