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(153.3 + ,154.5 + ,155.2 + ,156.9 + ,157 + ,157.4 + ,157.2 + ,157.5 + ,158 + ,158.5 + ,159 + ,159.3 + ,160 + ,160.8 + ,161.9 + ,162.5 + ,162.7 + ,162.8 + ,162.9 + ,163 + ,164 + ,164.7 + ,164.8 + ,164.9 + ,165 + ,165.8 + ,166.1 + ,167.2 + ,167.7 + ,168.3 + ,168.6 + ,168.9 + ,169.1 + ,169.5 + ,169.6 + ,169.7 + ,169.8 + ,170.4 + ,170.9 + ,171.9 + ,171.9 + ,172 + ,172 + ,172.4 + ,173 + ,173.7 + ,173.8 + ,173.8 + ,173.9 + ,174.6 + ,175 + ,175.9 + ,176 + ,175.1 + ,175.6 + ,175.9 + ,176.7 + ,176.1 + ,176.1 + ,176.2 + ,176.3 + ,177.8 + ,178.5 + ,179.4 + ,179.5 + ,179.6 + ,179.7 + ,179.7 + ,179.8 + ,179.9 + ,180.2 + ,180.4 + ,180.4 + ,181.3 + ,181.9 + ,182.5 + ,182.7 + ,183.1 + ,183.6 + ,183.7 + ,183.8 + ,183.9 + ,184.1 + ,184.4 + ,184.5 + ,185.9 + ,186.6 + ,187.6 + ,187.8 + ,187.9 + ,188 + ,188.3 + ,188.4 + ,188.5 + ,188.5 + ,188.6 + ,188.6 + ,189.4 + ,190 + ,191.9 + ,192.5 + ,193 + ,193.5 + ,193.9 + ,194.2 + ,194.9 + ,194.9 + ,194.9 + ,194.9 + ,195.5 + ,196 + ,196.2 + ,196.2 + ,196.2 + ,196.2 + ,197 + ,197.7 + ,198 + ,198.2 + ,198.5 + ,198.6 + ,199.5 + ,200 + ,201.3 + ,202.2 + ,202.9 + ,203.5 + ,203.5 + ,204 + ,204.1 + ,204.3 + ,204.5 + ,204.8 + ,205.1 + ,205.7 + ,206.5 + ,206.9 + ,207.1 + ,207.8 + ,208 + ,208.5 + ,208.6 + ,209 + ,209.1 + ,209.7 + ,209.8 + ,209.9 + ,210 + ,210.8 + ,211.4 + ,211.7 + ,212 + ,212.2 + ,212.4 + ,212.9 + ,213.4 + ,213.7 + ,214 + ,214.3 + ,214.8 + ,215 + ,215.9 + ,216.4 + ,216.9 + ,217.2 + ,217.5 + ,217.9 + ,218.1 + ,218.6 + ,218.9 + ,219.3 + ,220.4 + ,220.9 + ,221 + ,221.8 + ,222 + ,222.2 + ,222.5 + ,222.9 + ,223.1 + ,223.4 + ,224 + ,225.1 + ,225.5 + ,225.9 + ,226.3 + ,226.5 + ,227 + ,227.3 + ,227.8 + ,228.1 + ,228.4 + ,228.5 + ,228.8 + ,229 + ,229.1 + ,229.3 + ,229.6 + ,229.9 + ,230 + ,230.2 + ,230.8 + ,231 + ,231.7 + ,231.9 + ,233 + ,235.1 + ,236 + ,236.9 + ,237.1 + ,237.5 + ,238.2 + ,238.9 + ,239.1 + ,240 + ,240.2 + ,240.5 + ,240.7 + ,241.1 + ,241.4 + ,242.2 + ,242.9 + ,243.2 + ,243.9) > par20 = '' > par19 = '' > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '' > par12 = '' > par11 = '' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '' > par4 = '' > par3 = '' > par2 = '' > par1 = '12' > ylab = '' > xlab = '' > main = '' > #'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] 224 > (np <- floor(n / par1)) [1] 18 > 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] 19 19 19 19 19 19 19 19 18 18 18 18 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 153.3 160.0 165.0 169.8 173.9 176.3 180.4 184.5 188.6 194.9 198.6 204.8 [2,] 154.5 160.8 165.8 170.4 174.6 177.8 181.3 185.9 189.4 195.5 199.5 205.1 [3,] 155.2 161.9 166.1 170.9 175.0 178.5 181.9 186.6 190.0 196.0 200.0 205.7 [4,] 156.9 162.5 167.2 171.9 175.9 179.4 182.5 187.6 191.9 196.2 201.3 206.5 [5,] 157.0 162.7 167.7 171.9 176.0 179.5 182.7 187.8 192.5 196.2 202.2 206.9 [6,] 157.4 162.8 168.3 172.0 175.1 179.6 183.1 187.9 193.0 196.2 202.9 207.1 [7,] 157.2 162.9 168.6 172.0 175.6 179.7 183.6 188.0 193.5 196.2 203.5 207.8 [8,] 157.5 163.0 168.9 172.4 175.9 179.7 183.7 188.3 193.9 197.0 203.5 208.0 [9,] 158.0 164.0 169.1 173.0 176.7 179.8 183.8 188.4 194.2 197.7 204.0 208.5 [10,] 158.5 164.7 169.5 173.7 176.1 179.9 183.9 188.5 194.9 198.0 204.1 208.6 [11,] 159.0 164.8 169.6 173.8 176.1 180.2 184.1 188.5 194.9 198.2 204.3 209.0 [12,] 159.3 164.9 169.7 173.8 176.2 180.4 184.4 188.6 194.9 198.5 204.5 209.1 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [1,] 209.7 213.7 218.6 223.4 228.5 231.9 240.5 [2,] 209.8 214.0 218.9 224.0 228.8 233.0 240.7 [3,] 209.9 214.3 219.3 225.1 229.0 235.1 241.1 [4,] 210.0 214.8 220.4 225.5 229.1 236.0 241.4 [5,] 210.8 215.0 220.9 225.9 229.3 236.9 242.2 [6,] 211.4 215.9 221.0 226.3 229.6 237.1 242.9 [7,] 211.7 216.4 221.8 226.5 229.9 237.5 243.2 [8,] 212.0 216.9 222.0 227.0 230.0 238.2 243.9 [9,] 212.2 217.2 222.2 227.3 230.2 238.9 NA [10,] 212.4 217.5 222.5 227.8 230.8 239.1 NA [11,] 212.9 217.9 222.9 228.1 231.0 240.0 NA [12,] 213.4 218.1 223.1 228.4 231.7 240.2 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 1.2 0.8 0.8 0.6 0.7 1.5 0.9 1.4 0.8 0.6 0.9 0.3 0.1 [2,] 0.7 1.1 0.3 0.5 0.4 0.7 0.6 0.7 0.6 0.5 0.5 0.6 0.1 [3,] 1.7 0.6 1.1 1.0 0.9 0.9 0.6 1.0 1.9 0.2 1.3 0.8 0.1 [4,] 0.1 0.2 0.5 0.0 0.1 0.1 0.2 0.2 0.6 0.0 0.9 0.4 0.8 [5,] 0.4 0.1 0.6 0.1 -0.9 0.1 0.4 0.1 0.5 0.0 0.7 0.2 0.6 [6,] -0.2 0.1 0.3 0.0 0.5 0.1 0.5 0.1 0.5 0.0 0.6 0.7 0.3 [7,] 0.3 0.1 0.3 0.4 0.3 0.0 0.1 0.3 0.4 0.8 0.0 0.2 0.3 [8,] 0.5 1.0 0.2 0.6 0.8 0.1 0.1 0.1 0.3 0.7 0.5 0.5 0.2 [9,] 0.5 0.7 0.4 0.7 -0.6 0.1 0.1 0.1 0.7 0.3 0.1 0.1 0.2 [10,] 0.5 0.1 0.1 0.1 0.0 0.3 0.2 0.0 0.0 0.2 0.2 0.4 0.5 [11,] 0.3 0.1 0.1 0.0 0.1 0.2 0.3 0.1 0.0 0.3 0.2 0.1 0.5 [12,] 0.7 0.1 0.1 0.1 0.1 0.0 0.1 0.0 0.0 0.1 0.3 0.6 0.3 [,14] [,15] [,16] [,17] [,18] [,19] [1,] 0.3 0.3 0.6 0.3 1.1 0.2 [2,] 0.3 0.4 1.1 0.2 2.1 0.4 [3,] 0.5 1.1 0.4 0.1 0.9 0.3 [4,] 0.2 0.5 0.4 0.2 0.9 0.8 [5,] 0.9 0.1 0.4 0.3 0.2 0.7 [6,] 0.5 0.8 0.2 0.3 0.4 0.3 [7,] 0.5 0.2 0.5 0.1 0.7 0.7 [8,] 0.3 0.2 0.3 0.2 0.7 NA [9,] 0.3 0.3 0.5 0.6 0.2 NA [10,] 0.4 0.4 0.3 0.2 0.9 NA [11,] 0.2 0.2 0.3 0.7 0.2 NA [12,] 0.5 0.3 0.1 0.2 0.3 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/1n94i1257935572.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/23ebi1257935572.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/3kab81257935572.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/4fm8k1257935572.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] [1,] 153.30 154.50 155.20 156.90 157.00 157.40 157.20 157.50 158.00 158.50 [2,] 175.10 176.20 176.75 177.65 177.75 177.35 177.65 177.80 176.70 176.10 [3,] 194.90 195.50 196.00 196.20 196.20 196.20 196.20 197.00 195.95 196.45 [4,] 216.15 216.45 216.80 217.60 217.95 218.45 219.10 219.45 217.20 217.50 [5,] 240.50 240.70 241.10 241.40 242.20 242.90 243.20 243.90 238.90 239.10 [,11] [,12] [1,] 159.00 159.3 [2,] 176.10 176.2 [3,] 196.55 196.7 [4,] 217.90 218.1 [5,] 240.00 240.2 $n [1] 19 19 19 19 19 19 19 19 18 18 18 18 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 180.0203 180.9103 181.4828 181.7190 181.6284 181.3022 181.1753 181.9028 [2,] 209.7797 210.0897 210.5172 210.6810 210.7716 211.0978 211.2247 212.0972 [,9] [,10] [,11] [,12] [1,] 180.8674 181.0322 180.9833 181.0960 [2,] 211.0326 211.8678 212.1167 212.3040 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5bn761257935572.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,] 0.1 0.1 0.10 0.00 0.00 -0.2 0.00 0.1 0.1 0.0 0.0 0.0 [2,] 0.3 0.4 0.45 0.15 0.10 0.1 0.15 0.2 0.1 0.1 0.1 0.1 [3,] 0.7 0.5 0.90 0.20 0.30 0.3 0.30 0.3 0.3 0.2 0.2 0.1 [4,] 0.9 0.7 1.05 0.55 0.55 0.5 0.45 0.6 0.5 0.4 0.3 0.3 [5,] 1.5 1.1 1.90 0.90 0.90 0.8 0.80 1.0 0.7 0.5 0.5 0.6 $n [1] 19 19 19 19 19 19 19 18 18 18 18 18 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.4825139 0.3912569 0.6825139 0.05500926 0.1368854 0.1550093 0.1912569 [2,] 0.9174861 0.6087431 1.1174861 0.34499074 0.4631146 0.4449907 0.4087431 [,8] [,9] [,10] [,11] [,12] [1,] 0.1510362 0.1510362 0.08827713 0.1255181 0.02551809 [2,] 0.4489638 0.4489638 0.31172287 0.2744819 0.17448191 $out [1] 2.1 -0.9 -0.6 0.9 0.7 0.7 $group [1] 2 5 9 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(0.100000000000023, 0.300000000000011, 0.699999999999989, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6dawh1257935572.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] [1,] 153.30 160.00 165.00 169.80 173.90 177.80 180.40 185.90 188.60 194.90 [2,] 156.05 162.20 166.65 171.40 175.05 178.95 182.20 187.10 190.95 196.10 [3,] 157.30 162.85 168.45 172.00 175.90 179.65 183.35 187.95 193.25 196.20 [4,] 158.25 164.35 169.30 173.35 176.10 179.85 183.85 188.45 194.55 197.85 [5,] 159.30 164.90 169.70 173.80 176.70 180.40 184.40 188.60 194.90 198.50 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [1,] 198.60 204.80 209.70 213.70 218.60 223.40 228.50 231.90 240.50 [2,] 200.65 206.10 209.95 214.55 219.85 225.30 229.05 235.55 240.90 [3,] 203.20 207.45 211.55 216.15 221.40 226.40 229.75 237.30 241.80 [4,] 204.05 208.55 212.30 217.35 222.35 227.55 230.50 239.00 243.05 [5,] 204.50 209.10 213.40 218.10 223.10 228.40 231.70 240.20 243.90 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 156.2966 161.8694 167.2413 171.1106 175.4211 179.2395 182.5974 187.3343 [2,] 158.3034 163.8306 169.6587 172.8894 176.3789 180.0605 184.1026 188.5657 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 191.608 195.4018 201.6492 206.3325 210.4781 214.8729 220.2597 225.3738 [2,] 194.892 196.9982 204.7508 208.5675 212.6219 217.4271 222.5403 227.4262 [,17] [,18] [,19] [1,] 229.0886 235.7264 240.599 [2,] 230.4114 238.8736 243.001 $out [1] 176.3 184.5 $group [1] 6 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" NA Warning message: In bxp(list(stats = c(153.3, 156.05, 157.3, 158.25, 159.3, 160, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7txrr1257935572.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,] 195.6000 195.500 195.6250 [2,] 196.9409 195.975 196.6750 [3,] 197.6250 196.200 196.9875 [4,] 198.2553 196.500 197.8750 [5,] 199.0421 197.000 198.6250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 197.0255 195.9605 196.4402 [2,] 198.2245 196.4395 197.5348 $out [1] 194.9 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(195.6, 196.940935672515, 197.625, 198.255263157895, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1n94i1257935572.ps tmp/1n94i1257935572.png") > system("convert tmp/23ebi1257935572.ps tmp/23ebi1257935572.png") > system("convert tmp/3kab81257935572.ps tmp/3kab81257935572.png") > system("convert tmp/4fm8k1257935572.ps tmp/4fm8k1257935572.png") > system("convert tmp/5bn761257935572.ps tmp/5bn761257935572.png") > system("convert tmp/6dawh1257935572.ps tmp/6dawh1257935572.png") > system("convert tmp/7txrr1257935572.ps tmp/7txrr1257935572.png") > > > proc.time() user system elapsed 1.200 0.938 1.529