# NOT RUN {
# Generate simulated data
n <- 200
p <- 20
k <- 5
rho <- 0.4
SNR <- 10
cortype <- 1
seed <- 10
Tbeta <- rep(0, p)
Tbeta[1:k*floor(p/k):floor(p/k)] <- rep(1, k)
Data <- gen.data(n, p, k, rho, family = "gaussian", cortype = cortype, snr = SNR, seed = seed)
lm.bss <- bess(Data$x, Data$y, method = "sequential")
logLik(lm.bss)
logLik(lm.bss, best.model = FALSE)
# }
Run the code above in your browser using DataLab