R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(80,111,122,131,192,188,216,238,173,160,93,67,60,32,126,131,134,162,230,232,200,143,85,66,54,81,100,126,204,218,227,220,220,120,110,67,81,52,106,156,187,204,204,196,204,124,53,77,77,50,105,125,165,194,263,225,263,140,127,86,71,95,95,133,178,160,250,251,250,173,103,21,29,39,71,148,144,199,206,224,206,152,88,35,23,92,117,120,173,202,217,256,217,143,95,77,76,100,108,132,195,198,204,212,204,129,73,77,80,64,109,138,185,198,237,223,237,146,102,77,70,86,98,141,195,205,191,226,191,147,100,74,56,77,80,120,186,196,229,229,229,176,104,61,72,99,113,140,174,209,205,229,215,136,113,57,55,66,125,149,176,230,238,245,238,124,111,72,63,78,100,149,166,201,214,231,214,151,97,68,81,55,99,146,170,218,218,207,218,178,105,67,47,55,73,124,185,213,278,205,278,171,125,92,96,92,118,185,183,215,207,214,207,142,102,66,87,90,90,133,205,201,220,210,220,136,95,52,40,60,100,169,184,202,226,239,226,149,121,50) > 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] 240 > (np <- floor(n / par1)) [1] 20 > 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] 20 20 20 20 20 20 20 20 20 20 20 20 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 80 60 54 81 77 71 29 23 76 80 70 56 72 [2,] 111 32 81 52 50 95 39 92 100 64 86 77 99 [3,] 122 126 100 106 105 95 71 117 108 109 98 80 113 [4,] 131 131 126 156 125 133 148 120 132 138 141 120 140 [5,] 192 134 204 187 165 178 144 173 195 185 195 186 174 [6,] 188 162 218 204 194 160 199 202 198 198 205 196 209 [7,] 216 230 227 204 263 250 206 217 204 237 191 229 205 [8,] 238 232 220 196 225 251 224 256 212 223 226 229 229 [9,] 173 200 220 204 263 250 206 217 204 237 191 229 215 [10,] 160 143 120 124 140 173 152 143 129 146 147 176 136 [11,] 93 85 110 53 127 103 88 95 73 102 100 104 113 [12,] 67 66 67 77 86 21 35 77 77 77 74 61 57 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] 55 63 81 47 96 87 40 NA [2,] 66 78 55 55 92 90 60 NA [3,] 125 100 99 73 118 90 100 NA [4,] 149 149 146 124 185 133 169 NA [5,] 176 166 170 185 183 205 184 NA [6,] 230 201 218 213 215 201 202 NA [7,] 238 214 218 278 207 220 226 NA [8,] 245 231 207 205 214 210 239 NA [9,] 238 214 218 278 207 220 226 NA [10,] 124 151 178 171 142 136 149 NA [11,] 111 97 105 125 102 95 121 NA [12,] 72 68 67 92 66 52 50 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 31 -28 27 -29 -27 24 10 69 24 -16 16 21 27 [2,] 11 94 19 54 55 0 32 25 8 45 12 3 14 [3,] 9 5 26 50 20 38 77 3 24 29 43 40 27 [4,] 61 3 78 31 40 45 -4 53 63 47 54 66 34 [5,] -4 28 14 17 29 -18 55 29 3 13 10 10 35 [6,] 28 68 9 0 69 90 7 15 6 39 -14 33 -4 [7,] 22 2 -7 -8 -38 1 18 39 8 -14 35 0 24 [8,] -65 -32 0 8 38 -1 -18 -39 -8 14 -35 0 -14 [9,] -13 -57 -100 -80 -123 -77 -54 -74 -75 -91 -44 -53 -79 [10,] -67 -58 -10 -71 -13 -70 -64 -48 -56 -44 -47 -72 -23 [11,] -26 -19 -43 24 -41 -82 -53 -18 4 -25 -26 -43 -56 [12,] -7 -12 14 0 -15 8 -12 -1 3 -7 -18 11 -2 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] 11 15 -26 8 -4 3 20 NA [2,] 59 22 44 18 26 0 40 NA [3,] 24 49 47 51 67 43 69 NA [4,] 27 17 24 61 -2 72 15 NA [5,] 54 35 48 28 32 -4 18 NA [6,] 8 13 0 65 -8 19 24 NA [7,] 7 17 -11 -73 7 -10 13 NA [8,] -7 -17 11 73 -7 10 -13 NA [9,] -114 -63 -40 -107 -65 -84 -77 NA [10,] -13 -54 -73 -46 -40 -41 -28 NA [11,] -39 -29 -38 -33 -36 -43 -71 NA [12,] -9 13 -20 4 21 -12 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/wessaorg/rcomp/tmp/1w10c1321546304.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/wessaorg/rcomp/tmp/20szi1321546304.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/wessaorg/rcomp/tmp/31c9b1321546304.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/wessaorg/rcomp/tmp/44yx81321546304.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] [,11] [,12] [1,] 23.0 32.0 71.0 120.0 165.0 188.0 191.0 196.0 173.0 120.0 73.0 35 [2,] 54.5 55.0 96.5 128.5 171.5 197.0 206.5 213.0 205.0 136.0 94.0 59 [3,] 70.5 77.5 102.5 135.5 183.5 201.5 219.0 225.5 217.5 144.5 102.0 67 [4,] 80.0 92.0 115.0 148.5 189.5 211.0 233.5 235.0 233.0 156.0 110.5 77 [5,] 96.0 111.0 126.0 169.0 205.0 230.0 263.0 256.0 263.0 178.0 127.0 92 $n [1] 20 20 20 20 20 20 20 20 20 20 20 20 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 61.49088 64.42795 95.96397 128.434 177.1406 196.5538 209.4609 217.7274 [2,] 79.50912 90.57205 109.03603 142.566 189.8594 206.4462 228.5391 233.2726 [,9] [,10] [,11] [,12] [1,] 207.6076 137.434 96.17057 60.64062 [2,] 227.3924 151.566 107.82943 73.35938 $out [1] 185 134 144 162 160 278 278 53 21 $group [1] 4 5 5 6 6 7 9 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(23, 54.5, 70.5, 80, 96, 32, 55, 77.5, 92, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5pvet1321546304.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] [,12] [1,] -29 0.0 3.0 -4.0 -18.0 -14 -38.0 -39.0 -123.0 -73.0 -56.0 -20 [2,] -10 11.5 24.0 20.5 10.0 3 -9.0 -17.5 -87.5 -65.5 -43.0 -12 [3,] 13 23.5 39.0 42.5 23.0 14 4.5 -7.0 -76.0 -47.5 -37.0 -2 [4,] 24 44.5 49.5 61.0 33.5 36 17.5 9.0 -55.5 -34.0 -25.5 6 [5,] 69 94.0 77.0 78.0 55.0 69 39.0 38.0 -13.0 -10.0 -18.0 21 $n [1] 20 20 20 20 20 20 20 20 20 20 20 19 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.9878428 11.84114 29.99088 28.1914 14.69748 2.341142 -4.862417 [2,] 25.0121572 35.15886 48.00912 56.8086 31.30252 25.658858 13.862417 [,8] [,9] [,10] [,11] [,12] [1,] -16.362417 -87.30556 -58.62891 -43.18273 -8.524583 [2,] 2.362417 -64.69444 -36.37109 -30.81727 4.524583 $out [1] 90 -73 -65 73 24 -82 4 -71 $group [1] 6 7 8 8 11 11 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-29, -10, 13, 24, 69, 0, 11.5, 23.5, 44.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6xy191321546304.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] [,10] [,11] [,12] [1,] 67.0 32.0 54.0 52 50.0 21.0 29.0 23.0 73.0 64.0 70 56.0 [2,] 102.0 75.5 90.5 79 95.5 95.0 55.0 93.5 88.5 91.0 92 78.5 [3,] 145.5 132.5 123.0 140 133.5 146.5 146.0 131.5 130.5 142.0 144 148.0 [4,] 190.0 181.0 219.0 200 209.5 214.0 202.5 209.5 201.0 210.5 193 212.5 [5,] 238.0 232.0 227.0 204 263.0 251.0 224.0 256.0 212.0 237.0 226 229.0 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] 57 55.0 63.0 55.0 47.0 66.0 52.0 40 NA [2,] 106 91.5 87.5 90.0 82.5 99.0 90.0 80 NA [3,] 138 137.0 150.0 158.0 148.0 162.5 134.5 159 NA [4,] 207 234.0 207.5 212.5 209.0 207.0 207.5 214 NA [5,] 229 245.0 231.0 218.0 278.0 215.0 220.0 239 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 105.3626 84.38074 64.39029 84.81109 81.50383 92.2233 78.72426 [2,] 185.6374 180.61926 181.60971 195.18891 185.49617 200.7767 213.27574 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 78.59162 79.18799 87.49525 97.93322 86.8817 91.93322 72.00479 [2,] 184.40838 181.81201 196.50475 190.06678 209.1183 184.06678 201.99521 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] 95.26719 102.1269 90.3025 113.2405 80.90746 97.8817 NA [2,] 204.73281 213.8731 205.6975 211.7595 188.09254 220.1183 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" NA Warning message: In bxp(list(stats = c(67, 102, 145.5, 190, 238, 32, 75.5, 132.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7mn8x1321546304.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,] 64.900 67.00 67.3750 [2,] 86.900 89.75 87.9375 [3,] 143.400 140.00 142.0000 [4,] 210.575 209.50 211.0000 [5,] 225.600 225.50 223.5000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 86.991 85.38122 85.87037 [2,] 199.809 194.61878 198.12963 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(64.9, 86.9, 143.4, 210.575, 225.6, 67, 89.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1w10c1321546304.ps tmp/1w10c1321546304.png",intern=TRUE)) character(0) > try(system("convert tmp/20szi1321546304.ps tmp/20szi1321546304.png",intern=TRUE)) character(0) > try(system("convert tmp/31c9b1321546304.ps tmp/31c9b1321546304.png",intern=TRUE)) character(0) > try(system("convert tmp/44yx81321546304.ps tmp/44yx81321546304.png",intern=TRUE)) character(0) > try(system("convert tmp/5pvet1321546304.ps tmp/5pvet1321546304.png",intern=TRUE)) character(0) > try(system("convert tmp/6xy191321546304.ps tmp/6xy191321546304.png",intern=TRUE)) character(0) > try(system("convert tmp/7mn8x1321546304.ps tmp/7mn8x1321546304.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.585 0.438 2.050