  Recent technological advances have made it possible to simultaneously measure
multiple protein activities at the single cell level. With such data collected
under different stimulatory or inhibitory conditions, it is possible to infer
the causal relationships among proteins from single cell interventional data.
In this article we propose a Bayesian hierarchical modeling framework to infer
the signaling pathway based on the posterior distributions of parameters in the
model. Under this framework, we consider network sparsity and model the
existence of an association between two proteins both at the overall level
across all experiments and at each individual experimental level. This allows
us to infer the pairs of proteins that are associated with each other and their
causal relationships. We also explicitly consider both intrinsic noise and
measurement error. Markov chain Monte Carlo is implemented for statistical
inference. We demonstrate that this hierarchical modeling can effectively pool
information from different interventional experiments through simulation
studies and real data analysis.
