Я использую модель logit с использованием пакета zelig в R:
z.out <- zelig(trade961a ~ age962 + education962 + personal962 + economy962 + partisan962 + employment962 + union962 + home962 + market962 + race962 + income962, data=data96, model="logit")
Зависимая переменная trade961a является переменной дихотомического фактора. Все остальные переменные являются числовыми.
Меня интересует личный962, а именно:
> 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
Есть идеи, что я делаю не так?