Now handles NA values as expected in all cases (i.e trims them from groups)
Bug in initial cluster centre randomization fixed.