##### R-code for generating data used in Exercise 8.16

set.seed(1000)

#### Generate 500 observations of the IC data from the
#### IC process distribution t_4/sqrt(2)

y = rt(500,4)/sqrt(2)

#### Define the boundary q1, q2, ... q_p-1 where p=10

y1=sort(y)      # ordered y from the smallest to the largest

q1=(y1[50]+y1[51])/2
q2=(y1[100]+y1[101])/2
q3=(y1[150]+y1[151])/2
q4=(y1[200]+y1[201])/2
q5=(y1[250]+y1[251])/2
q6=(y1[300]+y1[301])/2
q7=(y1[350]+y1[351])/2
q8=(y1[400]+y1[401])/2
q9=(y1[450]+y1[451])/2

write.table(round(matrix(y,ncol=10),digits=2), 
            file="ex816.dat",col.names=F, row.names=F)