我正在对两种不同的癌症新治疗的结果进行倾向得分匹配分析,其中结果是二元的(无癌或非癌症(。在成功匹配后,我得到了配对的2x2列联表,用于我匹配的配对之间的结果,如下所示;
**Treatment 1**
Not-Cancer Free Cancer Free
**Treatment 2**. Not-Cancer Free 42 39
Cancer Free 53 50
为了获得我的比值比,我使用McNemars Exact Test,使用R包exact2x2(函数称为mcnemar.cexact(,并获得以下结果
Exact McNemar test (with central confidence intervals)
data: MatchedCasesTable
b = 39, c = 53, p-value = 0.175
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
0.4738071 1.1339142
sample estimates:
odds ratio
0.7358491
我知道p值和95%置信区间来自双侧测试方法,但我想对匹配配对的比值比进行单侧测试,因为我只对单侧测试的较低置信区间感兴趣。有没有任何方法可以修改测试或任何其他方法,这样我就可以得到一个单侧测试,也可以提供准确的较低比值比置信区间?提前感谢您!
您可能想要的是:
tt = data.table(tx1 = c(1,0,1,0), tx2 = c(0,0,1,1), n = c(53,42,50,39))[rep(1:.N,n)]
exact2x2(tt[,table(tx1,tx2)], paired=T, alternative = "greater")
输出
Exact McNemar-type test
data: tt[, table(tx1, tx2)]
b = 39, c = 53, p-value = 0.9413
alternative hypothesis: true odds ratio is greater than 1
95 percent confidence interval:
0.507341 Inf
sample estimates:
odds ratio
0.7358491