多元线性回归演示文稿_第1页
多元线性回归演示文稿_第2页
多元线性回归演示文稿_第3页
多元线性回归演示文稿_第4页
多元线性回归演示文稿_第5页
已阅读5页,还剩67页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

多元线性回归演示文稿现在是1页\一共有72页\编辑于星期五多元线性回归现在是2页\一共有72页\编辑于星期五第12章多元线性回归§12.1多元线性回归模型§12.2回归方程的拟合优度§12.3显著性检验§12.4多重共线性§12.5利用回归方程进行估计和预测(删去)§12.6变量选择与逐步回归(删去)§12.7虚拟自变量的回归现在是3页\一共有72页\编辑于星期五§12.1多元线性回归模型12.1.1多元回归模型与回归方程12.1.2估计的多元回归方程12.1.3参数的最小二乘估计现在是4页\一共有72页\编辑于星期五12.1多元线性回归模型一个因变量与两个及两个以上自变量的回归问题就是多元回归。12.1.1多元回归模型与回归方程设因变量y,k个自变量分别为x1,x2,…,xk,描述因变量y如何依赖自变量x1,x2,…,xk和误差项

的方程,称为多元回归模型(multipleregressionmodel)。多元回归模型一般形式为:其中,b0

,b1,b2

,,bk是参数

是被称为误差项的随机变量y是x1,,x2

,,xk

的线性函数加上误差项

包含在y里面但不能被k个自变量的线性关系所解释的变异性现在是5页\一共有72页\编辑于星期五12.1.1多元回归模型与回归方程(1).误差项ε是一个期望值为0的随机变量,即E()=0。即:(2).对于自变量x1,x2,…,xk的所有值,的方差2都相同(3).误差项ε是一个服从正态分布的随机变量,即ε~N(0,2),且相互独立。独立性意味着对于自变量x1,x2,…,xk的一组特定值所对应的ε与x1,x2,…,xk任意一组其他值所对应的ε不相关。正态性意味着对于给定的x1,x2,…,xk的值,因变量y也是一个服从正态分布的随机变量。现在是6页\一共有72页\编辑于星期五12.1.1多元回归模型与回归方程根据回归模型的假定有E(y)=0+1x1+2x2+…+

kxk,上式称为多元回归方程(multipleregressionequation),它描述了因变量y的期望值与自变量x1,x2,...,xk之间的关系。现在是7页\一共有72页\编辑于星期五12.1.1多元回归模型与回归方程二元线性回归模型(观察到的y)回归面0ix1yx2(x1,x2)}现在是8页\一共有72页\编辑于星期五12.1.2估计的多元回归方程现在是9页\一共有72页\编辑于星期五12.1.3参数的最小二乘估计2.求解各回归参数的标准方程如下1.使因变量的观察值与估计值之间的离差平方和达到最小来求得

。即现在是10页\一共有72页\编辑于星期五12.1.3参数的最小二乘估计【例12.1】继续沿用第11章中例11.6。一家大型商业银行在多个地区设有分行,其业务主要是进行基础设施建设、国家重点项目建设、固定资产投资等项目的贷款。近年来,该银行的贷款额平稳增长,但不良贷款额也有较大比例的提高,这给银行业务的发展带来较大压力。为弄清楚不良贷款形成的原因,抽取了该银行所属的25家分行2002年的有关业务数据。试建立不良贷款(y)与贷款余额(x1)、累计应收贷款(x2)、贷款项目个数(x3)和固定资产投资额(x4)的线性回归方程,并解释各回归系数的含义

用Excel进行回归

