> fsr.fast.include(x=ncaa2[,1:19],y=ncaa2$y,inc=c(12,3,5))
   var   pval  pvmax    Rsq
1   12 0.0000 0.0000 0.6042
2    3 0.0000 0.0000 0.6783
3    5 0.0000 0.0000 0.6874
4    2 0.0000 0.0000 0.7710
5    4 0.0043 0.0043 0.7914
6    7 0.0028 0.0043 0.8118
7   17 0.0539 0.0539 0.8198
8   15 0.0458 0.0539 0.8281
9    6 0.0976 0.0976 0.8336
10   9 0.0962 0.0976 0.8391
11   8 0.0281 0.0976 0.8484
12  10 0.2864 0.2864 0.8505
13  13 0.3163 0.3163 0.8524
14  18 0.2697 0.3163 0.8546
15  11 0.4953 0.4953 0.8555
16   1 0.6326 0.6326 0.8559
17  14 0.7056 0.7056 0.8562
18  19 0.8605 0.8605 0.8563
19  16 0.9032 0.9032 0.8563
$mod

Call:
lm(formula = y ~ x)

Coefficients:
(Intercept)          x12           x3           x5           x2           x4  
   -37.6733      -0.1770       0.2259       0.2916       3.3911       0.6966  
         x7  
    -2.5675  


$size
[1] 6

$x.ind
[1] 12  3  5  2  4  7

$alphahat.ER
[1] 0.01538462

$inc
[1] 12  3  5