相対危険度(対応のない場合)     Last modified: Apr 13, 2004

目的

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

使用法

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


・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI

Made with Macintosh