2023年回归大作业基于多元线性回归的期权价格预测模型.doc

基于多元线性回归的期权价格预测模型王某某（北京航空航天大学 计算机学院 北京 100191）作者简介：王某某，北京航空航天大学研究生 邮箱：摘 要：期权是国际市场成熟、普遍的金融衍生品，是金融市场极为重要的金融工具2023年2月9日，上海证券交易所正式推出了我国首支场内交易期权——上证50ETF期权，翻开了境内场内期权市场的新篇章50ETF期权上市以来，市场规模逐步扩大，其发展情况境外期权产品相同时期本文以此为研究背景，以“50ETF购12月1.95”这支期权为研究对象，以今日开盘价、收盘价、最高价、最低价、结算价、成交量、成交额、持仓量、涨停价和跌停价为解释变量，通过多元线性回归模型，预测该期权的明日收盘价本次研究以多元线性回归的全模型（模型1）为出发点，通过异方差检查、残差的独立性检查、误差的正太分布检查以及多重共线性检查，说明该模型不违反回归的基本假设条件进而通过主成分回归（模型4）和逐步回归（模型5）进行降维，结果表白因变量与解释变量之间存在强烈的线性相关关系，且主成分回归和逐步回归相比全模型有更好的预测能力关键词：期权价格 多元线性回归 50ETF 多重共线性 因子分析一、引言期权（option）是依据合约形态划分的一种衍生品，指赋予其购买方在规定期限内按买卖双方约定的价格（即协议价格或行权价格）购买或者出售一定数量某种金融资产（即标的资产）的权利的合约。

期权购买方为了获得这个权利，必须支付给期权出售方一定的费用，称为权利金或期权价格Error! Reference source not found.2023年2月9日，上海证券交易所正式推出了我国首支场内交易期权——上证50ETF，翻开了境内场内期权市场的新篇章期权是与期货并列的基础衍生产品，是金融市场极为重要的金融工具之一自50ETF上市以来，市场规模逐步扩大2023年2月日均合约成交面值为5.45亿元，12月就达成了47.69亿元，增长了7.75倍；2月日均合约成交量为2.33万张，12月就达成了19.81万张，增长了7.5倍；2月权利金总成交额为2.48亿元，12月就达成了35.98亿元，增长了13.51倍Error! Reference source not found.我国股票市场有上亿的个人投资者，是一个较为典型的散户市场Error! Reference source not found.相较于专业投资机构讲，散户缺少时间，精力以及专业分析，投资具有很大的投机行为对于这些投资者来说，期权价格的变动则是他们最为关注的问题，其变化直接影响到自身的收益在实际情况中，影响股票价格的因素很多，涉及到金融政策、利率政策以及国际市场等因素，其作用机制也相称复杂Error! Reference source not found.。

因此，对于期权价格预测的研究，则可以减少投资者的投资风险，及时调整投资结构，从而保障自身的收益本文选择“50ETF购12月1.95（期权代码：10000629）”这支期权作为研究对象，根据过去一个月内期权的交易数据，以今日开盘价、收盘价、最高价、最低价、结算价、成交量、成交额、持仓量、涨停价和跌停价为解释变量，通过多元线性回归模型，预测该期权的明日收盘价下文由如下几部分构成：第二部分介绍了本次研究的数据集，涉及数据来源、和数据字段；第三部分重点介绍了各个多元线性回归模型，涉及全模型及异方差检查，残差的独立性检查、误差的正太分布检查和多种共线性检查，在第4小节和第5小节分别采用主成分回归和逐步回归对模型加以改善；第四部分运用第三部分建立的各个模型对期权价格进行了预测；第五部分对本文研究进行了总结并未来的研究加以展望二、数据说明本次研究的数据来源于Wind资讯金融终端，从上面获取了“50ETF购12月1.95”这支期权自2023年10月24日至2023年11月24日（只包含工作日）共计24日的交易数据通过整理后得到最终的数据字段，见表1表 1 期权交易数据字段收盘价开盘价最高价最低价结算价成交额成交量持仓量涨停价跌停价 期权交易数据见附录1。

