Clustering of non-Gaussian data by variational Bayes for normal inverse Gaussian mixture models