A Poisson binomial-based statistical testing framework for comorbidity discovery across electronic health record datasets