#--- simulation ---
set.seed(1)
n <- 100; p <- 500
X <- matrix(rnorm(n=n*p),nrow=n,ncol=p)
beta <- rnorm(p)
prior <- beta + rnorm(p)
y <- X %*% beta
#--- print.transreg ---
object <- transreg(y=y,X=X,prior=prior)
object
Run the code above in your browser using DataLab