道草 Rで 全因子26データからMHRDを作る記述
保存したデータ
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")
【生起因子で26因子を調整して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 )
}
# * }