ivreg 还是 ivprobit 用于二进制变量？

问：

我估计的是参与培训对就业的因果效应，其中参与是二元变量，结果（就业）是二元的，城市培训的存在是二元的。R 中的 ivprobit 函数是否适合此非线性模型，或者我可以使用 ivreg？

我试过了，但我在 R 中收到错误消息ivprobit

Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : 
  contrasts can be applied only to factors with 2 or more levels

我将不胜感激 ivprobit 的信息或链接，我不确定在代码中在哪里包含城市和年份固定效果。谢谢。