版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、协整检验与误差修正模型一、长期均衡与协整分析一、长期均衡与协整分析Equilibrium and Cointegration1、问题的提出、问题的提出 经典回归模型(经典回归模型(classical regression modelclassical regression model)是建立在平稳)是建立在平稳数据变量基础上的,对于非平稳变量,不能使用经典回归模型,否则数据变量基础上的,对于非平稳变量,不能使用经典回归模型,否则会出现虚假回归等诸多问题。会出现虚假回归等诸多问题。由于许多经济变量是非平稳的,这就给经典的回归分析方法带来了很由于许多经济变量是非平稳的,这就给经典的回归分析方法带来
2、了很大限制。大限制。 但是,如果变量之间有着长期的稳定关系,即它们之间是协但是,如果变量之间有着长期的稳定关系,即它们之间是协整的(整的(cointegration)cointegration),则是可以使用经典回归模型方法建立,则是可以使用经典回归模型方法建立回归模型的。回归模型的。 例如,中国居民人均消费水平与人均例如,中国居民人均消费水平与人均GDPGDP变量的例子变量的例子, , 从经济从经济理论上说,人均理论上说,人均GDPGDP决定着居民人均消费水平,它们之间有着长决定着居民人均消费水平,它们之间有着长期的稳定关系,即它们之间是协整的。期的稳定关系,即它们之间是协整的。 经济理论指
3、出,某些经济变量间确实存在着长期均衡关系,经济理论指出,某些经济变量间确实存在着长期均衡关系,这种均衡关系意味着经济系统不存在破坏均衡的内在机制,这种均衡关系意味着经济系统不存在破坏均衡的内在机制,如果变量在某时期受到干扰后偏离其长期均衡点,则均衡机如果变量在某时期受到干扰后偏离其长期均衡点,则均衡机制将会在下一期进行调整以使其重新回到均衡状态。制将会在下一期进行调整以使其重新回到均衡状态。 假设X与Y间的长期“均衡关系”由式描述 2 2、长期均衡、长期均衡tttXY10该均衡关系意味着该均衡关系意味着: :给定X的一个值,Y相应的均衡值也随之确定为0+1X。 在在t-1期末,存在下述三种情形
4、之一:期末,存在下述三种情形之一: Y等于它的均衡值:等于它的均衡值:Yt-1= 0 0+ + 1 1Xt ; Y小于它的均衡值:小于它的均衡值:Yt-1 0 0+ + 1 1Xt ; 在时期在时期t,假设,假设X有一个变化量有一个变化量 Xt,如果变量,如果变量X与与Y在时期在时期t与与t-1末期仍满足它们间的长期均衡关系,即上末期仍满足它们间的长期均衡关系,即上述第一种情况,则述第一种情况,则Y的相应变化量为的相应变化量为: :tttvXY1vt=t-t-1 如果如果t-1期末,发生了上述第二种情况,即期末,发生了上述第二种情况,即Y的值的值小于其均衡值,则小于其均衡值,则t期末期末Y的变
5、化往往会比第一种的变化往往会比第一种情形下情形下Y的变化大一些;的变化大一些; 反之,如果反之,如果t-1期末期末Y的值大于其均衡值,则的值大于其均衡值,则t期末期末Y的变化往往会小于第一种情形下的的变化往往会小于第一种情形下的 Yt 。 可见,如果可见,如果Yt= = 0 0+ + 1 1Xt+ + t t正确地提示了正确地提示了X与与Y间间的长期稳定的的长期稳定的“均衡关系均衡关系”,则意味着,则意味着Y对其均衡对其均衡点的偏离从本质上说是点的偏离从本质上说是“临时性临时性”的。的。 一个重要的假设就是一个重要的假设就是: :随机扰动项随机扰动项 t t必须是平稳序必须是平稳序列。列。如果
6、如果 t t有随机性趋势(上升或下降),则会导致有随机性趋势(上升或下降),则会导致Y对其均衡点的任何偏离都会被长期累积下来而不对其均衡点的任何偏离都会被长期累积下来而不能被消除。能被消除。 式Yt= =0+1Xt+t中的随机扰动项也被称为非均衡误非均衡误差(差(disequilibrium error),它是变量X与Y的一个线性组合:tttXY10 如果如果X与与Y间的长期均衡关系正确,该式表述的非均间的长期均衡关系正确,该式表述的非均衡误差应是一平稳时间序列,并且具有零期望值,即衡误差应是一平稳时间序列,并且具有零期望值,即是具有是具有0均值的均值的I(0)序列。序列。 非稳定的时间序列,
7、它们的线性组合也可能成为平稳的。非稳定的时间序列,它们的线性组合也可能成为平稳的。称变量称变量X与与Y是协整的(是协整的(cointegrated)。)。3 3、协整、协整 如果序列如果序列XX1t1t,X,X2t2t, ,X,Xktkt 都是都是d d阶单整,存在向量阶单整,存在向量 =(=( 1 1, , 2 2, , , k k) ),使得,使得Z Zt t= = X XT T I(d-b)I(d-b), 其中,其中,b0b0,X=(XX=(X1t1t,X,X2t2t, ,X,Xktkt) )T T,则认为序列,则认为序列XX1t1t,X,X2t2t, ,X,Xktkt 是是(d,b)(
8、d,b)阶协整,记为阶协整,记为X XttCI(d,b)CI(d,b), 为为协整向量(协整向量(cointegrated vector)。)。 如果两个变量都是单整变量,只有当它们的单整阶数相同如果两个变量都是单整变量,只有当它们的单整阶数相同时,才可能协整;如果它们的单整阶数不相同,就不可能时,才可能协整;如果它们的单整阶数不相同,就不可能协整。协整。 3 3个以上的变量,如果具有不同的单整阶数,有可个以上的变量,如果具有不同的单整阶数,有可能经过线性组合构成低阶单整变量。能经过线性组合构成低阶单整变量。) 2(),2(),1 (IUIVIWttt)0()1 (IePcWQIbUaVPtt
9、tttt)1 ,1(,)1 ,2(,CIPWCIUVtttt (d,d)阶协整是一类非常重要的协整关系,它的经济)阶协整是一类非常重要的协整关系,它的经济意义在于:意义在于:两个变量,虽然它们具有各自的长期波两个变量,虽然它们具有各自的长期波动规律,但是如果它们是(动规律,但是如果它们是(d,dd,d)阶协整的,则它)阶协整的,则它们之间存在着一个长期稳定的比例关系。们之间存在着一个长期稳定的比例关系。 例如,中国例如,中国CPCCPC和和GDPPCGDPPC,它们各自都是,它们各自都是2 2阶单整,如果它们阶单整,如果它们是是(2,2)(2,2)阶协整,说明它们之间存在着一个长期稳定的比例关
10、阶协整,说明它们之间存在着一个长期稳定的比例关系,从计量经济学模型的意义上讲,建立如下居民人均消费系,从计量经济学模型的意义上讲,建立如下居民人均消费函数模型是合理的。函数模型是合理的。tttGDPPCCPC10 尽管两个时间序列是非平稳的,也可以用经典的尽管两个时间序列是非平稳的,也可以用经典的回归分析方法建立回归模型。回归分析方法建立回归模型。 从这里,我们已经初步认识到:从这里,我们已经初步认识到:检验变量之间的协检验变量之间的协整关系,在建立计量经济学模型中是非常重要的。整关系,在建立计量经济学模型中是非常重要的。 而且,从变量之间是否具有协整关系出发选择模而且,从变量之间是否具有协整
11、关系出发选择模型的变量,其数据基础是牢固的,其统计性质是优良型的变量,其数据基础是牢固的,其统计性质是优良的的。二、协整检验二、协整检验EG检验检验 1 1、两变量的、两变量的Engle-GrangerEngle-Granger检验检验 为了检验两变量Yt,Xt是否为协整,Engle和Granger于1987年提出两步检验法,也称为EG检验。 第一步,第一步,用OLS方法估计方程 Yt= =0+1Xt+t并计算非均衡误差,得到: tttttYYeXY10称为协整回归协整回归( (cointegrating)或静态回归静态回归( (static regression) )。 非均衡误差的单整性的
12、检验方法仍然是非均衡误差的单整性的检验方法仍然是DFDF检验或者检验或者ADFADF检验。检验。 需要注意是,这里的需要注意是,这里的DF或或ADF检验是针对协整回检验是针对协整回归计算出的误差项,而非真正的非均衡误差。归计算出的误差项,而非真正的非均衡误差。 而而OLS法采用了残差最小平方和原理,因此估计量法采用了残差最小平方和原理,因此估计量 是向下偏倚的,这样将导致拒绝零假设的机会比是向下偏倚的,这样将导致拒绝零假设的机会比实际情形大。实际情形大。 于是对于是对e et t平稳性检验的平稳性检验的DFDF与与ADFADF临界值应该比正常临界值应该比正常的的DFDF与与ADFADF临界值还
13、要小。临界值还要小。 MacKinnon(1991)通过模拟试验给出了协整检验的临界值。 例例 检验中国居民人均消费水平检验中国居民人均消费水平CPCCPC与人均国内生产总值与人均国内生产总值GDPPCGDPPC的协整关系。的协整关系。 已知CPC与GDPPC都是I(2)序列,已知它们的回归式 ttGDPPCCPC45831. 0764106.49R2=0.9981 对该式计算的残差序列作ADF检验,适当检验模型为: 31127. 249. 155. 1tttteeee (-4.47) (3.93) (3.05) LM(1)=0.00 LM(2)=0.00 t=-4.47-3.75=ADF0.
14、05,拒绝存在单位根的假设,残差项是平稳的。因此中国居民人均消费水平与人均中国居民人均消费水平与人均GDPGDP是是(2,2)(2,2)阶协整的,阶协整的,说明了该两变量间存在长期稳定的说明了该两变量间存在长期稳定的“均衡均衡”关系。关系。 2 2、多变量协整关系的检验、多变量协整关系的检验扩展的扩展的E-GE-G检验检验 多变量协整关系的检验要比双变量复杂一些,主要在于协整变量间可能存在多种稳定的线性组合协整变量间可能存在多种稳定的线性组合。 假设有4个I(1)变量Z、X、Y、W,它们有如下的长期均衡关系:tttttYXWZ3210非均衡误差项t应是I(0)序列: tttttYXWZ3210
15、 然而,如果Z与W,X与Y间分别存在长期均衡关系:tttvWZ110tttvYX210 则非均衡误差项v1t、v2t一定是稳定序列I(0)。于是它们的任意线性组合也是稳定的。例如tttttttYXWZvvv110021 由于vt象t一样,也是Z、X、Y、W四个变量的线性组合,由此vt 式也成为该四变量的另一稳定线性组合。 (1, -0,-1,-2,-3)是对应于t 式的协整向量,(1,-0-0,-1,1,-1)是对应于vt式的协整向量。 一定是I(0)序列。 检验程序:检验程序: 对于多变量的协整检验过程,基本与双变量情形相同对于多变量的协整检验过程,基本与双变量情形相同,即需检验变量是否具有
16、同阶单整性,以及是否存在稳定的线性组即需检验变量是否具有同阶单整性,以及是否存在稳定的线性组合合。 在检验是否存在稳定的线性组合时在检验是否存在稳定的线性组合时,需通过设置一个变量为被解释变量,其他变量为解释变量,进行OLS估计并检验残差序列是否平稳。 如果不平稳如果不平稳,则需更换被解释变量,进行同样的OLS估计及相应的残差项检验。 当所有的变量都被作为被解释变量检验之后,仍不当所有的变量都被作为被解释变量检验之后,仍不能得到平稳的残差项序列,则认为这些变量间不存能得到平稳的残差项序列,则认为这些变量间不存在(在(d,dd,d)阶协整。)阶协整。 检验残差项是否平稳的检验残差项是否平稳的DF
17、与与ADF检验临界值要比通常的检验临界值要比通常的DF与与ADF检验临界值小,而且该临界值还受到所检验的变量个数的影检验临界值小,而且该临界值还受到所检验的变量个数的影响。响。 MacKinnon(1991)通过模拟试验得到的不同变量协整检验的临界值。三、协整检验三、协整检验JJJJ检验检验( (教材教材6.4.3)6.4.3) JJ JJ检验的原理检验的原理 Johansen于于1988年,以及与年,以及与Juselius一起于一起于1990年年提出了一种用向量自回归模型进行检验的方法,提出了一种用向量自回归模型进行检验的方法,通常称为通常称为Johansen检验,或检验,或JJ检验,是一种
18、进行检验,是一种进行多重多重I(1)序列协整检验的较好方法。序列协整检验的较好方法。 没有移动平均项的向量自回归模型表示为:没有移动平均项的向量自回归模型表示为: tptpttyyy11tjtjpjtyy1ttjtpjjt11yyy差分Yt为为M个个I(1)过程构成的向量过程构成的向量 I(0)过程I(0)过程只有产生协整,才能保证新生误差是平稳过程 将将y的协整问题转变为讨论矩阵的协整问题转变为讨论矩阵的性质问题的性质问题ttjtpjjt11yyy 于是,将于是,将yt中的协整检验变成对矩阵中的协整检验变成对矩阵的分析问题。这就的分析问题。这就是是JJ检验的基本原理。检验的基本原理。 JJ
19、JJ检验的预备工作检验的预备工作 第一步:用第一步:用OLSOLS分别估计下式中的每一个方程,计分别估计下式中的每一个方程,计算残差,得到残差矩阵算残差,得到残差矩阵S S0 0,为一个,为一个(M(MT)T)阶矩阵。阶矩阵。 tjtpjjtyy1MTMMTT212222111211 第一步:用第一步:用OLSOLS分别估计下式中的每一个方程,计分别估计下式中的每一个方程,计算残差,得到残差矩阵算残差,得到残差矩阵S S1 1,也为一个,也为一个(M(MT)T)阶矩阵。阶矩阵。 tjtpjjtyy11 第三步:构造上述残差矩阵的积矩阵:第三步:构造上述残差矩阵的积矩阵: 00100SSRT10
20、101SSRT01110SSRT11111SSRT 第四步:计算有序特征值和特征向量。第四步:计算有序特征值和特征向量。 第五步:设定似然函数。第五步:设定似然函数。 JJ JJ检验之一检验之一特征值轨迹检验特征值轨迹检验 服从Johansen分布。被称为特征值轨迹统计量。 1, 2 , 1 , 0)1ln()(1MrTrMMrii ,一直检验下去,直到出现第一个不显著的,一直检验下去,直到出现第一个不显著的(Mr)为止,说明存在为止,说明存在r个协整向量。这个协整向量。这r个协整向量就个协整向量就是对应于最大的是对应于最大的r个特征值的经过正规化的特征向个特征值的经过正规化的特征向量。量。
21、JJ JJ检验之一检验之一最大特征值检验最大特征值检验 该统计量被称为最大特征值统计量。于是该检验被称为最大特征值检验。 )1ln() 1(rTr 由 Johansen和Juselius于1990年计算得到 Johansen分布临界值表。 JJJJ检验实例检验实例 GDP、CONSR、CONSP、INV取对数后为取对数后为I(1)序列。序列。即即lnGDP、lnCONSR、lnCONSP、lnINV。 对它们之间的协整关系进行检验。对它们之间的协整关系进行检验。两种方法的结论是一致的。两种方法的结论是一致的。如何处理高阶单整序列?如何处理高阶单整序列? 从理论上讲。从理论上讲。JJ JJ 检验
22、只适用于多个检验只适用于多个1 1阶单整序列。阶单整序列。 多个同阶高阶单整序列,差分为多个同阶高阶单整序列,差分为1 1阶后再检验,显然是可行阶后再检验,显然是可行的。但是意义发生变化。的。但是意义发生变化。 没有看到关于高阶多重协整检验的文献,难度太大。没有看到关于高阶多重协整检验的文献,难度太大。 能否先检验,然后建立均衡方程,通过对误差项的单位能否先检验,然后建立均衡方程,通过对误差项的单位根检验以判断发生何种协整?未见经典。根检验以判断发生何种协整?未见经典。如何选择截距和时间趋势项?如何选择截距和时间趋势项? 分别考虑分别考虑CE和和VAR中是否有截距和时间趋势项中是否有截距和时间
23、趋势项 作为假设作为假设 显著性检验显著性检验 重新检验重新检验 对协整关系检验结果无显著影响(检验统计量发生对协整关系检验结果无显著影响(检验统计量发生变化,但临界值同时发生变化)变化,但临界值同时发生变化)如何在多个协整关系中作出选择?如何在多个协整关系中作出选择? 一般选择对应于最大特征值的第一般选择对应于最大特征值的第1个协整关系个协整关系 从应用的目的出发选择从应用的目的出发选择四、误差修正模型四、误差修正模型Error Correction Model, ECM1 1、一般差分模型的问题、一般差分模型的问题 对于非稳定时间序列,可通过差分的方法将其化对于非稳定时间序列,可通过差分的
24、方法将其化为稳定序列,然后才可建立经典的回归分析模型。为稳定序列,然后才可建立经典的回归分析模型。tttXY10tttvXY11tttv模型只表达了模型只表达了X与与Y间的短期关系,间的短期关系,而没有揭示它们间的长期关系而没有揭示它们间的长期关系。关关于变量水平值的重要信息将被忽于变量水平值的重要信息将被忽略。略。误差项t不存在序列相关, t是一个一阶移动平均时间序列一阶移动平均时间序列,因而是序列相关的。是序列相关的。2 2、误差修正模型、误差修正模型 是一种具有特定形式的计量经济学模型,它的主是一种具有特定形式的计量经济学模型,它的主要形式是由要形式是由DavidsonDavidson、
25、 HendryHendry、SrbaSrba和和YeoYeo于于19781978年提出的,称为年提出的,称为DHSYDHSY模型。模型。tttXY10tttttYXXY11210tttttttttXYXYXXY12101111211011)1 ()1 ()(tttttXYXY)(11011由于现实经济中很少处在均衡点上,假设具有(1, 1)阶分布滞后形式 Y Y的变化决定于的变化决定于X X的变化以及前一时期的非均衡程的变化以及前一时期的非均衡程度度。 一阶误差修正模型一阶误差修正模型( (first-order error correction model) )的形式:的形式:ttttecm
26、XY11若若(t-1)(t-1)时刻时刻Y Y大于其长期均衡解大于其长期均衡解 0 0+ + 1 1X X,ecmecm为正,则为正,则(-(- ecm)ecm)为负,使得为负,使得 Y Yt t减少;减少;若若(t-1)(t-1)时刻时刻Y Y小于其长期均衡解小于其长期均衡解 0 0+ + 1 1X X ,ecmecm为负,则为负,则(-(- ecm)ecm)为正,使得为正,使得 Y Yt t增大。增大。体现了长期非均衡误差对短期变化的控制。体现了长期非均衡误差对短期变化的控制。 复杂的复杂的ECM形式,例如:形式,例如:tttttttYYXXXY2211231210tttttttXYXXY
27、Y)(110113112tttttttYZZXXY12211210tttttttZXYZXY)(12110111 误差修正模型的优点:误差修正模型的优点:如: a)一阶差分项的使用消除了变量可能存在的趋势因素,从而避免了虚假回归问题; b)一阶差分项的使用也消除模型可能存在的多重共线性问题; c)误差修正项的引入保证了变量水平值的信息没有被忽视; d)由于误差修正项本身的平稳性,使得该模型可以用经典的回归方法进行估计,尤其是模型中差分项可以使用通常的t检验与F检验来进行选取;等等。3 3、误差修正模型的建立、误差修正模型的建立 Granger 表述定理表述定理(Granger represen
28、taion theorem) Engle 与与 Granger 1987年提出年提出 如果变量如果变量X X与与Y Y是协整的,则它们间的短期非均衡是协整的,则它们间的短期非均衡关系总能由一个误差修正模型表述。关系总能由一个误差修正模型表述。tttecmXYlaggedY1),(模型中没有明确指出Y与X的滞后项数,可以是多阶滞后;由于一阶差分项是I(0)变量,因此模型中允许采用X的非滞后差分项Xt 。 建立误差修正模型建立误差修正模型,需要: 首先首先对变量进行协整分析,以发现变量之间的协整关系,即长期均衡关系,并以这种关系构成误差修正项。 然后然后建立短期模型,将误差修正项看作一个解释变量,
29、连同其它反映短期波动的解释变量一起,建立短期模型,即误差修正模型。Engle-Granger两步法两步法 第一步,第一步,进行协整回归(OLS法),检验变量间的协整关系,估计协整向量(长期均衡关系参数); 第二步第二步,若协整性存在,则以第一步求到的残差作为非均衡误差项加入到误差修正模型中,并用OLS法估计相应参数。 需要注意的是需要注意的是:在进行变量间的协整检验时,如有必要可在协整回归式中加入趋势项,这时,对残差项的稳定性检验就无须再设趋势项。 另外另外,第二步中变量差分滞后项的多少,可以残差项序列是否存在自相关性来判断,如果存在自相关,则应加入变量差分的滞后项。 用打开误差修正项括号的方
30、法直接估计误差修正模用打开误差修正项括号的方法直接估计误差修正模型。型。 一般不采用。一般不采用。 经济理论指出,居民消费支出是其实际收入的函数。 以中国国民核算中的居民消费支出经过居民消费价格指数缩减居民消费支出经过居民消费价格指数缩减得到中国居民实际消费支出时间序列(得到中国居民实际消费支出时间序列(C); 以支出法GDP对居民消费价格指数缩减近似地代表国民收对居民消费价格指数缩减近似地代表国民收入时间序列入时间序列(GDP)。 时间段为19782000(表) 例例 中国居民消费的误差修正模型 (1 1)对数据)对数据lnC与与lnGDP进行单整检验进行单整检验 容易验证lnC与lnGDP
31、是一阶单整的,它们适合的检验模型如下: (3.81)(-4.01) (2.66) (2.26) (2.54) LM(1)=0.38 LM(2)=0.67 LM(3)=2.34 LM(4)=2.46 首先,建立首先,建立lnC与与lnGDP的回归模型的回归模型(2)检验)检验lnC与与lnGDP的协整性,并建立长期均衡关系的协整性,并建立长期均衡关系 ttGDPCln923. 0047. 0ln (0.30) (57.48) R2=0.994 DW=0.744 发现有残关项有较强的一阶自相关性。考虑加入适当的滞后项,得lnC与lnGDP的分布滞后模型 11ln361. 0ln622. 0ln698. 0152. 0lnttttGDPCGDPC (1.63) (6.62) (4.92) (-2.17) R2=0.994 DW=1.92 LM(1)=0.00 LM(2)=2.31自相关性消除,因此可初步认为是lnC与lnGDP的长期稳定关系。 残差项的稳定性检验:残差项的稳定性检验:
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 医学教学中的临床技能竞赛设计
- 医学心理学与临床路径整合
- 医学影像技术术语与报告方法学解读
- 福建货架安全管理培训
- 2025年上半年护理核心度专项考核试题(含答案解析)
- 计算机机房突发事件应急专项预案
- 《有机化学》-第3章
- 就业指导与心理调适策略
- 小学安全教育演讲稿
- 预制构件厂培训室管理细则
- (市质检二检)福州市2024-2025学年高三年级第二次质量检测 历史试卷(含答案)
- OptiStruct结构分析与工程应用
- 2025中考数学复习专题:八类最值问题汇-总(瓜豆隐圆胡不归阿氏圆将军饮马逆等线费马点构造二次函数求最值)(原卷版)
- 柴油发电机施工方案
- 药物临床试验质量管理规范(GCP)
- 交通运输驾驶员安全承诺书
- 《建筑工程设计文件编制深度规定》(2022年版)
- 2024NEA水性气硅涂膏隔热保温墙体构造
- 物流外包与供应链管理课件
- 《热力发电厂》热力发电厂全面性热力系统
- 温病学--温病学课件
评论
0/150
提交评论