【doc】适用于一类半岛全站医学动态监测数据的稳健参数估计方法

　　半岛全站【doc】适用于一类医学动态监测数据的稳健参数估计方法doc

　　适用于一类医学动态监测数据的稳健参数估计方法260?中国卫生统计1999年10适用于一类医学动态监测数据的稳健参数估计方法西安交通大学生物医学工程研究所(7loo49)三—墨-郜索勋西安医科大学第二临床学院龚均,/张茹【撮要】目的介绍并分析适用于一类医学动态监测数据的稳健参效估计方法.方法结合分布假设和稳健统计给出了参数估计方法.首次提出一种选代求取形状参数稳健估计的算法.结果该方法对医学动态监涮中某生理量异常所占百分比值统计量的参数估计优于临床普遍采用的方法.结论通过对分布拟合,估计优效性和可实现性,估计稳健性等理论分析和模拟实验,表明该方法是一种高教,可靠的稳健参数估计方对于正偏态分布的总体,采用对数转换或其他趋势曲线…将其正态化.从而用样本均值和标准差来估计总体期望和方差的方法虽比较普遍,但却缺乏理论依据.且在很多情况下效果一般.Weibuu分布具有与正偏态样本拟合优度好,信息量丰富等特点.许多医学统计量服从或近似服从该分布】.一般的参数估计方法只考虑理想情况(样本独立同分布于理想模型)下的估计效果,但实际上由于总体不一定与假设严格相符,样本未必来自同一总体等原因口】,必须要考虑估计方法的稳健性.本文针对所举统计量,通过拟合实验进行webul1分布的假设检验,并且设计出一种合适的稳健参数估计方法.方法Weibul1分布的概率密度函数为为位置参数,0为形状参数.po为刻度参数.期望EX=r(1+)?,.+,方差DX=(r(1+2/)一r(1+1/))(r()表示Gatllma函数).当a=l时Weibul1分布化为指数分布.首先分析分布假设.数据来自1994～1999月西安医科大学第二临床学院消化内科对主述烧心,胃食道返流等症状的就诊者所做的24小时食管pH动态监测,共165例(其中男86倒.女79例,平均年龄4l岁),取每例数据pH4所占检测时间百分比值作为样本(目前临床将食道pH4判断为胃食道返流,且认为pH4占监测时间百分比值对判断胃食道返流具有最好的特异性】.样本频数图直观地表明样本呈明显正偏态——单边厚尾分布.绘在正态概率纸,wdbuu概率纸上的拟合图表明样本与正态分布的拟合程度差,与Weibul1分布的拟合程度较好.t检验,Kolmogorov检验,峰度——偏度联合检验等检验方法进一步肯定上述分析.另外由统计量得样本与Weibul1分布的拟合程度要优于与对数正态分布的拟合程度.样本的频效图尽管计算机数据采集技术可将读数错误,记录错误等减小到最低限度,但由于监测时间长,监测环境复杂等原因,样本中必定存在一定数量的异常值],对异常值分布我们所知甚少.同时考虑到前述分布假设及对其的猜测成分,使得简单应用样本均值和标准差的参数估计方法必将存在较大误差.我们结合次序统计量和稳健统计,在假设总体具有有限期望和方差的*国家八五科技攻关项目及国家科委资助项目,项目编号925—30—04Ao前提下,设计如下的稳健参数估计方法:对样本l,Y2…,YN从小到大排序得[1】,[2,N.计算对较大一侧单边截尾口后的样本一,二阶原点矩统计量二=__?-三NO,二:=1?(),令三=三2/},之后利用等式r(1+2/:)/rz(1+1/:):.(为与:有关的无偏因子)选代求取ao最后令三.:..ul(为子),CDY=N.(r(1+2/:)一P2(1+1/a(r(1+2/a)+(N一1)./2(1+1/a),得到三/r(1+1):,由)=:.+3)及由=(1).估计结果见下表1.稳健参数估计结果讨论百分比统计量的有效取值范围应为01,而261Weibul1分布的定义域为o.,两者并不一致.但考虑到由定义域变宽带来的误差很小,如当n=0.5,口=0.2,=0时,F(1)竺0.995.为方便起见,我们并未修正Weibul1分布的定义域.不同的参数估计标准会产生极大似然估计,线性无偏和线性不变估计,矩估计等不同估计方法】.极大似然估计需求解联立超越方程,计算和分析有很大不便;线性估计方法则在很大程度上依赖不稳健的极大值分布,另外当样本值为0时其极值变换函数=ln(f)无意义.相对而盲,矩估计法是容易实现,便于稳健分析且当样本容量较大时效果较好的估计方法.上述分析的模拟实验见表2(WeibuU分布样本由随机数发生器产生,分为150250和250015两组.不模拟实验结果(实际值为所列值乘10)总体服从W(,,时,【1)服从W(口,p/N,),当N—o..fl_将收敛到,即渐进无偏.另外,通过引入无偏因子可以保证n和等估计统计量渐进无偏.故上述估计统计量的MSE应与渐进方差相同.^1N对于矩法,EY=?三Yi的方差达到Rao-Cramer界/N,是最小方差无偏估计统计量;DY的估计统计量有c)1Y=CDYCDYI?(?一(?.萋)),D:y=计量有c(南.善一(南..曼t)),DzycDY21?.),Dayy南.看;(cDv因子)等,当总体服从We]bull时,三者的渐进方差比见图2,我们选择在n1时渐进方差最小的做为DX的估计统计量.二价矩估计统计量的方差比图总体服从w(1,1/^,O)时,f(z)=^?N/(N一计误差小于O.01,对其他较复杂的情况,模拟实验也1)!?xN-1?一～,7O,为均值统计量的概率密得到相似结论.另外,对本文所举统计量来说,样本容N100时该统计量可以95%的概率保证估量过大会导致收集时问的延长,收集手段的不一致262?等,反而影响到样本的平稳性和参数估计的效果.我们认为值的选择是合适的.从泛函角度分析,EY,DY等几个估计统计量的影响函数是广义有界的,而连续有界的影响函数对应对离群值不敏感及截断误差小的稳健估计统计量J.另由定义域知估计统计量=Y.】]稳健ian(x)和medianlx—median(x)l,Winsor【3也都属于稳健统计量,计算它们与无偏矩估计法对污染和未污染的Weibul1分布样本的参数估计的渐进方中国卫生统计1999年10月第16差,表明我们选择的单边截尾均值和单边截尾二阶矩的效果最好(表3).对Weibul1分布中形状参数的估计是对其他参数估计的基础.由于分析繁复,无偏因子与分布参数有关等原因,形状参数的稳健估计鲜有报道.我们在矩估计法基础上设计的迭代算法,利用了weibul1分布中形状参数的无偏因子仅与形状参数有关,形状参数的矩估计函数的单调性等特点,迭代求取形状参数的无偏稳健估计.由表3,与直接求取(不迭代)相比估计统计量的MSE要小一个数量级以上.不同总体下几种方法对:的渐进方差(实际值为所列值乘10一一)结论本方法适用于医学动态监测中某生理量异常所占百分比值的统计量,这类统计量一般具有明显的正偏态分布,多数样本集中于较小一侧等特点.除文中分析的统计量外,24 小时动态监测心电数据的ST 常累计时间占监测时间百分比值(该值有重要临床意义)也属于此类.由于较好地结合了分布假设与稳健 统计,本方法对此类统计量的参数估计优于临床普遍 采用的方法.理论分析,模拟实验和实际应用表明其 不失为一种高效,可靠的稳健参数估计方法.另外,随 着动态监测技术的日益成熟和临床统计手段的深入发 展,相信本方法会有更多的实际应用. 估计统计量的稳健性和效率是对立统一的关系, 我们的目的是在可实现基础上寻找最优折衷,虽然一 些具体的稳健方法已在临床上得到应用.但由于稳健 统计的发展相对较短,灵活性大,这类应用还不十分成 熟,我们所做的理论分析和模拟实验能为稳健统计提 供一定的临床应用基础.接下来我们可通过搜集更全 面的资料(包括正常组与确诊组的对照)来完善估计方 法.在得到进一步结论的基础上可涉及假设检验等与 置信概率有关的问题.从而针对临床比较系统地引入 稳健分析. TheRobmtParameterEstimationMethodAppiledtoaClz~sof MedicalDynamicMonitoringDataWangLei.ZhengChongxun HJiaoTongUniversityBiomedicalEngineeringInstitute(710049).Xi 【AbstractlObjectiveArobuperaInerestimatJonmethod apedtoacLassofmedicaldynamicmonito~ngdataisimmdueedand iyzedinthispaperMethodsThispat~eterestimationmethodisintro— ducedbybeingcombinedwithdistributionhypothesisaodrobuststatistics Theiteration 目d 口sI methodis0fhighcredibilityandeHi~iencyafterthe[henryarm3mdtrim—uhtioat 哪ofthegoodaess-of-[itdiscrlbution.meritaodworka~lJ.tyaod robustnessofestuni0nstatic. 【KeywordslRobustnessParameterestimationDynamic MonitoringGoodness-of-fitofdistribution 参考文献 1.上海第一医学院卫生统计学教研组编医学统计方法.上海:上海科 学拄术出版社,1979,160 周士楷威布尔分布在正偏态分布资料灿理中的应用中国卫生统计,1988.5(6):26 3F～kR.HampeLRobustEstimati~JohnWiley&Sonslne. 1986 4.FBa[di,etTechnics[ptsdimrMumitmlpH- metryinm叽:CDiWeDt Sla[USandrecommendationsGut,1987.28:l177-1188