现在是11页\一共有72页\编辑于星期五12.1.3参数的最小二乘估计现在是12页\一共有72页\编辑于星期五§12.2回归方程的拟合优度12.2.1多重判定系数12.2.2估计标准误差现在是13页\一共有72页\编辑于星期五12.2多重判定系数多元回归中因变量离差平方和的分解:SST=SSR+SSE多重判定系数(multiplecoefficientofdetermination)是多元回归中的回归平方和占总平方和的比例,它是度量多元回归方程拟合程度的一个统计量,反映了在因变量y的变差中被估计的回归方程所解释的比例。计算公式为现在是14页\一共有72页\编辑于星期五12.2多重判定系数注:由于自变量个数的增加,将影响到因变量中被估计回归方程中所解释的变差数量。当增加自变量时,会使预测误差变得比较小,从而减少残差平方和SSE,由于回归平方和SSR=SST-SSE,当SSE变小时,SSR会变大,从而R2也会变大。如果模型中增加一个自变量,即使这个自变量在统计上并不显著,R2也会变大,为避免这种情况,提出调整的多重判定系数(adjustedmultiplecoefficientofdetermination)计算公式为现在是15页\一共有72页\编辑于星期五12.2多重判定系数调整的多重判定系数的解释与R2类似,不同的是:(1).同时考虑了样本量和模型中的自变量的个数的影响,这就使得的值永远小于R2,而且的值不会由于模型中自变量个数的增加而越来越接近1。因此,在多元回归分析中,通常用调整的多重判定系数。(2).R2的平方根称为多重相关系数,也称为复相关系数,它度量了因变量同k个自变量的相关程度。现在是16页\一共有72页\编辑于星期五12.2.2估计标准误差多元回归分析中的估计标准误差也是对误差项的标准差的一个估计值,它是衡量多元回归方程的拟合优度方面也起着重要作用。计算公式为多元回归中对se的解释:由于se所估计的是预测误差的标准差,其含义是根据自变量x1,x2,…,xk来预测因变量y时的平均预测误差。现在是17页\一共有72页\编辑于星期五§12.3显著性检验12.3.1线性关系检验12.3.2回归系数检验和推断现在是18页\一共有72页\编辑于星期五12.3.1线性关系检验1.检验因变量与所有自变量之间的关系是否显著,也被称为总体显著性检验。2.检验方法是将回归平方和(SSR)同残差平方和(SSE)加以比较,应用F检验来分析二者之间的差别是否显著。如果是显著的,因变量与自变量之间存在线性关系如果不显著,因变量与自变量之间不存在线性关系现在是19页\一共有72页\编辑于星期五12.3.1线性关系检验第1步:提出假设H0:12k=0线性关系不显著H1:1,2,,k至少有一个不等于0第2步:计算检验统计量F第3步:作出统计决策。给定显著性水平和分子自由度k、分母自由度n-k-1找出临界值F,若F>F,拒绝H0;若F<F,则不拒绝H0。也可利用P值来判断。现在是20页\一共有72页\编辑于星期五12.3.1线性关系检验【例12.2】根据例12.1建立的回归方程,对回归方程线性关系的显著性进行检验(a=0.05)解:提出假设H0:1=2=3=4=0H1:1,2,3,4至少有一个不等于0计算检验统计量F作出统计决策。给定显著性水平=0.05和分子自由度4、分母自由度25-4-1=20找出临界值F=2.87,则F>F,所以拒绝H0,这意味着不良贷款与贷款余额、累计应收贷款、贷款项目个数和固定资产投资额之间的线性关系是显著的。现在是21页\一共有72页\编辑于星期五12.3.2回归系数检验和推断1.在回归方程通过线性关系检验后,就可以对各个回归系数有选择地进行一次或多次检验。但究竟要对哪几个回归系数进行检验,通常需要在建立模型之前作出决定。对回归系数检验的个数进行限制,以避免犯过多的第一类错误(弃真错误)2.对每一个自变量都要单独进行检验3.应用

t检验统计量现在是22页\一共有72页\编辑于星期五12.3.2回归系数检验和推断4.回归系数检验的具体步骤:提出假设H0:bi

=0(自变量xi

因变量y没有线性关系)H1:bi

0(自变量xi

因变量y有线性关系)计算检验的统计量

