肾综合征出血热与气象因子相关性分析及其预报研究

目的 根据连云港赣榆县1990～2002年的气象资料及肾综合征出血热(HFRS)的发病情况，探讨气象因子对HFRS高发区发病率的影响及其相关性，并建立基于气象因子的HFRS发生的预测预报模型。方法用数据处理软件SAS的相关性分析工具和多重线性回归模型。结果以湿度(x1)、日照(x2)、降水(x3)、平均温度(x4)和最低气温(x6)等气象因子为自变量，HFRS发病率为因变量，采用强迫引入法和逐步回归法，分别建立回归方程，y=8．766 0．148x1 0．015x2-0．028x3-0．567x4(①)，复相关系数(r)=0．45及y=1．742-0．054x6(②)，复相关系数(r)=0．44，经t检验，呈显著相关。结论方程②可作为HFRS的预测预报模型，同时进行了区间估计，证实可行。     

误差为平稳m—相依的线性模型的移动区段再抽样法

引　　言对于一般的线性回归模型Ｙ =Ｘβ ε,假定误差ε独立且服从正态分布 ,才能进行有效的参数区间估计及假设检验 ,当ε的分布未知时 ,若按通常的方法 ,参数的区间估计与假设检验就不稳健 ,1981年Ｆｒｅｅｄ ｍａｎ将Ｂｏｏｔｓｔｒａｐ法应用回归模型 ,提出参数及其方差的Ｂｏｏｔｓｔｒａｐ估计 ,1986年Ｗｕ又进一步加以改进 ,使方差的Ｂｏｏｔｓｔｒａｐ估计达到稳健 ,假设检验也能够顺利实现。然而 ,当误差ε或响应变量序列存在一定的相依关系时 ,上述方法就失效 ,而误差或响应变量序列为相依情形的线性模型是广泛存在的 ,比如目前…     

用ＭＣＭＣ方法研究统筹医疗保险的损失分布模型

在传统的统计模型分析中,仅限于反应变量或误差变量服从正态分布时才有较好的估计性质,当反应变量服从非正态分布时,多采用变量变换使其服从正态,且通常Box-cox变换（y=(1)(y^λ－1)/λ (2)iny λ≠０,（３）λ＝０参数λ的估计一般是采用迭人大似然法^[1]）但变换后资料常有服从正态。例如医疗保险的损失分布常为对数正态、Ｗeibull或Pareto分布^[2],其损失额的分布均是非负的偏态分布,且同一级别医院的医疗费用间存在一定相关即非独立,对这类非独立偏态分布数据,宜用Ｍarkov chain Monte Carlo(MCMC)法^[3]来处理,本文应ＭＣＭＣ法来探讨统筹医疗保险的损失分布模型。     

计数资料的统计分析模型

当反应变量为定性变量时 ,试验或观察结果称为计数资料 ( countdata)。计数资料在医学研究中很常见。例如发病率的调查结果 ,临床试验中疗效的评价等。定性变量的常见分布类型有二项分布、多项分布、Poisson分布、负二项分布等。表 1列出它们的概率函数、期望与方差。表 1　定性变量的常见分布分布类型概率函数期望 E(y)方差 V(y)二项分布B(p,n) p(y=1) =(rn) pr(1-p) n- r np np (1-p )多项分布 p(r1 ,r2 … ,rm) =n!r1 !r2 !… rm!p1 r1 p2 r2 … pmrk npi npi(1-p) iPoisson分布 p(y) =λyy!e-λ λλ负二项分布 p(y) =Γ (k+ y)Γ (k)…     

在广义线性混合模型框架下多结局测量的多变量meta分析

