two stage residual inclusion 程序

求教各位大咖，如何在stata实现two stage residual inclusion？小弟按照Terza（2017）的介绍撰写的代码不尽如人意，特别是无法计算标准差，标准差矩阵全是缺漏值。另外，若内生解释变量为二元变量，又如何使用2RIS？感激不尽。