t作出统计决策。给定显著性水平,并进行决策t>t,拒绝H0;t<t,不拒绝H0现在是23页\一共有72页\编辑于星期五12.3.2回归系数检验和推断【例12.3】根据例12.1建立的回归方程,对回归方程各系数的显著性进行检验(a=0.05)解:提出假设H0:i=0(i=1,2,3,4)H1:i≠0计算检验统计量,由excel可知,t1=3.84,t2=1.88,t3=0.17,t4=1.88作出统计决策。给定显著性水平=0.05和自由度25-4-1=20找出临界值t/2=2.1,则只有b1通过了检验。现在是24页\一共有72页\编辑于星期五12.3.2回归系数检验和推断5.回归系数在(1-)%置信水平下的置信区间为回归系数的抽样标准差现在是25页\一共有72页\编辑于星期五§12.4多重共线性12.4.1多重共线性及其所产生的问题12.4.2多重共线性的判别12.4.3多重共线性问题的处理现在是26页\一共有72页\编辑于星期五12.4多重共线性当回归模型中使用两个或两个以上的自变量时,这些自变量之间往往会提供多余的信息。也就是说,这些变量之间彼此相关。例如,在例12.1所建立的回归方程中,使用了4个变量,即贷款余额、累计应收贷款、贷款项目个数和固定资产投资额。虽然对预测不良贷款都有作用,但由于4个自变量之间本身存在相关关系,在预测中所用的信息就是重复的。从直观上看,贷款余额与累计应收贷款之间就有较高的相关关系,这两个变量所提供的信息就是重复的,或许只用一个自变量就可以了。现在是27页\一共有72页\编辑于星期五12.4.1多重共线性及其所产生的问题1.多重共线性(multicollinearity)是指在回归模型中,有两个或两个以上的自变量彼此相关。2.多重共线性带来的问题有(1).变量之间高度相关时,可能会使回归的结果造成混乱,甚至会把分析引入歧途(2).多重共线性可能对参数估计值的正负号产生影响,特别是各回归系数的正负号有可能同我们预期的正负号相反所以,当存在多重共线性时,对回归系数的解释将是危险的。现在是28页\一共有72页\编辑于星期五12.4.2多重共线性的判别1.检测多重共线性的最简单的一种办法是计算模型中各对自变量之间的相关系数,并对各相关系数进行显著性检验。如果有一个或多个相关系数是显著的,就表示模型中所用的自变量之间相关,因而存在着多重共线性问题2.如果出现下列情况,暗示存在多重共线性(1).模型中各对自变量之间显著相关。(2).当模型的线性关系检验(F检验)显著时,几乎所有回归系数的t检验却不显著。(3).回归系数的正负号与预期的相反。现在是29页\一共有72页\编辑于星期五12.4.2多重共线性的判别【例12.4】利用例12.1的数据,按上述方法判别所建立的回归方程是否存在多重共线性。贷款余额,应收贷款,贷款项目个数,固定资产投资额之间的相关矩阵现在是30页\一共有72页\编辑于星期五12.4.2多重共线性的判别计算出相关系数统计量:1.t(25-2)=2.07,所有统计量t>t(25-2)=2.07,所以均拒绝原假设,说明这4个自变量两两之间都有显著的相关关系现在是31页\一共有72页\编辑于星期五12.4.2多重共线性的判别2.由表Excel输出的结果可知,回归模型的线性关系显著(Significance-F=1.03539E-06<=0.05)。而回归系数检验时却有3个没有通过t检验(P-Value=0.075、0.86、0.067>=0.05)。这也暗示了模型中存在多重共线性3.固定资产投资额的回归系数为负号(-0.029),与预期的不一致由以上三点可以判断回归模型中存在多重共线性。现在是32页\一共有72页\编辑于星期五12.4.3多重共线性问题的处理多重共线性问题的解决办法:(1).将一个或多个相关的自变量从模型中剔除,使保留的自变量尽可能不相关(2).如果要在模型中保留所有的自变量,则应该:避免根据t统计量对单个参数进行检验;对因变量y值的推断(估计或预测)限定在自变量样本值的范围内。现在是33页\一共有72页\编辑于星期五12.4.3多重共线性问题的处理【例12.5】利用例12.1所建立的回归方程,对多重共线性问题进行处理。解:首先,考虑将一些相关的自变量从模型中剔除。从前表可以看出,贷款余额与贷款项目个数的相关系数最高,而且从定性角度看,贷款余额与应收贷款之间也有很强的相关关系。因此将贷款项目个数和累积应收贷款这两个自变量剔除,建立不良贷款与贷款余额和固定资产投资额的线性模型。从结果可以看出,线性关系和各回归系数在0.05的显著水平下是显著的,多重共线性问题不存在了。现在是34页\一共有72页\编辑于星期五12.4.3多重共线性问题的处理多重共线性问题带来的主要麻烦是对单个回归系数的解释和检验。在求因变量的置信区间和预测区间时一般不会受其影响,但必须保证用于估计或预测的自变量的值是在样本数据的范围之内。因此,如果仅仅是为了估计或预测,可以将所有自变量都保留在模型中。在建立多元线性回归模型时,不要试图引入更多的自变量,除非有必要。特别是在社会科学的研究中,由于所使用的大多数数据都是非试验性质的,因此,在某些情况下,得到的结果往往并不令人满意,但这不一定是选择的模型不适合,而是数据的质量不好,或者是由于引入的自变量不合适。现在是35页\一共有72页\编辑于星期五§12.7虚拟自变量的回归12.7.1在模型中引进虚拟变量12.7.2含有一个虚拟自变量的回归现在是36页\一共有72页\编辑于星期五12.7.1在模型中引进虚拟变量1.虚拟变量(dummyvariable)是指用数字代码表示的定性自变量2.当模型中使用虚拟自变量时,称为虚拟自变量的回归。3.虚拟变量的取值只能是0,14.回归分析中引入虚拟自变量的方法:(1).如果定性自变量只有两个水平时。比如,性别(男,女)现在是37页\一共有72页\编辑于星期五12.7.1在模型中引进虚拟变量(2).有两个以上水平的虚拟自变量。比如,贷款企业的类型(家电,医药,其他)一般而言,如果定性自变量有k个水平,需要引入k-1个虚拟变量。5.对于含有一个虚拟自变量的回归,采用下列形式的回归方程:现在是38页\一共有72页\编辑于星期五12.7.1在模型中引进虚拟变量【例12.8】为研究考试成绩与性别之间的关系,从某大学商学院随机抽取男女学生各8名,得到他们的市场营销学课程的考试成绩如下表对性别引入虚拟变量现在是39页\一共有72页\编辑于星期五12.7.2含有一个虚拟自变量的回归1.模型中只含有一个虚拟变量的回归:建立回归模型为:y=0+1x+e回归方程可写:E(y)=0+1x男(x=0):E(y)=0—男学生考试成绩的期望值女(x=1):E(y)=0+1—女学生考试成绩的期望值注意:当指定虚拟变量0—1时0总是代表与虚拟变量值0所对应的那个分类变量水平的平均值1总是代表与虚拟变量值1所对应的那个分类变量水平的平均响应与虚拟变量值0所对应的那个分类变量水平的平均值的差值,即平均值的差值=(0+1)-0=1现在是40页\一共有72页\编辑于星期五12.7.2含有一个虚拟自变量的回归解:散点图男女【例12.9】沿用例12.8。试建立考试成绩与性别之间的线性回归模型,并解释回归系数的含义。现在是41页\一共有72页\编辑于星期五12.7.2含有一个虚拟自变量的回归【例12.10】为研究工资水平与工作年限和性别之间的关系,在某行业中随机抽取10名职工,所得数据如下表现在是42页\一共有72页\编辑于星期五12.7.2含有一个虚拟自变量的回归2.模型中含有一个数值型自变量和一个虚拟变量回归引进虚拟变量时,回归方程可写:

