2013-07-01 76 views
5

私はy~x(ちょうど1つの独立変数と1つの独立変数)の回帰をしたいが、私はheteroskedasticityを持っている。 xの増加に伴い、yの変動性が増加します。それに対処するため、Rの"gls()"関数を使って加重最小二乗法を使用したいと考えています。 しかし、私はそれをどのように使用するのか分かりません。 関数の "重み"引数に分散関数を適用する必要があります。しかし、私はどれを選ぶべきか、どのように使うべきかはわかりません。あなたは私に助けてくれますか?加重最小二乗 - R

ありがとうございました。

答えて

5

ここでは、データのようなポアソン数の世話をする例を示します。ここでは、変動は平均と比例します(あなたのように見えます)。

fit = lm (y ~ x, data=dat,weights=(1/dat$x^2)) 

値を掛けるため、レシピェールを重量として使用します。分散は単位の二乗を持つため、ポアソン計数データの世話をするためにそれを正方形にします。

fit = lm (y ~ x, data=dat,weights=(1/dat$x)) 

単純にx値でスケーリングしてよりうまくいくかどうかを確認します。

+0

どのように効果があるか評価できます。 Bartlettテストから返される値は、 "weights"引数の有無と同じです。ここ は私のコードです: A2 = read.tableを( "total37.txt"、ヘッダ= TRUE) M1 = LM(RES〜ModeF、A2、重み= 1/A2 $ ModeF^2) 平方メートル=のLM( RES〜ModeF、A2) bartlett.test(残差(M1)〜A2 $ ModeF)) バートレット検定 データ:$ ModeF バートレットのK-二乗残差(M1)とa2 = 35.2706、DF = 11、P-値= 0.0002236 bartlett.test(残差(M2)〜A2 $ ModeF)) バートレット検定 データ:残差(M2)とa2 $ ModeF バートレットのK-乗= 35.2706、DF = 11、p値= 0.0002236 – user1671537

関連する問題