2017-02-27 15 views
0

私は現在、多重帰属データセットに含まれているキー前処理共変量を照合し、このデータに対して線形回帰を実行する必要がある論文を作成しようとしています。私は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が私のコンソールです。この問題に関する洞察は高く評価されます。

答えて

0

私は最近、多重帰属データフレームにzeligのロジットを使用してこのエラーが発生しました。私の解決策は、data =パラメータを適切に指定することでした。

lm.out.match <- zelig(dv ~ treat*year + treat + year + 
age + gender + ideo + married, data = mi(matched$imputations), model = "ls") 

私はzeligの古いバージョンは、このフォーマットを必要としなかったことを思い出すが、代わりに-があるとして「アメリア」クラスのオブジェクトを扱うように見えた:あなたはこれを試してみてください。

関連する問題