Bio Technical フォーラム

  • バイオ関連の実験をする上での、試薬、機器、プロトコールなどの情報交換の場です。
  • 新しいテーマで話を始める場合、質問をする場合は「新しいトピックを作る」から書き込みをしてください。
  • 質問に対して解答できる方は是非、書き込んで下さい。
  • このフォーラムにふさわしくないと管理人が判断した投稿は予告なく削除します。

新しいトピックを作る | トピック一覧 | 研究留学ネットに戻る

ひとつ前のフォーラム(readのみ)

このスレッドをはてなブックマークに追加このスレッドをはてなブックマークに追加

ANOVA の結果 トピック削除
No.4922-TOPIC - 2016/03/16 (水) 11:41:47 - おお
ちょっと結果の見かたについて教えてください
例えばRでOne way anovaをかけた結果

> summary(AnovaModel.1)
Df Sum Sq Mean Sq F value Pr(>F)
factor(var1) 2 3235 1617 6.367 0.0266 *
Residuals 7 1778 254
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1


でPost hocは
> pairwise.t.test(Dataset$var2, Dataset$var1, p.adj="bonferroni")

Pairwise comparisons using t tests with pooled SD

data: Dataset$var2 and Dataset$var1

1 2
2 1.000 -
3 0.043 0.088

P value adjustment method: bonferroni

> TukeyHSD(AnovaModel.1, "factor(var1)")
Tukey multiple comparisons of means
95% family-wise confidence level

Fit: aov(formula = var2 ~ factor(var1), data = Dataset, na.action = na.omit)

$`factor(var1)`
diff lwr upr p adj
2-1 6.217252 -32.108606 44.54311 0.8837720
3-1 39.490291 3.639733 75.34085 0.0333967
3-2 33.273039 -2.577519 69.12360 0.0668032

Multiple Comparisons of Means: Dunnett Contrasts


Fit: aov(formula = var2 ~ group.factor, data = Dataset)

Linear Hypotheses:
Estimate Std. Error t value Pr(>|t|)
2 - 1 == 0 6.217 13.014 0.478 0.849
3 - 1 == 0 39.490 12.173 3.244 0.025 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
(Adjusted p values reported -- single-step method)

質問はDunnettではSignif. code *をつけてくれているんですが、その他の検定方法ではP値が0.05以下でもSignif. code がついていません。これは有意差がないという結論ではないですよね。。。
 
- このトピックにメッセージを投稿する -



13件 ( 1 〜 13 )  前 | 次  1/ 1. /1


(無題) 削除/引用
No.4922-13 - 2016/03/25 (金) 12:01:19 - おお
>[Re:12] CDさんは書きました :
> R以外のソフトでの挙動はわからないですが、RでのTukey検定で有意差のついたP valueにマークが付かないのは仕様だと思います。

ありがとうございました。仕様ということですね。なにかこの仕様にしたほうがいい理由とかあったんでしょうか、それとも単にないということなんでしょうか、、、

(無題) 削除/引用
No.4922-12 - 2016/03/24 (木) 11:22:45 - CD
R以外のソフトでの挙動はわからないですが、RでのTukey検定で有意差のついたP valueにマークが付かないのは仕様だと思います。

(無題) 削除/引用
No.4922-11 - 2016/03/24 (木) 09:38:26 - おお
>[Re:10] 別の通りすがりさんは書きました :
> Rは勉強中で詳しくはないのですが
>
> >その他の検定方法ではP値が0.05以下でもSignif. code がついていません。
>
> そもそも、件のパッケージにSignif. codeを付与する機能がついていないということはないでしょうか?
>

ちょっとそれはわかりません、なので詳しい方がいればなぁとおもっています、ただ違うソフトでもなんか同じようなことが起こっているようなので、、、今手元で分析できませんが、だいぶ後になりますけど、そういう状況の説明をアップしたいと思います。

(無題) 削除/引用
No.4922-10 - 2016/03/22 (火) 09:44:09 - 別の通りすがり
Rは勉強中で詳しくはないのですが

>その他の検定方法ではP値が0.05以下でもSignif. code がついていません。

そもそも、件のパッケージにSignif. codeを付与する機能がついていないということはないでしょうか?

(無題) 削除/引用
No.4922-9 - 2016/03/21 (月) 00:32:55 - おお
ありがとうございました。もしかしたら統計専門のQAサイトなどで相談したほうがいいような気がしてきました。でもちょっと時間的余裕が無いので一週間後とかになりそうです。。。

(無題) 削除/引用
No.4922-8 - 2016/03/21 (月) 00:30:33 - 通りすがり
と思って閉手順を見直してみたんですが、1回目の有意水準は重要でも、p値は重要ではないですね。ANOVAだけに何かあるんでしょうか? よく分かりません。

(無題) 削除/引用
No.4922-7 - 2016/03/21 (月) 00:04:24 - 通りすがり
確かにTurkyのp adjでも星がついてないですね。。。こう考えてみましたがどうでしょう。
私の理解ではpost hoc testは閉手順で実施していると考えています。なので、2回目以降の検定は最初の検定より強く検出されないと(p値が小さくならないと)閉手順が成り立たない。
表題の場合では最初のANOVAのp=0.0266なので、post hoc testの調整済みp値は0.0266より小さくないと閉手順が成り立たない。すなわち、有意さがあるとはいえない。
私の使っている統計ソフトでやって見たらANOVA p=0.027、Tyrky p=0.037で星がついていたんで、違うかもしれませんが。。。

