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

統計学備忘録(R言語のメモ)

since2016 ときどきTEXのメモ

カイ自乗近似は不正確かもしれません

このコメントが出る場合は以下のようなときです 期待度数が 0になるセルがある場合、もしくは期待度数が 5未満になるセルが全体の 20% を超える場合 R: Fisher's Exact Test for Count Data that is if no cell has expected counts less than 1 and more th…

統計検定2級

以下のサイトに移行しました y2pt.com

決定係数(寄与率)

下記のサイトに移転しました y2pt.com

テキストファイルに検定結果を出力

下記のサイトに移転いたしました y2pt.com

九州理学療法士・作業療法士合同学会2017in宮崎で紹介したスクリプト

モンテカルロ法によるt検定の検定力分析 モンテカルロ法については以下を参照ください yoshida931.hatenablog.com tval <- numeric(length = 10000) #t値(空のベクトルを設定) cs <- 0 for(i in 1:10000){ x <- rnorm(n = 30, mean = 10, sd = 3) # xの乱…

独立性の検定

下記のサイトに移転しましたy2pt.com

統計学的に独立ということ

以前の独立という記事を削除して再投稿二つの事象A,Bが独立であるということ 一方の事象が起こるかどうかが、他方の事象の起こる確率に影響しないこと 2つの事象AとBがを満たす場合、事象AとBは統計的に独立である.このときAとBの余事象も独立であるまたAの…

二項分布からの最尤推定 (旧)

二項分布からの最尤推定 投稿日2017.6.13 更新日2017.10.4 理解してなかったのでもう一度勉強します二項分布の復習 確率1/2が最も尤もらしい例 右側の数値は同時確率 どう見てもP=0.5が最も尤もらしい確率です!最尤原理 最尤原理は、現在起きている事象の起…

配列 (array) の基本

配列 (array) の基本 投稿日2017.10.4リスト:異なるデータ構造をまとめて一つのオブジェクトとして保存.解析結果の多くはリスト形式. 配列:次元を自由に設定できるオブジェクト.RのサンプルHairEyeColorを通して勉強していきます > HairEyeColor , , Se…

Rのサンプル ~ データフレーム編

Rのサンプル ~ データフレーム編 投稿日2017.10.2 更新日2017.10.4 少し見やすくしましたこういうのが統計学の学習に役立つかも? サンプルの詳細を知りたいときはここ 間瀬茂 R 基本統計関数マニュアル https://cran.r-project.org/doc/contrib/manuals-jp…

Z検定

下記のサイトに移転しましたy2pt.com

対応のないt検定の検定力分析(事後の分析)

検定力(1-β):帰無仮説が「偽」であるとき(母集団に差があるとき)にサンプルから有意差を得る確率 効果量(effect size):標準化された平均値差(検定により色々あります...d family, r family) 標本数:研究者が決定しなければならない α(有意水準):慣例的…

2標本 平均の差の検定

下記のサイトに移転しましたy2pt.com

指数分布

下記のサイトに移転しましたy2pt.com

1標本t検定のまとめ

下記のサイトに移転しましたy2pt.comy2pt.com

確率点と累積分布

確率点と累積分布 投稿日2017.7.22 更新日2017.9.1lower.tailのdefaultはTRUE この意味が重要です.勉強のためにdefaultのTRUEも挿入しときます. lower.tail logical; if TRUE (default), probabilities are P[X ≤ x] otherwise, P[X > x]. Rのhelpより 正規…

ブログの再構成

自分が書いたブログなのに「見にくい&理解できない」 という悩みが募ってきました. よって、少しずつ再構成することにしました. はてな記法なるものに初挑戦! 見直したとき、せめて・・・ せめて書いた本人が理解できるような内容に改めたいと思っており…

モンティ・ホール問題

モンティ・ホール問題 投稿日2017.8.30モンティ・ホール問題、モンティ・ホール・ジレンマ Rを使って何とか解いてみようと思います. 3つのドアのうち1つのドアの後ろには高級車(車の扉)があり、残りの2つのドアにはヤギ(ヤギの扉)がいます。司会者は…

モンテカルロ法で円周率を求める

モンテカルロ法で円周率を求める 投稿日2017.8.28 モンテカルロ法 対象となるある現象がいくつかの確率変数の関数で表現できる場合(数値モデル)、各変数で仮定される確率分布に沿った標本値を大量に生成(乱数)して、その計算結果の分布を推定(結果の推…

平均と分散

以下のサイトに引っ越しました y2pt.com

二元配置分散分析(各水準の信頼区間推定)

二元配置分散分析(各水準の信頼区間推定) 投稿日2017.8.24 各要因における各水準の母平均の区間推定 繰り返し数 n n = 繰り返し数 要因1水準Aの有意水準αの信頼区間は、xai ± t α/2 (自由度) * 標準誤差 となります. 母分散が未知なので、標準誤差を残差…

標準偏回帰係数の求め方

下記のサイトに移転いたしました y2pt.com

トービット回帰直線

トービット回帰直線 投稿日2017.8.4 打ち切りデータの場合、つまり天井効果や床効果が生じているデータの場合には、トービット回帰直線で分析します. 参考書はもちろん豊田秀樹 (著, 編集);回帰分析入門 (Rで学ぶ最新データ解析) ,東京図書 ,2012 Rのサンプ…

残差分析

下記のサイトに移転しました y2pt.com

回帰係数の区間推定

下記のサイトに移転しました y2pt.com

ユールのQ

ユール ( Yule ) のQユールの関連係数、ユールの連関係数とも呼ばれています.関連が強いほど1または-1に近い値をとります。 Q= ( a*d - b*c ) / ( a*d + b*c ) = ( オッズ比 - 1 ) / ( オッズ比 + 1 )例)2つの質問 ( q1 , q2 ) に関する答え ( yes , no ) …

適合度の検定 (ピアソンの 𝜒^2 適合度検定)

下記のサイトに移転しました y2pt.com

Rstudioでの data.frame 表示の不具合

データフレームを表示したときに、ラベルが見えなくなりました RStudio-1.0.143 から RStudio-1.0.153 へバージョンアップすることで不具合は解消しました.https://www.rstudio.com/products/rstudio/download/preview/ Rstudioのサポートページの掲示板で…

ファイ係数

連関係数 クロス集計表における2つの変数間の関連性の程度を表す指標として,連関係数が提案されています. 連関係数<|0.2| 連関はほとんどない|0.2|≦連関係数<|0.4| やや連関がある|0.4|≦連関係数<|0.7| 強い連関がある|0.7|≦連関係数 かな…

イェーツの補正

イェーツの補正 / イェーツの連続修正(Yate's continuity correction) 離散型分布を連続型分布に近似させて統計的検定を行う際に使用する修正です.2×2分割表のデータに対して行われ、より正確な検定が可能になります.「連続的なカイ二乗分布(自由度1のχ…