2016-04-14 6 views
0

私はソフトウェアRXLStatで作業しています。私は片道ANOVAを実行しました(私のカテゴリ変数は3モーダル(1,2,3)で、私の応答変数は1〜10のスケールで定量的です)。RとXLStatを使用したAnovaのベータ係数

私はR and XLStatでこのANOVAを実行しました。Fフィッシャー、p値、係数推定、t値、stdエラーなどの出力はまったく同じです。

しかし、XLstatは標準出力係数(あまりにもベータ係数と呼ばれます)を出力します。まず、カテゴリ変数のベータ係数を計算することはできないと思っていたので、私は驚きました。私が読んだ参考文献によれば、それは意味がありません。

とにかく、私は見つけたユニークな式のおかげで、これらの係数を見つけようとしました:beta = estimate * sd(x)/sd(y)。 sd(x)はカテゴリ変数の標準偏差(sd(x)を計算するためにRを持つ数値変数として自動的に変換されます)は論理的に見え、sd(y)は私の応答変数の標準偏差です。

Rで得られた最初のベータIはXLstatと同じですが、2番目と3番目のベータは同じではありません。最初のものがRとXLStatと同じであるとすれば、Xlstatは数値変数のカテゴリ変数も変換すると思います(これは意味がありませんが、これは問題ではありません)。 さらに、XLStatが間違いを犯したかどうかを確かめるためにStatisticaでanovaを実行しましたが、ベータ係数の出力はXlstatと同じです...

これは私の質問です。 1つの方法でのベータ係数Anova?

次に、これらのベータ係数のカテゴリ変数への関連性についてお聞きしたいと思います。私の思考と私が読んだ出版物によると、それは意味がありません...

psとは対照的に、Xlstatは合計(ai)= 0です。ベータ係数の場合、XLStatは切片を削除します。私はこの事実が重要であると思うが、何とかわからない。

答えて

1

ANOVAのメトリック係数からベータ係数を得る式は、線形回帰の場合と同じです。係数は(カテゴリー変数のための)有意義な解釈を持たないが、標準化された係数はIVの相対的効果を異なるメトリクスと比較するのに有用である。

Rでは、モデルをフィットさせる前にscale()を使用してz-スコアに変換するか、lm()の代わりにlm.beta()を使用します。

XLStatで異なるベータ係数を取得する理由は明確ではありませんが、エラーでない場合は自由度と関係する可能性があります。 This exampleは、Rのlm.beta()関数をSASと比較し、同じ係数を取得します。

+0

ご返信ありがとうございます。 しかし、変数(レスポンス、カテゴリまたは両方)と関数lm.betaのスケーリングを含む多くのメソッドを試しましたが、XLStatとStatisticAと同じベータ係数を取得しません。 XlstatとStatisticaは私の以前のメッセージ 'beta = estimate * sd(x)/ sd(y)'に書いたものよりも別の式を使っていました。しかし、私はまだこの術式やlm.beta関数を使ってRで得た最初のベータ版はXLStatとStatisticaの最初のベータ版と同じだから、まだ驚いています... – Asselin

+0

さらに、私はすでにlmを比較したこの例を読んでいます。 RおよびSAS出力のβ関数しかし、それは定量的IVによる線形回帰に関係し、ANOVAにあるようにカテゴリー的ではない。 – Asselin

+0

ANOVAは連続変数の代わりにダミー変数を使用する線形回帰であり、係数はまったく同じ方法で計算されます。再現可能な実例がありますか? –

関連する問題