質問するログイン新規登録

Q&A

0回答

1979閲覧

Rのfixest,lfe,tidyverseでの重回帰分析での切片を確認したい

umimarine

総合スコア6

R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

統計

統計は、集団現象を数量で把握することです。また、調査で得られた性質や傾向を数量的に表したデータのことをいいます。

0グッド

0クリップ

投稿2022/01/06 05:26

0

0

Rのfixest,lfe,tidyverseというパッケージを用いて、
重回帰分析(最小二乗法)を行っています。
この際、固定効果モデル(feols)で推定を行っています。

こちらのドキュメント(https://lrberge.github.io/fixest/reference/feols.html)
を見る限り、summary()で(intercept)という切片の相関係数、標準偏差、p値
を求められているのですが、自分の手元ではこれが求まりません。

以下に手元のソースコードを記します。

R

1# パッケージ 2install.packages("tidyverse") 3install.packages("lfe") 4install.packages("fixest") 5library(tidyverse) 6library(lfe) 7library(fixest) 8# データセットアップ 9data <- read.csv(file.choose()) 10# 従属変数 11model1 <- feols(ESGScore ~ ROA+DEBT+PTNCUL+PTNWOM+PRODUCTAFF+RECALL+INDPPTN+AUDIT+CONPENSATION+NOMINATING+AUDITAFF+CSRSUS+CSRREPORT+SOIL+ENVAFF+RENEWABLE+CO2+ENERGY+INWATER+OUTWATER+WASTE+SALES+STOCK|name+Year+Category+JP,data) 12summary(model1) 13

結果は以下のようになっています。

R

1OLS estimation, Dep. Var.: ESGScore 2Observations: 561 3Fixed-effects: name: 187, Year: 3, Category: 7, JP: 2 4Standard-errors: Clustered (name) 5 Estimate Std. Error t value Pr(>|t|) 6ROA -0.020872 0.041197 -0.506628 6.1302e-01 7DEBT 0.041211 0.054006 0.763091 4.4638e-01 8...

interceptが得られていません。
本当はこうなるはずです。

Document

1# To estimate Two stage least squares, 2# insert a formula describing the endo. vars./instr. relation after a pipe: 3 4base = iris 5names(base) = c("y", "x1", "x2", "x3", "fe1") 6base$x_inst1 = 0.2 * base$x1 + 0.7 * base$x2 + rpois(150, 2) 7base$x_inst2 = 0.2 * base$x2 + 0.7 * base$x3 + rpois(150, 3) 8base$x_endo1 = 0.5 * base$y + 0.5 * base$x3 + rnorm(150, sd = 2) 9base$x_endo2 = 1.5 * base$y + 0.5 * base$x3 + 3 * base$x_inst1 + rnorm(150, sd = 5) 10 11# Using 2 controls, 1 endogenous var. and 1 instrument 12res_iv = feols(y ~ x1 + x2 | x_endo1 ~ x_inst1, base) 13 14# The second stage is the default 15summary(res_iv) 16#> TSLS estimation, Dep. Var.: y, Endo.: x_endo1, Instr.: x_inst1 17#> Second stage: Dep. Var.: y 18#> Observations: 150 19#> Standard-errors: IID 20#> Estimate Std. Error t value Pr(>|t|) 21#> (Intercept) -5.444083 82.90731 -0.065665 0.94773 22#> fit_x_endo1 4.156853 44.67892 0.093038 0.92600 23#> x1 0.533483 1.80938 0.294843 0.76853 24#> x2 -1.358091 19.67381 -0.069030 0.94506 25...

どのようにすれば良いでしょうか?

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだ回答がついていません

会員登録して回答してみよう

アカウントをお持ちの方は

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.29%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問