More updates to RxNseq function:
  - return each coefficient only if significant in anova() call
  - name coefficient based on model formulat in returned data.frame