目的

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

使用法

relative.risk(a, b, c, d)

引数


	
 
 
所見（結果）あり
 
所見（結果）なし
	
 
対象群（要因あり）
 
a
 
b
	
 
対照群（要因なし）
 
c
 
d


ソース

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

# 相対危険度（対応のない場合）とその信頼限界を求める
relative.risk <- function(   a,                      # 対象群・所見あり
                                b,                      # 対象群・所見なし
                                c,                      # 対照群・所見あり
                                d)                      # 対照群・所見なし
{
        cl <- function(x)
        {
                exp(log(rr)+c(1, -1)*qnorm(x)*sqrt(b/a/(a+b)+d/c/(c+d)))
        }
        rr <- a*(c+d)/c/(a+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.risk(76, 399, 129, 332)

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

$conf	# 信頼区間
              下側信頼限界値 上側信頼限界値
90%信頼区間        0.4624827      0.7069145
95%信頼区間        0.4440632      0.7362370
99%信頼区間        0.4101544      0.7971041
99.9%信頼区間      0.3740365      0.8740745