我正在使用 R 中的 zelig 包运行 logit 模型:
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 ...
所以我运行以下命令来模拟当person962为1(更好)时支持国际贸易的概率:
x.out <- setx(z.out, personal962=1)
我收到以下错误消息:
#Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] :
# incorrect number of dimensions
关于我做错了什么的任何想法?