基本的に、お仕事に関するネタは書かないようにしているのですが。
ふとジェネリックレビューというジェネリック医薬品の評価サイトを見つけ、どういう手法で評価しているのかなとサイト運営者のblogを見てみたのだけど肝心要のRによるコードが削除されてるというのがなんとも。とりあえずCohen's dを求めれば良いらしいということでWikipediaやEffect Size Calculator for T-Testや効果量と検定力分析入門(PDF)を参考にコード書いてみました。C屋なので、当然Cで記述します。
gcc -O2 cohen.c; ./a.out < etizolam-0.5.txt > result1.txtという感じでやっておりますが、リファレンス(クラビットのCmax)を使って計算するとRで得られたものと比べて傾向は似ているものの計算結果がなんとなく違うんですよね。どうも対数変換をサボっているためにこういう違いが出ている感じです。
とりあえず、ジェネリック医薬品のメーカー変更を行う際にどの程度のリスクがあるかというのを知る手がかりって何か無いのかなーとは前々から思っていて、こういうやり方があるならちょっと追ってみようかなと考えているところです。59.20kg(23:15)