关键词: Meta分析;Bayes方法;MCMC模拟;老年痴呆
摘 要:目的 研究我国老年痴呆患病率及其危险因素. 方法 运用MCMC模拟进行Bayes方法的Meta分析，算法的实现通过WinBUGS软件. 结果 高龄、女性和低文化水平是可能的危险因素. 结论 Meta分析的Bayes模型充分考虑了模型的不确定性，估计更可靠，更合理.
Keywords:Meta analysis;Bayesian method;MCMC simula-tion;dementia
Abstract:AIM To study dementia prevalence in China.METHODS Bayesian Meta analysis was made using Win-BUGS implementation of MCMC simulation.RESULTS Results show that old age，female and low education level are possible risk factors.CONCLUSION Estimation of Bayesian Meta analysis models are more reliable.
0 引言
老年痴呆（Dementia，简称D），其中以Alzheimer病（Alzheimer Dementia，简称AD）和血管性痴呆（Vascular Dementia，简称VD）为主，是严重危害老年人健康的疾病之一.痴呆，特别是AD，其病因仍不太清楚.近10来年间，国内对老年痴呆进行了大量的调查研究.通过对这些资料的综合分析，可以对我国老年痴呆的患病情况有一个全面的了解.
1 资料和方法
1.1 资料 我们收集了全国1986/1996年间10个地区老年痴呆的调查资料［1-10］ .
1.2 方法 Meta分析是将一些独立的研究综合在一起，对一个问题得出一个总的评价的一种定量分析方法.多数研究显示女性可能是痴呆的危险因素，虽然各研究同年龄段的患病率有差别，但都有着随年龄增长患病率增高的趋势，为了去除年龄的影响，我们拟合以性别为分组因素，年龄作为连续性自变量的Logistic回归，对这些资料进行Bayes方法的Meta分析.以rijk ，nijk ，pijk 分别记男性（女性）第j个年龄段的第k个研究的患病人数，调查人数，患病率.模型如下:rijk ～Bin（nijk ，pijk ）logit（p1jk ）=α+β・tjlogit（p2jk ）=α+β・tj +γ
我们关心的量为:ORfm =exp（γ）女性与男性比患病的OR值;p1j =logit-1 （α+β・tj ）不同年龄段男性的患病率;p2j =logit-1 （α+β・tj +γ）不同年龄段女性的患病率.
Bayes分析是根据Bayes公式用后验分布进行推断的，所以后验确定后，问题就转化为后验统计量的计算问题了，如计算后验均值，95%置信区间等.这是Bayes方法的一大难点.这里我们用Markov链Monte Carlo（MCMC）模拟［11］ 实现计算.通过MCMC模拟，可以构造一个Markov过程，它的不动分布就是指定的联合后验，即通过运行一个足够长的模拟，可以使抽样的分布充分接近这个不动分布.那么通过截取这个链的足够长之后的部分就可以得到我们所需要的后验的一个样本，用这样一个足够大的样本我们就可以代替后验总体进行推断了.
2 结果
对于上述模型的计算，我们采用MCMC模拟的方法实现，结果由专门用于MCMC模拟的Win-BUGS软件［12］ 得出.去除性别的影响，对于AD和VD，随着年龄的增长仍明显存在患病率增高的趋势;同时，考虑了年龄因素后，对于AD仍看出女性是危险因素（Tab1）.
表1 不同年龄段男女AD、VD患病率（%）与有关的OR值 略
3 讨论
Meta分析常常因为不相似的研究被不适当的结合在一起而产生出不可信的结论.在经典方法中，模型参数的不确定性常被忽略，或由点估计代替，或者为了计算的简便可行用近似分布代替，而Bayes方法不必局限于计算的方便而选用简单方便的模型.它可适用于任何参数后验分布的计算.特别适合于处理高维问题和多余参数问题.Bayes Meta分析模型充分考虑了模型参数的不确定性，得出的是一种“光滑”估计.在这里我们选用了真实分布 二项分布，而不是正态近似，对于小样本，特别是四格表中出现0时仍有效. 影响患病率的有关因素:①性别差异:不少研究指出，痴呆患病率有性别差异问题，我们综合分析的结果表明:AD以女性患病率为高，男女差异显著.而对于VD，差异不显著.在考虑了年龄的影响后，关系仍未改变.曾经提出的解释有女性受教育程度低，但我们没有收集到不同性别不同文化程度的患病率资料，差异的原因有待进一步研究;②年龄差异:尽管各研究提出不同的痴呆患病率，但均有患病率随年龄增加而上升的趋势.对于AD，增长速度更快，我们得出结论，不论男女，AD的患病率以每5a大于100%的速度递增，而VD以约50%的速度递增，高龄是痴呆的危险因素.
参考文献
［1］Chang JL，Zhang LL.An epidemiological study on senile de-mentia among1390elderly people in Haidian district，Beijing［J］.Zhonghua Liuxingbingxue Zazhi（Chin J Epidemiol），1998;19（1）:18-20.
［2］Chen CY，Chang L.An epidemiological study on visual demen-tia and senile dementia in Shihezi district，Xinjiang ［J］.Zhong-guo Shenjingmiangyixue He Shenjingbingxue Zazhi（Chin J Neuroimmunol Neurol），1998;5（1）:44-47.
［3］Lu SC，Yu HM.An epidemiological study on senile dementia in Dinghai district，Zhoushan ［J］.Zhonghua Jingshenke Zazhi（Chin J Psychiatry），1998;31（4）:225-227.
［4］Kuang XF，Huang XB.An epidemiological study on senile de-mentia in shangrao ［J］.Jiangxi Yiyao Zazhi（Jiangxi Med J），1993;28:213-216.
［5］Li ZJ，Dong ZF.An epidemiological study on senile dementia prevalence in Jianguomen Block，East district，Beijing ［J］.Zhonghua Laonianxue Zazhi（Chin J Gerontol），1999;19:3-4.
［6］Tang XT，Tian DZ.An epidemiological study on dementia of the aged Hui and Han residents in Niujie area of Beijing ［J］.Zhonghua Laonianyixue Zazhi（Chin J Geriatr），1998;17（5）:307-310.
［7］Tang MN，Guo YB.Epidemiology of senile dementia and Alzheimer disease in the rural area ［J］.Linchuang Jingshenyix-ue Zazhi（J Clin Psychol Med），1999;9（1）:20-22.
［8］Xue GH，Shao YC.An epidemiological study on senile dementia in Guangdong ［J］.Shiyong Yixue Zazhi（J Pract Med），1997;13:371-372.
［9］Yu JL，Feng RM.A prevalence study of senile dementia in lit-toral rural of Guangdong ［J］.Zhongguo Shenjingjingshenjibing Zazhi（Chin J Nerv Ment Dis），1998;24:30-32.
［10］ Zhang HR，Wen GQ.An epidemiological study on senile prophase and senile dementia in Hainan ［J］.Hainan Yixue（Hainan Med J），1996;3:207-208.
［11］Gelman A，Rubin DB.Markov chain Monte Carlo methods in biostatistics ［J］.Stat Methods Med Res，1996:5，339-355.［12］Spiegelhalter D，Thomas A，Best N.WinBUGS Version1.3［M］.Cambridge:MRC Biostatistics Unit，Institute of Public Health，1999.