1

私には2つの扶養家族がいて、どちらも2つの変数と相互に依存しています。これは、Rでモデル化できますか(必ず!)、方法がわかりません。誰かヒントがありますか?

明確に言えば:

次のモデルでデータをモデル化したい:

Y1=X1*coef1+X2*coef2
Y2=X1*coef2+X2*coef3

注:coef2は両方の行Xiに表示され、Yiはそれぞれ入力データと出力データです。

私はここまで来ました:

lm(Y1~X1+X2,mydata)

相互依存関係を含むモデルの2行目を追加するにはどうすればよいですか?

あなたの助けは大歓迎です!乾杯、バスティアーン

4

1 に答える 1

4

これを試して:

# sample data - true coefs are 2, 3, 4
set.seed(123)
n <- 35
DF <- data.frame(X1 = 1, X2 = 1:n, X3 = (1:n)^2)
DF <- transform(DF, Y1 = X1 * 2 + X2 * 3 + rnorm(n), 
          Y2 = X1 * 3 + X2 * 4 + rnorm(n))

# construct data frame for required model
DF2 <- with(DF, data.frame(y = c(Y1, Y2), 
x1 = c(X1, 0*X1),
x2 = c(X2, X1),
x3 = c(0*X2, X2)))

lm(y ~. - 1, DF2)

確かに、2、3、4の真のcoefsを回復することがわかります。

> lm(y ~. - 1, DF2)

Call:
lm(formula = y ~ . - 1, data = DF2)

Coefficients:
   x1     x2     x3  
2.084  2.997  4.007  
于 2010-12-18T09:14:33.653 に答える