Code by Jiang and Wild for the Auckland semiparametric approach to two-phase models
This is mostly functional now, but needs packaging updates and a few code corrections to run in contemporary R.
Longer term, I want to make it run off twophase() objects from the survey package