個別メモ

Latest update on 2012年3月5日 (月) at 10:54:46.

【第2125回】メール返信とか（2012年2月26日）

日曜は仕事メールへの返信が多かった。
久々に青木先生の掲示板を覗いたら，正確なウィルコクソンの順位和検定がcoinパッケージでできなかったと書かれていたので，少し調べてみた。拙著『Rによる保健医療データ解析演習』では，exactRankTestsパッケージのwilcox.exact関数を使えばできると記述していて，それは今でも間違いないのだが，exactRankTestsパッケージをロードすると，既に開発終了しているので，代わりにcoinパッケージをお薦めするというような意味のメッセージが表示されるので，できればcoinパッケージでやりたい。そこで調べてみると，coinパッケージのwilcox_test関数では，オプションとしてdistribution="exact"を付ければ正確な確率を計算してくれることがわかったので，青木先生の掲示板に書き込んでみた。
なお，これら２つのパッケージには微妙な違いがあって，exactRankTestパッケージのwilcox.exact関数ではチルダの右側にint型の変数もおけるけれども，coinパッケージのwilcox_test関数では，チルダの右側には2水準のfactorでなくてはエラーが出てしまう。as.factor()すればいいのだが，少し面倒だ。

▼前【2124】（雑務多し（2012年2月24～25日））　▲次【2126】（週明け（2012年2月27日））　●Top

△Read/Write COMMENTS