E(y)=0+1x1+2x2女(

x2=0):E(y|女性)=0+1x1男(x2=1):E(y|男性)=(0+2)+1x10的含义表示:女性职工的期望月工资收入(0+2)的含义表示:男性职工的期望月工资收入1含义表示:工作年限每增加1年,男性或女性工资的平均增加值2含义表示:男性职工的期望月工资收入与女性职工的期望月工资收入之间的差值(0+2)-0=2现在是43页\一共有72页\编辑于星期五习题选讲

【习题12.01】根据下面的数据用Excel进行回归,并对回归结果进行讨论,计算x1=200、x2=7时y的预测值。SUMMARYOUTPUT回归统计MultipleR0.459234RSquare0.210896AdjustedRSquare-0.01456标准误差13.34122观测值10现在是44页\一共有72页\编辑于星期五习题选讲

方差分析

dfSSMSFSignificanceF回归分析2332.9837166.49190.935410.436485残差71245.916177.988总计91578.9

Coefficients标准误差tStatP-valueIntercept25.028722.278631.123440.298298XVariable1-0.049710.105992-0.469040.653301XVariable21.9281691.472161.3097550.231624现在是45页\一共有72页\编辑于星期五习题选讲

【习题12.02】根据下面Excel输出的回归结果,说明模型中涉及多少个自变量?多少个观察值?写出回归方程,并根据F、se、R2及修正的的值对模型进行讨论。

