我正在尝试估计 中包含 5 个方程的非线性看似不相关的回归 (SUR) 模型R,并且我正在研究该包systemfit。一切都很顺利,直到需要对我的方程设置一些限制。使用包systemfit,该函数nlsystemfit()适用于非线性方程。restrict.matrix但它不允许使用该选项/参数nlsystemfit()(它适用于函数 中的线性方程systemfit())。
一个简化的例子是(我认为显示数据与这里无关):
EQ_1 <- Y1 ~ (c1 - x)*Q + c11*G11 + c12*G12 + c13*G13
EQ_2 <- Y2 ~ (c2 - x)*Q + c21*G11 + c22*G12 + c23*G13
EQ_3 <- Y3 ~ (c3 - x)*Q + c31*G11 + c32*G12 + c33*G13
start.values <- c(c1 = 0,c2= 0,c3 = 0,
c11 = 0,c12 = 0,c13 = 0,
c21 = 0,c22 = …Run Code Online (Sandbox Code Playgroud) r cran systemfit non-linear-regression seemingly-unrelated-regression