目的基于SAS PROC NLMIXED探讨在广义线性混合模型框架下实现多结局测量的多变量meta分析的方法。方法Stijnen等人提出一种基于广义线性混合效应模型(generalized linear mixed model,GLMM)的确切的研究内似然模型(exact within-study likelihood models,EWLM)meta分析技术,可以极方便地实现多结局测量的多变量meta分析。本文重新对Voko等人的荟萃分析数据进行分析,该数据旨在研究颈动脉狭窄患者采用颈动脉内膜切除术或保守治疗对联合终点卒中或死亡的影响。结果基于PROC NLMIXED拟合的多变量随机效应模型的meta分析显示,手术组术后第1个月内的事件发生率为6.31%,SE=0.0083,95%CI(0.0456,0.08063)。手术治疗一个月后的事件发生率为5.5%,SE=0.0024,95%CI(0.04953,0.05962);而对照组事件发生率为7.45%,SE=0.0058,95%CI(0.06231,0.08672)。组间比较比值比效应值OR=0.732,SE=0.0498,95%CI(0.6269,0.838),P=0.0003。结论多变量分析比单变量分析具有更大的潜在信息量,基于广义线性混合效应模型的确切的研究内似然模型meta分析技术可以成为多结局测量的优选方法。     

湖沼地区预测钉螺数量的模型研究

目的探讨并建立湖沼地区钉螺数量的预测模型,为钉螺的定量化研究提供方法学依据。方法在安徽省贵池区秋浦河沿岸随机数字表法随机选择滩地作为研究现场,根据植被类型分层随机抽样,以交叉复核随机抽检法（随机数字表法）调查钉螺,分别用差分GPS仪、T＆D Recorder for Windows、卷尺及烘干称重法收集高程、土壤温度（简称土温）和气温、植被高度、土壤湿度（简称土湿）和植被类型6个环境变量信息。对数据进行预处理,通过偏差量和Akaike信息准则比较不同误差分布和连接函数组合的广义线性模型,确定最佳模型结构,建立并验证预测模型效果。结果建模样本量为162框,变量之间存在着复杂的相关性,钉螺数与植被高度呈正相关（r=0．36）,与土湿呈负相关（r=-0.22）,气温与土温呈正相关（r=0．59）,土温与植被高度呈负相关（r=-0．36）,土湿与土温和气温均呈负相关（r=-0.34和-0．12）。广义线性模型的最佳结构是以gamma分布为误差分布、倒数为连接函数和均数平方为方差函数的模型结构。模型拟合结果显示高程、土湿、土温、植被类型和植被高度对于预测钉螺的数量有统计学意义,t值分别为-3．202,3．124,-1．989,2．668和-2．371,P值分别为0．00166,0．00214,0．04849,0．00846和0．01897,而气温的作用被土温取代没有进入模型。结论广义线性模型可用于建立钉螺的预测模型,为钉螺的定量化研究提供了广阔的研究前景。     

Logistic回归模型使用注意事项和结果表达

