相関係数は比べてもいいのか。

相関係数の比較の話。

相関係数について、今まで、比べてはならないものだと思い込んでいた。

たぶん学部の時にそう習ったのか、テスト前にネットサーフィンしているのか忘れたけど。

たしか、文系科目として英語と国語、理系科目として物理と化学。その２つの相関係数を比べても意味ない。的な話で妙に納得していたのを覚えている。

ただ、今回医師の先生から質問されて、相関係数は比べてはならない。って即答してしまったけど、なぜか違和感があった。

そこでいろいろ調べてみた。


結論から書くと、比べてもいい時もある。

例えば同じような指標を作ろうとしている時。

上の例で言うと英語と国語を比べている時と英語と古典を比べている時、文系科目としての相関が高いのはどっちかってのは比べる意味がある。

医学でいうと、OSとその他のサロゲートとの比較は意味があると思う。

ただ、相関係数にも弱点があって、それぞれの変数が正規分布してないと、例えば二峰性であったりすると、おかしな結果を出してしまうということ。

つまり、個人的には相関係数を比較するときはそれぞれの変数のヒストグラム、代表値を確認して正規分布だとわかってから使うべき。との解釈。

つまり予備検定としてはKruskalWallisして確認した後、相関の差の検定をするべきか。

医学の分野だと正の値しか取らないから基本歪んでるよね。

つまり、統計を正しく使える環境にない人はあまり相関係数の差に対しての言及はやめたほうがいいってことになるかな。

難しいところだけども。

追記

実際の検定はこちらから。毎度のことA先生のページ。

http://aoki2.si.gunma-u.ac.jp/R/diff_r.html

強制投入した場合の推定値の信頼区間

強制投入で解析したら推定値の信頼区間でないじゃないかという要望から。

confint（モデル式）

で出てきます。

青木先生の関数fisher exact testを使う時の注意点

http://aoki2.si.gunma-u.ac.jp/R/fisher.html
このページから。

source("http://aoki2.si.gunma-u.ac.jp/R/src/fisher.R", encoding="euc-jp")

をコピペしたらfisher exact testが出来る。

ｘに適当なデータを入れてやってみると

> fisher(x)

 以下にエラー asymptotic() :  関数 "printf" を見つけることができませんでした 

と出てくる事がある。

この場合はこのページ

http://aoki2.si.gunma-u.ac.jp/taygeta/statistics.cgi?mode=res&no=15678

の真ん中あたりで回答があり、

printf <- function(fmt, ...)
{
cat(sprintf(fmt, ...))
}

と打ち込めば使えるようになります。

何度か間違ってしまったのでまとめ。

all.R入れろって話ですかね。うーん。

あと、レトロな解析な場合は基本ｔ検定とカイ自乗検定は使わずにマンホイットニー（ウィルコクソン）とフィッシャーでいったほうがいいと思います。

英語

英語に対しての能力が伸び悩んでいる。

多分言葉にすると「スライドを使いながら話しているプレゼンテーションならいけるけど、ただ話しているだけの授業は内容を理解できない」レベル。

サブスクリプトを見たら内容は大体わかるけど、スクリプトなしでは聞き取れないレベルって感じか。

どうしたらええんやろなー。

もう一回基礎に戻って単語とかやり直すか・・・？

いい案ある人は教えていただければ幸いです。

エクセルでシートの変更他ショートカット

エクセルのショートカットで今日覚えた

ｃｔｒｌ＋PageUp（Down)でシートの変更

F12で名前をつけて保存

エクセルファイルでアンケートをまとめてもらって、一つ一つをCSVファイルにするときにすべてキーボードで出来るようになると意外に楽だった。

30シートほどしかなかったので、マクロ組んでVBAいじる程でもなかった。

強制投入法以外で変数候補を探索する。Rで。

まずRの画面で

library(MASS)

と入力し実行。（ダウンロードは必要ないと思います。）

