目的

相対危険度（対応のある場合）とその信頼限界を求める

使用法

relative.risk2(b, c)

引数

以下のような分割表における b と c


	
 
 
症例群（要因あり）
 
症例群（要因なし）
	
 
対照群（要因あり）
 
a
 
b
	
 
対照群（要因なし）
 
c
 
d


ソース

インストールは，以下の 1 行をコピーし，R コンソールにペーストする
source("http://aoki2.si.gunma-u.ac.jp/R/src/relative_risk2.R", encoding="euc-jp")

# 相対危険度（対応のある場合）とその信頼限界を求める
relative.risk2 <- function(  b,                      # 対照群要因あり，症例群要因なし
                                c)                      # 対照群要因なし，症例群要因あり
{
        cl <- function(x)
        {
                exp(log(rr)+c(1, -1)*qnorm(x)*sqrt(1/b+1/c))
        }
        rr <- c/b                                    # 相対危険度
        conf <- rbind(cl90=cl(0.05), cl95=cl(0.025), cl99=cl(0.005), cl999=cl(0.0005))
        colnames(conf) <- paste(c("下側","上側"), "信頼限界値", sep="")
        rownames(conf) <- paste(c(90, 95, 99, 99.9), "%信頼区間", sep="")
        list(rr=rr, conf=conf)
}


使用例

> relative.risk2(4, 25)

$rr	# 相対危険度
[1] 6.25

$conf	# 信頼区間
              下側信頼限界値 上側信頼限界値
90%信頼区間         2.577450       15.15548
95%信頼区間         2.175174       17.95834
99%信頼区間         1.561200       25.02081
99.9%信頼区間       1.062451       36.76640