R中匹配配对数据奇数比的单侧置信区间



我正在对两种不同的癌症新治疗的结果进行倾向得分匹配分析,其中结果是二元的(无癌或非癌症(。在成功匹配后,我得到了配对的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 

最新更新