Rによる条件付きロジスティック回帰分析と自作関数2016/11/212020/05/04Rによる条件付きロジスティック回帰分析と自作関数R言語入門ライター:高橋光太郎TwitterFacebook条件付きロジスティックモデルは、ロジスティックモデルを、対応のある2値データへ拡張したものです。R言語のsurvivalというパッケージには、条件付きロジスティック回帰分析を行うための関数、clogit()が入っています。ここでは、この関数clogit()の使い方を実例を交えて紹介していきます。解析結果をcsvファイルへの出力まで行ってくれる自作関数も作りましたので、そちらも是非お使いください。dat)として、使います。yには目的変数の列、xには説明変数の列を入れてください。説明変数の数だけ引数の中身も長くなっていきます。datには、解析用のデータが入ったデータフレームを入れてあげましょう。この関数を使うためには、survivalパッケージの読み込みを必ず行いましょう。見つけることができませんでしたというエラーに見舞われます。条件付きロジスティック回帰分析の自作関数プログラムlibrary(survival)#ライブラリの読み込み#datは目的変数列,説明変数列,条件列の順。#回帰分析の詳細な結果を得る。#出力結果がわかりやすいように列(行)に名前をつける。rbind(c(“目的変数”,“説明変数”,“条件”),names)write.table(names,file,append=T,quote=F,sep=“,”,row.names=F,col.names=F,fileEncoding=“CP932”)write.table(matrix(c(“”,colnames(result)),nrow=1),file,append=T,quote=F,sep=“,”,row.names=F,col.names=F,fileEncoding=“CP932”)write.table(result,file,append=T,quote=F,sep=“,”,row.names=T,col.names=F,fileEncoding=“CP932”)}以上の関数に、解析用のデータを体裁を整えて、引数として入れれば、条件付きロジスティック回帰分析が行われます。この関数では、説明変数は1つにしてあります。多変量で条件付きでロジスティックを重回帰させたい場合は、少々書き直して使っていただければと思います。以下に関数の使用例を載せます。match_logi_fun(dat,”条件付きロジスティック回帰分析結果.csv”)その結果、このように出力されます。条件付きロジスティック回帰分析では、条件付きオッズ比が推定されます。回)TwitterFacebookライター:高橋光太郎R言語入門POPULAR最弱オセロを初めて攻略した天才オセロ高校生。負け方を解説!G検定に落ちた人、合格した人。勉強法の違いはどこにある?G検定は文系でも合格できる。TrendCOMMENTコメントを残すコメントをキャンセルメールアドレスが公開されることはありません。コメント*お名前*メールアドレス*プライバシーポリシーはこちら上に表示された文字を入力してください。