Двухэтапный метод наименьших квадратов в R

Я хочу запустить двухэтапную пробитную регрессию наименьших квадратов в R. Кто-нибудь знает, как это сделать? Есть ли там какой-нибудь пакет? Я знаю, что это можно сделать с помощью Stata, поэтому я полагаю, что это можно сделать с помощью R.