WaldtestがRでstmazerと表示されたplmとresultでF統計量を取得しますか？

Question

私はアンバランスなショートパネルを使っています。 bankFull.xlsx

私が実際に望むのは、2つの側面の固定効果とStataで非常に簡単な堅牢なS.Eが報告されているだけです。私はオンラインチュートリアルに従いましたが、いつも何か問題がありました

# Adjust F statistic 
wald_results <- waldtest(FE1, vcov = cov1) 
Error in model.matrix.pFormula(formula, data, rhs = 1, model = model, : 
    NA in the individual index variable 

どのようにデータを調整しても問題ありません。それは私を夢中にさせる。コードが示すように

bankFull <- openxlsx::read.xlsx("bankFull.xlsx",1) 

attach(bankFull) 
library(plm) 

FE1 = plm( RoA ~ 
       log(1+degreeNW)+ 
       ln_assets+ 
       log(no_of_board_members/staffNo)+ 
       log(no_of_branch_covered_city)+ 
       log(operation_year)+ 
       `RoA-1`+ 
       log(staffNo), 
      data = bankFull, index = c("name","year"), 
      effect="twoways",na.action = na.omit, 
      model= "within") 

# robust S.E.----------- 
library(sandwich) 
library(lmtest) # waldtest; see also coeftest. 
library(stargazer) 

# Adjust standard errors 
cov1   <- vcovHC(FE1, type = "HC1") 
robust_se <- sqrt(diag(cov1)) 

# Adjust F statistic 
wald_results <- waldtest(FE1, vcov = cov1) 

# show results. how can I get the F value? 
stargazer(FE1, FE1, type = "text", 
      se  = list(NULL, robust_se), 
      omit.stat = "f") 

第二に、私は結果を実証するスターゲイザーを使用します。

は、ここに私のコードです。私はまた、調整されたF値を表に示す必要があります。パッケージには使用できるオプションがありますか？

Answer 1

編集：plmのCRANのバージョン1.6から4を使用したPLMの1.6から4

を解放CRANに応じて情報を更新する機能pwaldtestして、モデルの堅牢なFテストをサポートしています（この機能は、中Ftestと呼ばれていましたCRANのリリース前に改名されました）。

例：

data("Grunfeld", package = "plm") 
mod_fe <- plm(inv ~ value + capital, data = Grunfeld, model = "within") 
plm::pwaldtest(mod_fe, test = "F") 


# with robust vcov 
plm::pwaldtest(mod_fe, test = "F", vcov = vcovHC(mod_fe)) 
plm::pwaldtest(mod_fe, test = "F", vcov = function(x) vcovHC(x, type = "HC3")) 
summary(mod_fe, vcov = vcovHC) 

ロバスト値飼料に（ロバスト標準誤差、Tおよびp値を、F値関連p値）を使用引数 se、t、pおよびFのためのstargazerコマンドの単純にadd.linesをテストします（デフォルトでは、stargazerによって生成されたF統計を省略します）。 http://jakeruss.com/cheatsheets/stargazer.html（堅牢な標準エラー（Stataのロバストなオプションを複製する） "のセクションを参照してください）。