source('../Chapter6/chap_6_prob_7_data.R')

#
# This data comes in 2 batches:
#
inds_1 = 1:20
inds_2 = 21:40
n = 5

# From Table C:
#
A2 = 0.577
D3 = 0
D4 = 2.115
d2 = 2.326

# Each batch one at a time:
#
xbar_1 = mean(Batch_1$xbar)
rbar_1 = mean(Batch_1$R)
LCL_xbar_1 = xbar_1 - A2*rbar_1
UCL_xbar_1 = xbar_1 + A2*rbar_1
LCL_range_1 = D3*rbar_1
UCL_range_1 = D4*rbar_1
print(sprintf('xbar_1= %f; rbar_1= %f', xbar_1, rbar_1))

xbar_2 = mean(Batch_2$xbar)
rbar_2 = mean(Batch_2$R)
LCL_xbar_2 = xbar_2 - A2*rbar_2
UCL_xbar_2 = xbar_2 + A2*rbar_2
LCL_range_2 = D3*rbar_2
UCL_range_2 = D4*rbar_2
print(sprintf('xbar_2= %f; rbar_2= %f', xbar_2, rbar_2))

#postscript("../../WriteUp/Graphics/Chapter6/chap_6_prob_7_plot.eps", onefile=FALSE, horizontal=FALSE)

par(mfrow=c(2,1))

xbar_min = min( c( Batch_1$xbar, Batch_2$xbar, LCL_xbar_1, UCL_xbar_1, LCL_xbar_2, UCL_xbar_2 ) )
xbar_max = max( c( Batch_1$xbar, Batch_2$xbar, LCL_xbar_1, UCL_xbar_1, LCL_xbar_2, UCL_xbar_2 ) )

plot( c( Batch_1$xbar, Batch_2$xbar ), pch=19, ylim=c(xbar_min, xbar_max), col='black', ylab='x bars' )

lines( inds_1, rep( xbar_1, length(inds_1) ), col='black', lwd=2 )
lines( inds_1, rep( LCL_xbar_1, length(inds_1) ), col='red', lwd=2 )
lines( inds_1, rep( UCL_xbar_1, length(inds_1) ), col='red', lwd=2 )
lines( inds_2, rep( xbar_2, length(inds_2) ), col='black', lwd=2 )
lines( inds_2, rep( LCL_xbar_2, length(inds_2) ), col='red', lwd=2 )
lines( inds_2, rep( UCL_xbar_2, length(inds_2) ), col='red', lwd=2 )
grid()


range_min = min( c( Batch_1$R, Batch_2$R, LCL_range_1, UCL_range_1, LCL_range_2, UCL_range_2 ) )
range_max = max( c( Batch_1$R, Batch_2$R, LCL_range_1, UCL_range_1, LCL_range_2, UCL_range_2 ) )

plot( c( Batch_1$R, Batch_2$R ), pch=19, ylim=c(range_min, range_max), col='black', ylab='range' )
lines( inds_1, rep( rbar_1, length(inds_1) ), col='black', lwd=2 )
lines( inds_1, rep( LCL_range_1, length(inds_1) ), col='red', lwd=2 )
lines( inds_1, rep( UCL_range_1, length(inds_1) ), col='red', lwd=2 )
lines( inds_2, rep( rbar_2, length(inds_2) ), col='black', lwd=2 )
lines( inds_2, rep( LCL_range_2, length(inds_2) ), col='red', lwd=2 )
lines( inds_2, rep( UCL_range_2, length(inds_2) ), col='red', lwd=2 )
grid()

par(mfrow=c(1,1))

#dev.off()