Example call for NCAA Data.

> source("k:\\LSA.R.txt")
> source("K:\\lsa.linear.txt")
> lsa.linear(x=ncaa2[,1:19],y=ncaa2$y)->out.ncaa.lsa

> names(out.ncaa.lsa)
[1] "fit"       "st"        "mse"       "x.ind"     "coeff"     "intercept"


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

> out.ncaa.lsa$intercept
[1] -45.47996

> out.ncaa.lsa$coeff
 [1]  0.00000000  3.41431119  0.24059475  0.66570683  0.19436932  0.06292095
 [7] -1.99666379 -0.27302515  1.58174126  0.00000000  0.00000000  0.00000000
[13]  0.00000000  0.00000000  0.00000000  0.00000000  0.00000000  0.00000000
[19]  0.00000000

For comparison, the output from using the lars package was

> 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 

which is identical.