三、建模1 符号说明 各个变量及其符号说明见表2表 2 各个变量及其符号说明变量符号明日收盘价今日开盘价今日收盘价今日最高价今日最低价今日结算价今日成交额今日成交量今日持仓量今日涨停价今日跌停价2 解释变量与指标变量的散点图在建立模型之前，一方面运用MATLAB绘制各个解释变量与指标变量（明日收盘价）之间的散点图，观测各个解释变量与指标变量之间的关系，散点图结果见图1图 1 各个解释变量与指标变量（明日收盘价）的散点图通过图一中的散点图可以看出，明日收盘价与今日收盘价、今日开盘价、今日最高价、今日最低价、今日结算价、今日持仓量、今日涨停价以及今日跌停价之间有较为明显的线性关系；而与今日成交额以及今日成交量之间的线性关系并不明显3 全模型 通过散点图我们观测到指标变量与各个解释变量之间大体上为线性关系，所以考虑多元线性回归模型进行建模3.1 模型建立与求解一方面建立各个解释变量与指标变量的全模型，即考虑如下模型：式中，，，…，，都是与解释变量无关的未知参数，其中，，…，称为回归系数在MATLAB中使用regress函数即可求解此多元线性回归模型，求解结果见表3.表 3 MATLAB求解全模型结果回归系数回归系数的估计值回归系数置信区间0.1555-0.1938 0.50491.8288-0.9212 4.57870.0586-1.5671 1.68430.6199-0.4249 1.6647-1.0532-1.9703 -0.1362-1.2208-3.8759 1.43420.0000-0.0000 0.0000-0.0000-0.0001 0.0000-0.0000-0.0000 0.00000.3233-1.0317 1.678300 0 = 0.9489 = 26.7989 < 0.0001 = 0.0001从上表容易看出模型（1）成立，此时指标变量与解释变量满足的回归模型为： （模型1）在MATLAB中画出模型1的残差分布，见图2。

图 2 模型1的残差分布除第12个数据外其余残差的置信区间均包含零点，因此第12个点应当视为异常点，将其剔除后重新计算，得到模型2： （模型2）模型2的 = 0.9736， = 49.1655， < 0.0001， = 0.0001但同时我们发现，各个回归系数的置信区间仍然均包含零点画出模型2的残差分布，发现又出现了异常点，按照上面的方法剔除异常点直至没有异常点产生最终一共去掉了5个异常点（分别是2，12，15，22以及23），最终得到模型3： （模型3）模型3的 = 0.9969， = 281.6079， < 0.0001， < 0.0001可见模型2相称显著，各个回归系数的置信区间均不包含零点；但是去除的数据点有5个（总共23个），去除的异常点比较多，容易去除有用信息，所以我们只把模型3作为后期预测的一个对比模型，在之后的分析中，我们只考虑去除第12个异常点的数据接下来进行异方差检查，残差的独立性检查以及多重共线性检查。

3.2 异方差检查图 3 残差对预测值的散点图通过SPSS画出残差对预测值的散点图，见图3根据图3可以认为数据基本不存在异方差性，即认为随机误差具有相同的方差3.3 残差的独立性检查通过SPSS的Durbin-watson检查得到DW值为2.545，模型2中解释变量数量k（涉及常数项）为7，样本数量n 为22，查DW分布表得到相应的0.510， 2.015由于，认为随机误差之间无自相关Error! Reference source not found.3.4 误差的正态分布检查通过SPSS画出正态P-P图以及回归标准化残差频率分布直方图，见图4图 4 误差的正太分布检查图从上图可以认为误差服从均值为0的正态分布3.5 多重共线性检查多元线性回归模型的基本假设中规定设计矩阵X的列向量之间不存在密切的线性关系Error! Reference source not found.，假如存在，就称它们之间存在多重共线性表 4 多重共线性诊断结果通过SPSS进行多重共线性诊断，结果如下表从上表可以看出绝大部分变量的VIF值远远大于10，说明严重存在多重共线性，针对多重共线性的问题，接下来用主成分分析加以解决。

4 主成分回归运用MATLAB对这个十个解释变量进行主成分分析，分析过程如下1）对原始数据进行标准化解决即每个元素的值减去该列的均值，然后除以这一列的方差2）计算相关系数矩阵R，结果见表表 5 解释变量的相关系数1.000.920.980.961.000.610.50-0.880.920.920.921.000.900.960.910.430.31-0.880.990.990.980.901.000.940.980.680.58-0.880.910.910.960.960.941.000.950.460.34-0.840.960.961.000.910.980.951.000.640.54-0.880.910.910.610.430.680.460.641.000.99-0.540.460.460.500.310.580.340.540.991.00-0.450.340.34-0.88-0.88-0.88-0.84-0.88-0.54-0.451.00-0.86-0.860.920.990.910.960.910.460.34-0.861.001.000.920.990.910.960.910.460.34-0.861.001.003）计算特性值和累计奉献率见表6表 6 各个解释变量的特性根和累计奉献率变量特性根累计奉献率8.103081.03%1.501196.04%0.187897.92%0.155399.47%0.024099.71%0.017399.89%0.007499.96%0.002899.99%0.0013100.0%0.00001。