Let i index gene target-indication pairs (gi, ti), i = 1, …, N. Let yi ∈ {0, 1} be 1 if pair i is found in at least one US/EU approved drug and 0 otherwise. Let X be an N × d design matrix where d is the number of non-genetic predictors with ith row xi′. yi∼indepBernoulli(logit-1(α+ηi+xi′β))i=1,…,N where ηi={∑j=0pγjEC(gi,ti)jthereexists(g,t)∈Asuchthatg=gi0otherwise