2017/01/17
2020/05/04
Rでフィッシャーの正確確率検定 そのまま使える自作関数例
独立性の検定を行うときに、データ数が少ないときや0になるセルがあるときは、カイ二乗検定のカイ二乗分布への近似が不正確になる場合があります。そういうときに使うのが、フィッシャーの正確確率検定です。Rでのカイ二乗検定はこちら→Rで独立性のカイ二乗検定 そのまま使える自作関数
今回は、R言語でフィッシャーの正確確率検定を行なって、その結果をcsvファイルに出力する自作関数を作りましたので、お役立ていただければと思います。
fisher_fun <- function(data,filename)
{
ans <- fisher.test(data)
p_value <- ans[1]
odds_low <- ans$conf.int[1]
odds_high <- ans$conf.int[2]
oddsrate <- ans[3]
N <- data[1]+data[2]+data[3]+data[4]
result = cbind(data,p_value,odds_low,odds_high,oddsratio,N)
result[2:nrow(result),c(3:7)]=“”
write.table(matrix(c(“”,colnames(result)),nrow=1),filename,append=T,quote=F,sep=“,”,row.names=F,col.names=F)
write.table(result,filename,append=T,quote=F,sep=“,”,row.names=T,col.names=F)
write.table(“”,filename,append=T,quote=F,sep=“,”,row.names=F,col.names=F)
}
fisher_fun(データ,出力ファイル名)
として使うようにしてください。データは、分割表の形にして入れるようにしてください。
出力結果は、
元データ、p値、オッズ値の95%信頼区間、オッズ比、データの大きさ
の順番になります。
R言語について、さらに学びたい方はR言語入門をご覧ください。
最新投稿記事
-
【独占】コインバンクの卜部代表、既存案件への影響限定的ー新型コロナで
2020年7月3日
-
【独占】データサイエンスコンペ参加者、３年後に１万人超目指すーNishikaの松田CTOインタビュー
2020年7月2日
-
【独占】純真学園大学・村中光副学長「G検定、全受講者に取得してもらいたい」ー「人工知能学」新設
2020年7月1日
-
OKI、９月下旬に「AIエッジ・コンテスト」ー「AE2100」で課題解決力競う
2020年6月30日
-
G検定に数学力はどれくらい必要？文系でも受かるのか？
2020年6月26日
-
G検定に落ちた人、合格した人。勉強法の違いはどこにある？
2020年6月25日
-
G検定の難易度・合格ライン・内容・評判を徹底調査！
2020年6月19日
-
コニカミノルタ、パートナーと画像IoT開発推進、実践でAI人材育成ー高槻サイトで
2020年6月19日
-
「第１回全国高専ディープラーニングコンテスト」本選を８月22日に開催ーJDLA
2020年6月18日
-
G検定、94%が「誰かに勧めたい」―JDLA、1400人アンケート調査
2020年6月16日
週間ランキング
G検定の対策方法まとめ。本、講座、過去問など徹底調査。 2020年6月6日 
GAFA３社が採用内定の理系日本人へ「年収6000万円」のオファー、「AI人材の争奪… 2020年5月1日 
G検定に落ちた人、合格した人。勉強法の違いはどこにある？ 2020年6月25日 
G検定に数学力はどれくらい必要？文系でも受かるのか？ 2020年6月26日 
OKI、９月下旬に「AIエッジ・コンテスト」ー「AE2100」で課題解決力競う 2020年6月30日 
G検定の難易度・合格ライン・内容・評判を徹底調査！ 2020年6月19日 
【独占】データサイエンスコンペ参加者、３年後に１万人超目指すーNishikaの松田CTOインタビュー 2020年7月2日 
【独占】純真学園大学・村中光副学長「G検定、全受講者に取得してもらいたい」ー「人工知能学」新設 2020年7月1日 
E資格合格談、「コード読解力ついた」ーVSNのデータサイエンティスト・原さん 2020年6月12日 
落合陽一氏がアプリ「音声文字変換」を熱く語るーツイート動画再生15万回超 2020年5月14日
この記事を読んだ方はこちらの記事も読まれています
- 2016年9月13日 Rでデータ集計に便利なコマンド集と列ごと集計の関数例
- 2017年5月15日 RでYouden Indexを用いたROC曲線の最適カットオフ
- 2016年11月22日 Rでカイ二乗分布のグラフを描く方法
- 2017年2月6日 Rで散布図をプロットしてpngで保存する自作関数
- 2016年11月21日 Rによる条件付きロジスティック回帰分析と自作関数
- 2017年4月19日 Rでサイコロの和のシミュレートをしてグラフ化する
- 2016年11月17日 Mac版のRでマルチバイトエラーが出た時は…
- 2016年11月3日 Rでt分布のグラフを描く方法【dt()の使い方】
- 2016年11月2日 Rでベクトルにその要素が無いなら追加する関数
- 2016年7月22日 R言語の便利なコマンド集【第2.5回】
- 2016年7月22日 R言語で線形モデルによる回帰分析【第3回】
- 2016年7月22日 Rでロジスティック回帰分析 そのまま使える自作関数例
- 2016年11月24日 Rで新しくパッケージをインストールするには？
- 2016年7月27日 R言語で線形モデルによる回帰分析の自作関数例
- 2016年7月22日 実行するだけでR言語入門出来るプログラム
- 2016年7月22日 Rでグラフをプロットして保存する【第5回】
POPULAR
- G検定の対策方法まとめ。本、講座、過去問など徹底調査。
- GAFA３社が採用内定の理系日本人へ「年収6000万円」のオファー、「AI人材の争奪…
- G検定に落ちた人、合格した人。勉強法の違いはどこにある？
- G検定に数学力はどれくらい必要？文系でも受かるのか？
- OKI、９月下旬に「AIエッジ・コンテスト」ー「AE2100」で課題解決力競う
- G検定の難易度・合格ライン・内容・評判を徹底調査！
- 【独占】データサイエンスコンペ参加者、３年後に１万人超目指すーNishikaの松田CTOインタビュー
- 【独占】純真学園大学・村中光副学長「G検定、全受講者に取得してもらいたい」ー「人工知能学」新設
- E資格合格談、「コード読解力ついた」ーVSNのデータサイエンティスト・原さん
- 落合陽一氏がアプリ「音声文字変換」を熱く語るーツイート動画再生15万回超
COMMENT
コメントを残す
-
foo.bar 2017.1.18 9:44 AM
細かいことですが，odds rate ではなく odds ratio です
-
Take 2017.1.23 12:56 PM
いつも研究で使用させていただいています。
ひとつ質問なんですがfisher_funに入れる型はどのような形にすればよろしいでしょうか。 -
IMIN 2017.1.23 9:22 PM
いつもありがとうございます。
引数には解析用データを分割表の形にして、入れるようにしてください。以下に使用例を示します。data <- matrix(c(12,7,3,2,15,2),2,2,) #分割表形式の解析用データ fn = “result.csv” #出力ファイル名。 fisher_fun(data,fn) #関数の実行。
-
IMIN 2017.1.23 9:34 PM
ご指摘ありがとうございます。訂正させていただきました。
foo.bar 2017.1.18 9:44 AM
細かいことですが，odds rate ではなく odds ratio です
Take 2017.1.23 12:56 PM
いつも研究で使用させていただいています。
ひとつ質問なんですがfisher_funに入れる型はどのような形にすればよろしいでしょうか。
IMIN 2017.1.23 9:22 PM
いつもありがとうございます。
引数には解析用データを分割表の形にして、入れるようにしてください。以下に使用例を示します。
data <- matrix(c(12,7,3,2,15,2),2,2,) #分割表形式の解析用データ fn = “result.csv” #出力ファイル名。 fisher_fun(data,fn) #関数の実行。
IMIN 2017.1.23 9:34 PM
ご指摘ありがとうございます。訂正させていただきました。