保存したデータ
　csv　■◇　ただしい症例定義・・　.csv　；　　 全メニュー
　cludeなRDデータ名：crd
------------　
　data = kanzi　  記述　　
dd<- NULL  
dd<- kanzi
     #  dd= read.csv(file.choose())　　#　dd：csv　新データ
　　 #　空白削除
　　d26 <-complete.cases(dd)
　　d <- d[d26,]　　　　　　　　　#　d　：Rに保存された全因子データ

【cludeなRD】曝露数、部分リスクをも計算

   crd<- NULL
　t_ec<- NULL
　hidari<-NULL
     migi<-NULL
　 for　(x in 1:27) {　　　　　　　 　 #　x　　　全因子を調整対象とし、、、、　　　　　　　
　　a1 <- sum( d[,1] *d[x]　)　　　　#　発生有
　　b1 <- sum( d[,1] *(1-d[x]) )　　　#　
   　c1 <- sum( (1-d[,1]) *d[x] )　　 　 # 発生なし
　　d1 <- sum( (1-d[,1]) *(1-d[x]) )
                                e_count <- a1+c1　　 # 曝露数
 　n1 <-a1+b1+c1+d1　
 　  k11<-a1+c1 
      k10<- b1+d1
　　　　　rd1 <- a1/k11-b1/k10
                           　hidari<- c(hidari,a1/k11) 　　# リスク差　左
                              migi <- c(migi,b0/k10)
crd<-c(crd,rd1)　
                           t_ec<-c(t_ec,e_count)    # 曝露数ベクトル
}
    ----------
 書き出し　
　例　write.csv( データフレーム名, "ファイル名.csv")
　例　write.csv(data4, "data4.csv")
　　write.csv(hidari　, "hidari.csv")
　　write.csv(migi　, "migi.csv")
　　　　
　【生起因子で２６因子を調整してMHRD】
　　　　# 26因子データから　tで調整したMHOR を計算
　　　y <-13　　#　固定する　*　yについて全因子を指定すると全因子×全因子のMHRD
mhrd　<-NULL
tmhrd<-NULL
for　(x in 1:27) {　　　　　　 #　x　　　全因子を調整対象とし、、、、　　　　　　　
# for　(y in 1:27) {　　　　　 #　y　13 ; t 因子固定*　してt 有無についてtableを　　
　a1 <- sum( d[,1] *d[x]*d[y]　)　　　#　y　t あり
　b1 <- sum( d[,1] *(1-d[x]) *d[y])　　#　xの2×2を 固定yの tで層化
　c1 <- sum( (1-d[,1]) *d[x] *d[y]) 　
　d1 <- sum( (1-d[,1]) *(1-d[x])*d[y] )
　　a0<- sum( d[,1] *d[x]*(1-d[y])　)　　#　y t なし　
　　b0 <- sum( d[,1] *(1-d[x]) *(1-d[y] ) )
　　c0 <- sum( (1-d[,1]) *d[x] *(1-d[y]) ) 　
　　d0 <- sum( (1-d[,1]) *(1-d[x])*(1-d[y]) )
#　式　；　　　　　　　ありなしによる、　の mhrd；；
　　n1 <-a1+b1+c1+d1
　　 n0 <-a0+b0+c0+d0　
　 　k11<-a1+c1
k10<- b1+d1
k01<-a0+c0
k00<- b0+d0
rd1 <- a1/k11-b1/k10
rd0 <-a0/k01 -b0/k00
wei0 <- k01*k00/(k01+k00)
wei1 <- k11*k10/(k11+k10)
　mhrd<- (wei0*rd1+wei1*rd0)/(wei0+wei1)
mhrd
tmhrd<-c(tmhrd,mhrd)
　# print( tmhrd )　
}
# * }