where N total = N c + N t + N b is the total number of observations from all the data sets. The gradient and Hessian of R are given in Appendix S2 of Supporting Information 1. We plug in the estimates (??), (??), and (??) to compute cov( θctb ).