现在是46页\一共有72页\编辑于星期五习题选讲

【习题12.03】根据两个自变量得到的多元回归方程为现在是47页\一共有72页\编辑于星期五习题选讲

【习题12.04】一家电器销售公司的管理人员认为,每月的销售额是广告费用的函数,并想通过广告费用对月销售额作出估计。下面是近8个月的销售额与广告费用数据。要求:(1).用电视广告费用作自变量,月销售额作因变量,建立估计的回归方程。(2).用电视广告费用和报纸广告费用作自变量,月销售额作因变量,建立估计的回归方程。(3).上述(1)和(2)所建立的估计方程,电视广告费用的系数是否相同?对其回归系数分别进行解释。(4).根据问题(2)所建立的估计方程,在销售收入的总变差中,被估计的回归方程所解释的比例是多少?(5).根据问题(2)所建立的估计方程,检验回归系数是否显著(a=0.05)?现在是48页\一共有72页\编辑于星期五习题选讲

现在是49页\一共有72页\编辑于星期五习题选讲

【习题12.05】某农场通过试验取得早稻收获量与春季降雨量和春季温度的数据如下。要求:(1).试确定早稻收获量对春季降雨量和春季温度的二元线性回归方程。(2).解释回归系数的实际意义。(3).根据你的判断,模型中是否存在多重共线性?

现在是50页\一共有72页\编辑于星期五习题选讲

【习题12.12】为分析某行业中的薪水有无性别歧视,从该行业中随机抽取15名员工,有关的数据如下。要求:用Excel进行回归,并对结果进行分析。现在是51页\一共有72页\编辑于星期五习题选讲

【习题12.06】一家房地产评估公司想对某城市的房地产销售价格y与地产的评估价值x1、房产的评估价值x2和使用面积x3建立一个模型,以便对销售价格作出合理预测。为此,收集了20栋住宅的房地产评估数据如下。用Excel进行回归,回答下面的问题:(1).写出估计的多元回归方程。(2).在销售价格的总变差中,被估计的回归方程所解释的比例是多少?(3).检验回归方程的线性关系是否显著(a=0.05)。(4).检验各回归系数是否显著

(a=0.05)?现在是52页\一共有72页\编辑于星期五习题选讲

【习题12.07】根据11.4题中的数据,回答下面的问题:(1).a=0.01的水平下,检验二元回归模型线性关系的显著性。(2).在a=0.05的水平下,检验回归系数b1的显著性,你认为x1应该从模型中剔除吗?(3).在a=0.05的水平下,检验回归系数b2的显著性,你认为x2应该从模型中剔除吗?现在是53页\一共有72页\编辑于星期五习题选讲

【习题12.08】根据下面的数据回答下面的问题:(1).计算y与x1之间的相关系数,有无证据表明二者之间存在线性关系(a=0.05)?(2).计算y与x2之间的相关系数,有无证据表明二者之间存在线性关系(a=0.05)?(3).根据上面的结论,你认为E(y)=b0+b1x1+b2x2对预测y是否有用?(4).用Excel进行回归,并对模型进行检验,所得的结论与(3)是否相同?(5)计算x1与x2之间的相关系数,所得结果意味着什么?现在是54页\一共有72页\编辑于星期五习题选讲

