OR的置信区间及如何由置信区间求解标准差.pdf

OR的置信区间及如何由置信区间求解标准差

OR的置信区间及如何由置信区间求解标准差

董圣杰(Dongsj)

02-08-2014

关于OR、lnOR及置信区间的计算

暴露⾮暴露合计处理组

59a）28b）87n1）对照组

8c）

33d）

41n2）

合计67m1）61m2）128N）

1、处理组的暴露⽐值=a/b，对照组的暴露⽐值=c/d，两组的优势⽐OR=abcd

ad

bc

2、OR置信区间的计算：主要有计算⽅法，即Miettinen法和Woolf法，其中Woolf法

是Meta分析中采⽤的⽅法；先简单介绍两种⽅法。

(1)Miettinen法：计算中要利⽤四个表的卡⽅值，计算如下：

(2)Woolf法：该法采⽤了对数转换，利⽤了正态分布的进⾏计算，总体OR的α的置信区间：

59338.69

828

OR??()()()()()

2

226.07

adbcn

abcdacbdχ-?++++()

/2

2

1.961126.07

8.69

3.39,19.93ZOR

αχ

±±

标准正态分布的分位数

[]

/2lnORexplnORZSEα±?ln111111110.456

5928833

OR

SEabcd=+++=+++=()()

expln8.691.960.4563.56,21.24±?

如何采⽤Stata计算OR及其置信区间命令很简单：

Woolf法计算的结果

知道置信区间如何求解标准差（误）

经典统计学置信区间的计算，是通过构建枢轴量，根据枢轴量的分布获得的。

置信区间的求解

1、寻找⼀个样本X1,X2,…,Xn和θ的函数W=W(X1,X2,…,Xn;θ)，使W的分布不依赖于以及其他未知参数，称具有这种性质的函数W为枢轴量；?2、对于给定的置信⽔平1-α，定出两个常数a,b

使得

从求解得到与之等价的θ的不等式，其中和都是统计量，那么就是θ的⼀个置信⽔平为1-α的置信区间。

(){}12,,...;1nPaWXXXbθα

=-()12,,...;naWXXXbθθθθ(),θθ12(),,,nXXXθ12(),,,nXXXθ

如果知道不同统计量的置信区间，求解标准差的过程是相似的，以lnOR为例，设LCI为置信区间的下限，UCI为上限，则

以前⾯的例题为例，lnOR=ln(8.69)，其95%的CI为(1.27，3.06)，则

知道置信区间如何求解标准差（误）

/2ln/2ln/2lnln/2

lnln22OROR

OR

OR

LCIORZSEUCIORZSEUCILCIZSEUCILCISEZαααα=-?=+?-=??-=

ln3.061.271.790.4566

21.963.92

OR

SE-===?原始数据计算⽽得

置信区间计算⽽得

关于对数的⼀点说明

如果a的x次⽅等于N（a0，且a不等于1），那么数x叫做以a为底N的对数（logarithm），记作x=logaN，其中，a叫做对数的底数，N叫做真数，x叫做“以a为底N的对数”。通俗的讲就是求幂。对数可

以⼤⼤的简化运算，特别的在统计学中对数化后，往往可以是数据具有更好的统计性质，更有利于分析和计算。在多数统计软件采⽤的是⾃然对数(Naturallogarithm)，即底数a为⽆理数

e(e=2.71828…)，软件的对应的即log()运算函数。

以BCG数据为例，采⽤R语⾔计算

BCG例中有原始数据，在此旨在演⽰如何计算，并以原始数据为参照进⾏对⽐。数据如下：

数据中的lc及uc为or置信区间的上下限，如果利⽤置信区间及效应值进⾏meta分析时，需要注意应该采⽤对数转换后的数据（⾄少在R、stata软件中是这样）。R、stata进⾏运算时其默认的效应值是

对数转换后的，即lnOR及lnOR的标准误（这⼀点尤为重要）。

RCode

dat=read.csv(“c:/bcg.csv”)

attach(dat)

lnlci=log(lci)

lnuci=log(uci)

se=(lnuci-lnlci)/3.92

lnor=log(or)

library(metafor)

res=rma(yi=lnor,sei=se,method=“FE”)

r