R version 2.7.1 (2008-06-23) 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(100.58,118.48,79.58,81.97,127.13,120.76,120.26,74.9,67.59,87.73,102.87,144.94,110.48,96.34,100.43,90.88,128.28,101.21,73.76,73.64,66.4,57.34,113.59,123.53,102.87,102.99,95.8,98.43,102.65,129.55,100.37,101.93,101.94,93.87,100.91,92.64,101.67,88.67,129.86,98.07,166.45,176.52,82.07,92.18,95.02,84.69,103.01,107.9,204.13,101.99,119.23,95.65,160.95,111.06,150.41,94.79,160.34,104.08,101.07,111.5,136.9,141.71,153.98,134.27,124.71,72.89,101.2,73.28,174.05,111.9,97.06,105.23,109.13,84.04,118.82,90.84,144.28,110.16,86.09,59.87,108.97,94.93,87.36,143.52,108.7,121.13,210.25,110.2,161.46,99.41,132.72,174.29,69.93,83.43,127.53,187.58) > 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] 96 > (np <- floor(n / par1)) [1] 8 > 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] 8 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 100.58 110.48 102.87 101.67 204.13 136.90 109.13 108.70 NA [2,] 118.48 96.34 102.99 88.67 101.99 141.71 84.04 121.13 NA [3,] 79.58 100.43 95.80 129.86 119.23 153.98 118.82 210.25 NA [4,] 81.97 90.88 98.43 98.07 95.65 134.27 90.84 110.20 NA [5,] 127.13 128.28 102.65 166.45 160.95 124.71 144.28 161.46 NA [6,] 120.76 101.21 129.55 176.52 111.06 72.89 110.16 99.41 NA [7,] 120.26 73.76 100.37 82.07 150.41 101.20 86.09 132.72 NA [8,] 74.90 73.64 101.93 92.18 94.79 73.28 59.87 174.29 NA [9,] 67.59 66.40 101.94 95.02 160.34 174.05 108.97 69.93 NA [10,] 87.73 57.34 93.87 84.69 104.08 111.90 94.93 83.43 NA [11,] 102.87 113.59 100.91 103.01 101.07 97.06 87.36 127.53 NA [12,] 144.94 123.53 92.64 107.90 111.50 105.23 143.52 187.58 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 17.90 -14.14 0.12 -13.00 -102.14 4.81 -25.09 12.43 NA [2,] -38.90 4.09 -7.19 41.19 17.24 12.27 34.78 89.12 NA [3,] 2.39 -9.55 2.63 -31.79 -23.58 -19.71 -27.98 -100.05 NA [4,] 45.16 37.40 4.22 68.38 65.30 -9.56 53.44 51.26 NA [5,] -6.37 -27.07 26.90 10.07 -49.89 -51.82 -34.12 -62.05 NA [6,] -0.50 -27.45 -29.18 -94.45 39.35 28.31 -24.07 33.31 NA [7,] -45.36 -0.12 1.56 10.11 -55.62 -27.92 -26.22 41.57 NA [8,] -7.31 -7.24 0.01 2.84 65.55 100.77 49.10 -104.36 NA [9,] 20.14 -9.06 -8.07 -10.33 -56.26 -62.15 -14.04 13.50 NA [10,] 15.14 56.25 7.04 18.32 -3.01 -14.84 -7.57 44.10 NA [11,] 42.07 9.94 -8.27 4.89 10.43 8.17 56.16 60.05 NA [12,] -34.46 -20.66 9.03 96.23 25.40 3.90 -34.82 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/1o2xb1219002530.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/23m4g1219002530.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/31oto1219002530.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/4pci21219002530.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] [1,] 100.580 84.040 79.580 81.970 102.650 72.890 73.760 59.87 66.400 [2,] 102.270 92.505 98.115 90.860 125.920 100.310 84.080 73.46 68.760 [3,] 108.915 102.490 119.025 96.860 136.280 110.610 100.785 83.54 98.480 [4,] 123.690 119.805 141.920 104.315 161.205 125.155 126.490 98.36 134.655 [5,] 136.900 141.710 153.980 110.200 166.450 129.550 150.410 101.93 174.050 [,10] [,11] [,12] [1,] 83.430 87.360 92.640 [2,] 84.060 98.985 106.565 [3,] 90.800 101.970 117.515 [4,] 99.505 108.300 144.230 [5,] 111.900 113.590 187.580 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 96.94948 87.23983 94.5549 89.34384 116.5693 96.73123 77.09417 [2,] 120.88052 117.74017 143.4951 104.37616 155.9907 124.48877 124.47583 [,8] [,9] [,10] [,11] [,12] [1,] 69.6305 61.67011 82.1722 96.7665 96.47479 [2,] 97.4495 135.28989 99.4278 107.1735 138.55521 $out [1] 204.13 210.25 134.27 176.52 174.29 57.34 127.53 $group [1] 1 3 4 6 8 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(100.58, 102.27, 108.915, 123.69, 136.9, 84.04, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5utr71219002530.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,] -25.090 -38.900 -31.790 -9.56 -62.050 -94.450 -55.620 -7.310 -62.150 [2,] -19.615 -1.550 -29.885 20.81 -50.855 -28.315 -36.640 -7.275 -35.150 [3,] -6.440 14.755 -21.645 48.21 -30.595 -12.285 -13.170 1.425 -9.695 [4,] 8.620 37.985 -3.580 59.37 1.850 30.810 5.835 57.325 2.715 [5,] 17.900 89.120 2.630 68.38 26.900 39.350 41.570 100.770 20.140 [,10] [,11] [,12] [1,] -14.84 -8.270 -34.820 [2,] -5.29 6.530 -27.560 [3,] 11.09 10.185 3.900 [4,] 31.21 49.115 17.215 [5,] 56.25 60.050 25.400 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -22.212476 -7.329819 -36.33935 26.66983 -60.036770 -45.31307 -36.89714 [2,] 9.332476 36.839819 -6.95065 69.75017 -1.153230 20.74307 10.55714 [,8] [,9] [,10] [,11] [,12] [1,] -34.66149 -30.84693 -9.299424 -13.60359 -22.83891 [2,] 37.51149 11.45693 31.479424 33.97359 30.63891 $out [1] -102.14 -100.05 -104.36 96.23 $group [1] 1 3 8 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-25.09, -19.615, -6.44, 8.62, 17.9, -38.9, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6szmc1219002530.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] [,8] [,9] [1,] 67.590 57.340 92.640 82.070 94.790 72.890 59.870 69.930 NA [2,] 80.775 73.700 97.115 90.425 101.530 99.130 86.725 104.055 NA [3,] 101.725 98.385 101.420 99.870 111.280 118.305 101.950 124.330 NA [4,] 120.510 112.035 102.760 118.880 155.375 139.305 114.490 167.875 NA [5,] 144.940 128.280 102.990 129.860 204.130 174.050 144.280 210.250 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 83.6016 80.90015 98.84528 86.89148 86.72093 99.98091 89.2862 [2,] 119.8484 115.86985 103.99472 112.84852 135.83907 136.62909 114.6138 [,8] [,9] [1,] 95.22127 NA [2,] 153.43873 NA $out [1] 129.55 166.45 176.52 $group [1] 3 4 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(67.59, 80.775, 101.725, 120.51, 144.94, 57.34, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7crym1219002530.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,] 89.74625 83.5400 85.0625 [2,] 102.10688 97.6700 95.8500 [3,] 106.38938 102.2300 105.5062 [4,] 123.90063 114.0625 114.7200 [5,] 139.48875 136.2800 125.5538 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 96.4491 94.75327 96.89952 [2,] 116.3297 109.70673 114.11298 $out [1] 143.8012 $group [1] 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(89.74625, 102.106875, 106.389375, 123.900625, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1o2xb1219002530.ps tmp/1o2xb1219002530.png") > system("convert tmp/23m4g1219002530.ps tmp/23m4g1219002530.png") > system("convert tmp/31oto1219002530.ps tmp/31oto1219002530.png") > system("convert tmp/4pci21219002530.ps tmp/4pci21219002530.png") > system("convert tmp/5utr71219002530.ps tmp/5utr71219002530.png") > system("convert tmp/6szmc1219002530.ps tmp/6szmc1219002530.png") > system("convert tmp/7crym1219002530.ps tmp/7crym1219002530.png") > > > proc.time() user system elapsed 1.424 0.952 1.689