Hierarchical models extend traditional multiple regression methods for exploring main effects and interactions in an epidemiological dataset by regressing the first-level coefficients on external data65-67. External information can include simple pathway indicator variables68, genomic annotation or pathway ontologies69, functional assays70, in silico predictions of function or evolutionary conservation71, or simulation of pathway kinetics72,73.