目的:利用时空分布模型和集成嵌套拉普拉斯近似算法研究上海市肺结核的时空特征及其影响因素，为区域化肺结核疫情防控措施的制定提供理论依据。方法:基于中国疾病预防控制信息系统结核病管理信息系统下载的2013-2020年上海市户籍肺结核病例，采用层次贝叶斯模型对肺结核病例数据进行拟合，识别肺结核的时空变化特征，从时空角度，评估肺结核发病风险与人口、经济、卫生服务等因素的关系。结果:2013-2020年，上海市共报告肺结核病例29 281例，年均报告发病率为25.224/10万。2013-2020年的发病率呈先上升后下降的趋势，2014年的报告发病率最高（27.991/10万）。上海市肺结核的发生具有空间聚集性，通过肺结核报告发病率的热点分析和风险分析发现，上海市肺结核的高风险区域为郊区，低风险区域为中心城区。肺结核的发病风险与地区人均生产总值（ RR=0.48）、每万人医疗机构床位数（ RR=0.56）、归一化植被指数（ RR=0.50）、夜间灯光指数（ RR=0.80）有关。 结论:在上海市中心城区肺结核防控工作的稳步推进下，需要格外关注郊区的防控。在后续的研究中，应进一步探索肺结核发病与归一化植被指数和夜间灯光指数之间的联系。

目的:研究大气NO 2短期暴露与心力衰竭住院风险的关联。 方法:基于China-PEACE心力衰竭研究，将2015年1月1日至12月31日来自全国62个城市92家医院的117 364例心力衰竭住院患者纳入研究，并获取相应时间段的日均NO 2暴露水平、温度和相对湿度等信息。采用广义可加模型和贝叶斯层次模型分析NO 2短期暴露（滞后 0~3 d）与心力衰竭住院风险的关联，并比较不同亚组间NO 2短期暴露与心力衰竭住院风险的关联。 结果:研究对象年龄为（70.32±12.22）岁，2015年1月1日至12月31日，62个城市NO 2全年365日日均暴露浓度的中位数为26.4 μg/m 3，最小值为2.33 μg/m 3，最大值为150.25 μg/m 3。随着NO 2暴露浓度每升高10 μg/m 3，心力衰竭住院风险增加2.2%（ OR=1.022，95% CI：1.012~1.031）。在住院日前0~1、0~2和0~3 d，NO 2移动平均暴露水平与心力衰竭住院风险的关联均具有统计学意义［ OR值分别为1.020（95% CI：1.009~1.030）、1.016（95% CI：1.004~1.028）和1.013（95% CI：1.001~1.026）］。不同年龄、不同区域、不同季节间NO 2短期暴露水平与心力衰竭住院风险的关联均具有统计学意义（均 P<0.05），但各组间差异无统计学意义（均 P>0.05）。 结论:NO 2短期暴露水平升高可增加心力衰竭住院风险。

目的 探讨中国62个城市中大气CO短期暴露对心力衰竭住院的影响.方法 基于China-PEACE心力衰竭回顾研究,纳入2015年1月1日至2015年12月31日62个城市的心力衰竭住院患者,同时获取62个城市每日CO暴露浓度和气象监测数据等资料.采用广义相加模型和贝叶斯层次模型分析大气CO浓度与心力衰竭住院风险的关系,根据年龄、地区和季节对患者进行亚组分析.结果 本研究共纳入了 117 364名心力衰竭患者,在2015年1月1日至12月31日,日均大气CO暴露浓度为1.07 mg/m3.CO浓度每升高1 mg/m3,暴露当天心力衰竭住院的OR值为1.041(95％CI:1.015～1.068).在住院前0～1 d,CO移动平均暴露水平与心力衰竭住院风险的关联有统计学意义(OR=1.032,95％CI:1.001～1.064).亚组分析结果显示,不同季节、不同年龄的亚组人群中,大气CO暴露水平升高均与心力衰竭住院风险的关联有统计学意义.在北方地区CO暴露水平升高可以增加心力衰竭住院风险,而在南方地区,住院日前1 d和前2d对应的CO暴露水平与心力衰竭住院风险的统计学关联相反[OR值分别为1.102(95％C/:1.054～1.153),0.949(95％CI:0.907～0.995)].结论 大气CO短期暴露水平升高可导致心力衰竭住院风险升高.

