R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-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(426.113 + ,383.703 + ,232.444 + ,70.939 + ,226.731 + ,947.293 + ,611.281 + ,158.047 + ,33.999 + ,37.028 + ,388.3 + ,506.652 + ,392.25 + ,180.818 + ,198.296 + ,217.465 + ,275.562 + ,1030.944 + ,57.47 + ,136.452 + ,556.277 + ,213.361 + ,274.482 + ,220.553 + ,236.71 + ,260.642 + ,2763.544 + ,213.923 + ,169.861 + ,403.064 + ,449.594 + ,406.167 + ,206.893 + ,156.187 + ,257.102 + ,62.156 + ,662.883 + ,251.422 + ,171.328 + ,350.089 + ,221.588 + ,4.813 + ,183.186 + ,190.379 + ,223.166 + ,232.669 + ,356.725 + ,109.215 + ,475.834 + ,315.955 + ,694.87 + ,8.95 + ,278.741 + ,308.16 + ,207.533 + ,192.797 + ,601.162 + ,289.714 + ,293.671 + ,386.688 + ,699.645 + ,85.094 + ,131.812 + ,645.285 + ,197.549 + ,308.174 + ,86.58 + ,242.205 + ,238.502 + ,187.881 + ,140.321 + ,440.31 + ,421.403 + ,218.761 + ,1305.923 + ,137.55 + ,262.517 + ,348.821 + ,150.034 + ,64.016 + ,261.596 + ,259.7 + ,171.26 + ,203.077 + ,249.148 + ,211.655 + ,252.64 + ,438.555 + ,239.89 + ,401.915 + ,216.886 + ,184.641 + ,380.155 + ,653.641 + ,313.906 + ,366.936 + ,236.302 + ,229.641 + ,235.577 + ,103.898 + ,263.906 + ,241.171 + ,216.548 + ,295.281 + ,193.299 + ,204.386 + ,257.567 + ,136.813 + ,240.755 + ,59.609 + ,213.511 + ,380.531 + ,242.344 + ,250.407 + ,183.613 + ,191.835 + ,266.793 + ,246.542 + ,330.563 + ,403.556 + ,208.108 + ,324.04 + ,308.532 + ,199.297 + ,200.156 + ,262.875 + ,287.069 + ,190.157 + ,199.746 + ,265.777 + ,435.956 + ,72.844 + ,756.46 + ,206.771 + ,4202.361 + ,401.422 + ,216.046 + ,39.047 + ,441.437) > par1 = '1' > par1 <- '1' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Wessa P., (2012), Standard Deviation Plot (v1.0.2) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_sdplot.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > par1 <- as.numeric(par1) > (n <- length(x)) [1] 139 > (np <- floor(n / par1)) [1] 139 > arr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + if (j == par1) j = 0 + } > ari [1] 139 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 426.113 383.703 232.444 70.939 226.731 947.293 611.281 158.047 33.999 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [1,] 37.028 388.3 506.652 392.25 180.818 198.296 217.465 275.562 1030.944 57.47 [,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] 136.452 556.277 213.361 274.482 220.553 236.71 260.642 2763.544 213.923 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 169.861 403.064 449.594 406.167 206.893 156.187 257.102 62.156 662.883 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [1,] 251.422 171.328 350.089 221.588 4.813 183.186 190.379 223.166 232.669 [,47] [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [1,] 356.725 109.215 475.834 315.955 694.87 8.95 278.741 308.16 207.533 [,56] [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [1,] 192.797 601.162 289.714 293.671 386.688 699.645 85.094 131.812 645.285 [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 197.549 308.174 86.58 242.205 238.502 187.881 140.321 440.31 421.403 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [1,] 218.761 1305.923 137.55 262.517 348.821 150.034 64.016 261.596 259.7 [,83] [,84] [,85] [,86] [,87] [,88] [,89] [,90] [,91] [1,] 171.26 203.077 249.148 211.655 252.64 438.555 239.89 401.915 216.886 [,92] [,93] [,94] [,95] [,96] [,97] [,98] [,99] [,100] [1,] 184.641 380.155 653.641 313.906 366.936 236.302 229.641 235.577 103.898 [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [,109] [1,] 263.906 241.171 216.548 295.281 193.299 204.386 257.567 136.813 240.755 [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 59.609 213.511 380.531 242.344 250.407 183.613 191.835 266.793 246.542 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [1,] 330.563 403.556 208.108 324.04 308.532 199.297 200.156 262.875 287.069 [,128] [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [1,] 190.157 199.746 265.777 435.956 72.844 756.46 206.771 4202.361 401.422 [,137] [,138] [,139] [,140] [1,] 216.046 39.047 441.437 NA > arr.sd <- array(NA,dim=par1) > arr.range <- array(NA,dim=par1) > arr.iqr <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) + arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) + arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) + } > overall.sd <- sqrt(var(x)) > overall.range <- max(x) - min(x) > overall.iqr <- quantile(x,0.75) - quantile(x,0.25) > postscript(file="/var/wessaorg/rcomp/tmp/1p1411354379682.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.sd,0) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/22tvu1354379682.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.range,0) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/33uc01354379682.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.iqr,0) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/44eew1354379682.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] [1,] 4.8130 [2,] 193.0480 [3,] 241.1710 [4,] 361.8305 [5,] 611.2810 $n [1] 139 $conf [,1] [1,] 218.5518 [2,] 263.7902 $out [1] 947.293 1030.944 2763.544 662.883 694.870 699.645 645.285 1305.923 [9] 653.641 756.460 4202.361 $group [1] 1 1 1 1 1 1 1 1 1 1 1 $names [1] "1" > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5ktbo1354379682.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,] 426.113 383.703 232.444 70.939 226.731 947.293 611.281 158.047 33.999 [2,] 426.113 383.703 232.444 70.939 226.731 947.293 611.281 158.047 33.999 [3,] 426.113 383.703 232.444 70.939 226.731 947.293 611.281 158.047 33.999 [4,] 426.113 383.703 232.444 70.939 226.731 947.293 611.281 158.047 33.999 [5,] 426.113 383.703 232.444 70.939 226.731 947.293 611.281 158.047 33.999 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [1,] 37.028 388.3 506.652 392.25 180.818 198.296 217.465 275.562 1030.944 57.47 [2,] 37.028 388.3 506.652 392.25 180.818 198.296 217.465 275.562 1030.944 57.47 [3,] 37.028 388.3 506.652 392.25 180.818 198.296 217.465 275.562 1030.944 57.47 [4,] 37.028 388.3 506.652 392.25 180.818 198.296 217.465 275.562 1030.944 57.47 [5,] 37.028 388.3 506.652 392.25 180.818 198.296 217.465 275.562 1030.944 57.47 [,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] 136.452 556.277 213.361 274.482 220.553 236.71 260.642 2763.544 213.923 [2,] 136.452 556.277 213.361 274.482 220.553 236.71 260.642 2763.544 213.923 [3,] 136.452 556.277 213.361 274.482 220.553 236.71 260.642 2763.544 213.923 [4,] 136.452 556.277 213.361 274.482 220.553 236.71 260.642 2763.544 213.923 [5,] 136.452 556.277 213.361 274.482 220.553 236.71 260.642 2763.544 213.923 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 169.861 403.064 449.594 406.167 206.893 156.187 257.102 62.156 662.883 [2,] 169.861 403.064 449.594 406.167 206.893 156.187 257.102 62.156 662.883 [3,] 169.861 403.064 449.594 406.167 206.893 156.187 257.102 62.156 662.883 [4,] 169.861 403.064 449.594 406.167 206.893 156.187 257.102 62.156 662.883 [5,] 169.861 403.064 449.594 406.167 206.893 156.187 257.102 62.156 662.883 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [1,] 251.422 171.328 350.089 221.588 4.813 183.186 190.379 223.166 232.669 [2,] 251.422 171.328 350.089 221.588 4.813 183.186 190.379 223.166 232.669 [3,] 251.422 171.328 350.089 221.588 4.813 183.186 190.379 223.166 232.669 [4,] 251.422 171.328 350.089 221.588 4.813 183.186 190.379 223.166 232.669 [5,] 251.422 171.328 350.089 221.588 4.813 183.186 190.379 223.166 232.669 [,47] [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [1,] 356.725 109.215 475.834 315.955 694.87 8.95 278.741 308.16 207.533 [2,] 356.725 109.215 475.834 315.955 694.87 8.95 278.741 308.16 207.533 [3,] 356.725 109.215 475.834 315.955 694.87 8.95 278.741 308.16 207.533 [4,] 356.725 109.215 475.834 315.955 694.87 8.95 278.741 308.16 207.533 [5,] 356.725 109.215 475.834 315.955 694.87 8.95 278.741 308.16 207.533 [,56] [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [1,] 192.797 601.162 289.714 293.671 386.688 699.645 85.094 131.812 645.285 [2,] 192.797 601.162 289.714 293.671 386.688 699.645 85.094 131.812 645.285 [3,] 192.797 601.162 289.714 293.671 386.688 699.645 85.094 131.812 645.285 [4,] 192.797 601.162 289.714 293.671 386.688 699.645 85.094 131.812 645.285 [5,] 192.797 601.162 289.714 293.671 386.688 699.645 85.094 131.812 645.285 [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 197.549 308.174 86.58 242.205 238.502 187.881 140.321 440.31 421.403 [2,] 197.549 308.174 86.58 242.205 238.502 187.881 140.321 440.31 421.403 [3,] 197.549 308.174 86.58 242.205 238.502 187.881 140.321 440.31 421.403 [4,] 197.549 308.174 86.58 242.205 238.502 187.881 140.321 440.31 421.403 [5,] 197.549 308.174 86.58 242.205 238.502 187.881 140.321 440.31 421.403 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [1,] 218.761 1305.923 137.55 262.517 348.821 150.034 64.016 261.596 259.7 [2,] 218.761 1305.923 137.55 262.517 348.821 150.034 64.016 261.596 259.7 [3,] 218.761 1305.923 137.55 262.517 348.821 150.034 64.016 261.596 259.7 [4,] 218.761 1305.923 137.55 262.517 348.821 150.034 64.016 261.596 259.7 [5,] 218.761 1305.923 137.55 262.517 348.821 150.034 64.016 261.596 259.7 [,83] [,84] [,85] [,86] [,87] [,88] [,89] [,90] [,91] [1,] 171.26 203.077 249.148 211.655 252.64 438.555 239.89 401.915 216.886 [2,] 171.26 203.077 249.148 211.655 252.64 438.555 239.89 401.915 216.886 [3,] 171.26 203.077 249.148 211.655 252.64 438.555 239.89 401.915 216.886 [4,] 171.26 203.077 249.148 211.655 252.64 438.555 239.89 401.915 216.886 [5,] 171.26 203.077 249.148 211.655 252.64 438.555 239.89 401.915 216.886 [,92] [,93] [,94] [,95] [,96] [,97] [,98] [,99] [,100] [1,] 184.641 380.155 653.641 313.906 366.936 236.302 229.641 235.577 103.898 [2,] 184.641 380.155 653.641 313.906 366.936 236.302 229.641 235.577 103.898 [3,] 184.641 380.155 653.641 313.906 366.936 236.302 229.641 235.577 103.898 [4,] 184.641 380.155 653.641 313.906 366.936 236.302 229.641 235.577 103.898 [5,] 184.641 380.155 653.641 313.906 366.936 236.302 229.641 235.577 103.898 [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [,109] [1,] 263.906 241.171 216.548 295.281 193.299 204.386 257.567 136.813 240.755 [2,] 263.906 241.171 216.548 295.281 193.299 204.386 257.567 136.813 240.755 [3,] 263.906 241.171 216.548 295.281 193.299 204.386 257.567 136.813 240.755 [4,] 263.906 241.171 216.548 295.281 193.299 204.386 257.567 136.813 240.755 [5,] 263.906 241.171 216.548 295.281 193.299 204.386 257.567 136.813 240.755 [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 59.609 213.511 380.531 242.344 250.407 183.613 191.835 266.793 246.542 [2,] 59.609 213.511 380.531 242.344 250.407 183.613 191.835 266.793 246.542 [3,] 59.609 213.511 380.531 242.344 250.407 183.613 191.835 266.793 246.542 [4,] 59.609 213.511 380.531 242.344 250.407 183.613 191.835 266.793 246.542 [5,] 59.609 213.511 380.531 242.344 250.407 183.613 191.835 266.793 246.542 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [1,] 330.563 403.556 208.108 324.04 308.532 199.297 200.156 262.875 287.069 [2,] 330.563 403.556 208.108 324.04 308.532 199.297 200.156 262.875 287.069 [3,] 330.563 403.556 208.108 324.04 308.532 199.297 200.156 262.875 287.069 [4,] 330.563 403.556 208.108 324.04 308.532 199.297 200.156 262.875 287.069 [5,] 330.563 403.556 208.108 324.04 308.532 199.297 200.156 262.875 287.069 [,128] [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [1,] 190.157 199.746 265.777 435.956 72.844 756.46 206.771 4202.361 401.422 [2,] 190.157 199.746 265.777 435.956 72.844 756.46 206.771 4202.361 401.422 [3,] 190.157 199.746 265.777 435.956 72.844 756.46 206.771 4202.361 401.422 [4,] 190.157 199.746 265.777 435.956 72.844 756.46 206.771 4202.361 401.422 [5,] 190.157 199.746 265.777 435.956 72.844 756.46 206.771 4202.361 401.422 [,137] [,138] [,139] [,140] [1,] 216.046 39.047 441.437 NA [2,] 216.046 39.047 441.437 NA [3,] 216.046 39.047 441.437 NA [4,] 216.046 39.047 441.437 NA [5,] 216.046 39.047 441.437 NA $n [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [38] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [75] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [112] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 426.113 383.703 232.444 70.939 226.731 947.293 611.281 158.047 33.999 [2,] 426.113 383.703 232.444 70.939 226.731 947.293 611.281 158.047 33.999 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [1,] 37.028 388.3 506.652 392.25 180.818 198.296 217.465 275.562 1030.944 57.47 [2,] 37.028 388.3 506.652 392.25 180.818 198.296 217.465 275.562 1030.944 57.47 [,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] 136.452 556.277 213.361 274.482 220.553 236.71 260.642 2763.544 213.923 [2,] 136.452 556.277 213.361 274.482 220.553 236.71 260.642 2763.544 213.923 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 169.861 403.064 449.594 406.167 206.893 156.187 257.102 62.156 662.883 [2,] 169.861 403.064 449.594 406.167 206.893 156.187 257.102 62.156 662.883 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [1,] 251.422 171.328 350.089 221.588 4.813 183.186 190.379 223.166 232.669 [2,] 251.422 171.328 350.089 221.588 4.813 183.186 190.379 223.166 232.669 [,47] [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [1,] 356.725 109.215 475.834 315.955 694.87 8.95 278.741 308.16 207.533 [2,] 356.725 109.215 475.834 315.955 694.87 8.95 278.741 308.16 207.533 [,56] [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [1,] 192.797 601.162 289.714 293.671 386.688 699.645 85.094 131.812 645.285 [2,] 192.797 601.162 289.714 293.671 386.688 699.645 85.094 131.812 645.285 [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 197.549 308.174 86.58 242.205 238.502 187.881 140.321 440.31 421.403 [2,] 197.549 308.174 86.58 242.205 238.502 187.881 140.321 440.31 421.403 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [1,] 218.761 1305.923 137.55 262.517 348.821 150.034 64.016 261.596 259.7 [2,] 218.761 1305.923 137.55 262.517 348.821 150.034 64.016 261.596 259.7 [,83] [,84] [,85] [,86] [,87] [,88] [,89] [,90] [,91] [1,] 171.26 203.077 249.148 211.655 252.64 438.555 239.89 401.915 216.886 [2,] 171.26 203.077 249.148 211.655 252.64 438.555 239.89 401.915 216.886 [,92] [,93] [,94] [,95] [,96] [,97] [,98] [,99] [,100] [1,] 184.641 380.155 653.641 313.906 366.936 236.302 229.641 235.577 103.898 [2,] 184.641 380.155 653.641 313.906 366.936 236.302 229.641 235.577 103.898 [,101] [,102] [,103] [,104] [,105] [,106] [,107] [,108] [,109] [1,] 263.906 241.171 216.548 295.281 193.299 204.386 257.567 136.813 240.755 [2,] 263.906 241.171 216.548 295.281 193.299 204.386 257.567 136.813 240.755 [,110] [,111] [,112] [,113] [,114] [,115] [,116] [,117] [,118] [1,] 59.609 213.511 380.531 242.344 250.407 183.613 191.835 266.793 246.542 [2,] 59.609 213.511 380.531 242.344 250.407 183.613 191.835 266.793 246.542 [,119] [,120] [,121] [,122] [,123] [,124] [,125] [,126] [,127] [1,] 330.563 403.556 208.108 324.04 308.532 199.297 200.156 262.875 287.069 [2,] 330.563 403.556 208.108 324.04 308.532 199.297 200.156 262.875 287.069 [,128] [,129] [,130] [,131] [,132] [,133] [,134] [,135] [,136] [1,] 190.157 199.746 265.777 435.956 72.844 756.46 206.771 4202.361 401.422 [2,] 190.157 199.746 265.777 435.956 72.844 756.46 206.771 4202.361 401.422 [,137] [,138] [,139] [,140] [1,] 216.046 39.047 441.437 NA [2,] 216.046 39.047 441.437 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" [13] "13" "14" "15" "16" "17" "18" "19" "20" "21" "22" "23" "24" [25] "25" "26" "27" "28" "29" "30" "31" "32" "33" "34" "35" "36" [37] "37" "38" "39" "40" "41" "42" "43" "44" "45" "46" "47" "48" [49] "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60" [61] "61" "62" "63" "64" "65" "66" "67" "68" "69" "70" "71" "72" [73] "73" "74" "75" "76" "77" "78" "79" "80" "81" "82" "83" "84" [85] "85" "86" "87" "88" "89" "90" "91" "92" "93" "94" "95" "96" [97] "97" "98" "99" "100" "101" "102" "103" "104" "105" "106" "107" "108" [109] "109" "110" "111" "112" "113" "114" "115" "116" "117" "118" "119" "120" [121] "121" "122" "123" "124" "125" "126" "127" "128" "129" "130" "131" "132" [133] "133" "134" "135" "136" "137" "138" "139" NA > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6w63k1354379682.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) > names(z) <- list('S.D.','Range','IQR') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 434.8749 4197.548 168.7825 [2,] 434.8749 4197.548 168.7825 [3,] 434.8749 4197.548 168.7825 [4,] 434.8749 4197.548 168.7825 [5,] 434.8749 4197.548 168.7825 $n [1] 1 1 1 $conf [,1] [,2] [,3] [1,] 434.8749 4197.548 168.7825 [2,] 434.8749 4197.548 168.7825 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" > dev.off() null device 1 > > try(system("convert tmp/1p1411354379682.ps tmp/1p1411354379682.png",intern=TRUE)) character(0) > try(system("convert tmp/22tvu1354379682.ps tmp/22tvu1354379682.png",intern=TRUE)) character(0) > try(system("convert tmp/33uc01354379682.ps tmp/33uc01354379682.png",intern=TRUE)) character(0) > try(system("convert tmp/44eew1354379682.ps tmp/44eew1354379682.png",intern=TRUE)) character(0) > try(system("convert tmp/5ktbo1354379682.ps tmp/5ktbo1354379682.png",intern=TRUE)) character(0) > try(system("convert tmp/6w63k1354379682.ps tmp/6w63k1354379682.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 5.113 1.123 6.298