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

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

ユーザ用ツール

サイト用ツール


00.初めての医療統計rとezr:19.rでrepeatedmeasuresanova


差分

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

この比較画面へのリンク

両方とも前のリビジョン前のリビジョン
次のリビジョン
前のリビジョン
00.初めての医療統計rとezr:19.rでrepeatedmeasuresanova [2022/03/23] – [EZRでezr-repeated-measures-anova] adash33300.初めての医療統計rとezr:19.rでrepeatedmeasuresanova [2023/01/18] (現在) – [参考文献] adash333
行 9: 行 9:
 →分かりやすい! →分かりやすい!
  
 +
 +===== ソースコード =====
 +https://colab.research.google.com/drive/1XSBY6hiXNXOWUUmw_X3dSVCqDYwzq_3G?usp=sharing
 +220323_RepeatedMeasuredANOVA.ipynb
  
 ===== 対応のある3群以上の平均値の比較。正規分布に従い、等分散であり、球面性がある?場合のみ、Repeated Measures ANOVA ===== ===== 対応のある3群以上の平均値の比較。正規分布に従い、等分散であり、球面性がある?場合のみ、Repeated Measures ANOVA =====
行 20: 行 24:
 Rむけのデータ表を作る Rむけのデータ表を作る
 →pivot_longerとの違いは? →pivot_longerとの違いは?
 +
 +===== 正規分布の確認 =====
 +今回は、HbA1cの経過を時系列変化で有意に下がっているかどうかを検定します。
 +
 +{{ :00.初めての医療統計rとezr:book2.xlsx |}}
 +
 +{{:00.初めての医療統計rとezr:pasted:20220323-143106.png}}
 +
 +{{ :00.初めての医療統計rとezr:book2.csv |}}
 +
 +反復測定分散分析はパラメトリック検定
 +  -正規分布
 +  -等分散性
 +の確認が必要。
 +
 +各群が正規分布しているかどうかを、Kolmogorov-Smirnov検定する
 +
 +<code>
 +library(tidyverse)
 +d <- tribble(
 +  ~ID,~age,~sex,~DM,~HbA1c00M,~HbA1c01M,~HbA1c02M,~HbA1c03M,
 +  1,56,"M",1,7.3,7.4,7.3,7.4,
 +  2,77,"F",1,8.3,8.4,8.3,8.2,
 +  3,68,"M",1,7.1,7,7,7,
 +  4,81,"F",1,7.6,7.4,7.3,7.4,
 +  5,42,"M",1,10.1,8.4,7.3,6.5,
 +  6,65,"M",1,6.1,6.2,6.2,6.3,
 +  7,68,"F",1,7.8,7.4,7.9,8.4
 +)
 +d
 +</code>
 +
 +{{:00.初めての医療統計rとezr:pasted:20220323-150024.png}}
 +
 +
 +<code>
 +# データの取り込み
 +library(tidyverse)
 +d <- read_csv("book2.csv"
 +d
 +
 +# 正規性の検定 Kolmogorov-Smirnov test
 +# p>=0.05 であれば、正規分布に従っているとみなすことができる
 +vx <- d$HbA1c00M
 +ks.test(x=vx,y="pnorm",mean=mean(vx),sd=sd(vx))
 +</code>
 +
 +{{:00.初めての医療統計rとezr:pasted:20220323-145949.png}}
 +
 +https://data-science.gr.jp/implementation/ist_r_kolmogorov_smirnov_test.html
 +Rによるコルモゴロフ・スミルノフ検定
 +
 +
 +===== モークリーの球面性の検定 =====
 +Mauchly Tests for Sphericity
 +等分散性の確認。
 +
 +
 +
 +
 +https://www.r-bloggers.com/2021/04/repeated-measures-of-anova-in-r-complete-tutorial/
 +Repeated Measures of ANOVA in R Complete Tutorial
 +Posted by finnstats
 +
 +
 +
 +https://yusuke-ujitoko.hatenablog.com/entry/2018/11/29/000314
 +RでMauchly’s sphericity testを行う
 +2018-11-29
 +
 +https://www.datanovia.com/en/lessons/mauchlys-test-of-sphericity-in-r
 +Mauchly’s Test of Sphericity in R
 +
 +https://bellcurve.jp/statistics/glossary/2194.html
 +Mauchlyの球面性検定
 +Mauchly's sphericity test
  
  
行 78: 行 158:
 Rで一元配置分散分析(対応あり) Rで一元配置分散分析(対応あり)
 投稿日時: 2012-08-13 投稿日時: 2012-08-13
 +
 +http://babybear.site/anova.html
 +分散分析(analysis of variance: ANOVA)
 +ANOVAとは
 +
 +https://mikuhatsune.hatenadiary.com/entry/20160530/1464615103
 +時系列データにt 検定を行う
 +20160530
  
  


00.初めての医療統計rとezr/19.rでrepeatedmeasuresanova.1648031085.txt.gz · 最終更新: 2022/03/23 by adash333

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki