Customizable number of inner knots. Extract AIC from REML output. Changed ntraits() to get_ntraits() for the sake of good naming practice.