何度もやってるのにどうも覚えられないのでここにまとめてメモ書き。自分の研究は、質問紙などの順序尺度データだったり、正規性のないデータだったりを使うことが多いので、相関分析をする場合には必然的にピアソンの積率相関係数ではなく、順位相関係数を使うことが多い。基本的にはスピアマンを使っているが、ケンドールが報告されているのをあまり見かけないのでいつかドヤ顔で使ってみたいという厨二心もある。順位相関係数を出すだけならすぐにできる。Rだとこう(ダジャレではない)。#ケンドール無相関検定とか信頼区間の算出をする場合はcor.test関数を使う(引数が2つ必要なことに注意)。値を計算することができはいでてきt、あ、あれ?調べてみるとcor.test関数では、ピアソンの積率相関の場合のみにしか信頼区間を出してくれないらしい。http://www.inside-r.org/r-doc/stats/cor.test次はこれを試してみた。corrgramパッケージのcorrgram関数。信頼区間が一気出しできるpanel.confオプションが便利。https://cran.r-project.org/web/packages/corrgram/corrgram.pdfだがここにも落とし穴。母相関係数の信頼区間を一気出しするためにcorrgramを使ったら、出てきた相関係数(スピアマンの順位相関)が、psychパッケージのpairs.panels(method=”spearman”)で出した相関係数と一致せず、おかしいなあと思って調べてみた。2015年7月31日マニュアルを読みなおしてみると、どうもcorrgramで信頼区間を出すpanel.confオプションは、オプションで順位相関を選んだ場合(cor.method=”spearmman/kendall”)機能しないそうだ。結局は中でcor.testが動いているので、順位相関では信頼区間が出せないいうことだそうだ。【解決策1】langtestを使う。スピアマンでもケンドールでも使える。簡単かつ手っ取り早い。http://langtest.jp/shiny/cor/【解決策2】スピアマンの場合、DescToolsパッケージのSpearmanRho関数を使う。0.8506335この結果はlangtestの結果と一致する。【解決策3】これもスピアマンの場合。RVAideMemoireパッケージのspearman.ci関数を使う。https://cran.r-project.org/web/packages/RVAideMemoire/RVAideMemoire.pdfこれは上の2つとは違い、スピアマン順位相関のブートストラップ信頼区間を出してくれる。estimates:rho0.6544704【解決策4】これはケンドールの場合のみ。NSM3パッケージのkendall.ci関数を使う。https://cran.r-project.org/web/packages/NSM3/NSM3.pdfどうもこの本の計算方法に基づいて信頼区間を算出するようだ。0.838ブートストラップを使用する場合は、Bの部分で抽出回数を指定する。種類があって、どれを使うのがいいかはよくわからないし、もう全部langtestでOKじゃなかろうかと思ったりもするが、手数は多いに越したことはないだろう。【追記1】2016.03.13このようなご指摘があったのでここで紹介させていただきます。より,スピアマンの順位相関係数と信頼区間が求まる(表示はピアソンの積率相関係数だが)。#15にて、このブログの内容について発表しました。以下は発表スライドです。【追記3】2016.03.27スピアマンの順位相関係数の信頼区間について、以下のページをご紹介いただきました。RjpWikicor.test関数とrank関数の組み合わせで算出した信頼区間と、ブートストラップで算出した信頼区間との間には大きな差はないようだ。川口勇作この夏に開催される関連学会にて筆頭発表者としていくつか発表をさせていただくことになっているので,ここにまとめていきたいと思います。あと筆頭発表者ではないものもまとめさせていただいてます。全国研究大会にて発表させていただきます(私の発表は2日目です)。試み:英語を学習する大学生を対象として」です。草薙邦広さん(名古屋大学大学院)との共同研究です。いう尺度を作成して,その妥当性,信頼性を検証するという内容です。発表スライドはこちらです。同日に(というより発表の直後に)石井雄隆さん(早稲田大学大学院),石井卓巳さん(筑波大学大学院),草薙さん,阿部大輔さん,福田純也さん(いずれも名古屋大学大学院)との公募シンポジウムにも登壇します。Project―」というタイトルです。ライティングプロセス記録・分析ソフト「WritingMaetriX」を用いて記録したデータで大規模なライティングプロセスコーパスを作成しようという試みについての内容です。(JASELE)第40回徳島研究大会にて発表させていただきます(発表はこれまた2日目です。こちらは「エッセイライティングにおける増加語数の時系列推移傾向はエッセイ評価を予測するか―線形回帰モデルおよびポアソン分布へのフィッティングを用いて―」というやたら長いタイトルのものです。こちらは室田大介さん,後藤亜希さん(いずれも名古屋大学大学院)との共同研究です。記録したライティングプロセスを線形回帰モデル(最初から最後まで一貫して語数が伸び続けるライティングプロセス)とポアソン分布(望ましいとされるライティングプロセスに形が似ている)にあてはめてみて,それぞれへの当てはまりのよさがエッセイの評定と相関するかを見た研究です。発表スライドはこちらです。おける文法的慎重性尺度の開発」もあります。こちらは,「文法的慎重性」という学習者固有の要因を測定する質問紙の作成を行い,その妥当性や信頼性を検証したものです。関心のある方はぜひお越しください。部会長の福田純也先生がご自身のブログでご紹介くださっていたので,リブログさせていただきました。みなさまの温かいご支援をよろしくお願いいたします。川口勇作お久しぶりです。この間LET中部での発表が終わりましたのでそのご報告です。発表スライドはこちらです。配布資料はこちらで公開しています。https://sites.google.com/site/kwsk3939/home/closetなぜかSlideShareをWordpressに貼り付けると大きいサイズになってしまうんだけどなんでだろう。質疑ではたくさんのご質問・ご助言をいただいてありがたかったです。発表にお越しいただいた皆様にこの場をお借りしてお礼申し上げます。先輩にもたくさん助けていただきました。早く独り立ちできるようにがんばります。発表する度に自分の勉強不足を痛感して嫌になってきますが,これからもめげずにやっていきたいと思います。今年の発表も残すところあと一つ。12月7日のNagoya.Rでの入門者講習です。http://atnd.org/events/43634R入門者が講師の入門者講習という斬新な企画・・ぜひお越しください。寒くなってきました,もう完全に秋ですね。いうことで学会・研究会の情報を。秋田(秋田県)メソ研には初めての参加になるのでワクワクしています。東北に行くのも初めてなので,そちらも楽しみです。