差分

このページの2つのバージョン間の差分を表示します。

--- 00.初めての医療統計rとezr:14.rでone-wayanova [2021/06/14] – [等分散性がない場合または気にしない場合はWelch Test(Welchの分散分析)] adash333
+++ 00.初めての医療統計rとezr:14.rでone-wayanova [2021/06/14] (現在) – [one-way ANOVAリンク] adash333
@@ 行 70: / 行 70: @@
 (引用ここまで)
-===== 等分散性がない場合または気にしない場合はWelch Test(Welchの分散分析) =====
+http://stack.sakura.ne.jp/stat/index.php?%E7%90%83%E9%9D%A2%E6%80%A7
+球面性　sphericity
+> ３回以上繰り返しのある分散分析を行うには、球面性の仮定が保たれてないといけない。
+===== 等分散性がない場合または気にしない場合はWelch Test(Welchの分散分析) Welch ANOVA =====
+https://www.statology.org/welchs-anova-in-r/
+ APRIL 12, 2021 BY ZACH
+How to Perform Welch’s ANOVA in R (Step-by-Step)
+https://www.reddit.com/r/statistics/comments/ay7zsf/welchs_anova_post_hoc_in_r/
+Posted byu/Half_Flatworm
+years ago
+Welch's ANOVA post hoc in R
+→ R 4.0.0では、install.packages("biostat") ができず。
 https://stackoverflow.com/questions/28587498/post-hoc-tests-for-one-way-anova-with-welchs-correction-in-r
 Post-hoc tests for one-way ANOVA with Welch's correction in R
 Asked 6 years, 3 months ago
-Active 5 years, 11 months ago
-Viewed 10k times
+https://rdrr.io/github/GegznaV/BioStat/man/posthoc_anova.html
+posthoc_anova: *!!* Post-hoc tests for ANOVA and Welch ANOVA
@@ 行 106: / 行 125: @@
+===== Welch ANOVAでpost hoc test =====
+Welch ANOVAを用いる（oneway.test(y ~ group, data=data, var.equal=FALSE)）のは、『正規分布』であるが、『等分散ではない』とき。（等分散でも用いてもよい。ぶっちゃけ、等分散は調べないで、最初からone way ANOVAではなく、Welch ANOVAで計算する方法もあり。）
+https://statisticsbyjim.com/anova/welchs-anova-compared-to-classic-one-way-anova/
+Benefits of Welch’s ANOVA Compared to the Classic One-Way ANOVA
+By Jim Frost
+> The Games-Howell post hoc test, like Welch’s analysis of variance, does not require the groups to have equal standard deviations. Conversely, Tukey’s method does require equal standard deviations.
+http://www.f.kpu-m.ac.jp/y/math/Med_Stat/files/stat4h.pdf
+の、スライドp25
+(引用ここから)
+{{:00.初めての医療統計rとezr:pasted:20210614-150612.png}}
+(引用ここまで)
+http://aoki2.si.gunma-u.ac.jp/R/oneway-ANOVA.html
+一元配置分散分析と多重比較　　　　　Last modified: Apr 09, 2015
+<code>
+# インストールは，以下の 1 行をコピーし，R コンソールにペーストする
+source("http://aoki2.si.gunma-u.ac.jp/R/src/tukey.R", encoding="euc-jp")
+</code>
+関数のインストールが必要。
+<code>
+source("http://aoki2.si.gunma-u.ac.jp/R/src/tukey.R", encoding="euc-jp")
+</code>
+使用例
+<code>
+x <- c(
+, 206, 164, 190, 194, 203,
+, 221, 197, 185,
+, 265, 197, 220, 212, 281,
+, 276, 237, 254, 230
+  )
+g <- rep(paste("A", 1:4, sep=""), c(6, 4, 6, 5))
+data <- list(x,g)
+oneway.test(x~g, data=data, var.equal=FALSE)
+</code>
+{{:00.初めての医療統計rとezr:pasted:20210614-153017.png}}
+p<0.05 なので、Welch ANOVAで4群の平均値に有意差あり。
+→ Games-Howellで、post hoc Testを行う。
+<code>
+# tukey()関数を用いるためには、source("http://aoki2.si.gunma-u.ac.jp/R/src/tukey.R", encoding="euc-jp")　が必要
+tukey(x, g, method="Games-Howell")
+</code>
+{{:00.初めての医療統計rとezr:pasted:20210614-153329.png}}
+→　しかし、今回は、全ての組み合わせのp > 0.05 なので、どの2群にも有意差は無し。