Robust methods for analysis of 3-way compositional data in R