We outline a fairly simple method to obtain in Stata Murphy–Topel-corrected variances for a two-step estimation of a heckprobit model with endogeneity in the main equation. The procedure uses predict’s score option and the powerful matrix tool accum in Stata and builds on previous works by Hardin (2002, Stata Journal 2: 253–266) and Hole (2006, Stata Journal 6: 521–529).