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(93.2 + ,96 + ,95.2 + ,77.1 + ,70.9 + ,64.8 + ,70.1 + ,77.3 + ,79.5 + ,100.6 + ,100.7 + ,107.1 + ,95.9 + ,82.8 + ,83.3 + ,80 + ,80.4 + ,67.5 + ,75.7 + ,71.1 + ,89.3 + ,101.1 + ,105.2 + ,114.1 + ,96.3 + ,84.4 + ,91.2 + ,81.9 + ,80.5 + ,70.4 + ,74.8 + ,75.9 + ,86.3 + ,98.7 + ,100.9 + ,113.8 + ,89.8 + ,84.4 + ,87.2 + ,85.6 + ,72 + ,69.2 + ,77.5 + ,78.1 + ,94.3 + ,97.7 + ,100.2 + ,116.4 + ,97.1 + ,93 + ,96 + ,80.5 + ,76.1 + ,69.9 + ,73.6 + ,92.6 + ,94.2 + ,93.5 + ,108.5 + ,109.4 + ,105.1 + ,92.5 + ,97.1 + ,81.4 + ,79.1 + ,72.1 + ,78.7 + ,87.1 + ,91.4 + ,109.9 + ,116.3 + ,113 + ,100 + ,84.8 + ,94.3 + ,87.1 + ,90.3 + ,72.4 + ,84.9 + ,92.7 + ,92.2 + ,114.9 + ,112.5 + ,118.3 + ,106 + ,91.2 + ,96.6 + ,96.3 + ,88.2 + ,70.2 + ,86.5 + ,88.2 + ,102.8 + ,119.1 + ,119.2 + ,125.1 + ,106.1 + ,102.1 + ,105.2 + ,101 + ,84.3 + ,87.5 + ,92.7 + ,94.4 + ,113 + ,113.9 + ,122.9 + ,132.7 + ,106.9 + ,96.6 + ,127.3 + ,98.2 + ,100.2 + ,89.4 + ,95.3 + ,104.2 + ,106.4 + ,116.2 + ,135.9 + ,134 + ,104.6 + ,107.1 + ,123.5 + ,98.8 + ,98.6 + ,90.6 + ,89.1 + ,105.2 + ,114 + ,122.1 + ,138 + ,142.2 + ,116.4 + ,112.6 + ,123.8 + ,103.6 + ,113.9 + ,98.6 + ,95 + ,116 + ,113.9 + ,127.5 + ,131.4 + ,145.9 + ,131.5 + ,131 + ,130.5 + ,118.9 + ,114.3 + ,85.7 + ,104.6 + ,105.1 + ,117.3 + ,142.5 + ,140 + ,159.8 + ,131.2 + ,125.4 + ,126.5 + ,119.4 + ,113.5 + ,98.7 + ,114.5 + ,113.8 + ,133.1 + ,143.4 + ,137.3 + ,165.2 + ,126.9 + ,124 + ,135.7 + ,130 + ,109.4 + ,117.8 + ,120.3 + ,121 + ,132.3 + ,142.9 + ,147.4 + ,175.9 + ,132.6 + ,123.7 + ,153.3 + ,134 + ,119.6 + ,116.2 + ,118.6 + ,130.7 + ,129.3 + ,144.4 + ,163.2 + ,179.4 + ,128.1 + ,138.4 + ,152.7 + ,120 + ,140.5 + ,116.2 + ,121.4 + ,127.8 + ,143.6 + ,157.6 + ,166.2 + ,182.3 + ,153.1 + ,147.6 + ,157.7 + ,137.2 + ,151.5 + ,98.7 + ,145.8 + ,151.7 + ,129.4 + ,174.1 + ,197 + ,193.9 + ,164.1 + ,142.8 + ,157.9 + ,159.2 + ,162.2 + ,123.1 + ,130 + ,150.1 + ,169.4 + ,179.7 + ,182.1 + ,194.3 + ,161.4 + ,169.4 + ,168.8 + ,158.1 + ,158.5 + ,135.3 + ,149.3 + ,143.4 + ,142.2 + ,188.4 + ,166.2 + ,199.2 + ,182.7 + ,145.2 + ,182.1 + ,158.7 + ,141.6 + ,132.6 + ,139.6 + ,147 + ,166.6 + ,157 + ,180.4 + ,210.2 + ,159.8 + ,157.8 + ,168.2 + ,158.4 + ,152 + ,142.2 + ,137.2 + ,152.6 + ,166.8 + ,165.6 + ,198.6 + ,201.5 + ,170.7 + ,164.4 + ,179.7 + ,157 + ,168 + ,139.3 + ,138.6 + ,153.4 + ,138.9 + ,172.1 + ,198.4 + ,217.8 + ,173.7 + ,153.8 + ,175.6 + ,147.1 + ,160.3 + ,135.2 + ,148.8 + ,151 + ,148.2 + ,182.2 + ,189.2 + ,183.1 + ,170 + ,158.4 + ,176.1 + ,156.2 + ,153.2 + ,117.9 + ,149.8 + ,156.6 + ,166.7 + ,156.8 + ,158.6 + ,210.8 + ,203.6 + ,175.2 + ,168.7 + ,155.9 + ,147.3 + ,137 + ,141.1 + ,167.4 + ,160.2 + ,191.9 + ,174.4 + ,208.2 + ,159.4 + ,161.1 + ,172.1 + ,158.4 + ,114.6 + ,159.6 + ,159.7 + ,159.4 + ,160.7 + ,165.5 + ,205 + ,205.2 + ,141.6 + ,148.1 + ,184.9 + ,132.5 + ,137.3 + ,135.5 + ,121.7 + ,166.1 + ,146.8 + ,162.8 + ,186.8 + ,185.5 + ,151.5 + ,158.1 + ,143 + ,151.2 + ,147.6 + ,130.7 + ,137.5 + ,146.1 + ,133.6 + ,167.9 + ,181.9 + ,202 + ,166.5 + ,151.3 + ,146.2 + ,148.3 + ,144.7 + ,123.6 + ,151.6 + ,133.9 + ,137.4 + ,181.6 + ,182 + ,190 + ,161.2 + ,155.5 + ,141.9 + ,164.6 + ,136.2 + ,126.8 + ,152.5 + ,126.6 + ,150.1 + ,186.3 + ,147.5 + ,200.4 + ,177.2 + ,127.4 + ,177.1 + ,154.4 + ,135.2 + ,126.4 + ,147.3 + ,140.6 + ,152.3 + ,151.2 + ,172.2 + ,215.3 + ,154.1 + ,159.3 + ,160.4 + ,151.9 + ,148.4 + ,139.6 + ,148.2 + ,153.5 + ,145.1 + ,183.7 + ,210.5 + ,203.3 + ,153.3 + ,144.3 + ,169.6 + ,143.7 + ,160.1 + ,135.6 + ,141.8 + ,159.9 + ,145.7 + ,183.5 + ,198.2 + ,186.8 + ,172 + ,150.6 + ,163.3 + ,153.7 + ,152.9 + ,135.5 + ,148.5 + ,148.4 + ,133.6 + ,194.1 + ,208.6 + ,197.3 + ,164.4 + ,148.1 + ,152 + ,144.1 + ,155 + ,124.5 + ,153 + ,146 + ,138 + ,190 + ,192 + ,192 + ,147 + ,133 + ,163 + ,150 + ,129 + ,131 + ,145 + ,137 + ,138 + ,168 + ,176 + ,188 + ,139 + ,143 + ,150 + ,154 + ,137 + ,129 + ,128 + ,140 + ,143 + ,151 + ,177 + ,184) > par1 = '37' > #'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] 456 > (np <- floor(n / par1)) [1] 12 > 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] 13 13 13 13 13 13 13 13 13 13 13 13 12 12 12 12 12 12 12 12 12 12 12 12 12 [26] 12 12 12 12 12 12 12 12 12 12 12 12 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 93.2 84.4 94.3 98.2 114.3 116.2 130.0 152.6 166.7 162.8 147.5 186.8 [2,] 96.0 87.2 87.1 100.2 85.7 118.6 150.1 166.8 156.8 186.8 200.4 172.0 [3,] 95.2 85.6 90.3 89.4 104.6 130.7 169.4 165.6 158.6 185.5 177.2 150.6 [4,] 77.1 72.0 72.4 95.3 105.1 129.3 179.7 198.6 210.8 151.5 127.4 163.3 [5,] 70.9 69.2 84.9 104.2 117.3 144.4 182.1 201.5 203.6 158.1 177.1 153.7 [6,] 64.8 77.5 92.7 106.4 142.5 163.2 194.3 170.7 175.2 143.0 154.4 152.9 [7,] 70.1 78.1 92.2 116.2 140.0 179.4 161.4 164.4 168.7 151.2 135.2 135.5 [8,] 77.3 94.3 114.9 135.9 159.8 128.1 169.4 179.7 155.9 147.6 126.4 148.5 [9,] 79.5 97.7 112.5 134.0 131.2 138.4 168.8 157.0 147.3 130.7 147.3 148.4 [10,] 100.6 100.2 118.3 104.6 125.4 152.7 158.1 168.0 137.0 137.5 140.6 133.6 [11,] 100.7 116.4 106.0 107.1 126.5 120.0 158.5 139.3 141.1 146.1 152.3 194.1 [12,] 107.1 97.1 91.2 123.5 119.4 140.5 135.3 138.6 167.4 133.6 151.2 208.6 [13,] 95.9 93.0 96.6 98.8 113.5 116.2 149.3 153.4 160.2 167.9 172.2 197.3 [14,] 82.8 96.0 96.3 98.6 98.7 121.4 143.4 138.9 191.9 181.9 215.3 164.4 [15,] 83.3 80.5 88.2 90.6 114.5 127.8 142.2 172.1 174.4 202.0 154.1 148.1 [16,] 80.0 76.1 70.2 89.1 113.8 143.6 188.4 198.4 208.2 166.5 159.3 152.0 [17,] 80.4 69.9 86.5 105.2 133.1 157.6 166.2 217.8 159.4 151.3 160.4 144.1 [18,] 67.5 73.6 88.2 114.0 143.4 166.2 199.2 173.7 161.1 146.2 151.9 155.0 [19,] 75.7 92.6 102.8 122.1 137.3 182.3 182.7 153.8 172.1 148.3 148.4 124.5 [20,] 71.1 94.2 119.1 138.0 165.2 153.1 145.2 175.6 158.4 144.7 139.6 153.0 [21,] 89.3 93.5 119.2 142.2 126.9 147.6 182.1 147.1 114.6 123.6 148.2 146.0 [22,] 101.1 108.5 125.1 116.4 124.0 157.7 158.7 160.3 159.6 151.6 153.5 138.0 [23,] 105.2 109.4 106.1 112.6 135.7 137.2 141.6 135.2 159.7 133.9 145.1 190.0 [24,] 114.1 105.1 102.1 123.8 130.0 151.5 132.6 148.8 159.4 137.4 183.7 192.0 [25,] 96.3 92.5 105.2 103.6 109.4 98.7 139.6 151.0 160.7 181.6 210.5 192.0 [26,] 84.4 97.1 101.0 113.9 117.8 145.8 147.0 148.2 165.5 182.0 203.3 147.0 [27,] 91.2 81.4 84.3 98.6 120.3 151.7 166.6 182.2 205.0 190.0 153.3 133.0 [28,] 81.9 79.1 87.5 95.0 121.0 129.4 157.0 189.2 205.2 161.2 144.3 163.0 [29,] 80.5 72.1 92.7 116.0 132.3 174.1 180.4 183.1 141.6 155.5 169.6 150.0 [30,] 70.4 78.7 94.4 113.9 142.9 197.0 210.2 170.0 148.1 141.9 143.7 129.0 [31,] 74.8 87.1 113.0 127.5 147.4 193.9 159.8 158.4 184.9 164.6 160.1 131.0 [32,] 75.9 91.4 113.9 131.4 175.9 164.1 157.8 176.1 132.5 136.2 135.6 145.0 [33,] 86.3 109.9 122.9 145.9 132.6 142.8 168.2 156.2 137.3 126.8 141.8 137.0 [34,] 98.7 116.3 132.7 131.5 123.7 157.9 158.4 153.2 135.5 152.5 159.9 138.0 [35,] 100.9 113.0 106.9 131.0 153.3 159.2 152.0 117.9 121.7 126.6 145.7 168.0 [36,] 113.8 100.0 96.6 130.5 134.0 162.2 142.2 149.8 166.1 150.1 183.5 176.0 [37,] 89.8 84.8 127.3 118.9 119.6 123.1 137.2 156.6 146.8 186.3 198.2 188.0 [,13] [1,] 139 [2,] 143 [3,] 150 [4,] 154 [5,] 137 [6,] 129 [7,] 128 [8,] 140 [9,] 143 [10,] 151 [11,] 177 [12,] 184 [13,] NA [14,] NA [15,] NA [16,] NA [17,] NA [18,] NA [19,] NA [20,] NA [21,] NA [22,] NA [23,] NA [24,] NA [25,] NA [26,] NA [27,] NA [28,] NA [29,] NA [30,] NA [31,] NA [32,] NA [33,] NA [34,] NA [35,] NA [36,] NA [37,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 2.8 2.8 -7.2 2.0 -28.6 2.4 20.1 14.2 -9.9 24.0 52.9 -14.8 [2,] -0.8 -1.6 3.2 -10.8 18.9 12.1 19.3 -1.2 1.8 -1.3 -23.2 -21.4 [3,] -18.1 -13.6 -17.9 5.9 0.5 -1.4 10.3 33.0 52.2 -34.0 -49.8 12.7 [4,] -6.2 -2.8 12.5 8.9 12.2 15.1 2.4 2.9 -7.2 6.6 49.7 -9.6 [5,] -6.1 8.3 7.8 2.2 25.2 18.8 12.2 -30.8 -28.4 -15.1 -22.7 -0.8 [6,] 5.3 0.6 -0.5 9.8 -2.5 16.2 -32.9 -6.3 -6.5 8.2 -19.2 -17.4 [7,] 7.2 16.2 22.7 19.7 19.8 -51.3 8.0 15.3 -12.8 -3.6 -8.8 13.0 [8,] 2.2 3.4 -2.4 -1.9 -28.6 10.3 -0.6 -22.7 -8.6 -16.9 20.9 -0.1 [9,] 21.1 2.5 5.8 -29.4 -5.8 14.3 -10.7 11.0 -10.3 6.8 -6.7 -14.8 [10,] 0.1 16.2 -12.3 2.5 1.1 -32.7 0.4 -28.7 4.1 8.6 11.7 60.5 [11,] 6.4 -19.3 -14.8 16.4 -7.1 20.5 -23.2 -0.7 26.3 -12.5 -1.1 14.5 [12,] -11.2 -4.1 5.4 -24.7 -5.9 -24.3 14.0 14.8 -7.2 34.3 21.0 -11.3 [13,] -13.1 3.0 -0.3 -0.2 -14.8 5.2 -5.9 -14.5 31.7 14.0 43.1 -32.9 [14,] 0.5 -15.5 -8.1 -8.0 15.8 6.4 -1.2 33.2 -17.5 20.1 -61.2 -16.3 [15,] -3.3 -4.4 -18.0 -1.5 -0.7 15.8 46.2 26.3 33.8 -35.5 5.2 3.9 [16,] 0.4 -6.2 16.3 16.1 19.3 14.0 -22.2 19.4 -48.8 -15.2 1.1 -7.9 [17,] -12.9 3.7 1.7 8.8 10.3 8.6 33.0 -44.1 1.7 -5.1 -8.5 10.9 [18,] 8.2 19.0 14.6 8.1 -6.1 16.1 -16.5 -19.9 11.0 2.1 -3.5 -30.5 [19,] -4.6 1.6 16.3 15.9 27.9 -29.2 -37.5 21.8 -13.7 -3.6 -8.8 28.5 [20,] 18.2 -0.7 0.1 4.2 -38.3 -5.5 36.9 -28.5 -43.8 -21.1 8.6 -7.0 [21,] 11.8 15.0 5.9 -25.8 -2.9 10.1 -23.4 13.2 45.0 28.0 5.3 -8.0 [22,] 4.1 0.9 -19.0 -3.8 11.7 -20.5 -17.1 -25.1 0.1 -17.7 -8.4 52.0 [23,] 8.9 -4.3 -4.0 11.2 -5.7 14.3 -9.0 13.6 -0.3 3.5 38.6 2.0 [24,] -17.8 -12.6 3.1 -20.2 -20.6 -52.8 7.0 2.2 1.3 44.2 26.8 0.0 [25,] -11.9 4.6 -4.2 10.3 8.4 47.1 7.4 -2.8 4.8 0.4 -7.2 -45.0 [26,] 6.8 -15.7 -16.7 -15.3 2.5 5.9 19.6 34.0 39.5 8.0 -50.0 -14.0 [27,] -9.3 -2.3 3.2 -3.6 0.7 -22.3 -9.6 7.0 0.2 -28.8 -9.0 30.0 [28,] -1.4 -7.0 5.2 21.0 11.3 44.7 23.4 -6.1 -63.6 -5.7 25.3 -13.0 [29,] -10.1 6.6 1.7 -2.1 10.6 22.9 29.8 -13.1 6.5 -13.6 -25.9 -21.0 [30,] 4.4 8.4 18.6 13.6 4.5 -3.1 -50.4 -11.6 36.8 22.7 16.4 2.0 [31,] 1.1 4.3 0.9 3.9 28.5 -29.8 -2.0 17.7 -52.4 -28.4 -24.5 14.0 [32,] 10.4 18.5 9.0 14.5 -43.3 -21.3 10.4 -19.9 4.8 -9.4 6.2 -8.0 [33,] 12.4 6.4 9.8 -14.4 -8.9 15.1 -9.8 -3.0 -1.8 25.7 18.1 1.0 [34,] 2.2 -3.3 -25.8 -0.5 29.6 1.3 -6.4 -35.3 -13.8 -25.9 -14.2 30.0 [35,] 12.9 -13.0 -10.3 -0.5 -19.3 3.0 -9.8 31.9 44.4 23.5 37.8 8.0 [36,] -24.0 -15.2 30.7 -11.6 -14.4 -39.1 -5.0 6.8 -19.3 36.2 14.7 12.0 [37,] -5.4 9.5 -29.1 -4.6 -3.4 6.9 15.4 10.1 16.0 -38.8 -11.4 -49.0 [,13] [1,] 4 [2,] 7 [3,] 4 [4,] -17 [5,] -8 [6,] -1 [7,] 12 [8,] 3 [9,] 8 [10,] 26 [11,] 7 [12,] NA [13,] NA [14,] NA [15,] NA [16,] NA [17,] NA [18,] NA [19,] NA [20,] NA [21,] NA [22,] NA [23,] NA [24,] NA [25,] NA [26,] NA [27,] NA [28,] NA [29,] NA [30,] NA [31,] NA [32,] NA [33,] NA [34,] NA [35,] NA [36,] NA [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/rcomp/tmp/13zix1272298206.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/rcomp/tmp/2w9hi1272298206.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/rcomp/tmp/3w9hi1272298206.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/rcomp/tmp/4w9hi1272298206.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] [,10] [,11] [,12] [1,] 84.4 85.7 85.6 72.0 69.2 64.8 70.1 94.3 112.5 100.2 100.7 91.2 [2,] 98.2 96.0 95.2 95.3 104.2 106.4 116.2 126.4 130.7 118.3 116.4 119.4 [3,] 130.0 143.0 150.0 129.3 144.4 143.0 135.5 140.0 138.4 137.0 139.3 135.3 [4,] 152.6 166.8 165.6 163.3 177.1 163.2 161.4 155.9 147.3 151.0 152.3 151.2 [5,] 186.8 200.4 185.5 210.8 203.6 194.3 179.4 179.7 168.8 168.0 194.1 184.0 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] 93.00 82.80 80.5 70.20 69.90 67.50 75.70 94.20 89.30 101.1 105.20 [2,] 97.70 97.45 89.4 84.55 95.85 101.10 112.45 128.55 116.90 120.2 111.00 [3,] 132.75 130.15 135.0 147.80 147.70 149.05 142.80 144.95 134.55 144.8 135.45 [4,] 164.05 173.15 163.1 177.45 159.90 163.65 162.95 155.75 147.35 158.2 143.35 [5,] 197.30 215.30 202.0 208.20 217.80 199.20 182.70 175.60 182.10 160.3 190.00 [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [1,] 102.10 92.50 84.40 81.40 79.10 72.10 70.40 74.80 75.90 109.90 [2,] 118.95 101.15 107.45 94.90 91.25 104.35 104.15 120.25 122.65 124.85 [3,] 135.00 124.50 146.40 142.35 136.85 145.80 142.40 152.90 135.90 137.15 [4,] 155.45 171.15 156.85 174.40 162.10 171.85 159.05 162.35 160.95 144.35 [5,] 192.00 210.50 203.30 205.00 205.20 183.10 210.20 193.90 176.10 168.20 [,34] [,35] [,36] [,37] [1,] 98.70 100.90 96.60 84.80 [2,] 127.60 115.45 122.15 119.25 [3,] 136.75 128.80 146.00 132.25 [4,] 155.55 152.65 164.15 171.45 [5,] 159.90 168.00 183.50 198.20 $n [1] 13 13 13 13 13 13 13 13 13 13 13 13 12 12 12 12 12 12 12 12 12 12 12 12 12 [26] 12 12 12 12 12 12 12 12 12 12 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 106.1612 111.9745 119.1498 99.5015 112.4543 118.1095 115.6928 127.0727 [2,] 153.8388 174.0255 180.8502 159.0985 176.3457 167.8905 155.3072 152.9273 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 131.1257 122.6704 123.5681 121.3648 102.4873 95.62272 101.3849 105.4277 [2,] 145.6743 151.3296 155.0319 149.2352 163.0127 164.67728 168.6151 190.1723 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 118.4864 120.5205 119.7666 132.5439 120.6616 127.4679 120.6949 118.3521 [2,] 176.9136 177.5795 165.8334 157.3561 148.4384 162.1321 150.2051 151.6479 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 92.57253 123.8683 106.0895 104.5348 115.0128 117.3597 133.6979 118.4311 [2,] 156.42747 168.9317 178.6105 169.1652 176.5872 167.4403 172.1021 153.3689 [,33] [,34] [,35] [,36] [,37] [1,] 128.2559 124.0018 111.8328 126.8435 108.4412 [2,] 146.0441 149.4982 145.7672 165.1565 156.0588 $out [1] 77.3 79.5 97.7 208.6 71.1 86.3 $group [1] 8 9 9 12 20 33 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" Warning message: In bxp(list(stats = c(84.4, 98.2, 130, 152.6, 186.8, 85.7, 96, 143, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/560z31272298206.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] [,10] [,11] [,12] [1,] -28.6 -10.8 -49.8 -17.0 -30.8 -19.2 -12.8 -22.7 -29.4 -12.3 -23.2 -24.70 [2,] -7.2 -1.6 -17.9 -6.2 -15.1 -6.5 -3.6 -8.6 -10.3 0.1 -12.5 -11.25 [3,] 2.8 -0.8 0.5 2.9 -0.8 -1.0 12.0 -0.6 2.5 2.5 -0.7 -5.00 [4,] 14.2 7.0 10.3 12.2 8.3 5.3 16.2 3.0 8.0 11.7 14.5 14.40 [5,] 24.0 19.3 52.2 15.1 25.2 16.2 22.7 10.3 21.1 26.0 26.3 34.30 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] -32.90 -17.5 -35.50 -48.80 -12.90 -30.5 -37.50 -43.8 -25.80 -25.10 -9.00 [2,] -13.80 -15.9 -3.85 -11.55 -6.80 -11.3 -11.25 -24.8 -5.45 -18.35 -4.15 [3,] -0.25 -4.6 1.60 0.75 2.70 5.1 -1.00 -3.1 8.00 -6.10 2.75 [4,] 9.60 11.1 21.05 16.20 9.55 12.8 19.05 6.4 14.10 2.50 12.40 [5,] 43.10 33.2 46.20 19.40 33.00 19.0 28.50 36.9 28.00 11.70 14.30 [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [1,] -52.80 -11.9 -50.0 -22.30 -13.00 -25.90 -11.60 -52.40 -43.30 -14.40 -35.30 [2,] -19.00 -5.7 -15.5 -9.45 -6.55 -13.35 -0.55 -26.45 -14.65 -5.95 -20.00 [3,] 0.65 2.5 4.2 -2.95 1.90 -0.20 6.45 1.00 5.50 3.70 -4.85 [4,] 5.05 7.9 13.8 1.95 22.20 8.60 17.50 9.15 10.40 13.75 1.75 [5,] 26.80 10.3 39.5 7.00 44.70 29.80 36.80 28.50 18.50 25.70 30.00 [,35] [,36] [,37] [1,] -19.30 -39.10 -49.00 [2,] -10.05 -17.25 -20.25 [3,] 5.50 -8.30 -4.00 [4,] 27.70 13.35 9.80 [5,] 44.40 36.20 16.00 $n [1] 13 13 13 13 13 13 13 13 13 13 13 12 12 12 12 12 12 12 12 12 12 12 12 12 12 [26] 12 12 12 12 12 12 12 12 12 12 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -6.577762 -4.568633 -11.85761 -5.163122 -11.054188 -6.170915 3.323380 [2,] 12.177762 2.968633 12.85761 10.963122 9.454188 4.170915 20.676620 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] -5.683273 -5.519301 -2.583273 -12.53176 -16.699137 -10.92290 -16.914881 [2,] 4.483273 10.519301 7.583273 11.13176 6.699137 10.42290 7.714881 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] -9.757057 -11.90696 -4.757345 -5.892172 -14.82003 -17.33053 -0.9168862 [2,] 12.957057 13.40696 10.157345 16.092172 12.82003 11.13053 16.9168862 [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] -15.609825 -4.798566 -10.31937 -3.703051 -9.163927 -8.149617 -11.21307 [2,] 3.409825 10.298566 11.61937 8.703051 17.563927 2.249617 15.01307 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [1,] -10.211542 -1.782726 -15.2374 -5.925473 -5.285302 -14.770321 -11.71803 [2,] 9.811542 14.682726 17.2374 16.925473 12.685302 5.070321 22.71803 [,36] [,37] [1,] -22.256865 -17.706007 [2,] 5.656865 9.706007 $out [1] 52.9 -23.2 -21.4 49.7 -32.9 -51.3 -28.6 20.9 -32.7 -28.7 60.5 -61.2 [13] -44.1 45.0 52.0 38.6 44.2 47.1 -45.0 -28.8 30.0 -63.6 -50.4 $group [1] 1 2 2 4 6 7 8 8 10 10 10 14 17 21 22 23 24 25 25 27 27 28 30 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" Warning message: In bxp(list(stats = c(-28.6, -7.2, 2.800, 14.2000000000000, 24, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/660z31272298206.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] [,10] [,11] [,12] [1,] 64.8 69.2 70.2 89.1 98.7 98.7 130.0 117.9 114.6 123.6 126.4 124.5 [2,] 77.1 79.1 90.3 100.2 117.3 129.3 145.2 151.0 147.3 141.9 145.1 144.1 [3,] 84.4 92.5 96.6 113.9 126.5 145.8 158.5 160.3 159.7 151.3 153.3 152.0 [4,] 96.3 97.7 113.0 127.5 137.3 159.2 169.4 175.6 172.1 166.5 172.2 172.0 [5,] 114.1 116.4 132.7 145.9 165.2 197.0 199.2 201.5 208.2 202.0 210.5 208.6 [,13] [1,] 128.0 [2,] 138.0 [3,] 143.0 [4,] 152.5 [5,] 154.0 $n [1] 37 37 37 37 37 37 37 37 37 37 37 37 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 79.4128 87.66864 90.70367 106.8088 121.305 138.0335 152.2140 153.9101 [2,] 89.3872 97.33136 102.49633 120.9912 131.695 153.5665 164.7860 166.6899 [,9] [,10] [,11] [,12] [,13] [1,] 153.2582 144.9101 146.2608 144.7530 136.3865 [2,] 166.1418 157.6899 160.3392 159.2470 149.6135 $out [1] 85.7 175.9 210.2 217.8 210.8 215.3 177.0 184.0 $group [1] 5 5 7 8 9 11 13 13 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" NA Warning message: In bxp(list(stats = c(64.8, 77.1, 84.4, 96.3, 114.1, 69.2, 79.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7z9yo1272298206.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,] 129.6923 124.5 124.300 [2,] 134.3083 135.0 131.300 [3,] 136.6667 138.4 134.650 [4,] 137.8750 144.8 140.025 [5,] 142.0667 152.9 144.750 $n [1] 37 37 37 $conf [,1] [,2] [,3] [1,] 135.7402 135.8544 132.3837 [2,] 137.5931 140.9456 136.9163 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > try(system("convert tmp/13zix1272298206.ps tmp/13zix1272298206.png",intern=TRUE)) character(0) > try(system("convert tmp/2w9hi1272298206.ps tmp/2w9hi1272298206.png",intern=TRUE)) character(0) > try(system("convert tmp/3w9hi1272298206.ps tmp/3w9hi1272298206.png",intern=TRUE)) character(0) > try(system("convert tmp/4w9hi1272298206.ps tmp/4w9hi1272298206.png",intern=TRUE)) character(0) > try(system("convert tmp/560z31272298206.ps tmp/560z31272298206.png",intern=TRUE)) character(0) > try(system("convert tmp/660z31272298206.ps tmp/660z31272298206.png",intern=TRUE)) character(0) > try(system("convert tmp/7z9yo1272298206.ps tmp/7z9yo1272298206.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.870 1.510 2.693