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

#
# Break this data into 3 batches:
#
inds_1 = 1:20
Batch_1 = DF[inds_1, ]
inds_2 = 21:40
Batch_2 = DF[inds_2, ]
inds_3 = 41:60
Batch_3 = DF[41:60, ]

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))

xbar_3 = mean(Batch_3$xbar)
rbar_3 = mean(Batch_3$R)
LCL_xbar_3 = xbar_3 - A2*rbar_3
UCL_xbar_3 = xbar_3 + A2*rbar_3
LCL_range_3 = D3*rbar_3
UCL_range_3 = D4*rbar_3
print(sprintf('xbar_3= %f; rbar_3= %f', xbar_3, rbar_3))


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

par(mfrow=c(2,1))

xbar_min = min( c( DF$xbar, LCL_xbar_1, UCL_xbar_1, LCL_xbar_2, UCL_xbar_2, LCL_xbar_3, UCL_xbar_3 ) )
xbar_max = max( c( DF$xbar, LCL_xbar_1, UCL_xbar_1, LCL_xbar_2, UCL_xbar_2, LCL_xbar_3, UCL_xbar_3 ) )

plot( DF$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 )
lines( inds_3, rep( xbar_3, length(inds_3) ), col='black', lwd=2 )
lines( inds_3, rep( LCL_xbar_3, length(inds_3) ), col='red', lwd=2 )
lines( inds_3, rep( UCL_xbar_3, length(inds_3) ), col='red', lwd=2 )
grid()


range_min = min( c( DF$R, LCL_range_1, UCL_range_1, LCL_range_2, UCL_range_2, LCL_range_3, UCL_range_3 ) )
range_max = max( c( DF$R, LCL_range_1, UCL_range_1, LCL_range_2, UCL_range_2, LCL_range_3, UCL_range_3 ) )

plot( DF$R, pch=19, ylim=c(range_min, range_max), col='black', ylab='Rs' )
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 )
lines( inds_3, rep( rbar_3, length(inds_3) ), col='black', lwd=2 )
lines( inds_3, rep( LCL_range_3, length(inds_3) ), col='red', lwd=2 )
lines( inds_3, rep( UCL_range_3, length(inds_3) ), col='red', lwd=2 )
grid()

par(mfrow=c(1,1))

#dev.off()