Adds replace=T to sampling in model fitting to deal with low number of obs.