先ほどのデータを解析 source("http://aoki2.si.gunma-u.ac.jp/R/src/ROC.R", encoding="euc-jp") ROC0(acs,survive) AUCが最後に表示される ... $c.index [1] 0.6589782 ここによるとあんまり尤もらしい曲線じゃないけどね！ Epiパッケージを使うともっとす…

事象がおきたかどうかについて、ある値が高いか低いかでわけるとき R ではどうするか データを用意 読み込み d2 <- read.table(pipe("pbpaste"), header=TRUE) subsetを使って切り分け、そのグループをt検定（あるいはt.test(LHratio~ACS,d2) だけでもいいっ…

Prismで普段解析していますがGUIなだけに自動化しにくい．．． 特にグループ分けした場合それぞれを別のカラムに入力しなければいけないのがネック。 例えばgroup0-3まであった場合↓のようにexcelのマクロを使って整形しておかなければならない そこでRの出…