Variational Bayesian data fusion of multi-class discrete observations with applications to cooperative human-robot estimation