added ddply's to calculate bimodality by system and taxa