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(10284.5,12792,12823.61538,13845.66667,15335.63636,11188.5,13633.25,12298.46667,15353.63636,12696.15385,12213.93333,13683.72727,11214.14286,13950.23077,11179.13333,11801.875,11188.82353,16456.27273,11110.0625,16530.69231,10038.41176,11681.25,11148.88235,8631,9386.444444,9764.736842,12043.75,12948.06667,10987.125,11648.3125,10633.35294,10219.3,9037.6,10296.31579,11705.41176,10681.94444,9362.947368,11306.35294,10984.45,10062.61905,8118.583333,8867.48,8346.72,8529.307692,10697.18182,8591.84,8695.607143,8125.571429,7009.758621,7883.466667,7527.645161,6763.758621,6682.333333,7855.681818,6738.88,7895.434783,6361.884615,6935.956522,8344.454545,9107.944444) > 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,] 10284.50 11214.14 9386.444 9362.947 7009.759 NA [2,] 12792.00 13950.23 9764.737 11306.353 7883.467 NA [3,] 12823.62 11179.13 12043.750 10984.450 7527.645 NA [4,] 13845.67 11801.88 12948.067 10062.619 6763.759 NA [5,] 15335.64 11188.82 10987.125 8118.583 6682.333 NA [6,] 11188.50 16456.27 11648.312 8867.480 7855.682 NA [7,] 13633.25 11110.06 10633.353 8346.720 6738.880 NA [8,] 12298.47 16530.69 10219.300 8529.308 7895.435 NA [9,] 15353.64 10038.41 9037.600 10697.182 6361.885 NA [10,] 12696.15 11681.25 10296.316 8591.840 6935.957 NA [11,] 12213.93 11148.88 11705.412 8695.607 8344.455 NA [12,] 13683.73 8631.00 10681.944 8125.571 9107.944 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 2507.50000 2736.0879 378.2924 1943.4056 873.70805 NA [2,] 31.61538 -2771.0974 2279.0132 -321.9029 -355.82151 NA [3,] 1022.05129 622.7417 904.3167 -921.8310 -763.88654 NA [4,] 1489.96969 -613.0515 -1960.9417 -1944.0357 -81.42529 NA [5,] -4147.13636 5267.4492 661.1875 748.8967 1173.34849 NA [6,] 2444.75000 -5346.2102 -1014.9596 -520.7600 -1116.80182 NA [7,] -1334.78333 5420.6298 -414.0529 182.5877 1156.55478 NA [8,] 3055.16969 -6492.2805 -1181.7000 2167.8741 -1533.55017 NA [9,] -2657.48251 1642.8382 1258.7158 -2105.3418 574.07191 NA [10,] -482.22052 -532.3677 1409.0960 103.7671 1408.49802 NA [11,] 1469.79394 -2517.8823 -1023.4673 -570.0357 763.48990 NA [12,] -2469.58441 755.4444 -1318.9971 -1115.8128 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/12tl51257803233.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/224qm1257803233.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/3mv281257803233.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/4vn5i1257803233.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] [1,] 9362.947 7883.467 10984.45 6763.759 6682.333 7855.682 6738.88 [2,] 9362.947 9764.737 10984.45 10062.619 8118.583 8867.480 8346.72 [3,] 9386.444 11306.353 11179.13 11801.875 10987.125 11188.500 10633.35 [4,] 10284.500 12792.000 12043.75 12948.067 11188.824 11648.312 11110.06 [5,] 11214.143 13950.231 12823.62 13845.667 15335.636 11648.312 13633.25 [,8] [,9] [,10] [,11] [,12] [1,] 7895.435 9037.60 6935.957 8344.455 8125.571 [2,] 8529.308 9037.60 8591.840 8695.607 8631.000 [3,] 10219.300 10038.41 10296.316 11148.882 9107.944 [4,] 12298.467 10697.18 11681.250 11705.412 10681.944 [5,] 16530.692 10697.18 12696.154 12213.933 13683.727 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 8735.278 9167.296 10430.63 9763.025 8817.701 9223.57 8680.782 [2,] 10037.611 13445.409 11927.63 13840.725 13156.549 13153.43 12585.924 [,8] [,9] [,10] [,11] [,12] [1,] 7556.022 8865.755 8113.346 9022.162 7658.752 [2,] 12882.578 11211.068 12479.285 13275.603 10557.137 $out [1] 7009.759 7527.645 16456.273 15353.636 6361.885 $group [1] 1 3 6 9 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(9362.947368, 9362.947368, 9386.444444, 10284.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5kabc1257803233.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] [1,] 378.2924 -355.82151 -921.8310 -1960.94167 661.1875 -1116.802 -1334.7833 [2,] 873.7080 -355.82151 -763.8865 -1944.03572 661.1875 -1116.802 -414.0529 [3,] 1943.4056 -321.90294 622.7417 -613.05147 748.8967 -1014.960 182.5877 [4,] 2507.5000 31.61538 904.3167 -81.42529 1173.3485 -520.760 1156.5548 [5,] 2736.0879 31.61538 1022.0513 1489.96969 1173.3485 -520.760 1156.5548 [,8] [,9] [,10] [,11] [,12] [1,] -6492.281 -2657.4825 -532.3677 -2517.8823 -2469.5844 [2,] -1533.550 -2105.3418 -482.2205 -1023.4673 -1894.2907 [3,] -1181.700 574.0719 103.7671 -570.0357 -1217.4049 [4,] 2167.874 1258.7158 1408.4980 763.4899 -180.1842 [5,] 3055.170 1642.8382 1409.0960 1469.7939 755.4444 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 788.9723 -595.66487 -556.0065 -1929.1673 387.005 -1436.1212 -927.1998 [2,] 3097.8389 -48.14101 1801.4899 703.0644 1110.788 -593.7979 1292.3752 [,8] [,9] [,10] [,11] [,12] [1,] -3797.117 -1802.963 -1232.210 -1832.6952 -2571.5491 [2,] 1433.717 2951.107 1439.744 692.6238 136.7392 $out [1] -2771.097 2279.013 -4147.136 5267.449 2444.750 -5346.210 5420.630 $group [1] 2 2 5 5 6 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(378.292398000000, 873.708046, 1943.405572, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6dbs51257803233.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,] 10284.50 8631.00 9037.600 8118.583 6361.885 NA [2,] 12256.20 11129.47 9992.018 8438.014 6751.319 NA [3,] 12807.81 11201.48 10657.649 8781.544 7268.702 NA [4,] 13764.70 12876.05 11676.862 10379.900 7889.451 NA [5,] 15353.64 13950.23 12948.067 11306.353 9107.944 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 12119.77 10404.86 9889.18 7895.836 6749.593 NA [2,] 13495.84 11998.11 11426.12 9667.251 7787.811 NA $out [1] 16456.27 16530.69 $group [1] 2 2 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(10284.5, 12256.2, 12807.80769, 13764.69697, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7dijd1257803233.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,] 9451.559 9107.944 9653.703 [2,] 10069.245 10128.856 9776.057 [3,] 10442.079 10810.239 10168.527 [4,] 11089.519 11183.817 10846.128 [5,] 11203.249 11801.875 11514.100 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 9976.726 10329.06 9680.46 [2,] 10907.433 11291.41 10656.59 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(9451.5586586, 10069.2453023, 10442.0790684, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/12tl51257803233.ps tmp/12tl51257803233.png") > system("convert tmp/224qm1257803233.ps tmp/224qm1257803233.png") > system("convert tmp/3mv281257803233.ps tmp/3mv281257803233.png") > system("convert tmp/4vn5i1257803233.ps tmp/4vn5i1257803233.png") > system("convert tmp/5kabc1257803233.ps tmp/5kabc1257803233.png") > system("convert tmp/6dbs51257803233.ps tmp/6dbs51257803233.png") > system("convert tmp/7dijd1257803233.ps tmp/7dijd1257803233.png") > > > proc.time() user system elapsed 1.874 1.479 2.219