Our choice of p = 2 is supported by WAIC [33] [34]. Predictors in X were top-level MeSH category, target class, estimated time the target has been under development, and RVIS score [17]. Details are provided in S5 Text. Priors were α∼N(μa,σa2); βj∼iidN(0,σb2)j=1,…,d γj∼iidN(0,σg2)k=0,1,2