## 1. load results from Klein (2015a)
data(klein15a)
M <- klein15a$model.list
## 2. extract variables
X <- do.call(rbind.data.frame, M$X)
eta <- c(klein15a$coefs$eta, rep(0, length(M$X)-length(M$W)))
X <- cbind(X,eta)
y <- unlist(M$R)
## 3. apply KHB method
khb(X=X, y=y, z="eta")
Run the code above in your browser using DataLab