理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

分散分析

エラーバー付きのグラフ ggplot2

使用するパッケージ install.packages("ggplot2") library(ggplot2) データセットを作成します #乱数を使用のため、実施毎にグラフは異なります dat <- c(c(rnorm(15)),c(rnorm(15,3,1.2)),c(rnorm(15,10,2)),c(rnorm(15,17,1.5))) pre_post <- c(rep("前", …

対応のある一元配置分散分析(繰り返しのない二元配置分散分析)、多重比較

投稿日2018.7.21, 最終更新日2019.6.4 以下のデータセットを使用します x1 <- c(93, 89, 115, 90, 75) x2 <- c(121, 136, 121, 161, 125) x3 <- c(101, 115, 118, 122, 106) 一要因被験者内計画で対応のある一元配置分散分析です また、要因1=ID, 要因2=検…

級内相関係数 Case1

以前掲載していたのですが、以下の文献を参考に書き直しました。 主にRのプログラム (パッケージは使用してません) を忘れないようにメモしておきます。 理論は以下の文献を参照してください。 1) SHROUT, Patrick E.; FLEISS, Joseph L. Intraclass correla…

共分散分析を考えるためのグラフ

いつも勉強させていただいおります、下記サイトの共分散分析より グラフの描き方の備忘録です 我楽多頓陳館--雑学と統計学の館 以下のようなグラフになります データのサンプルは上記サイトから引用してdataに入れました # まずはこれを描いて考える・・・!…

分散分析の基本

最終更新日:2018.3.5 まだ理解できていない.なので書き直し・・・ 一元配置分散分析 言葉の整理 要因(factor), 因子(factor):実験結果に影響を与える要素.それぞれの分野で使い分ける場合もあるので注意.このブログでは要因と因子の区別をせず「要因」で統…

分散分析(繰り返し vs 反復)

月元敬;対応あり, 繰り返し, 反復: 心理データ解析基礎において区別すべき 3 つの 「R」. 岐阜大学教育学部研究報告. 人文科学, 2015, 64.1: 67-73. 「繰り返し」と「反復」は区別なく混同している参考書も多いようですが、今回は上記の文献をもとに違いを整…

二元配置分散分析 (2)

二元配置分散分析 (2) 投稿日2017.8.22 分析後の解釈について勉強します 二元配置分散分析(1)より 分散分析表 平方和の分解∑(f)∑(g)∑(n) ( xij - xa )^2 = ∑(f) g*n*( xai, - xa ) ^2 Aの主効果 + ∑(g) f*n*( xa-j - xa ) Bの主効果 + ∑(f)∑(g)n* ( xaij -…

二元配置分散分析 (1)

二元配置分散分析 投稿日2017.3.8更新日2017.8.17 主効果:因子の効果がその水準によって異なるとき主効果があるという 交互作用:別の因子がどのようなものかによって効果が異なるとき因子間に交互作用があるという 要因効果:主効果と交互作用の総称 2要因…

なぜt検定を繰り返してはいけないのか?

なぜt検定を繰り返してはいけないのか投稿日2017.2.24修正日2017.7.9 sample は「石村卓夫;分散分析のはなし,東京図書,1992,p137」 A<-c(12,14,16)B<-c(13,15,17)C<-c(15,17,19)D<-c(17,19,21) x<-c(rep("A",3),rep("B",3),rep("C",3),rep("D",3))y<-c(A,B,…

一元配置分散分析 (対応なし) F値の算出方法

更新日2017.3.13更新日2017.6.22 一元配置分散分析とは、一つの因子に基づく“3つ以上の母平均の差の検定”です 因子Factor:実験結果に影響を与えると考えられる要因水準Lebel:因子をいくつかに分ける基準 (水準=群) 例題)治療因子に基づく、患者20名の…

Rで簡単(反復測定による一元配置分散分析、データの対応あり、繰り返しのない二元配置分散分析)

方法は「Rで簡単(一元配置分散分析:対応なし)」と同じです 最下部の式が少し違うだけです サンプル 人 順 得点野原 1回目 12野原 2回目 8野原 3回目 9風間 1回目 8風間 2回目 5風間 3回目 8桜田 1回目 12桜田 2回目 10桜田 3回目 10佐藤 1回目 6佐藤 2回…

Rで簡単(一元配置分散分析:データの対応なし、繰り返し5回)

サンプルをコピーして(エクセルからでもOK)、下の2式を実行するのみ サンプル 治療方法 点数A 15A 9A 18A 14A 18B 13B 8B 8B 12B 7C 10C 6C 11C 7C 12D 10D 7D 3D 5D 7 x<-read.table("clipboard",header=T) 注)x$治療方法 が数字のときにはfactorに変換…