観測数最初の相違

Question

とPLMを使用しているとき、私は、以下のいずれかのようなデータセットでplmを使用してパネルデータとの回帰を実行した後、簡単な問題があります：あなたが見ることができるように

dataset <- data.frame(id = rep(c(1,2,3,4,5), 2), 
         time = rep(c(0,1), each = 5), 
         group = rep(c(0,1,0,0,1), 2), 
         Y = runif(10,0,1)) 
model <-plm(Y ~ time*group, method = 'fd', effect = 'twoways', data = dataset, 
      index = c('id', 'time')) 
summary(model) 

stargazer(model) 

を、両方のモデルsummaryで表示され、stargazerで表示される表は、私の観測数が10であると言います。しかし、最初の違いで時間要素を取り除いたので、N = 5と言うのは正しいとは限りませんか？

Answer 1

あなたの所見は正しいですか？しかし、あなたのコードはあなたがしたいものではありません（最初の差分モデル）。

あなたが最初の差分モデルをしたい場合は、modelに引数methodを切り替える（と、それは最初の差分モデルの意味がありませんので、引数effectを削除）：要約出力で

model <-plm(Y ~ time*group, model = 'fd', data = dataset, 
      index = c('id', 'time')) 
summary(model) 

## Oneway (individual) effect First-Difference Model 
## 
## Call: 
## plm(formula = Y ~ time * group, data = dataset, model = "fd", 
##  index = c("id", "time")) 
## 
## Balanced Panel: n = 5, T = 2, N = 10 
## Observations used in estimation: 5 
## 
## Residuals: 
##  Min. 1st Qu.  Median 3rd Qu.  Max. 
## -0.3067240 -0.0012185 0.0012185 0.1367080 0.1700160 
## [...] 

、あなたがすることができます元のデータの観測数（N = 10）とFDモデルで使用された観測数（5）を参照してください。