The regsubsets does work. I also notice the bestglm, will see how it work too. Thank you. -- View this message in context: http://r.789695.n4.nabble.com/input-matrix-for-leaps-algorithm-tp2954453p2956866.html Sent from the R help mailing list archive at Nabble.com.