《基于回归森林模型的蛋白质序列二硫键连接模式的预测方法.pdf》由会员分享，可在线阅读，更多相关《基于回归森林模型的蛋白质序列二硫键连接模式的预测方法.pdf（13页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN104063632A43申请公布日20140924CN104063632A21申请号201410303084622申请日20140627G06F19/1620110171申请人南京理工大学地址210000江苏省南京市孝陵卫200号72发明人李阳於东军胡俊沈红斌杨静宇74专利代理机构南京理工大学专利中心32203代理人朱显国54发明名称基于回归森林模型的蛋白质序列二硫键连接模式的预测方法57摘要本发明公开一种基于回归森林模型的蛋白质序列二硫键连接模式的预测方法，包括以下步骤步骤1、通过多视角特征提取和特征组合，得到蛋白质序列中每个半胱氨酸残基对的特征向量；步骤2、对于待预测的蛋。
2、白质序列信息和训练数据集合，生成待预测的蛋白质序列信息和训练数据集合中所有半胱氨酸残基对的特征向量，从而分别构成训练样本集合与待预测样本集合；步骤3、采用随机森林算法学习二硫键样本在特征空间中的分布规律，生成随机森林回归模型；步骤4、利用随机森林回归模型对待预测样本集合的特征向量进行预测，得到每个半胱氨酸残基对形成二硫键的倾向值，得分最高的二硫键连接模式即为最终预测得到的蛋白质序列中的二硫键连接模式。51INTCL权利要求书2页说明书9页附图1页19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书9页附图1页10申请公布号CN104063632ACN104063632A1/2页。
3、21一种基于回归森林模型的蛋白质序列二硫键连接模式的预测方法，其特征在于，包括以下步骤步骤1、特征提取，即基于输入的蛋白质序列信息，进行多视角特征提取和特征组合，得到蛋白质序列中每个半胱氨酸残基对最终的特征向量；步骤2、对于待预测的蛋白质序列信息和训练数据集合，采用上述步骤1的方法，生成待预测的蛋白质序列信息和训练数据集合中所有半胱氨酸残基对的特征向量，从而分别构成训练样本集合与待预测样本集合；步骤3、特征学习采用随机森林算法学习二硫键样本在特征空间中的分布规律，生成随机森林回归模型，其中，模型的决策树个数采用大于等于300的整数，每次分裂时随机选择的特征个数由前述步骤1所得的特征向量的维度决。
4、定；步骤4、二硫键连接模式预测利用前述步骤3的随机森林回归模型对待预测样本集合中的特征向量进行预测，得到每个半胱氨酸残基对形成二硫键的倾向值，记为PPCP；令P为所有可能的二硫键连接模式的数量，则第I个二硫键连接模式的得分SI用下述公式得到其中B代表第I个二硫键连接模式中二硫键的个数，PPCPJ是第I个二硫键连接模式中第J个半胱氨酸残基对的倾向值如果第I个二硫键的连接模式的得分最高，则该第I个模式就是最终预测得到的结果即该得分最高的二硫键连接模式I为最终预测得到的蛋白质序列中的二硫键连接模式。2根据权利要求1所述的基于回归森林模型的蛋白质序列二硫键连接模式的预测方法，其特征在于，所述步骤1中特。
5、征提取的具体实现包括以下步骤分别使用PSIBLAST算法抽取蛋白质序列的进化信息，利用PSIPRED算法抽取蛋白质序列的二级结构信息，再利用滑动窗口与特征串行组合方式从前述进化信息和二级结构信息中提取每个半胱氨酸残基对的多视角特征；然后计算蛋白质序列中半胱氨酸残基对在蛋白质序列中的距离特征信息DOC，利用PSIBLAST算法抽取蛋白质序列中每个半胱氨酸残基对的协变异信息CM，利用MODELLER同源建模软件获取半胱氨酸残基对的三维结构距离信息PDTCR；最后将上述从进化信息和二级结构信息中提取每个半胱氨酸残基对的多视角特征与前述距离特征信息、协变异信息、三维结构距离信息进行特征组合，得到蛋白质。
6、序列中每个半胱氨酸残基对最终的特征向量。3根据权利要求1所述的基于回归森林模型的蛋白质序列二硫键连接模式的预测方法，其特征在于，前述蛋白质序列中每个半胱氨酸残基对最终的特征向量，其特征维度为权利要求书CN104063632A2/2页3601维。4根据权利要求1所述的基于回归森林模型的蛋白质序列二硫键连接模式的预测方法，其特征在于，前述随机森林回归模型的每次分裂时随机选择的特征个数为24。5根据权利要求1所述的基于回归森林模型的蛋白质序列二硫键连接模式的预测方法，其特征在于，前述随机森林回归模型的决策树个数取值为500。权利要求书CN104063632A1/9页4基于回归森林模型的蛋白质序列二硫。
7、键连接模式的预测方法技术领域0001本发明涉及生物信息学蛋白质序列中二硫键预测技术领域，具体而言涉及一种基于回归森林模型的蛋白质序列二硫键连接模式的预测方法。背景技术0002二硫键是最重要的蛋白质结构特性之一。二硫键是在蛋白质多肽链中两个半胱氨酸残基之间形成的主共价键，它们可以在肽链的链间或者链内形成。二硫键在蛋白质折叠方式以及稳定性方面有着非常重要的作用。因此，预测蛋白质里面的半胱氨酸残基组成二硫键的方式在预测蛋白质结构和功能上起着举足轻重的作用。0003目前有很多种预测二硫键的方法，比如，DISULFIND法ACERONI,APASSERINI,AVULLOETAL,“DISULFINDA。
8、DISULFIDEBONDINGSTATEANDCYSTEINECONNECTIVITYPREDICTIONSERVER,”NUCLEICACIDSRESEARCH,VOL34,NOSUPPL2,PPW177W181,2006、PAIRWISESVM法CHTSAI,BJCHEN,CHCHANETAL,“IMPROVINGDISULFIDECONNECTIVITYPREDICTIONWITHSEQUENTIALDISTANCEBETWEENOXIDIZEDCYSTEINES,”BIOINFORMATICS,VOL21,NO24,PP44164419,2005、SS_SVR法JSONG,ZYUAN。
9、,HTANETAL,“PREDICTINGDISULFIDECONNECTIVITYFROMPROTEINSEQUENCEUSINGMULTIPLESEQUENCEFEATUREVECTORSANDSECONDARYSTRUCTURE,”BIOINFORMATICS,VOL23,NO23,PP31473154,2007、FS_SVR法LZHU,JYANG,JNSONGETAL,“IMPROVINGTHEACCURACYOFPREDICTINGDISULDECONNECTIVITYBYFEATURESELECTION,”JCOMPUTCHEM,VOL31,NO7,PP147885,MAY,201。
10、0、DBCP法HHLIN,ANDLYTSENG,“DBCPAWEBSERVERFORDISULDEBONDINGCONNECTIVITYPATTERNPREDICTIONWITHOUTTHEPRIORKNOWLEDGEOFTHEBONDINGSTATEOFCYSTEINES,”NUCLEICACIDSRESEARCH,VOL38,NOSUPPL2,PPW503W507,2010、DISLOCATE法DISLOCATEMIPICOVCSAVOJARDO,PFARISELLI,MALHAMDOOSHETAL,“IMPROVINGTHEPREDICTIONOFDISULDEBONDSINEUKARY。
11、OTESWITHMACHINELEARNINGMETHODSANDPROTEINSUBCELLULARLOCALIZATION,”BIOINFORMATICS,VOL27,NO16,PP222430,AUG15,2011、DMC法CSAVOJARDO,PFARISELLI,PLMARTELLIETAL,“PREDICTIONOFDISULFIDECONNECTIVITYINPROTEINSWITHMACHINELEARNINGMETHODSANDCORRELATEDMUTATIONS,”BMCBIOINFORMATICS,VOL14,NOSUPPL1,PPS10,2013、DIANNA法FFE。
12、RR,ANDPCLOTE,“DIANNA11ANEXTENSIONOFTHEDIANNAWEBSERVERFORTERNARYCYSTEINECLASSICATION,”NUCLEICACIDSRESEARCH,VOL34,NOSUPPL2,PPW182W185,2006等等。0004前述方法大体可以分为三类二硫键成键状态预测；已知半胱氨酸是否形成二硫键的信息来预测二硫键的组成模式；对上述两者都进行预测。最近，利用机器学习的方法预测蛋白质二硫键模式受到越来越多人的关注，从实验的结果来看，运用先进的机器学习算法可以很大程度上提高预测性能，选取有效的特征和利用强大的机器学说明书CN10406363。
13、2A2/9页5习算法成为可以显著影响影响二硫键连接模式预测性能的两个重要方面。综上所述，现有技术中虽然已经在预测二硫键模式方面取得了一定的发展，不过在预测精度方面仍有提升的空间。发明内容0005本发明目的在于提供一种基于回归森林模型的蛋白质序列二硫键连接模式的预测方法，提高预测性能。0006为达成上述目的，本发明所采用的的技术方案如下0007一种基于回归森林模型的蛋白质序列二硫键连接模式的预测方法，其实现包括以下步骤0008步骤1、特征提取，基于输入的蛋白质序列信息，进行多视角特征提取和特征组合，即0009分别使用PSIBLAST算法抽取蛋白质序列的进化信息，利用PSIPRED算法抽取蛋白质序。
14、列的二级结构信息，再利用滑动窗口与特征串行组合方式从前述进化信息和二级结构信息中提取每个半胱氨酸残基对的多视角特征；0010然后计算蛋白质序列中半胱氨酸残基对在蛋白质序列中的距离特征信息DOC，利用PSIBLAST算法抽取蛋白质序列中每个半胱氨酸残基对的协变异信息CM，利用MODELLER同源建模软件获取半胱氨酸残基对的三维结构距离信息PDTCR；0011最后将上述从进化信息和二级结构信息中提取每个半胱氨酸残基对的多视角特征与前述距离特征信息、协变异信息、三维结构距离信息进行特征组合，得到蛋白质序列中每个半胱氨酸残基对最终的特征向量；0012步骤2、对于待预测的蛋白质序列信息和训练数据集合，采。
15、用上述步骤1的方法，生成待预测的蛋白质序列信息和训练数据集合中所有半胱氨酸残基对的特征向量，从而构成训练样本集合与待预测样本集合；0013步骤3、特征学习0014采用随机森林算法RANDOMFORESTS,RF学习二硫键样本在特征空间中的分布规律，生成随机森林回归模型，其中，模型的决策树个数NTREE采用大于等于300的整数，每次分裂时随机选择的特征个数MTRY由前述步骤1所得的特征向量的维度决定；0015步骤4、二硫键连接模式预测0016利用前述步骤3的随机森林回归模型对待预测样本集合中的特征向量进行预测，得到每个半胱氨酸残基对形成二硫键的倾向值，记为PPCP；0017令P为所有可能的二硫键。
16、连接模式的数量，则第I个二硫键连接模式的得分SI用下述公式得到00180019其中B代表第I个二硫键连接模式中二硫键的个数，PPCPJ是第I个二硫键连接模式中第J个半胱氨酸残基对的倾向值0020如果第I个二硫键的连接模式的得分最高，则该第I个模式就是最终预测得到的结果说明书CN104063632A3/9页600210022即该得分最高的二硫键连接模式I为最终预测得到的蛋白质序列中的二硫键连接模式。0023由以上本发明的技术方案可知，本发明的有益效果在于00241、提高预测速度与精度PDTCR特征直接计算了两个半胱氨酸之间的距离，而这个距离值有跟二硫键的成键有着直接的关系，所以可以显著提高预测精。
17、度；00252、在随机森林中运用的随机策略已经证明要优于许多其他的分类器，并且随机森林还有抗过拟合的特性，随机森林的训练和预测过程相比支持向量回归更加快速，训练时参数获取也更加容易，预测精度也有一定的提高。附图说明0026图1为本发明一实施方式基于回归森林模型的蛋白质序列二硫键连接模式的预测方法的实现流程示意图。具体实施方式0027为了更了解本发明的技术内容，特举具体实施例并配合所附图式说明如下。0028如图1所示，根据本发明的较优实施例，一种基于回归森林模型的蛋白质序列二硫键连接模式的预测方法，其实现包括以下步骤0029步骤1、特征提取，基于输入的蛋白质序列信息，进行多视角特征提取和特征组合。
18、，即0030分别使用PSIBLAST算法抽取蛋白质序列的进化信息，利用PSIPRED算法抽取蛋白质序列的二级结构信息，再利用滑动窗口与特征串行组合方式从前述进化信息和二级结构信息中提取每个半胱氨酸残基对的多视角特征；0031然后计算蛋白质序列中半胱氨酸残基对在蛋白质序列中的距离特征信息DOC，利用PSIBLAST算法抽取蛋白质序列中每个半胱氨酸残基对的协变异信息CM，利用MODELLER同源建模软件获取半胱氨酸残基对的三维结构距离信息PDTCR；0032最后将上述从进化信息和二级结构信息中提取每个半胱氨酸残基对的多视角特征与前述距离特征信息、协变异信息、三维结构距离信息进行特征组合，得到蛋白质。
19、序列中每个半胱氨酸残基对最终的特征向量；0033步骤2、对于待预测的蛋白质序列信息和训练数据集合，采用上述步骤1的方法，生成待预测的蛋白质序列信息和训练数据集合中所有半胱氨酸残基对的特征向量，从而构成训练样本集合与待预测样本集合；0034步骤3、特征学习0035采用随机森林算法RANDOMFORESTS,RF学习二硫键样本在特征空间中的分布规律，生成随机森林回归模型，其中，模型的决策树个数NTREE采用大于等于300的整数，每次分裂时随机选择的特征个数MTRY由前述步骤1所得的特征向量的维度决定；0036步骤4、二硫键连接模式预测0037利用前述步骤3的随机森林回归模型对待预测样本集合中的特征。
20、向量进行预测，说明书CN104063632A4/9页7得到每个半胱氨酸残基对形成二硫键的倾向值，记为PPCP；0038令P为所有可能的二硫键连接模式的数量，则第I个二硫键连接模式的得分SI用下述公式得到00390040其中B代表第I个二硫键连接模式中二硫键的个数，PPCPJ是第I个二硫键连接模式中第J个半胱氨酸残基对的倾向值0041如果第I个二硫键的连接模式的得分最高，则该第I个模式就是最终预测得到的结果00420043即该得分最高的二硫键连接模式I为最终预测得到的蛋白质序列中的二硫键连接模式。0044作为可选的实施方式，前述步骤1的实现过程中，其中的蛋白质序列的进化信息即位置特异性得分矩阵、。
21、蛋白质序列的二级结构信息、半胱氨酸对在序列中的距离特征DOC、半胱氨酸残基对的协变异信息CM、以及半胱氨酸残基对的三维结构距离信息PDTCR，其具体的提取和组合方式如下。0045A蛋白质序列的进化信息，即位置特异性得分矩阵0046对于某一给定的含有L个氨基酸残基的特定蛋白质P，其原始位置特异性得分矩阵PSSM首先由PSIBLAST程序用默认的EVALUE参数从SWISSPROT数据库该数据库为现有技术中的公共数据库经过三次迭代得到的，该矩阵为一个N行20列矩阵，表达如下00470048其中A、R、V表示20种常见氨基酸残基，OI,J表示蛋白质的第I个氨基酸残基在进化过程中突变成20种常见氨基酸。
22、残基中的第J种氨基酸残基的可能性。0049然后，对原始位置特异性得分矩阵中的每一个元素X，利用下述函数公式2归一化。00500051归一化后的位置特异性得分矩阵，表达如下0052说明书CN104063632A5/9页80053再次，对于蛋白质序列P中的一个半胱氨酸CII为半胱氨酸残基在序列中的位置，以位置特异性得分矩阵中的第I行为中心，使用一个宽度为W的窗口，该窗口内所有元素构成第I个氨基酸残基的PSSM特征矩阵，表达如下00540055最后，将上述特征矩阵4按行优先的组合方式形成维数为20W的特征向量，取滑动窗口大小为13，即W13，则半胱氨酸残基CI的遗传信息可以用下面特征向量表示0056。
23、0057对于一个半胱氨酸残基对，分别计算两个半胱氨酸残基的遗传信息特征向量即前述表达5所示的特征向量，并串行组合，即可得到一个半胱氨酸残基对的遗传信息特征向量。0058B、蛋白质序列的二级结构特征PSS0059对于一个含有L个氨基酸残基的蛋白质P，首先通过PSIPRED算法生成蛋白质序列P的二级结构矩阵，该矩阵大小为L3，表达如下00600061其中，C、H、E表示蛋白质的三种二级结构COIL、HELIX、STRAND，SI,1表示蛋白质说明书CN104063632A6/9页9P中半胱氨酸残基CI的二级结构是COIL的概率，SI,2表示蛋白质中半胱氨酸残基CI的二级结构是HELIX的概率，SI。
24、,3表示蛋白质中半胱氨酸残基CI的二级结构是STRAND的概率；然后，使用上述同样的宽度为W的滑动窗口与按行优先的组合方式，一个半胱氨酸残基可以由一个维数为3W的特征向量表示00620063对于一个半胱氨酸残基对，分别计算两个半胱氨酸残基的二级结构信息特征向量即前述表达7所示的特征向量，并串行组合，即可得到一个半胱氨酸残基对的二级结构信息特征向量。0064C、半胱氨酸对在序列中的距离特征DOC0065根据TSAI等人首先提出的半胱氨酸对的距离特征DOC研究，本实施例中采用下述公式来表示一个半胱氨酸对的距离特征0066DOCI,J|IJ|0067其中，I和J分别是两个半胱氨酸残基在蛋白质序列中的。
25、位置。0068然后，再使用DOCLOG作为归一化方式，得到距离特征DOC的最终表达00690070D、协变异特征CM0071如果半胱氨酸残基对存在某种较强的联系，那么这个半胱氨酸残基对就可能独立于进化过程，鉴于此，借鉴现有技术中RUBINSTEIN和FISER等人通过协变异分析进行二硫键预测研究，本实施例中采用其所提出的打分模式，对每个半胱氨酸残基对进行打分，这样每个半胱氨酸残基对都获得一个介于0,1的得分。0072关于协变异的具体打分模式，如RUBINSTEIN和FISER所著“PREDICTINGDISULDEBONDCONNECTIVITYINPROTEINSBYCORRELATEDMU。
26、TATIONSANALYSISJBIOINFORMATICS,2008,244498504”0073E、半胱氨酸残基对三维结构距离特征PDTCR0074使用MODELLER同源建模软件从蛋白质序列得到该蛋白质的三维结构；然后再利用得到的结构中的坐标值计算出两个半胱氨酸的距离。半胱氨酸的间距越小，它们形成一个二硫键的可能性就越大。0075最终，每个半胱氨酸残基对的特征向量将由上面提到的五个类型的特征组成，是位置特异性得分矩阵PSSM520维，5202260、二级结构特征PSS78维，78239、距离特征DOC1维、协变异特征CM1维、三维结构距离特征PDTCR1维，最终得到的特征向量，其维度为5。
27、2078111601维。0076如图1所示，对于待预测的蛋白质序列信息和训练数据集合，均采用上述步骤1的方法特征提取和特征串行组合，生成待预测的蛋白质序列信息和训练数据集合中所有半胱氨酸残基对的特征向量，从而分别构成训练样本集合与待预测样本集合。0077然后，在步骤3的特征学习中，采用随机森林算法RANDOMFORESTS,RF学习二硫键样本在特征空间中的分布规律，生成随机森林回归模型，其中，模型的决策树个数NTREE采用大于等于300的整数，每次分裂时随机选择的特征个数MTRY由前述步骤1说明书CN104063632A7/9页10所得的特征向量的维度决定。0078本实施例中，前述决策树个数N。
28、TREE优选500。每次分裂时随机选择的特征个数MTRY取值为其中601即前述特征向量的特征维度。0079接下来在步骤4中，采用前述步骤3的随机森林回归模型对待预测样本集合中的特征向量进行预测，得到每个半胱氨酸残基对形成二硫键的倾向值，记为PPCP；0080令P为所有可能的二硫键连接模式的数量，则第I个二硫键连接模式的得分SI用下述公式得到00810082其中B代表第I个二硫键连接模式中二硫键的个数，PPCPJ是第I个二硫键连接模式中第J个半胱氨酸残基对的倾向值0083如果第I个二硫键的连接模式的得分最高，则该第I个模式就是最终预测得到的结果00840085即该得分最高的二硫键连接模式I为最终。
29、预测得到的蛋白质序列中的二硫键连接模式。0086如图1所示，以图1中的一对半胱氨酸残基对作为示例，该两个半胱氨酸的位置分别是3和N1，如前述步骤1所述的，获取此半胱氨酸残基对的PSSM、PSS、CM、DOC还有PDTCR特征后，进行特征组合可得到一个601维的一个特征向量。0087再用已经训练完毕的随机森林回归模型预测这个特征向量，得到该个半胱氨酸残基对形成二硫键的倾向值。0088如前述步骤4的处理，取得分最高的二硫键的连接模式记为最终的预测结果。0089如图1所示，其中所有可能的二硫键模式包括三种C1C2,C3C4，C1C3,C2C4，C1C4,C2C3，按照图1所示，其得分别为0090S1。
30、030407；0091S2010203；0092S3000505。0093由此可见，第一个模式的得分最高，即为最终所预测的二硫键连接模式。0094本实施例中，将PDTCR特征与PSSM、PSS、CM和DOC特征相结合，在SP39数据集和PDBCYSR数据集的基础上，进行了对比试验只用不包含PDTCR特征的四个特征即PSSM、PSS、CM和DOC特征通过交叉验证，同样地，在前述两个数据集SP39数据集和PDBCYSR数据集上用四个原始特征即PSSM、PSS、CM和DOC特征加上PDTCR特征再进行交叉验证，然后得到预测性能的结果对比。0095其中，前述SP39数据集如FARISELLIP,CAS。
31、ADIOR“PREDICTIONOFDISULDECONNECTIVITYINPROTEINS”JBIOINFORMATICS,2001,1710957640096前述PDBCYSR数据集是根据SAVOJARDOC,FARISELLIP,ALHAMDOOSHM,ETAL“IMPROVINGTHEPREDICTIONOFDISULDEBONDSINEUKARYOTESWITHMACHINELEARNINGMETHODSANDPROTEINSUBCELLULARLOCALIZATION”JBIOINFORMATICS,2011,271622说明书CN104063632A108/9页112430所提。
32、出的PDBCYS数据集，本实施例中，基于此PDBCYS数据集把二硫键个数小于2或者大于5的序列都移除，剩下的263条蛋白质序列形成了PDBCYSR数据集。0097由于现有的预测方法多数是基于上述两个数据集，因此本实施例中，在SP39数据集和PDBCYSR两个数据集上，分别采用4轮交叉验证和20轮交叉验证，表41显示了分别使用随机森林RF本发明提出的方法和支持向量回归SVR在有无PDTCR特征的情况下的实验结果。0098表41分别用支持向量回归和随机森林在有无PDTCR特征的情况下的实验结果00990100AWITHOUTPDTCR四个传统特征,PSSM,PSS,CM,和DOC；0101BWIT。
33、HPDTCR四个传统特征加上PDTCR特征；0102C通过才用4轮交叉验证得到的结果。0103D通过采用20轮交叉验证得到的结果。0104其中0105QC用于衡量在所有的二硫键中可以正确预测出来的二硫键的百分比，它的定义如下01060107其中NC是正确预测出来的二硫键个数，TC是测试数据集中所有的二硫键个数。0108QP用于衡量所有蛋白质中可以正确预测出来的蛋白质的百分比。只有一个蛋白质中所有的二硫键都被预测出来并且预测出来的二硫键数量要跟真是的数量相同，这样才可以判断这个蛋白质被正确预测出来。QP的定义如下01090110其中NP是正确预测出来的蛋白质的数量，而TP是数据集中所有蛋白质的数。
34、量。0111表41中的数据表明，结合新的PDTCR特征，预测性能可以得到明显的提升。在二硫键个数B2、3、4、5不同的情况下，QP和QC的值在用了PDTCR特征之后相比不用新特征几乎都有提高。在SP39数据集跟PDBCYSR数据集上无论用支持向量回归亦或随机森林，结果都是如此。我们以总的QP和QC为例，在SP39数据集上加上了PDTCR特征之后，用SVR方法两者分别提高了36跟28，用RF方法两者分别提高了16和19。而在PDBCYSR数据集上，用SVR方法两者分别提高了80和73，用RF则分别提高了61和44。表2中的数据足以说明们新加特征的有效性。预测性能得到很大提升的原因是说明书CN10。
35、4063632A119/9页12PDTCR特征直接计算了两个半胱氨酸之间的距离，而这个距离值有跟二硫键的成键有着直接的关系。明显，在同样的特征下，使用随机森林比使用SVR可以获得更高预测的精度。0112而且，在训练随机森林模型时，随机森林还进行了特征选择，从而降低了维度并消除冗余，提高预测性能同时减少预测时间。0113虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明。本发明所属技术领域中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。说明书CN104063632A121/1页13图1说明书附图CN104063632A13。