0

相対リスクの標準誤差をコーディングしようとしています。多くの場合、信頼限界のためにSEとNaNの無限大を取得しています。有限の標準誤差を得るには、パラメーターの制限を指定する必要があるかどうか疑問に思いました。これが私のコードです

    set.seed(111) 
    k<-63
    n<-60
    xt<-NULL
    xc<-NULL

    pt<-rbeta(k,3,3)   # treatment effect
    pc<-rbeta(k,1,3)   # control effect
    true_RR<-pt/pc   # True Relative Risk

    for(i in 1:k)
    {
    xt<-cbind(xt,rbinom(300,n,pt[i]))

    xc<-cbind(xc,rbinom(300,n,pc[i]))
     }
    a<-xt
    b<-n-a
    c<-xc
    d<-n-c

    RR<-xt/xc

    for (i in 1:300){

    for (j in 1:63) 
    {
    if (xc[i,j]==0)
    {
     RR[i,j]<-0
     }}}
     se_rate<-sqrt((b/(a*(a+b)))+(d/(c*(c+d))))

     l_cl<-RR*exp(-1.96*se_rate)
     u_cl<-RR*exp(1.96*se_rate)
4

1 に答える 1

2

あなたが得てInfいる理由は、0で割っているためです。特に、設定します

xc = cbind(xc,rbinom(300,n,pc[i]))
c = xc

そして後で計算します：

d/(c*(c+d))

cしたがって、ゼロに等しく設定される可能性はゼロではありません。代わりにすべきことは統計の問題であり、適切なアドバイスを提供するのに十分な情報がありません.

于 2012-08-10T17:47:22.597 に答える