##### R-code for generating data in Exercise 8.15

#### Generate 100 phase II observations for online process mean monitoring.
#### The first 50 observations are from the distribution t_4/sqrt(2), and
#### the remaining 50 observations are from t_4/sqrt(2)+1.

set.seed(100)

x=matrix(c(rt(50,4)/sqrt(2),rt(50,4)/sqrt(2)+1),ncol=10,byrow=T)

write.table(round(x,digits=2),"exercise815.dat",row.names=F,col.names=F)