morの解析ブログ

解析疫学、リスクにまつわるメモや計算

「推定」のまわりをさぐる.教科書では「解析はMHにより行う、因子が多ければ重回帰を用いる」という風で詳しい例は少ない.独自(のつもり)な思いつきで具体に試行.
 数理を用いるべきアセスメントにも切り込む.

Rで計算 cORを並べてみたい

・xtableを作ったのに加えて、cORをそろえて表示する
cor<-NULL
 for( i in 1:8)
  { 
  cor<-cbind(cor ,round(intetab[i,1]*intetab[i,4]/intetab[i,2]/intetab[i,3],digits=2) ,"")
   }     # ORの値と、空白が 交互に並ぶのがミソ
  # ORは小数点2までとした


・tableとくっつける
  tab4_cor<-rbind(abcd,cor) # 2×2表の下に cORが表示される
 

 いつのまにか、文字になっている.Inf の文字のせいか、数値に戻しがたい.
 ビュアーで指標を並べてみるならこれでよしとする.

×

非ログインユーザーとして返信する