A Poisson-multivariate normal hierarchical model for measuring microbial conditional independence networks from metagenomic count data