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. 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(10032.80,10204.59,10001.60,10411.75,10673.38,10539.51,10723.78,10682.06,10283.19,10377.18,10486.64,10545.38,10554.27,10532.54,10324.31,10695.25,10827.81,10872.48,10971.19,11145.65,11234.68,11333.88,10997.97,11036.89,11257.35,11533.59,11963.12,12185.15,12377.62,12512.89,12631.48,12268.53,12754.80,13407.75,13480.21,13673.28,13239.71,13557.69,13901.28,13200.58,13406.97,12538.12,12419.57,12193.88,12656.63,12812.48,12056.67,11322.38,11530.75,11114.08) > 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] 50 > (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 4 4 4 4 4 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] 10032.80 10554.27 11257.35 13239.71 11530.75 [2,] 10204.59 10532.54 11533.59 13557.69 11114.08 [3,] 10001.60 10324.31 11963.12 13901.28 NA [4,] 10411.75 10695.25 12185.15 13200.58 NA [5,] 10673.38 10827.81 12377.62 13406.97 NA [6,] 10539.51 10872.48 12512.89 12538.12 NA [7,] 10723.78 10971.19 12631.48 12419.57 NA [8,] 10682.06 11145.65 12268.53 12193.88 NA [9,] 10283.19 11234.68 12754.80 12656.63 NA [10,] 10377.18 11333.88 13407.75 12812.48 NA [11,] 10486.64 10997.97 13480.21 12056.67 NA [12,] 10545.38 11036.89 13673.28 11322.38 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 171.79 -21.73 276.24 317.98 -416.67 [2,] -202.99 -208.23 429.53 343.59 NA [3,] 410.15 370.94 222.03 -700.70 NA [4,] 261.63 132.56 192.47 206.39 NA [5,] -133.87 44.67 135.27 -868.85 NA [6,] 184.27 98.71 118.59 -118.55 NA [7,] -41.72 174.46 -362.95 -225.69 NA [8,] -398.87 89.03 486.27 462.75 NA [9,] 93.99 99.20 652.95 155.85 NA [10,] 109.46 -335.91 72.46 -755.81 NA [11,] 58.74 38.92 193.07 -734.29 NA [12,] 8.89 220.46 -433.57 208.37 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/1gn781229337207.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/2orsl1229337207.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/32xn91229337207.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/43jpt1229337207.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,] 10032.80 10204.59 10001.60 10411.75 10673.38 10539.51 10723.78 10682.06 [2,] 10554.27 10532.54 10162.95 10553.50 10750.59 10705.99 10847.49 10913.85 [3,] 11257.35 11114.08 11143.72 11440.20 11602.72 11692.68 11695.38 11669.76 [4,] 11530.75 11533.59 12932.20 12692.86 12892.30 12525.51 12525.52 12231.20 [5,] 11530.75 11533.59 13901.28 13200.58 13406.97 12538.12 12631.48 12268.53 [,9] [,10] [,11] [,12] [1,] 10283.19 10377.18 10486.64 10545.38 [2,] 10758.94 10855.53 10742.31 10791.13 [3,] 11945.65 12073.18 11527.32 11179.63 [4,] 12705.72 13110.11 12768.44 12497.83 [5,] 12754.80 13407.75 13480.21 13673.28 $n [1] 5 5 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 10567.37 10406.74 8956.011 9750.102 9910.772 10255.27 10369.73 10629.06 [2,] 11947.33 11821.42 13331.419 13130.298 13294.658 13130.10 13021.03 12710.47 [,9] [,10] [,11] [,12] [1,] 10407.70 10292.06 9926.673 9831.346 [2,] 13483.61 13854.30 13127.967 12527.924 $out [1] 13239.71 13557.69 $group [1] 1 2 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(10032.8, 10554.27, 11257.35, 11530.75, 11530.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5nc8c1229337207.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] [1,] -416.67 -208.23 -700.700 132.560 -868.85 -118.55 -362.950 -398.87 93.990 [2,] -21.73 -205.61 -239.335 162.515 -501.36 -9.92 -294.320 -154.92 96.595 [3,] 171.79 70.30 296.485 199.430 -44.60 108.65 -133.705 275.89 127.525 [4,] 276.24 386.56 390.545 234.010 89.97 151.43 66.370 474.51 404.400 [5,] 317.98 429.53 410.150 261.630 135.27 184.27 174.460 486.27 652.950 [,10] [,11] [,12] [1,] -755.810 -734.290 -433.570 [2,] -545.860 -347.685 -212.340 [3,] -131.725 48.830 108.630 [4,] 90.960 125.905 214.415 [5,] 109.460 193.070 220.460 $n [1] 5 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -38.75485 -397.5143 -201.1202 142.9490 -511.7507 -18.8165 -418.6501 [2,] 382.33485 538.1143 794.0902 255.9110 422.5507 236.1165 151.2401 [,8] [,9] [,10] [,11] [,12] [1,] -221.3597 -115.6410 -634.8128 -325.3061 -228.5064 [2,] 773.1397 370.6910 371.3628 422.9661 445.7665 $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(-416.67, -21.7299999999996, 171.790000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6tvub1229337207.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,] 10001.60 10324.31 11257.35 11322.38 11114.08 [2,] 10243.89 10624.76 12074.14 12306.72 11114.08 [3,] 10449.19 10921.83 12445.26 12734.56 11322.42 [4,] 10609.38 11091.27 13081.27 13323.34 11530.75 [5,] 10723.78 11333.88 13673.28 13901.28 11530.75 $n [1] 12 12 12 12 2 $conf [,1] [,2] [,3] [,4] [,5] [1,] 10282.49 10709.06 11985.89 12270.87 10856.90 [2,] 10615.90 11134.61 12904.62 13198.24 11787.93 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(10001.6, 10243.89, 10449.195, 10609.38, 10723.78, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7u1m41229337207.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,] 11322.98 11114.08 11033.07 [2,] 11560.05 11218.49 11378.85 [3,] 11633.83 11565.02 11631.25 [4,] 11743.85 11694.03 11701.51 [5,] 11982.82 12073.18 12028.00 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 11550.00 11348.12 11484.08 [2,] 11717.66 11781.91 11778.41 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(11322.976, 11560.05375, 11633.8325, 11743.84875, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1gn781229337207.ps tmp/1gn781229337207.png") > system("convert tmp/2orsl1229337207.ps tmp/2orsl1229337207.png") > system("convert tmp/32xn91229337207.ps tmp/32xn91229337207.png") > system("convert tmp/43jpt1229337207.ps tmp/43jpt1229337207.png") > system("convert tmp/5nc8c1229337207.ps tmp/5nc8c1229337207.png") > system("convert tmp/6tvub1229337207.ps tmp/6tvub1229337207.png") > system("convert tmp/7u1m41229337207.ps tmp/7u1m41229337207.png") > > > proc.time() user system elapsed 1.209 0.940 1.582