目的 探讨双相情感障碍抑郁发作(简称双相抑郁)的中医证候要素、证型分布及证候特点.方法 收集首都医科大学附属北京安定医院门诊293例双相抑郁患者的基本信息、四诊信息,以及汉密尔顿抑郁量表、杨氏躁狂量表评分.保留出现率＞12％的四诊信息,使用R语言dist函数,以欧氏距离的方法,计算四诊信息样本之间间距.再以离差平方和法,通过hclust函数对四诊信息进行层次聚类.5位研究者组成工作组,根据聚类结果确定双相抑郁证素的病性及病位,计算患者证素因子分,并使用PC算法构建证素贝叶斯网络模型,工作组参考模型结果对贝叶斯网络中的证素进行组合,进而提取常见证型.再将患者的各证素因子分相加后按照证型内证素个数进行标准化,得到患者证型因子分,按照证型因子分对每个患者进行证型的判定,再由工作组复核,对证型与临床不符的患者进行更正,获得患者最终证型分布.结果 双相抑郁中医症状频次大于12％的四诊信息有77种,出现频率较高的前15个四诊信息分别为脉滑、神疲乏力、脉弦、多思、常独处、烦躁易怒、口舌干燥、心悸、多梦、无故担忧、胸闷、舌苔薄白、善忘、善太息、纳差.四诊信息可聚为11类证候要素,病性类证素涉及火、气虚、血虚、气逆、阴虚、湿、热、郁火、痰;病位类证素涉及心、肝、脾、胃、肾、膀胱经、胆.通过贝叶斯网络模型构建与专家组意见得出双相抑郁的常见证候类型有6种,其中心胃蕴热型占比最高,为27.99％(82例),其次为心脾两虚(55例,18.77％)、心肾不交(49例,16.72％)、肝郁血虚(42例,14.33％)、心气虚(37例,12.63％)、肝胆湿热(28例,9.56％).结论 双相抑郁的中医病性类证素以火、热类属性为主,病位证素以心、肝、脾居多,证型以心胃蕴热型、心脾两虚占比最高.

为调查山谷型水库的鱼类群落营养结构及其食物资源利用情况,基于碳(δ13C)、氮(δ15N)稳定同位素技术研究了湖北省洈水水库2020年夏季与秋季鱼类群落结构特征,并运用贝叶斯混合模型评价4种潜在碳源(颗粒有机物、陆生植物、周丛藻类和有机碎屑)对不同食性鱼类的贡献.结果表明:洈水水库鱼类群落的营养结构存在季节差异,鱼类群落的营养层次(NR)、基础食物来源(CR)、生态位总空间(TA)、群落的整体密度(NND)、聚集均匀程度(SDNND)和核心生态位空间(SEAc)等参数为夏季高于秋季,表明夏季食物源多样性高、营养冗余程度低、食物网结构更加复杂而稳定.洈水水库夏季和秋季鱼类均依赖外源性碳源,但两个季节的利用方式存在显著差异(P＜0.05).整体上,陆生植物在两个季节对洈水水库鱼类的碳源贡献最大,但夏季周丛藻类为次重要碳源,秋季有机碎屑为次重要碳源,颗粒有机物(POM)在两个季节均不重要.夏季和秋季对外源碳的依赖程度最高的鱼类均为黄尾鲴(碎屑食性鱼类),对外源碳的依赖程度最低的鱼类为鳙(滤食性鱼类)和团头鲂(植食性鱼类).在夏季和秋季共有的9种鱼类中,肉食性鱼类(鳜、蒙古鲌)、杂食性鱼类(鲤、鲫)与碎屑食性鱼类(黄尾鲴)营养级变化显著.夏季鱼类营养级范围为2.5-3.5,平均3.01,而秋季鱼类营养级范围为2.9-3.8,平均3.32.基于稳定同位素技术分析不同食性鱼类的营养级发现,秋季鱼类的营养级普遍高于夏季(P＜0.05).研究揭示洈水水库鱼类夏季和秋季的碳源利用存在显著差别,鱼类的营养级秋季普遍高于夏季,这主要是由于洈水水库秋季的低水位不利于陆生植物进入水生食物网,鱼类食物资源竞争加剧、捕食作用加强所致.基于贝叶斯混合模型分析鱼类食物资源的利用情况表明:植食性鱼类和杂食性鱼类的鱼苗投放时间应选在早春(2-3月),碎屑食性鱼类鱼苗的投放时间在秋、冬季更加合适,这将为研究山谷型水库鱼类群落的营养结构及渔业结构调整方案提供科学依据.

目的 基于实时再生数建立COVID-19预测模型并开发新冠疫情可视化平台.方法 利用上海市卫生健康委员会报告的COVID-19每日感染者数,将2022年3月1日—4月29日作为模型拟合阶段,通过贝叶斯分层模型估算实时再生数的变化情况并预测疫情发展,将4月30日—5月6日作为预测性能评估阶段,使用随机化概率分布变换和分级概率评分评估模型预测结果.使用Python、R等开源软件开发可视化平台.结果 上海市疫情的实时再生数于4月5日首次降至1以下,模型预测本轮疫情将于6月中下旬实现"动态清零",模型预测结果的概率分布和实际观察值在统计学上具有一致性(P=0.12)且集中度较好(382.53±105.29).可视化平台可实现上海市疫情变化趋势描述和疫情预测两大功能,已于2022年4月19日开放访问.结论 该模型预测性能良好,可视化平台填补了国内对COVID-19疫情数据深层次挖掘的空白,具有较好的扩展价值和应用前景.

