R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(-32.0928134170003,-128.744158053173,-275.904809333092,471.760903738039,-69.1654637163138,346.218578823237,163.216457137027,-198.373745010757,204.401273185691,-175.940320227463,17.4783262653075,223.539436263922,51.4668945111487,260.738288102879,204.977677230572,-232.867546347215,-217.221406054262,387.746049688947,63.5945545423544,364.638829260013,96.4516787351593,-29.7694719663539,286.722082027268,211.416619028106,-185.782838773910,-141.320849733312,238.502069839504,209.511569776478,285.264826114983,20.084745860321,-487.431942052935,107.320923964348,34.3269970752997,141.696685920852,-133.438481853825,-229.255640647698,69.0555896372578,-139.443214460286,80.0405548902168,-140.866735807813,415.577107396144,151.513504277031,-24.7248794807892,8.4258084117008,272.69433043438,313.164475266341,179.140854480554,-495.203375413695) > 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] 48 > (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] 4 4 4 4 4 4 4 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] -32.09281 51.46689 -185.78284 69.055590 NA [2,] -128.74416 260.73829 -141.32085 -139.443214 NA [3,] -275.90481 204.97768 238.50207 80.040555 NA [4,] 471.76090 -232.86755 209.51157 -140.866736 NA [5,] -69.16546 -217.22141 285.26483 415.577107 NA [6,] 346.21858 387.74605 20.08475 151.513504 NA [7,] 163.21646 63.59455 -487.43194 -24.724879 NA [8,] -198.37375 364.63883 107.32092 8.425808 NA [9,] 204.40127 96.45168 34.32700 272.694330 NA [10,] -175.94032 -29.76947 141.69669 313.164475 NA [11,] 17.47833 286.72208 -133.43848 179.140854 NA [12,] 223.53944 211.41662 -229.25564 -495.203375 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] -96.65134 209.27139 44.46199 -208.49880 NA [2,] -147.16065 -55.76061 379.82292 219.48377 NA [3,] 747.66571 -437.84522 -28.99050 -220.90729 NA [4,] -540.92637 15.64614 75.75326 556.44384 NA [5,] 415.38404 604.96746 -265.18008 -264.06360 NA [6,] -183.00212 -324.15150 -507.51669 -176.23838 NA [7,] -361.59020 301.04427 594.75287 33.15069 NA [8,] 402.77502 -268.18715 -72.99393 264.26852 NA [9,] -380.34159 -126.22115 107.36969 40.47014 NA [10,] 193.41865 316.49155 -275.13517 -134.02362 NA [11,] 206.06111 -75.30546 -95.81716 -674.34423 NA [12,] -172.07254 -397.19946 298.31123 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/1pv2g1259942648.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/2shwb1259942648.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/3tdli1259942648.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/479yb1259942648.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,] -185.78284 -141.32085 -275.90481 -232.86755 -217.2214 20.08475 -487.43194 [2,] -108.93783 -140.38203 -97.93213 -186.86714 -143.1934 85.79913 -256.07841 [3,] 9.68704 -134.09369 142.50912 34.32242 108.0497 248.86604 19.43484 [4,] 60.26124 65.99707 221.73987 340.63624 350.4210 366.98231 113.40551 [5,] 69.05559 260.73829 238.50207 471.76090 415.5771 387.74605 163.21646 [,8] [,9] [,10] [,11] [,12] [1,] -198.37375 34.32700 -175.94032 -133.43848 -495.20338 [2,] -94.97397 65.38934 -102.85490 -57.98008 -362.22951 [3,] 57.87337 150.42648 55.96361 98.30959 -8.91951 [4,] 235.97988 238.54780 227.43058 232.93147 217.47803 [5,] 364.63883 272.69433 313.16448 286.72208 223.53944 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -123.9802 -297.1332 -110.0318 -382.4053 -281.9057 26.73132 -272.4575 [2,] 143.3543 28.9458 395.0500 451.0501 498.0051 471.00076 311.3271 [,8] [,9] [,10] [,11] [,12] [1,] -203.5802 13.63129 -204.9619 -131.5105 -466.8885 [2,] 319.3269 287.22166 316.8891 328.1297 449.0494 $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(-185.78283877391, -108.937826095455, 9.6870405470742, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5358u1259942648.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,] -208.49880 -147.16065 -437.8452 -540.9264 -265.18008 -507.5167 -361.5902 [2,] -152.57507 -101.46063 -329.3763 -262.6401 -264.62184 -415.8341 -164.2198 [3,] -26.09468 81.86158 -124.9489 45.6997 75.66022 -253.5768 167.0975 [4,] 126.86669 299.65334 359.3376 316.0985 510.17575 -179.6203 447.8986 [5,] 209.27139 379.82292 747.6657 556.4438 604.96746 -176.2384 594.7529 [,8] [,9] [,10] [,11] [,12] [1,] -268.1872 -380.34159 -275.13517 -674.34423 -397.19946 [2,] -170.5905 -253.28137 -204.57939 -385.08069 -284.63600 [3,] 95.6373 -42.87550 29.69751 -85.56131 -172.07254 [4,] 333.5218 73.91992 254.95510 65.37782 63.11934 [5,] 402.7750 107.36969 316.49155 206.06111 298.31123 $n [1] 4 4 4 4 4 4 4 4 4 4 4 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -246.8537 -235.0185 -669.0328 -411.5038 -536.4299 -440.18574 -316.476 [2,] 194.6643 398.7416 419.1351 502.9032 687.7503 -66.96788 650.671 [,8] [,9] [,10] [,11] [,12] [1,] -302.6114 -301.3645 -333.3347 -441.4235 -489.2996 [2,] 493.8860 215.6135 392.7298 270.3009 145.1546 $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(-208.498804097544, -152.575074366858, -26.0946777977874, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/698hz1259942648.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,] -275.904809 -232.86755 -487.43194 -495.20338 NA [2,] -152.342239 10.84871 -163.55184 -82.08405 NA [3,] -7.307244 150.71468 27.20587 74.54807 NA [4,] 213.970355 273.73019 175.60413 225.91759 NA [5,] 471.760904 387.74605 285.26483 415.57711 NA $n [1] 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [1,] -174.3849 30.81267 -127.4854 -65.93354 NA [2,] 159.7704 270.61668 181.8972 215.02969 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-275.904809333092, -152.342239140318, -7.3072435758464, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7cfv51259942648.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,] -72.37574 -134.09369 -85.64308 [2,] -30.76539 14.56094 -16.63822 [3,] 66.39540 56.91849 61.65694 [4,] 95.54473 125.27940 104.01911 [5,] 226.39072 248.86604 237.62838 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 8.784505 6.419054 6.624326 [2,] 124.006291 107.417919 116.689559 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-72.3757401923413, -30.7653877732995, 66.3953982023351, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1pv2g1259942648.ps tmp/1pv2g1259942648.png") > system("convert tmp/2shwb1259942648.ps tmp/2shwb1259942648.png") > system("convert tmp/3tdli1259942648.ps tmp/3tdli1259942648.png") > system("convert tmp/479yb1259942648.ps tmp/479yb1259942648.png") > system("convert tmp/5358u1259942648.ps tmp/5358u1259942648.png") > system("convert tmp/698hz1259942648.ps tmp/698hz1259942648.png") > system("convert tmp/7cfv51259942648.ps tmp/7cfv51259942648.png") > > > proc.time() user system elapsed 1.221 0.960 1.463