0

私は現在、mitools5 つの帰属データセットを分析するために使用しています (私は Amelia を使用しました)。分析の単位は国年です。各データセットには 3600 を超える観測値があります。すべての分析で対象となる従属変数は、ラグ カウント変数です。

データの性質 (パネルとカウント) を考慮して、pglm package. 次のコードを使用して、推定された 5 つのデータセットすべてで pglm モデルを推定し、これらの推定値を mitools と組み合わせました。

setwd("C:/Desktop")

library(mitools)
library(pglm)

#The data.dir code below establishes where the files are located
data.dir <- "C:/Desktop"

#The imputationList function combines all 5 imputed gtot datasets into one object
#called allgtot
allcountdata <- imputationList(lapply(list.files(data.dir, pattern = "panelcount.\\.dta", full=TRUE),read.dta, warn.missing.labels=FALSE))

#To see if the combining was successful
allcountdata

#Lets see if all of the column names are still there
colnames(allcountdata)


count_mitools<-pglm(count_lag ~     
                        HROsec
                        + hrofilled
                        + physint
                        + democracy
                        + log(PopTotal)
                        + GDPpercapita
                        + cell
                        + gini
                        + polity2
                        + CivilWar
                        + milper
                        + elp
                        + factor(year),
                        allcountdata,
                        family="negbin",
                        model="within", 
                        print.level=0,
                        method="nr",
                        index=c("cowcode", "year"))

ただし、次のエラーが表示されます。

Error in is.finite(x[[i]]) : 
  default method not implemented for type 'list'

次のように、リスト オブジェクトをデータフレームに強制しようとします。

as.data.frame(allcountdata)

しかし、R は私がそうすることができないと述べています。

pglm の見積もりを機能させるにはどうすればよいですか?

pglm を使用できない場合、次のいずれか (できれば両方) のパネル カウント データ回帰分析に適した代替モデル (および対応する R パッケージ) は何ですか?

a) mitools 分析と互換性がある
b) 推定のためのリスト オブジェクトの受け入れ

私が作業している 5 つの帰属データセットすべてを含む圧縮 zip ファイルへのリンクを提供しています。

4

0 に答える 0