私は現在、多重帰属データセットに含まれているキー前処理共変量を照合し、このデータに対して線形回帰を実行する必要がある論文を作成しようとしています。私はamelia
から得た多重帰属データに対して線形回帰を実行するコードをzelig
に書くことができました。しかし、cem
を使用して重要な治療前共存症との一致を試みると、モデルを実行するのにzelig
を取得できません。私は最後のエラーは意味し、なぜ私は多重代入を使用した後に設定一致したデータの回帰を実行することはできません何を理解していないのですAmelia + CEM + Zeligエラー
> # Run Amelia, 5 imputations!
> complete <- amelia(clean, m = 5, ts = "wave", cs = "country", idvars = c("idnum", "strata", "weight"))
-- Imputation 1 --
1 2 3 4 5 6 7 8 9 10 11 12
-- Imputation 2 --
1 2 3 4 5 6 7 8 9 10 11
-- Imputation 3 --
1 2 3 4 5 6 7 8 9 10 11 12
-- Imputation 4 --
1 2 3 4 5 6 7 8 9 10 11 12
-- Imputation 5 --
1 2 3 4 5 6 7 8 9 10 11 12
>
> # Estimate naive and controlled linear models with MI data sets using Zelig
>
> lm.out.naive <- zelig(dv ~ treat*year + treat + year, data = complete, model = "ls")
> lm.out.control <- zelig(dv~ treat*year + treat + year + age + gender + ideo + married, data = complete, model = "ls")
> # Estimate CEM linear model
>
> # First, generate matched data set
> imputed <- complete$imputations[1:15]
> matched <- cem("treat", drop = c("pais", "wave", "year", "idnum", "strata"),
+ datalist = imputed, data = clean)
>
> # Then, run regression
> lm.out.match <- zelig(dv ~ treat*year + treat + year + age + gender + ideo + married, data = matched, model = "ls")
Error in resolve_vars(new_groups, tbl_vars(.data)) :
unknown variable to group by : by
:ここに出力ROMが私のコンソールです。この問題に関する洞察は高く評価されます。