1

基本的に、多くの線形回帰を実行し、適合モデルごとに推定係数と標準誤差の結果を以前に作成したデータ フレームに格納するスクリプトの効率を上げようとしていますresults。

したがって、results回帰係数を格納する前に、データ フレームには必要なディメンションが既に含まれています。

また、各i-th 回帰について、次のようにします。

mod.fit <- plm(y ~ x1 + x2, index="group", sample)

そして、私は実行します:

  results[i,1] <- summary(m.fit)$coefficients[1,1]
  results[i,2] <- summary(m.fit)$coefficients[2,1]
  results[i,3] <- summary(m.fit)$coefficients[1,2]
  results[i,4] <- summary(m.fit)$coefficients[2,2]

上記のストレージステップを高速化する方法はありますか?

.

4

1 に答える 1

6

マトリックスインデックスを使用できます。

results[i,1:4] <- summary(m.fit)$coefficients[matrix(c(1,2,1,2,1,1,2,2),ncol=2)]

幅が4列しかない場合は、左側をresults削除できます。1:4

代わりに

results[i,] <- summary(m.fit)$coefficients[1:2,1:2]

Rは行列を列の最初の順序で格納するため、機能するはずです。

クラスに定義されている場合は、coef()アクセサではなくアクセサを使用することをお勧めします...$coefficientssummary.plm

于 2013-01-05T19:46:06.243 に答える