Model selection with rho
Fix derivative and test