目的 分析2013 2014年中国各地区发生人感染H7N9流感病毒事件的风险大小.方法 采用Mofan's I统计指数的方法来分析2013-2014年中国各地区人感染H7N9流感病毒的年均发病率的空间自相关特性,并用贝叶斯层次模型来预测2013 2014年中国各地区发生人感染H7N9流感病毒事件的风险大小.结果 2013年中国各地区人感染H7N9流感病毒的发病率为0.02/10万～1.44/10万,2014年发病率为0.03/10万～2.32/10万.2013-2014年中国各地区人感染H7N9流感病毒的年均发病率存在空间正相关性,因为它们的Moran's I统计指数值均＞0,均P＜0.05,空间点模式为聚集分布.2013年中国各地区发生人感染H7N9流感病毒事件的相对风险值为7.24～ 354.80,2014年的相对风险值为6.13～ 318.60.2014年中国各地区发生人感染H7N9流感病毒事件的风险相对高的地区数量相对于2013年有所减少,暴发人感染H7N9流感病毒风险的最大值和最小值也有所下降.结论 2013-2014年中国各地区发生人感染H7N9流感病毒事件的风险呈降低趋势.

目的 介绍适用于比例Meta分析的二项式-正态层次模型及其贝叶斯方法实现.方法 阐述二项式-正态层次模型和正态-正态层次模型,分别在此两模型框架下,选择随机效应模型对两个文献数据重新分析.采用贝叶斯和频率学方法拟合二项式-正态层次模型;对数据未经或经logit转换后,采用倒方差法等经典频率学方法拟合正态-正态层次模型.结果 在二项式-正态层次模型框架下,采用贝叶斯方法获得合并比例的点估计及95％可信区间分别为0.057(0.039,0.077)和0.799(0.693,0.897),研究间方差分别为0.488和0.919;采用频率学方法获得合并比例的点估计及95％置信区间分别为0.056(0.041,0.078)和0.798(0.692,0.875),研究间方差分别为0.400和0.589.在正态-正态层次模型框架下,经logit转换后获得合并比例的点估计及95％置信区间分别为0.073(0.057,0.094)和0.754(0.661,0.827),研究间方差分别为0.170和0.301;直接合并原始数据获得合并比例的点估计及95％置信区间分别为0.049(0.032,0.065)和0.804(0.719,0.888),研究间方差分别为0.001和0.018.结论 不同模型可以获得不同的结果,NNHM可能低估研究间方差.基于二项式-正态层次模型框架下贝叶斯方法更适用于比例的Meta分析.

目的 鉴于国内研究者熟悉的贝叶斯统计软件WinBUGS已停止更新,OpenBUGS已见介绍,本文介绍国外新近出现的贝叶斯统计软件Stan的基本原理与实际使用.方法 由于Stan安装具有其自身特点且较为复杂,本文首先介绍Windows平台下Stan安装与使用方法,接着介绍Stan语言及运行步骤,最后以线性-正态层次模型为例介绍Stan的应用.结果 Stan作为一种新型贝叶斯统计软件,使用了全新的汉密尔顿蒙托卡罗(Hamiltonian Monte Carlo,HMC)抽样方法,可以处理GIBBS抽样难以收敛的复杂多元层次模型.Stan使用的概率模型定义语言,相比较BUGS,逻辑清晰,程序更易理解.结论 Stan功能强大,应用方便,是贝叶斯分析的有力工具.

目的 介绍单臂试验连续型数据的Meta分析模型、贝叶斯方法及实现.方法 阐述正态-正态层次模型,基于该模型框架,以贝叶斯方法拟合随机效应模型,对效应参数μ和异质性参数τ分别选择不同的先验,使用R软件的bayesmeta包对两个文献数据重新分析.结果 在正态-正态层次模型框架下,基于不同的先验信息,贝叶斯Meta分析结果为:数据1参数μ的点估计及95％CI分别为-4.26(-6.97,-1.92)和-4.50(-9.27,-0.53),参数τ点估计及95％CI分别为1.51(0.41,2.75)和2.28(0.00,6.57);数据2参数μ的点估计及95％ CI分别为-4.07(-5.54,-2.71)和-4.12(-5.96,-2.46),参数τ点估计及95％CI分别为1.54(0.78,2.48)和1.81(0.74,3.51).结论 不同的先验可能影响参数估计值.基于NNHM框架下的贝叶斯方法适用于单臂试验连续型数据的Meta分析.Bayesmeta包以其简单、快速、准确、可重量性算法等可以用于实现贝叶斯随机效应模型Meta分析.