I think you need to set up a cut-off of Cp and then get the "good" values of Cp from adjr$Cp. -- View this message in context: http://r.789695.n4.nabble.com/best-model-cp-mallow-tp2550015p2550283.html Sent from the R help mailing list archive at Nabble.com.