(無題) 削除/引用
No.4922-6 - 2016/03/19 (土) 23:46:57 - おお
私も基本手計算なんですけどね。と言いますのはpaired, one sampleなどのテストもありANOVAかけるのがやりにくいものもあるので。

実は、これRではなくって他のソフトで見つかったんですけど、RとかのほうがポピュラーでOutputとかわかりやすいかなと思ったしだいです。

補正に関しては> pairwise.t.test(Dataset$var2, Dataset$var1, p.adj="bonferroni")とありますのでかかっているんじゃないかと思いますが、かかってないとしますと、他の方法(TukeyHSDをトピ上げたとき入れてますが)は自分で補正しないといけないんですかね。。。TukeyHSDもアウトプットはp.adjと表記されています。ちなみにTukeyHSDでもp.adj 0.0333967でSignif. codes *がついてません。

あとちょっとわかってないのですがPost hocのばあいAnovaの計算値を使い結果が通常のbonferroni 補正と変わるってことはないんでしょうか、何かちょっと数値が違ってるような気がしてます、、、

(無題) 削除/引用
No.4922-5 - 2016/03/19 (土) 22:13:41 - 通りすがり
Rは知らないのですが。。。

p=0.043ほ補正済みのp値なのですか? 私が使用している統計ソフトではp値は補正なしで、有意水準が補正されます。例えば今回t testを2回行っているとしたら、p=0.043(補正なし)と有意水準0.05/2=0.025を比較するので、有意さは無しとなります。
ttestとボンフェローニ補正だけなので、私ならとりあえずエクセルで確認します。

(無題) 削除/引用
No.4922-4 - 2016/03/18 (金) 22:28:21 - おお
ありがとうございました。それは存じ上げているのですが、何が起こっているのか不思議なんです。bonferroniでやろうとしていますが、bonferroniの結果の1と3の比較でp adj=0.043で0.05より小さいのですが、Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1のなかの*がついていません。なぜでしょうというかという質問です。p adj=0.043で有意差なしということなんでしょうか。

統計処理前に決めておくべきでは 削除/引用
No.4922-3 - 2016/03/18 (金) 15:07:04 - PD
どの多重比較で検定するかは、どういうサンプルかによってあらかじめ決めておいて統計処理するもので、全部やってみてこれで有意差がでたからというものではないと思いますが…

http://www.ibaraki-kodomo.com/toukei/posthoc.html
このサイトに多重比較法の比較があるので参考になるかもしれません。

(無題) 削除/引用
No.4922-2 - 2016/03/16 (水) 14:36:16 - おお
ありゃりゃみにくいですね。。。あとで修正しに来ます。。。

ANOVA の結果 削除/引用
No.4922-1 - 2016/03/16 (水) 11:41:47 - おお
ちょっと結果の見かたについて教えてください
例えばRでOne way anovaをかけた結果

> summary(AnovaModel.1)
Df Sum Sq Mean Sq F value Pr(>F)
factor(var1) 2 3235 1617 6.367 0.0266 *
Residuals 7 1778 254
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1


でPost hocは
> pairwise.t.test(Dataset$var2, Dataset$var1, p.adj="bonferroni")

Pairwise comparisons using t tests with pooled SD

data: Dataset$var2 and Dataset$var1

1 2
2 1.000 -
3 0.043 0.088

P value adjustment method: bonferroni

> TukeyHSD(AnovaModel.1, "factor(var1)")
Tukey multiple comparisons of means
95% family-wise confidence level

Fit: aov(formula = var2 ~ factor(var1), data = Dataset, na.action = na.omit)

$`factor(var1)`
diff lwr upr p adj
2-1 6.217252 -32.108606 44.54311 0.8837720
3-1 39.490291 3.639733 75.34085 0.0333967
3-2 33.273039 -2.577519 69.12360 0.0668032

Multiple Comparisons of Means: Dunnett Contrasts


Fit: aov(formula = var2 ~ group.factor, data = Dataset)

Linear Hypotheses:
Estimate Std. Error t value Pr(>|t|)
2 - 1 == 0 6.217 13.014 0.478 0.849
3 - 1 == 0 39.490 12.173 3.244 0.025 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
(Adjusted p values reported -- single-step method)

質問はDunnettではSignif. code *をつけてくれているんですが、その他の検定方法ではP値が0.05以下でもSignif. code がついていません。これは有意差がないという結論ではないですよね。。。

13件 ( 1 〜 13 )  前 | 次  1/ 1. /1


パスワードを入力してチェックした記事を チェックした記事を

このトピックにメッセージを投稿する
名前 
メール   アドレス非公開
   タイトル 
本文      
設定  クッキーを保存(次回の入力の手間を省けます)
上に上げない(トピックの一覧で一番上に移動させません)
解決(問題が解決した際にチェックしてください)
暗証  半角英数字8-12文字の暗証番号を入れると、あとで削除、修正ができます。
送信 

〔使い方〕
  • 「アドレス非公開」をチェックすれば、自分のメールアドレスを公開しないで他の方からメールを受け取れます。
  • 問題が解決した際には、解決ボタンをチェックして解決した旨のコメントをつけてください。これは、初めにトピックを作った人と管理人のみが可能です。
  • 半角カタカナ、機種依存文字(全角ローマ数字、○の中の数字等)は文字化けの原因となりますので使わないでください。