morの解析ブログ

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

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

Rで計算 2×2表でみたい... 配置するだけの作り方

2×2表を並べ替えて作る.
 1行目はab、2行目はcdのように並べる.
 intetab ;因子ごとの発生数dfから作る.
------------------

ad<-NULL
ab<-NULL
 for (i in 1:8)
{
ad<-cbind( intetab[i,1],intetab[i,2])
ab<-cbind(ab,ad)
}
ab
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15]
[1,] 139 0 11 128 95 44 131 8 104 35 98 41 132 7 123
[,16]
[1,] 16
--------------
cdは、
cdad<-NULL
> cd<-NULL
>
> for (i in 1:8)
+ {
+ cdad<-cbind( intetab[i,3],intetab[i,4])
+ cd<-cbind(cd,cdad)
+ }
> cd
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15]
[1,] 0 119 15 104 84 35 105 14 71 48 76 43 81 38 88
[,16]
[1,] 31
-----------------
2行を合わせて・・
 

           Rstudioで見る限りは、ズレない.
 カラム名をつけておく.
colnames(abcd) <-c("y","","wat","","tya","","mesi","","tori","","sake","","tam","","potesara","")

×

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