R/define.grps.R

In GRIN2: Genomic Random Interval (GRIN)

# define.grps: Function defines the lesion groups to be included in the Kruskal-Wallis test.
# Lesion groups with fewer subjects than specified by min.grp.size are excluded from the test.
# At least two eligible groups must remain for the Kruskal-Wallis test; otherwise one.KW.pvalue returns NA.

define.grps=function(grps,               # Lesion groups
                     min.grp.size=NULL)  # minimum number of subjects in the lesion group to be included in the KW test.
{
  grp.tbl=table(grps)
  small.grp=(grp.tbl<min.grp.size)
  exclude.grp=grps%in%(names(grp.tbl)[small.grp])
  grps[exclude.grp]="EXCLUDE"
  return(grps)
}