# NOT RUN {
# load HDtweedie library
library(HDtweedie)
# load auto data set
data(auto)
# fit the lasso
m0 <- HDtweedie(x=auto$x,y=auto$y,p=1.5)
# print out results
print(m0)
# define group index
group1 <- c(rep(1,5),rep(2,7),rep(3,4),rep(4:14,each=3),15:21)
# fit the grouped lasso
m1 <- HDtweedie(x=auto$x,y=auto$y,group=group1,p=1.5)
# print out results
print(m1)
# }
Run the code above in your browser using DataLab