# The data for Problem 1:

Y.values<-c(975.51,1734.02,2177.37,2965.93,1250.45,2821.97,0.00,2322.25,0.00,
972.20,4352.78,2546.13,2186.04,0.00,869.86,1281.62,3557.40,1234.57,741.99)
W.values <-c(1,1,1,1,1,1,0,1,0,1,1,1,1,0,1,1,1,1,1)


# The age and birthweight data for Problem 3:

y <- c(2289,3393,3271,2648,2897,3327,2970,2535,3138,3573)/1000
 x1 <- c(15,17,18,15,16,19,17,16,18,19)
 X <- as.matrix( cbind(rep(1, length(x1)), x1) )