R version 2.10.1 (2009-12-14) 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(52347,52407,50570,50442,50590,50040,50476,50268,50595,48708,48547,48196,48375,47915,46462,46132,46308,46532,46817,46824,46263,45992,46404,46995,48102,48719,47912,48430,50141,50608,51005,51857,52513,52406,53634,55165,57294,58026,56701,58706,60103,61153,62395,63850,64534,65765,66954,65741,65474,60687,59227,59373,59995,59532,59696,59507,60210,58782,59372,58827,60481,59508,56565,56201,56193,56431,56316,55316,54795,53310,51848,50618,52026,50120,46825,46374,45441,45392,45032,44302,42880,42101,41886,41415,43228,41633,39375,38603,37847,36881,36700,36477,35684,35896,37109,37612,39570,39518,37970,38343,37966,38942,39304,39438,38999,38110,40024,41050,42239,42313,41159,42067,42515,43554,45018,45797,46749,47291,48800,50566,54884,54002,51813,52751,54461,55364,56900,57795) > 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] 128 > (np <- floor(n / par1)) [1] 10 > 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] 11 11 11 11 11 11 11 11 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 52347 48375 48102 57294 65474 60481 52026 43228 39570 42239 54884 [2,] 52407 47915 48719 58026 60687 59508 50120 41633 39518 42313 54002 [3,] 50570 46462 47912 56701 59227 56565 46825 39375 37970 41159 51813 [4,] 50442 46132 48430 58706 59373 56201 46374 38603 38343 42067 52751 [5,] 50590 46308 50141 60103 59995 56193 45441 37847 37966 42515 54461 [6,] 50040 46532 50608 61153 59532 56431 45392 36881 38942 43554 55364 [7,] 50476 46817 51005 62395 59696 56316 45032 36700 39304 45018 56900 [8,] 50268 46824 51857 63850 59507 55316 44302 36477 39438 45797 57795 [9,] 50595 46263 52513 64534 60210 54795 42880 35684 38999 46749 NA [10,] 48708 45992 52406 65765 58782 53310 42101 35896 38110 47291 NA [11,] 48547 46404 53634 66954 59372 51848 41886 37109 40024 48800 NA [12,] 48196 46995 55165 65741 58827 50618 41415 37612 41050 50566 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 60 -460 617 732 -4787 -973 -1906 -1595 -52 74 -882 [2,] -1837 -1453 -807 -1325 -1460 -2943 -3295 -2258 -1548 -1154 -2189 [3,] -128 -330 518 2005 146 -364 -451 -772 373 908 938 [4,] 148 176 1711 1397 622 -8 -933 -756 -377 448 1710 [5,] -550 224 467 1050 -463 238 -49 -966 976 1039 903 [6,] 436 285 397 1242 164 -115 -360 -181 362 1464 1536 [7,] -208 7 852 1455 -189 -1000 -730 -223 134 779 895 [8,] 327 -561 656 684 703 -521 -1422 -793 -439 952 NA [9,] -1887 -271 -107 1231 -1428 -1485 -779 212 -889 542 NA [10,] -161 412 1228 1189 590 -1462 -215 1213 1914 1509 NA [11,] -351 591 1531 -1213 -545 -1230 -471 503 1026 1766 NA [12,] 179 1107 2129 -267 1654 1408 1813 1958 1189 4318 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/yougetitorg/rcomp/tmp/10pli1294083260.ps",horizontal=F,onefile=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/yougetitorg/rcomp/tmp/215s31294083260.ps",horizontal=F,onefile=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/yougetitorg/rcomp/tmp/349kt1294083260.ps",horizontal=F,onefile=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/yougetitorg/rcomp/tmp/4pp8u1294083260.ps",horizontal=F,onefile=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,] 39570 39518 37970.0 38343.0 37847 36881.0 36700 36477.0 35684 35896.0 [2,] 45665 45114 43810.5 44099.5 43978 44473.0 45025 45049.5 42880 42101.0 [3,] 52026 50120 47912.0 48430.0 50141 50040.0 50476 50268.0 48672 47999.5 [4,] 56089 56014 54189.0 54476.0 55327 55897.5 56608 56555.5 54795 53310.0 [5,] 65474 60687 59227.0 59373.0 60103 61153.0 62395 63850.0 64534 65765.0 [,11] [,12] [1,] 37109.0 37612 [2,] 41886.0 41415 [3,] 48673.5 49381 [4,] 53634.0 55165 [5,] 66954.0 65741 $n [1] 11 11 11 11 11 11 11 11 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 47060.13 44927.37 42967.81 43486.76 44734.47 44597.51 44958 44786.68 [2,] 56991.87 55312.63 52856.19 53373.24 55547.53 55482.49 55994 55749.32 [,9] [,10] [,11] [,12] [1,] 42718.79 42399.04 42803.73 42510.95 [2,] 54625.21 53599.96 54543.27 56251.05 $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(39570, 45665, 52026, 56089, 65474, 39518, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/5lhg31294083260.ps",horizontal=F,onefile=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] [,11] [1,] -1906 -3295.0 -772 -933.0 -966.0 -360.0 -1000.0 -1422 -1887 -1462.0 -1230 [2,] -1284 -2223.5 -347 -192.5 -256.0 24.5 -215.5 -561 -1428 -161.0 -545 [3,] -460 -1548.0 146 176.0 238.0 362.0 7.0 -56 -525 889.5 76 [4,] 67 -1389.0 713 1009.5 939.5 839.0 815.5 684 212 1228.0 1026 [5,] 732 -807.0 2005 1711.0 1050.0 1536.0 1455.0 952 1231 1914.0 1766 [,12] [1,] 179 [2,] 1107 [3,] 1531 [4,] 1958 [5,] 2129 $n [1] 11 11 11 11 11 11 11 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1103.6001 -1945.546 -358.9712 -396.6183 -331.5218 -26.01796 -484.1560 [2,] 183.6001 -1150.454 650.9712 748.6183 807.5218 750.01796 498.1560 [,8] [,9] [,10] [,11] [,12] [1,] -678.0516 -1344.4094 195.5002 -708.9342 1105.806 [2,] 566.0516 294.4094 1583.4998 860.9342 1956.194 $out [1] -4787 -267 4318 $group [1] 1 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1906, -1284, -460, 67, 732, -3295, -2223.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/684mc1294083260.ps",horizontal=F,onefile=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,] 48196.0 45992.0 47912.0 56701.0 58782.0 50618.0 41415.0 35684.0 37966.0 [2,] 49374.0 46285.5 48574.5 58366.0 59299.5 54052.5 42490.5 36588.5 38226.5 [3,] 50459.0 46497.0 50806.5 61774.0 59519.5 56197.0 45212.0 37360.5 39151.5 [4,] 50592.5 46909.5 52459.5 65137.5 60102.5 56498.0 46599.5 38989.0 39544.0 [5,] 52407.0 46995.0 55165.0 66954.0 60687.0 59508.0 52026.0 41633.0 41050.0 [,10] [,11] [1,] 41159 51813.0 [2,] 42276 53376.5 [3,] 44286 54672.5 [4,] 47020 56132.0 [5,] 50566 57795.0 $n [1] 12 12 12 12 12 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 49903.23 46212.39 49034.53 58685.47 59153.25 55081.59 43337.86 36265.62 [2,] 51014.77 46781.61 52578.47 64862.53 59885.75 57312.41 47086.14 38455.38 [,9] [,10] [,11] [1,] 38550.58 42122.23 53133.24 [2,] 39752.42 46449.77 56211.76 $out [1] 48375 47915 65474 60481 43228 $group [1] 2 2 5 6 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(48196, 49374, 50459, 50592.5, 52407, 45992, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/7wnz01294083260.ps",horizontal=F,onefile=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,] 48598.09 47912.0 48078.88 [2,] 49044.64 48551.0 48697.12 [3,] 49475.67 49710.5 49470.12 [4,] 50049.55 50204.5 50683.25 [5,] 51274.55 52026.0 50877.00 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 49017.33 48956.33 48564.24 [2,] 49934.02 50464.67 50376.01 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(48598.0909090909, 49044.6363636364, 49475.6727272727, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/10pli1294083260.ps tmp/10pli1294083260.png",intern=TRUE)) character(0) > try(system("convert tmp/215s31294083260.ps tmp/215s31294083260.png",intern=TRUE)) character(0) > try(system("convert tmp/349kt1294083260.ps tmp/349kt1294083260.png",intern=TRUE)) character(0) > try(system("convert tmp/4pp8u1294083260.ps tmp/4pp8u1294083260.png",intern=TRUE)) character(0) > try(system("convert tmp/5lhg31294083260.ps tmp/5lhg31294083260.png",intern=TRUE)) character(0) > try(system("convert tmp/684mc1294083260.ps tmp/684mc1294083260.png",intern=TRUE)) character(0) > try(system("convert tmp/7wnz01294083260.ps tmp/7wnz01294083260.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.500 1.220 2.994