Bayesian hierarchical modeling of means and covariances of gene expression data within families