x <- c(97.51,96.65,95.91,95.86,95.7,95.57,95.57,95.57,94.87,95.07,95.13,95.48,95.38,95.38,95.48,95.77,94.78,92.51,92.17,91.75,90.43,90.55,90.37,90.4,90.41,90.41,90.41,89.77,89.77,89.77,89.37,89.81,89.07,89.84,89.73,90.02,88.39,90.13,90.13,90.37,89.73,89.73,89.73,89.73,89.6,89.63,86.42,86.8,86.51,86.41,86.39,86.62,85.85,87.36,87.28,87.35,87.35,87.35,87.38,88.17,88.37,87.44,87.44,87.47,87.47,87.48,87.11,87.11,86.26,86.28,86.28,86.28)
par1 = '12'
par1 <- '12'
#'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))
(np <- floor(n / par1))
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
arr
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/1p1841354526010.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()
postscript(file="/var/wessaorg/rcomp/tmp/2duox1354526010.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()
postscript(file="/var/wessaorg/rcomp/tmp/3adh21354526010.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()
postscript(file="/var/wessaorg/rcomp/tmp/4vbez1354526010.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'))
dev.off()
postscript(file="/var/wessaorg/rcomp/tmp/53cvz1354526010.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'))
dev.off()
postscript(file="/var/wessaorg/rcomp/tmp/6xnh01354526010.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'))
dev.off()

try(system("convert tmp/1p1841354526010.ps tmp/1p1841354526010.png",intern=TRUE))
try(system("convert tmp/2duox1354526010.ps tmp/2duox1354526010.png",intern=TRUE))
try(system("convert tmp/3adh21354526010.ps tmp/3adh21354526010.png",intern=TRUE))
try(system("convert tmp/4vbez1354526010.ps tmp/4vbez1354526010.png",intern=TRUE))
try(system("convert tmp/53cvz1354526010.ps tmp/53cvz1354526010.png",intern=TRUE))
try(system("convert tmp/6xnh01354526010.ps tmp/6xnh01354526010.png",intern=TRUE))