没有合适的资源?快使用搜索试试~ 我知道了~
数据回归-保序回归的算法及应用.pdf
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 160 浏览量
2022-06-24
23:44:08
上传
评论
收藏 990KB PDF 举报
温馨提示
试读
33页
数据回归-保序回归的算法及应用.pdf
资源推荐
资源详情
资源评论
保序[口1归的算法及麻用
第一章绪论
1.1保序回归的背景介绍
在药物剂量反应中,我们感兴趣的是,随着药物剂量的增加,疗效或副作用
通常会呈现一定的趋势,比如:剂量越高,疗效越好;剂量越高,毒性越大等等。
临床试验研究的主要目标之一就是评估药物在不同剂量水平下的毒性,并且
建议一个对病人既安全又有效的剂量,即最大耐受剂量(Maximum
Tolerated
Dose),简记为MTD。固定价递增的剂量水平,假定随着剂量的增加,药物的毒
性概率是非减的。预先给定毒性靶水平,所谓靶水平是可接受的最大毒性概率,
MTD被定义为毒性概率不超过毒性靶水平的最高剂量水平。连续招募病人进行试
验,直到确定最佳估计的MTD。
临床试验通常假定随着剂量的增加,药效和毒性也是增加的。
然而基于每
个剂量水平下病人毒性反应的比率估计不同,剂量水平下的毒性概率可能不是剂
量水平的非减函数。于是我们可以采用保序回归方法(参见Ayer
et
a1.(1955),
Barlow
et
a1.(1972),Robertson
et
a1.(1988))。一般地,设有价递增的剂
量水平西≤吐≤…≤农。令Ⅳ,和R,分别表示在剂量水平d,下的试验人数和出现
中毒反应的人数。对于任一低于Z的剂量4(,.≤f)和任一高于Z的剂量以(s≥f),
记
则毒性概率B的保序回归估计为
∑q
所.I',=}
∑M
j=r
A
2罂嘴P,,,声
这就是保序回归的背景及其实际应用,其计算方法主要是PAVA(Pool
Adjecent
Violators
Algorithm)法,该算法是Ayer
etal.(1955)提出的。
江西师范大学硕上学位论文
1.2保序回归的发展现状
约束条件下的统计推断已经成为统计分析的一个重要领域,保序回归是约束
条件下的统计推断的典型问题。保序回归的研究从20世纪中叶就已经开始,到
六十年代,保序回归得到了广泛的重视。Bartholoremew,Barlow等对此进行了
研究,并且于1972年联合出版了保序回归的第一本专著《The
theory
and
application
of
Isotonic
Regression},使得保序回归问题得到了进一步的重
视和发展,成为热门话题。
保序回归是在约束条件下的一种回归,它的一个很重要的计算方法是PAVA
算法。它是以全序为例给出的算法,是Ayer
eral.(1955)提出的。对于其他的
半序,也在此基础上给出了类似的算法。比如,对于伞型半序,Shi(1988a)提出
了一个基于PAVA算法求保序回归的方法。与上述算法有关的子程序,可以在
gran(1980),Geng
and
Shi(1990)找到。另外,Robertson
et
a1.(1988)给出了
保序回归的最大最小公式;Michael
J.Best
and
Ni
lotpal
Chakravarti(1990)
也给出了MLS(Minimum
Lower
sets)等算法并证明其所得解为保序回归的解。
在保序回归和最大似然估计这方面的工作中,对二项分布、几何分布、泊松
分布和r分布等都进行了详细的讨论,这些可以在Robertson
et
a1.(1988)中
找到,并总结出:指数分布族下,全序约束下的最大似然估计就是无约束的最大
似然估计的保序回归。Shi(1993)也对此进行了研究,Shi在文中还对保序回归
的优良性进行了讨论。刑务强(2004)对保序回归解的大样本性质给出了详细证
明。关于指数分布,也有讨论,但结果很不理想(Vijayasree
and
Singh,1991)。
关于正态均值和标准差的比也有相关讨论(Chou
and
Owen,1991)。
多维保序回归是由Sasabuchi
et
a1.(1983)提出的,Shi
and
Geng(1992)
给出了相关的一般算法,同时还给出了一个子程序。Dykstra
and
Robertson(1982),
Robertson
et
a1.(1988),
Burdakov,Grimval
1,and
Hussian(2004)讨论了一类多变量全序约束下的估计问题,类似一维保序回归,
对一类指数分布族,求多维保序回归的最小解与求参数在约束下的MLE是一致
的。
在现有的研究成果中,还对广义保序回归进行了讨论。在正态分布下,求广
义保序回归的最小解与求带约束的正态均值的MLE是等价的,更为一般的结果可
2
保序pI归的算法及戍用
以在Kudo
and
Choi(1975),Wynn(1975),Nomakuchi(1983),Shapiro(1985,1988),
Raubertas
and
Lee(1986),Shi(1988,b)中找到。
保序回归问题也越来越多的被引入到二次的凸规划问题中,如Best
and
Chakravarti(1990)and
Best,Chakravarti,and
Ubhaya(2000)。近年来,应
用保序回归解决数学规划问题在Ahuja
and
Orlin(2001),and
Hansohm(2007)中
也可以找到。
以上提到的保序回归主要为经典的保序回归,即厶保序回归,在Shi(1993)
中最后给出了保序回归的一个自然推广:工。保序回归,文中提到,针对厶保序
回归是样本中位数的函数,而样本中位数比样本均值具有稳健性(Stigler,
1973),因此设想厶保序回归比三,保序回归更具有稳健性。
本文中我们就是要讨论厶保序回归,针对厶保序回归给出PAVA算法,并详
细的证明了该迭代算法收敛到的值就是对应最优化问题的解,类似厶保序回归,
发现厶保序回归与约束条件下双边指数分布的最大似然估计等价,并且证明其参
数估计值具有相合性和渐近正态性;由厶保序回归和厶保序回归的PAVA算法思
想及该两种保序回归与最大似然估计的关系,本文在最后针对带约束的一类最大
似然估计的求解,给出了一个比较直接的迭代方法。
1.3本文的工作意义与工作安排
6‘
基于已有的相关成果,本文在厶保序回归的基础上进一步推广,工作意义
是:给出厶保序回归的PAVA算法和详细的证明过程,其所得解比三:保序回归解
更具有稳健性。类似£,保序回归,给出厶保序回归与最大似然估计的关系,在
密度函数不可导的情况下简化最大似然估计值的求解过程,并且能够求出全序条
件下极大似然的估计值,同样也证明出其所求解具有相合性和渐近正态性。由第
二章与第三章的结论我们知道,厶保序回归与厶保序回归的算法都是利用PAVA
思想,因此在第四章中直接将PAVA与最大似然估计的求解联系起来,针对一类
全序约束下最大似然估计给出了更直接的求解方法。
具体工作安排如下:
1、厶保序回归与最大似然估计。主要是介绍厶保序回归的定义、算法和定
理,并且在一类指数分布族下,带全序约束的最大似然估计是无序约束的最大似
江两师范大学硕士学位论文
然估计的保序回归,特别地,针对厶保序回归解的优良性和大样本性质也做了
相关的介绍。因此根据保序回归与最大似然估计的关系,可以求出许多全序约束
模型下的参数估计和非参数估计,最后将保序回归应用到I型区间删失下生存函
数的估计问题中。
2、厶保序回归的算法及其应用。首先给出厶保序回归的概念及思想,然后
给出厶保序回归的算法并证明该迭代算法收敛到的值就是对应最优问题的解。最
后再给出厶保序回归问题的应用,我们发现全序约束下的双边指数分布的最大似
然估计等价于求厶保序回归的问题,所得到的参数估计值同样也满足相合性和渐
近正态性。
3、PAVA算法的应用及总结。正态分布下带约束的总体均值的MLE为无约束
总体均值的MLE的厶保序回归,双边指数分布下带约束的总体中位数的MLE为
无约束总体中位数的MLE的厶保序回归,而这两种保序回归的算法都是利用PAVA
思想,因此直接将PAVA与最大似然估计的求解联系起来,针对带约束的一类最
大似然估计的求解,给出了一个比较直接的迭代方法并给予了相应的证明。
4
剩余32页未读,继续阅读
资源评论
programxh
- 粉丝: 17
- 资源: 1万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- Python毕业设计基于Django的网易云数据分析可视化大屏系统的设计与实现+使用说明+全部资料(优秀项目).zip
- Python毕业设计基于Django的网易云数据分析可视化大屏系统的设计与实现+使用说明+全部资料(高分项目).zip
- JavaScript实现的鼠标手势
- 人工兔优化算法ARO MATLAB源码, 应用案例为函数极值求解以及优化svm进行分类,代码注释详细,可结合自身需求进行应用
- 人才网站设计-asp.net+sql-(系统源码)
- asp.net+sql人才网站设计-含系统源码
- C#应用的用户配置窗体方案
- python实现绘制爱心图形的代码
- JAVAWEB项目-校园订餐系统项目源码.zip
- flink-1.19.0-bin-scala-2.12.tgz flink-1.16.3-bin-scala-2.12.tgz
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功