4

私texregはいくつかの変量効果モデル (を使用して推定plm) の結果を表で報告するために使用します。

Hausman 検定の p 値 (各モデルを対応する固定効果と比較する) を、によって報告される適合度測定に追加するにはどうすればよいtexregですか? より一般的には、追加の適合度測定値を で報告するにはどうすればよいtexregですか?

4

2 に答える 2

5

というモデル オブジェクトがあると仮定するとmodel、次の 3 つの手順に従う必要があります。

1)オブジェクトtr <- extract(model)を作成するために実行しtexregます。

2) この texreg オブジェクトを操作します。例：

tr@gof <- c(tr@gof, 0.5)
tr@gof.names <- c(tr@gof.names, "new row")
tr@gof.decimal <- c(tr@gof.decimal, TRUE)

3) 元のtexregコマンドを実行しますがtr、元の の代わりにオブジェクトを含めますmodel。例: screenreg(tr).

デフォルトで texreg パッケージに含まれる抽出メソッドに GOF 測定を含める必要があると思われる場合は、R-Forge の texreg フォーラムで提案を行う必要があります。

于 2013-11-11T23:19:14.287 に答える
1

Hausman 検定の p 値 (およびその他の追加の GOF 測定値) を texreg テーブルに追加する方法を見つけました。このメソッドは GOF 測定値を自動的に取得するため、手動で行う必要はありません。ご不明な点がございましたらお問い合わせください。

「random1」と「fixed1」がモデルであるとします

最初にHausmanテストをオブジェクトにします

ph1 <- phtest(random1, fixed1)

次に、モデルの係数を抽出します

sf1 <- summary(fixed1)
names1 <- rownames(sf1$coef)
co1 <- sf1$coef[, 1]
se1 <- sf1$coef[, 2]
pval1 <- sf1$coef[, 4]

次に、選択した GOF を抽出します

pval.ftest1 <- (summary(fixed1)$fstatistic)$p.value
rs1 <- sf1$r.squared
adj1 <- sf1$adj.r.squared
n1 <- nobs(fixed1)
phtest1 <- fp.value.fixed1

テーブルに含める GOF を選択する

gof1 <- c(pval.ftest1, rs1, adj1, n1, phtest1)
gof.names1 <- c("p\\ (F-test)","R$^2$", "Adj.\\ R$^2$", "Num.\\ Obs.", "p\\ (Hausman\\ test)")

GOF の小数点以下の桁数を選択する

decimal.places <- c(TRUE, TRUE, TRUE, FALSE, TRUE)

texreg オブジェクトを作成する

tr1 <- createTexreg(coef.names = names1, coef = co1, se = se1, pvalues = pval1, gof.names = gof.names1, gof = gof1, gof.decimal = decimal.places)

screenreg ですべてが正常に見えることを確認する

screenreg(tr1)

必要に応じて、すべてのモデル (fixed2、fixed3...& random2、random3...) に対して同じことを繰り返します。

texreg または htmlreg 関数を使用する前に、すべてのモデルですべてが正常に見えることを確認してください

screenreg(list(tr1, tr2, tr3, tr4, tr5))
于 2016-12-14T21:37:09.237 に答える