r - R: クラス "c('pmg', 'panelmodel')" のオブジェクトの二重クラスター標準エラーを取得する方法はありませんか?

Question

Fama-Macbeth 回帰を推定しています。このサイトからコードを取得しました

fpmg <- pmg(Mumbo~Jumbo, test, index=c("year","firmid"))
summary(fpmg)     
Mean Groups model   
Call:
pmg(formula = Mumbo ~ Jumbo, data = superfdf, index = c("day","Firm")) 

Residuals
 Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
-0.142200 -0.006930  0.000000  0.000000  0.006093  0.142900
Coefficients
               Estimate  Std. Error z-value Pr(>|z|)
(Intercept) -3.0114e-03  3.7080e-03 -0.8121   0.4167
Jumbo        4.9434e-05  3.4309e-04  0.1441   0.8854
Total Sum of Squares: 1.6915
Residual Sum of Squares: 0.86425
Multiple R-squared: 0.48908`

fpmg を推定した後、二重クラスタリングでロバストな SE を推定します。

vcovDC <- function(x, ...){
vcovHC(x, cluster="group", ...) + vcovHC(x, cluster="time", ...) - 
    vcovHC(x, method="white1", ...)}
coeftest(fpmg, vcov=function(x) vcovHC(x, cluster="group", type="HC1"))

次のエラーが表示されます。

Error in UseMethod("estfun") : 
  no applicable method for 'estfun' applied to an object of class "c('pmg', 'panelmodel')"

このエラーを解決する方法を提案してください。

更新: 「multiwayvcov」パッケージも試しましたが、同じエラーが表示されます。これらのパッケージ (Sandwich、multiwayvcov など) ではオブジェクト クラスが許可されていないようです。Rは本質的に私のすべての労力を無駄にするようで、私は行き止まりにぶつかりました。上記をPythonで行う方法を見つけました（コードを意味します）が、それについての知識はありません。

Rで問題を解決する方法はありませんか?

---