data(amitriptyline)
Y <- amitriptyline[ , 1:2]
X1 <- amitriptyline[ , 4:7]
X2 <- amitriptyline[ , 3]
u <- u.eppls(X1, X2, Y)
u
m <- eppls(X1, X2, Y, 2)
pred.res <- pred.eppls(m, X1[1, ], X2[1])
pred.res
Run the code above in your browser using DataLab