の 5 つの方程式を使用して、非線形の一見無関係な回帰 (SUR) モデルを推定しようとしRていますが、パッケージを調べていましたsystemfit。私の方程式にいくつかの制限を設定する必要があるまで、すべてがうまくいきます。パッケージを使用すると、非線形方程式に対してsystemfit機能する関数になります。しかし、それが許可されていないnlsystemfit()オプション/パラメーター(関数内の線形方程式で機能します)。restrict.matrixnlsystemfit()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 = 0,c23 = 0,
c31 = 0,c32 = 0,c33 = 0)
model <- list(EQ_1 ,EQ_2 ,EQ_3)
model.sur <- nlsystemfit(method = "SUR",
eqns = model,
startvals = start.values,
data = as.data.frame(dat))
推定はこれまでのところ完全に機能しています。しかし今、次の制約を設定する必要があります。
Rest_1 <- c11 + c12 + c13 = 0
Rest_2 <- c21 + c22 + c23 = 0
Rest_3 <- c31 + c32 + c33 = 0
Rest_4 <- c1 + c2 + c3 = -1
明らかに、ここでのモデルは 3 つの方程式で線形ですが、これはアイデアを単純化しようとしているからです。しかし、現在のモデルには 5 つの非線形方程式とさらに多くのパラメーターがあります。
Rで制限のある非線形SUR推定を実行する方法について、誰でも私を案内できますか?
よろしくお願いします。