Ｌｏｇｉｓｔｉｃ回归模型是流行病学研究中用途较广的多因素分析方法。然而 ,有时人们因不参很好地理解模型的假设条件而误用 ,或因不能明确地表达和解释结果造成读者的误解。Ｌｏｇｉｓｔｉｃ回归模型简介设Ｄ为两项随机变量 (ｂｉｎｏｍｉａｌｒａｎｄｏｍｖａｒｉａｂｌｅ) ,Ｄ =1表示在特定的观察期间内 (ｔ0 到ｔ)发病 ,Ｄ =0为未发病。可能和疾病发生有关联的自变量有 χ1,χ2 ,… ,χｋ。在ｔ0 时测量自变量。疾病发生的条件概率 (ｐｒ)为 :ｐｒ(Ｄ =1| χ1,χ2 ,… ,χｋ) =Ｐ(ｘ)对Ｐ(ｘ)作“ｌｏｇｉｔ”转换变换 :ｌｎ[Ｐ(ｘ) /…     

尘肺患者生存年限的COX比例风险模型分析

目的探讨影响尘肺病患者生存年限的因素。方法以尘肺患者的死亡风险函数为应变量，生存年限为暴露时间变量，发病工龄、工种、是否并发肺结核、Ⅱ期晋升为Ⅲ期的年限及尘肺种类为方程自变量，运用STATA统计分析软件进行尘肺病例生存年限的COX比例风险模型分析。结果COX比例风险模型方程为：λ（t）=λ0（ε）exp（-0．0249×fb＋0．0344×gz＋0．2059×tb-0．0215×nx2—0．0280×cf），矽肺相对于煤工尘肺的COX比例风险模型方程为：λ（t）=λ0（t）exp（0．5289×cf）。由COX模型方程可得出各自变量的相对危险度（RR）值分别为1．025、1．035、1．229、1．022、1．697，提示，发病工龄对患者死亡的影响呈递减的形式，发病工龄每减少1年，尘肺病例在某时刻死亡的RR增加1．025倍；井下采矿工尘肺病例在某时刻死亡的RR值是非井下采矿工尘肺病例的1．035倍；尘肺并发肺结核在某时刻死亡的RR值是单纯尘肺的1．229倍；Ⅱ期尘肺晋升为Ⅲ期的年限每减少1年，尘肺病例在某时刻死亡的RR值增加1．022倍；矽肺病例在某时刻死亡的RR值是煤工尘肺的1．697倍。结论影响尘肺病例生存年限的因素有发病工龄、工种、是否并发肺结核、Ⅱ期尘肺晋升为Ⅲ期的年限及尘肺种类。     

中国不同地区群体遗传结构差异及调整策略研究

目的 描述中国不同地区群体遗传结构特征,探索并评价不同分析方案控制队列样本群体遗传结构混杂因素的效果。方法 通过中国慢性病前瞻性研究（CKB）队列10个地区4 500例样本的全基因组关联研究数据,通过主成分分析提取样本第一、二主成分,绘制主成分二维图,并与样本地区来源相比较,分析我国不同地区样本的遗传结构特征。以CKB队列数据为基础,生成存在遗传结构差异、亲缘关系等队列样本特征的模拟数据集,探索并评价不同分析策略对膨胀因子（λ）的控制效果。结果 我国不同地区人群存在显著的群体遗传结构差异,人群遗传结构主成分分布与项目地区的地理分布基本一致,第一主成分对应不同地区的纬度,第二主成分对应不同地区的经度。生成的模拟数据集,直接进行关联分析假阳性率较高（λ=1.16）,即使调整遗传结构主成分或根据地区进行亚组分析仍无法有效控制λ（λ>1.05）;使用混合线性模型引入亲属关系矩阵作为随机效应量后,无论是否进一步调整遗传结构主成分,λ均得到有效控制（λ=0.99）。结论 我国不同地区人群遗传结构存在较大差异,在分子流行病学研究中需要谨慎处理群体遗传结构造成的研究偏倚;针对大队列数据遗传结构复杂、亲缘关系广泛等特征,需要使用混合线性模型进行关联分析。     

用改善的自回归预测方法预测门诊人数

时间序列资料,建立自回归模型进行预测时,为了提高预测的准确性,要求数据列{x(k)}(k=1,…,n)为光滑离散函数,其主要条件是:ε>0,k0,当k>k0时,有x(k)/∑k-1i=1x(i)<ε〔1〕。即,当数据列{x(k)}与{y(k)}满足y(k)/∑k-1i=1y(i)0,y(k)=ln(x(k) x(2k) 1),则数据列{y(k)}比数据列{x(k)}光滑。即利用反双曲正弦函数变换提高了数据列的光滑程度。[证明]令g(x)=x-x2 1.ln(x x2 1),则g(0)=0。当x>0…     

简单催化模型拟合方法的探讨

用微分拟合法并借鉴灰色均值生成理论来拟合简单催化模型 ,取得了较好的效果。1　原理和公式简单催化模型微分方程为 :dydt=r(k- y) (1)y代表在给定年龄组中有阳性病史的比率 ,r为在单位时间内每个个体的有效接触数 ,t为时间 ,k代表所有经过有效接触且试验结果为阳性的比率。满足 t=0 ,y=0的解为 :y=k(1- e- rt) (2 )根据微分拟合法 ,取 dy=Δ y,即 y的累减值 ,dt=Δt,即年龄组中值的累减值。令 I=Δ yΔ t代入 (1)式则有 x=rk- ky,可用最小二乘法求解参数。2　举例及效果评价某地百日咳曾患率与年龄关系资料 ,试作催化模型拟合。表 1　某…     

影响三氯乙烯药疹样皮炎患者住院时间多因素分析

目的 探讨影响三氯乙烯药疹样皮炎患者住院时间的因素.方法 对120例三氯乙烯药疹样皮炎患者的临床资料进行回顾性分析,用多元线性回归分析筛选影响患者住院时间的因素并建立多元线性回归方程.结果 120例患者的住院天数回归方程为:y=54.132+1.027x1 +0.014x2-14.687x3-12.851x4.该模型的拟合系数为0.429,调整后的拟合系数为0.407.模型拟合情况一般.结论 发病到入本院时间、住院期间糖皮质激素使用总量、性别、病情分型为影响三氯乙烯药疹样皮炎患者住院时间的因素,其中前两者为正相关.     

线性结构模型及其医学应用

目的探：讨线性结构模型理论及其在医学中的应用。方法：通过建立线性结构模型，对医学实例进行分析。结果：由线性结构模型对医学实例的分析获得了较为满意的结果。结论：在生物医学领域中，当变量间具有较强的相关性时，以及当代医学模式主张“多因-多果”病因学观点的情况下，线性结构模型是一种值得推广应用的统计分析方法。     

赣东北地区父母健康素养与中学生健康危险行为关联

  目的  探究父母健康素养与中学生健康危险行为的关联,为提高父母健康素养和减少中学生健康危险行为提供理论依据。  方法  采取多阶段分层随机抽样的方法,于2020年6—9月对江西省东北区域1 993名中学生及其父母进行电子化问卷调查,使用Spearman相关分析父母健康素养与中学生5类不健康行为的相关性。使用广义线性模型分析父母健康素养评分与中学生5类健康危险行为的关联。  结果  广义线性模型分析显示,调整家庭经济状况等变量后父母健康素养与中学生健康危险行为关联,科学健康观(β=-0.05,95%CI=-0.09~-0.01)、健康信息素养(β=-0.06,95%CI=-0.08~-0.05)均与中学生物质滥用与暴力行为关联;健康信息素养与中学生健康饮食行为(β=-0.03,95%CI=-0.04~-0.01)、中学生运动与早餐(β=-0.03,95%CI=-0.05~-0.02)关联;慢性病防治素养(β=0.05,95%CI=0.00~0.09)、健康信息素养(β=-0.02,95%CI=-0.03~-0.00)与中学生视屏行为关联;科学健康观(β=-0.05,95%CI=-0.09~-0.01)、基本医疗素养(β=-0.08,95%CI=-0.13~-0.04)、健康信息素养(β=0.02,95%CI=0.00~0.03)均与中学生交通安全行为关联(P值均 < 0.05)。  结论  父母健康素养与中学生健康危险行为之间存在一定的关联,需要在队列研究中明确关联的方向和强度。     

COX模型在三硝基甲苯中毒与工龄、工种关系中的应用

目的探讨三硝基甲苯作业工人工种、工龄与职业性慢性三硝基甲苯中毒发生的关系。方法选择2005年某化工厂职工健康监护体检资料，使用STATA7．0统计软件进行COX比例风险模型分析。结果进行Log-rank定性检验及COX比例风险模型求证，以接尘工龄作为暴露时间变量t，COX比例风险模型为λ（t）=λ0（t）exp^（1.354588）（工种）。结论工种间的相对危险度比，RR=3．8752，一线作业工人的慢性三硝基甲苯中毒发生的RR是二线工作的工人的3．87倍，应用COX比例风险模型分析同类资料有指导意义。     

血肌酐测定方法的比较

目的：BeckmanCX7全自动生化仪在检测血肌酐时电极法与化学法的性能。方法：分别采用两种方法测定血肌酐值进行精密度试验,以化学法为比较法,电极法为实验方法,对两种方法的相关性作比较分析,同时按NCCLS线性评价方法进行线性评价试验。结果：(1)精确度：电极法变异系数为0．7％～2．5％．化学法变异系数为1．79％～5．63％;(2)相关性：y=0．983x-0．098 r=0．995;(3)线性：电极法和化学法测定血肌酐的线性范围分别为0．2250umol／L及0．2000umol／L。结论：两种方法的精密度、线性均符合临床要求。     

影响患者住院天数超长的因素分析

目的研究影响某院患者住院天数超长的因素,为缩短患者住院的天数,提高医院管理水平提供依据。方法通过提取某院2012年11月至2013年10月住院天数超过90天的患者的病历资料,共220例,应用多元线性回归模型方法,对住院天数超长的影响因素进行分析。结果性别、手术情况、转科情况、病例分型、入院病情和院感等变量,在回归方程拟合过程中被剔除,最终引入回归方程的是年龄、付款方式和病种3个自变量。模型的决定系数R。=0．071,方差分析F=5．487,P〈0．01,说明以上变量所拟合的方程有统计学意义,各自变量对住院天数的贡献依次为病种〉年龄〉付款方式。结论病种、年龄和付款方式是患者住院天数超长的影响因素,加强安全教育,减少损伤事故的发生,提高医疗质量的管理水平,缩短患者住院的天数,可减少患者的负担。     

重复测量线性混合模型中方差-协方差结构的选择

线性混合效应模型在医学重复测量资料的数据分析与设计中广泛应用。本文围绕医学实例阐述了重复测量线性混合效应模型在对不同方差 -协方差结构建模时协方差结构的选择 ,并对结果进行分析 ,为重复测量线性混合效应模型在医学研究中的应用提供了方便     

1990～1999年我国居民糖尿病死亡率分析及预测

目的分析 1990～ 1999年我国城市和农村居民糖尿病死亡率的变化并进行趋势预测。方法采用 1990～ 1999年《全国卫生统计年报资料》,对我国 1990～ 1999年糖尿病死亡率进行流行病学分析 ,并利用灰色动态模型 GM(1,1) ,预测我国 2 0 0 0～ 2 0 0 5年城市和农村糖尿病死亡率趋势。结果 1999年 ,糖尿病死亡率城市为 15 .37/10万 ,农村为 5 .13/10万 ;分别较 1990年上升 1.89倍和 1.71倍 ;女性高于男性 ,城市和农村男女之比分别为 0 .6∶ 1和 0 .7～ 0 .8∶ 1;所建立的灰色模型预测方程城市为 Xt=- 10 6 .4 817(1- e0 .0 74 58t) y0 ( t- 1 ) (t=1,2 ...,N) ;农村为 :Xt=- 33.6 0 5 3(1- e0 .0 7537t) y0 ( t- 1 )(t=1,2 ,...,N) ;预测到 2 0 0 5年我国城市和农村糖尿病死亡率将分别达到 2 5 .2 4 /10万和 8.15 /10万 ,精度检验证明 ,城市和农村预测模型拟合良好 ,外推测值是可信的。结论 1我国人群糖尿病死亡率上升趋势明显 ,提示有关部门应加强防治力度。 2灰色模型 GM(1,1) ,可作为糖尿病死亡率近期预测使用。     

选择"最优"多元线性回归模型的信息熵法

对同一组线性相关数据因变量Y与自变量X,有多种方法和多项指标用于选择“最优”多元线性回归模型。本文将熵(信息熵)〔1〕的概念引入线性回归模型。一个回归模型方程可视作一个系统,将系统的元序或混乱程度〔2〕与模型复杂性和拟合优度相联系,从而提出选择“最优”多元线性回归模型的信息熵法。原理与方法熵是一个系统无序或混乱程度的度量。将一个回归模型视作一个系统,信息熵最小的系统为不确定(无序)程度最小,混乱程度最低的系统,即欲得“最优”回归模型。因变量Y与自变量Xi(i=1,2,…,m,m为全模型自变量数目)呈线性相关。对于给定的准…