「統計学関連なんでもあり」の過去ログ--- 047

No.21624　多重性の調整　　【赤岳】　2015/04/25(Sat) 17:52

線形回帰モデルlmを用いて，Drug1～4を下記のように検定しました。
この結果は多重性が調整されていないと思うんですが，lmでは指定できないんでしょうか。Dunnettで調整する方法はどうすればよいでしょうか。

> LM<-lm(value~Drug, data=Dataset)
(一部省略）
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 791.500 9.869 80.198 < 2e-16 ***
Drug[T.Drug2] 70.400 13.957 5.044 1.31e-05 ***
Drug[T.Drug3] 134.500 13.957 9.636 1.65e-11 ***
Drug[T.Drug4] 151.500 13.957 10.854 6.64e-13 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 31.21 on 36 degrees of freedom
Multiple R-squared: 0.8024, Adjusted R-squared: 0.786
F-statistic: 48.74 on 3 and 36 DF, p-value: 9.252e-13
ご教示よろしくお願いします。

No.21625　Re: 多重性の調整　　【青木繁伸】　2015/04/25(Sat) 20:28

どのようなデータなんでしょうか。
重回帰分析出来るデータですか？

No.21628　Re: 多重性の調整　　【赤岳】　2015/04/25(Sat) 22:13

青木先生，いつもご教示ありがとうございます。
Dataは，次のとおりです。ダミー変数を使えば重回帰分析は可能かと思います。
よろしくお願いします。
Drug gender value
Drug1 m 774
Drug1 m 776
Drug1 f 816
Drug1 f 808
Drug1 f 773
Drug1 f 707
Drug1 m 774
Drug1 m 819
Drug1 f 852
Drug1 f 816
Drug2 f 848
Drug2 f 875
Drug2 m 858
Drug2 f 823
Drug2 f 858
Drug2 m 891
Drug2 m 843
Drug2 m 870
Drug2 m 872
Drug2 m 881
Drug3 f 925
Drug3 f 917
Drug3 m 912
Drug3 f 912
Drug3 m 949
Drug3 m 908
Drug3 f 908
Drug3 f 989
Drug3 f 931
Drug3 m 909
Drug4 f 919
Drug4 m 969
Drug4 m 924
Drug4 m 902
Drug4 f 880
Drug4 f 956
Drug4 m 964
Drug4 f 948
Drug4 m 983
Drug4 f 985

No.21629　Re: 多重性の調整　　【青木繁伸】　2015/04/25(Sat) 22:50

独立なデータですね。
一般的に，多変量解析で相関係数行列や偏回帰係数などの多重性は考えないことが普通だと思います。
このデータの場合，重回帰分析を行う必要があるんでしょうか？直接ダネットの検定を行えば良いのではないでしょうか？

No.21630　Re: 多重性の調整　　【赤岳】　2015/04/26(Sun) 00:30

青木先生，ご回答ありがとうございます。
そうですね。独立変数を増やしたかったので，lmでないとと思いましたが，multcompのglhtを使えば，複数の独立変数を組み込んで直接ダネットの検定を行うことも確認できました。ありがとうございました。