この宣言でステップワイズ法が使えるようになります。

デフォルトでstepという関数があるけど個人的にちょっと結果がみにくいと思っている。

計算自体は

stepAIC(coxph(Surv(time,outcome)~factor1+factor2+...,data=解析したいデータセットの名前),direction="backward")

で行えます。（環境によっては改行が入るかもしれませんが1行です）

coxph()の部分は考えるモデルで強制投入の時に作ったものと同じ。もちろんコックスハザード以外も使えますが、医学論文のための手法としてはこれかロジスティックでしょう。

あと

cox.analysis<-coxph(Surv(time,outcome)~factor1+factor2+...,data=data)

などとしていただいて、

stepAIC(cox.analysis,direction="backward")

としていただいても構いません。

後ろのdirectionの部分はbackward,forward,bothの3種類を指定できます。

減らしていく方法、増やしていく方法、どちらもやっていく方法です。

医学論文では基本的にいわゆるp値でしかみないけど、AICって指標を用いることでどの要因が効いているのかを評価できる方法もあります。

AICってのは日本人の赤池先生が作って世界中で使われている指標です。

ジムに行こうか検討している。 ポーンを３ｆへ。

題名で一人チェスやってみるっていうのどう？




学部生の頃からの趣味というかやめられないことにサウナってのがある。

それはたしか減量する為に行き始めたんだけど、やっぱ汗かくって気持ちいいやんか。

出た後、何も考えずに寝て、目が覚めたら、え？１時間くらいしか寝てないの？もう１日分寝たきがすんねんけど！みたいな感じになるのが好き。

で、東京に来てからその感覚を忘れていたというか毎日いけるところに銭湯すらなくて、ちょっとチャリこいだらあったんだけど、冬は冷えてしまう距離、夏はまた汗をかいてしまう距離みたいなところにあったのでイケてなかった。

ただ、今は車もあるし、それなりに安い温泉が近くに行って、正直毎日5,600円それにかけるくらいは稼いでるし、またサウナがルーチンワークになりつつある。

ただ、毎日毎日行っていると、まあ１万５０００円くらいそれに使ってることになるわけで、それなら、なんかジムとか行ったらええやん。最悪サウナだけになってもいいし、それで運動はじめるなら、それはそれでいいやん。みたいな感じに思ってきた。

で、近くにジムがあるのかなって見たら、家と勤務先の間にあるみたいで、まあ月額も１万円弱とかで、温泉なら２日に１回行くなら元が取れるレベル。

今のコトロ出張とかなかったら基本的に行ってしまうところがあるから、まあいっかなーって感じだけど、なぜか踏ん切りがつかない。

それは、多分、一回岡山にいた時に入ったけど、結局忙しくて全然イカなかったってのがあるんだと思う。

一応朝は空いてないけど、夜は２３時まで開いているし、大体デューティーのタスクは毎日７時に一度終わるような状態なので、７時に中抜けして９時半くらいに再出勤で１２時くらいまでいるというような形も可能と言えば可能。（というか最近いつもそんな感じ）

それを考えると、やらない理由がないのだけど。

とりあえず週末の出張が終われば、大きな出張は２ヶ月くらいないし、まあ行ってみるかーって感じで。

多分環境と、そこで友達（？）ができるとやるだろうし、取得しているコーチの資格の講座で学んだ、メニューの建て方とかを実践して、セルフコーチングみたいなこともやれたらなとか思っている。

そんな時間あるの？って感じだけど・・・まあ入会金もそんなに高くないし、無理ならやめたら一家ーって感じで。

正直６歳からずっと部活のある生活をしてきたので、運動ないってのがなんか考えられない。
とりあえず入ってやってみるかーって感じで。

というかフットサルとかのチームがあればいいんだけど、ここ静岡やん？

なんかサッカー大国やん？

よっていけなさそうやん？

っていう被害妄想をどうにかしたいけどね。うん。

っていう結局何も決まらない日記。