Subset selection for linear mixed models