r - probit with clustered standard error - Stata コマンドと同等

Question

Stataに次のprobitコマンドがあり、R で同等のコードを探します。

probit mediation viol ethniccomp  lncrisisdur  lncapratio  lnten_mean durable_avg neighbors totaldem_nbrhd geostr medprev jointdem if newcrisis==1, cluster(crisno)

係数の推定結果は再現できますが、修正された標準誤差 (クラスター化されている) は再現できません。

probit.3.1_1 <- glm(mediation ~           viol+ethniccomp+lncrisisdur+lncapratio+lnten_mean+durable_avg+neighbors+
                    totaldem_nbrhd+geostr+medprev+jointdem,
                    data=as.data.frame(basedata[basedata$newcrisis==1,]), family=binomial (link=probit)) 

私は基本的に R で Stata オプションに相当するものを探していますcluster(crisno)。

私はこの返信を見ましたが、私が知る限り、提案された解決策はロジットのみを参照しており、プロビットは参照していません。