【习题12.09】下面是随机抽取的15家大型商场销售的同类产品的有关数据(单位:元)。要求:(1).计算y与x1、y与x2之间的相关系数,是否有证据表明销售价格与购进价格、销售价格与销售费用之间存在线性关系?(2).根据上述结果,你认为用购进价格和销售费用来预测销售价格是否有用?(3).用Excel进行回归,并检验模型的线性关系是否显著(a=0.05)。(4).解释判定系数R2,所得结论与问题(2)中是否一致?(5).计算x1与x2之间的相关系数,所得结果意味着什么?(6).模型中是否存在多重共线性?你对模型有何建议?现在是55页\一共有72页\编辑于星期五习题选讲

【习题12.10】设因变量为y,一个数值型自变量x1和一个具有两个水平(水平1和水平2)的分类型自变量。要求:(1).写出因变量y关于自变量x1和分类自变量的多元回归方程。(2).对应于分类自变量水平1的y的期望值是多少?(3).对应于分类自变量水平2的y的期望值是多少?现在是56页\一共有72页\编辑于星期五习题选讲

【习题12.11】一家货物运输公司想研究运输费用与货物类型的关系,并建立运输费用与货物类型的回归模型,以此对运输费用作出预测。该运输公司所运输的货物分为两种类型:易碎品和非易碎品。下面给出了15个路程大致相同、而货物类型不同的运输费用数据。要求:(1).写出运输费用与货物类型之间的线性方程。(2).对模型中的回归系数进行解释。(3).检验模型的线性关系是否显著(a=0.05)?现在是57页\一共有72页\编辑于星期五§12.5利用回归方程进行估计和预测删去不讲现在是58页\一共有72页\编辑于星期五§12.5利用回归方程进行估计和预测STATISTICA输出的不良贷款的置信区间【例12.6】根据例12.1的数据,取x1=100、x2=10、x3=15、x4=60,建立不良贷款的95%的置信区间和预测区间。现在是59页\一共有72页\编辑于星期五§12.5利用回归方程进行估计和预测STATISTICA输出的不良贷款的预测区间现在是60页\一共有72页\编辑于星期五§12.5利用回归方程进行估计和预测STATISTICA输出的不良贷款的置信区间【例12.6】根据例12.1的数据,取x1=100、x2=10、x3=15、x4=60,建立不良贷款的95%的置信区间和预测区间。现在是61页\一共有72页\编辑于星期五§12.5利用回归方程进行估计和预测STATISTICA输出的不良贷款的预测区间现在是62页\一共有72页\编辑于星期五§12.6变量选择与逐步回归12.6.1变量选择过程12.6.2向前选择12.6.3向后剔除12.6.4逐步回归现在是63页\一共有72页\编辑于星期五变量选择过程1.在建立回归模型时,对自变量进行筛选2.选择自变量的原则是对统计量进行显著性检验(1).将一个或一个以上的自变量引入到回归模型中时,是否使得残差平方和(SSE)有显著地减少。如果增加一个自变量使SSE的减少是显著的,则说明有必要将这个自变量引入回归模型,否则,就没有必要将这个自变量引入回归模型(2).确定引入自变量是否使SSE有显著减少的方法,就是使用F统计量的值作为一个标准,以此来确定是在模型中增加一个自变量,还是从模型中剔除一个自变量3.变量选择的方法主要有:向前选择、向后剔除、逐步回归、最优子集等现在是64页\一共有72页\编辑于星期五§12.6.2向前选择

1.从模型中没有自变量开始2.对k个自变量分别拟合对因变量y的一元线性回归模型,共有k个,然后找出F统计量的值最高的模型及其自变量(P值最小的),并将其首先引入模型3.分别拟合引入模型外的k-1个自变量的线性回归模型,如此反复进行,直至模型外的自变量均无统计显著性为止现在是65页\一共有72页\编辑于星期五§12.6.3向后剔除

1.先对因变量拟合包括所有k个自变量的回归模型。然后考察p(p<k)个去掉一个自变量的模型(这

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论