Rでzeligパッケージを使用してロジットモデルを実行しています:
z.out <- zelig(trade961a ~ age962 + education962 + personal962 + economy962 + partisan962 + employment962 + union962 + home962 + market962 + race962 + income962, data=data96, model="logit")
従属変数 trade961a は二分因子変数です。他のすべての変数は数値です。
私はpersonal962に興味があります。
> str(personal962)
num [1:1714] 3 3 1 1 1 2 3 2 2 1 ...
したがって、次のコマンドを実行して、personal962 が 1 (より良い) の場合に国際貿易をサポートする確率をシミュレートします。
x.out <- setx(z.out, personal962=1)
次のエラー メッセージが表示されます。
#Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] :
# incorrect number of dimensions
私が間違っていることについてのアイデアはありますか?