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

x = seq(0.1,1,0.15)
n = length(x)
y = matrix(0,20,n)

set.seed(100)

### generate 10 IC profiles from the model y=2/(1+x)+e, where e~N(0,0.25^2)

for(i in 1:10){
   y[i,] = 2/(1+x)+rnorm(n,0,0.1)
}

### generate 10 OC profiles from the model y=1.5/(1+x^2)+e, where e~N(0,0.25^2)

for(i in 11:20){
   y[i,] = 1.5/(1+x^2)+rnorm(n,0,0.1)
}

write.table(round(y,digits=3),
            "exercise107.dat",row.names=F,col.names=F)