SASとRでの使い方メモ。

	R DescTools::BinomDiffCI method = TYPE	SAS PROC FREQ RISKDIFF (CL = TYPE)
Agresti-Caffo	"ac"	AC
Brown, Li's Jeffreys	"blj"	-
Farrington and Manning	-	FM
Haldane	"hal"	-
Hauck-Anderson	"ha"	HA
Jeffreys-Perks	"jp"	-
Mee	"mee"	MN (CORRECT= NO)
Miettinen and Nurminen	"mn"	MN
Newcombe (Corrected)	"scorecc"	NEWCOMBE (CORRECT) \| SCORE (CORRECT) \| WILSON (CORRECT)
Newcombe	"score"	NEWCOMBE \| SCORE \| WILSON
Wald (Corrected)	"waldcc"	WALD (CORRECT)
Wald	"wald"	WALD

Rでの実行

c00 <- 50; c01 <- 75; N0 <- c00+c01;
c10 <- 40; c11 <- 90; N1 <- c10+c11;
 
library(DescTools)
100 * BinomDiffCI(c00, N0, c10, N1, conf.level=.95, method="ac")
100 * BinomDiffCI(c00, N0, c10, N1, conf.level=.95, method="ha")
100 * BinomDiffCI(c00, N0, c10, N1, conf.level=.95, method="mee")
100 * BinomDiffCI(c00, N0, c10, N1, conf.level=.95, method="mn")
100 * BinomDiffCI(c00, N0, c10, N1, conf.level=.95, method="score")
100 * BinomDiffCI(c00, N0, c10, N1, conf.level=.95, method="scorecc")
100 * BinomDiffCI(c00, N0, c10, N1, conf.level=.95, method="wald")

SASでの実行

data ADS;
    length FACTOR YN CNT 8.;
    input FACTOR YN CNT @@;
    cards;
0 0 50  0 1 75  1 0 40  1 1 90
;
run;
  
proc freq data=ADS;
    table FACTOR * YN / riskdiff(cl=(wald ha fm ac mn newcombe));
    table FACTOR * YN / riskdiff(cl=(wald(correct) mn(correct=no) newcombe(correct)));
    weight CNT / zeros;
run;