Example call for NCAA data.

> lasso.adapt.bic2(x=ncaa2[,1:19],y=ncaa2$y)->out.ncaa2

> names(out.ncaa2)
[1] "fit"       "st"        "mse"       "x.ind"     "coeff"     "intercept"
[7] "object"    "bic2"      "step.bic2"

> out.ncaa2$x.ind
[1] 2 3 4 5 6 7 8 9

> out.ncaa2$intercept
[1] -45.47996

> out.ncaa2$coeff
         x1          x2          x3          x4          x5          x6 
 0.00000000  3.41431119  0.24059475  0.66570683  0.19436932  0.06292095 
         x7          x8          x9         x10         x11         x12 
-1.99666379 -0.27302515  1.58174126  0.00000000  0.00000000  0.00000000 
        x13         x14         x15         x16         x17         x18 
 0.00000000  0.00000000  0.00000000  0.00000000  0.00000000  0.00000000 
        x19 
 0.00000000 

> out.ncaa2$bic2
Intercept                                                                       
 519.5384  174.2636  162.2225  140.4111  139.7968  136.3603  132.0811  134.2287 
                                                                                
 129.4769  130.5416  134.5146  137.8518  142.3620  144.4098  147.6826  148.6982 
                                        
 151.4980  155.8465  160.3501  164.8659 

> plot(out.ncaa2$bic2,xlab="Model size",ylab="BIC2",ylim=c(100,200))
> title("Adaptive LASSO for NCAA Data\nBIC2 vs Model Size")