サルでもわかるWEBプログラミング

フリーソフトのみでホームページ作成

ユーザ用ツール

サイト用ツール


サイドバー

目次

サルでもわかる機械学習

このページへのアクセス
今日: 1 / 昨日: 0
総計: 6

sidebar

00.初めての医療統計rとezr:14.rでone-wayanova


文書の過去の版を表示しています。


14.Rでone-wayANOVA

参考ページ

http://www.sthda.com/english/wiki/one-way-anova-test-in-r#visualize-your-data-and-compute-one-way-anova-in-r
One-Way ANOVA Test in R

http://mizumot.com/handbook/?page_id=219
Rを使った分析(一元配置分散分析)
→多重比較の計算方法あり

library(tidyverse)

# Compute the analysis of variance
res.aov <- aov(weight ~ group, data = my_data)
# Summary of the analysis
summary(res.aov)

上記でp<0.05であった場合、3群の平均値に有意差があると言えて、以下の多重比較に進むことができる

attach(data)

#Tukey
TukeyHSD(res.aov) 
TukeyHSD(aov(weight~group)) 

#Holm
pairwise.t.test(res.aov)
pairwise.t.test(weight, group)

# bonferroni ← 迷ったらこれ
pairwise.t.test(res.aov, p.adjust.method="bonferroni")
pairwise.t.test(weight, group, p.adjust.method="bonferroni")

one-way ANOVAリンク

https://foxglovetree.wiki.fc2.com/m/wiki/Programming.R.pairwise.t.test
Kohji 2020-04-19

attach()を利用している

https://data-science.gr.jp/implementation/ist_r_oneway_anova.html
anova関数とaov関数

https://stats.biopapyrus.jp/stats/anova.html
anova関数とaov関数

https://qiita.com/FukuharaYohei/items/f6bbc10675eaa29ba18a
@FukuharaYohei
2018年01月21日に更新
【R入門】R言語の基本:一元配置分散分析(対応なし)

https://www.stats-guild.com/analytics/15634
1元配置分散分析(One-way ANOVA)
公開日:2020/03/18 最終更新日:2020/06/17

(引用ここから)
1元配置分散分析は、独立したサンプルのt検定と同じように、量的従属変数が正規分布にしたがう「正規性」と、各水準のばらつきが同程度である「等分散性」の仮定を持ち、これらを満たすことができない場合は別の検定手法を適用します。正規性を満たさない場合は、クラスカル-ウォリス検定 Kruskal-Wallis testを適用し、等分散性を満たさない場合はウェルチ検定 Welch test がよく用いられます。
(引用ここまで)

等分散性がない場合または気にしない場合はWelch Test(Welchの分散分析) Welch ANOVA

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

https://www.yodosha.co.jp/jikkenigaku/statistics/q3.html
Q3 3つ以上の群の差を調べるにはどうしたらよいですか?t検定は使えないのですか?
荻島創一:東北大学東北メディカル・メガバンク機構医療情報ICT部門
※基本編 4章 実験の目的に合った検定の選び方・実験計画Q36を抜粋(2012.10.02掲載)

https://www.ai.u-hyogo.ac.jp/~arima/lectures/JT-9.pdf
統計解析法の選定
post hoc

https://rstudio-pubs-static.s3.amazonaws.com/393500_6ae9e0f0ebb54d2f9c42c5e355c3668a.html
第7章 分散分析 1 of 2 Liu 2018/6/1

http://aoki2.si.gunma-u.ac.jp/lecture/mb-arc/arc043/12857.html
No.12857 「welchの方法によるt検定」と「welchの分散分析」  【ホンダ】 2010/06/17(Thu) 22:58

https://oku.edu.mie-u.ac.jp/~okumura/stat/ttest.html
t 検定

https://www.researchgate.net/post/How_do_I_choose_a_post_hoc_test_when_equal_variances_are_not_assumed_in_SPSS
Asked 17th Feb, 2016
Teresa Diaz Calvo
University of East Anglia
How do I choose a post hoc test when equal variances are not assumed in SPSS?

attachとdetach

http://cse.naro.affrc.go.jp/takezawa/r-tips/r/46.html

データフレームの要素を参照するときは『データフレーム名$要素のグループ』と $ を用いて参照する。

# data$score を表示したいとき
attach(data)
score

これで、data$scoreが表示される。
解除するときは、以下のようにする

# dataのattachを解除するとき
detach(data)

リンク


00.初めての医療統計rとezr/14.rでone-wayanova.1623674209.txt.gz · 最終更新: 2021/06/14 by adash333