《应用R进行QSAR的尝试.pptx》由会员分享，可在线阅读，更多相关《应用R进行QSAR的尝试.pptx（21页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、马斌：R在QSAR中的应用1关于QSAR的背景 定量构效关系(Quantitative Structure-Activity Relationship，简称QSAR)是指利用理论计算和统计分析工具来研究系列化合物结构与其效应之间的定量关系，即借助结构参数构建数学模型来描述化合物结构与活性之间的关系。第1页/共21页马斌：R在QSAR中的应用2环境科学中进行QSAR的目的由于测定化合物的各种毒性需要花费大量的人力、物力和财力，人们不可能对众多化学品进行逐一测定而利用QSAR 模型即可对化学品的生物毒性和环境行为进行预测，并筛选出具有潜在危害的化学品，这在环境科学研究中无疑是一件极具意义的工作 第
2、2页/共21页马斌：R在QSAR中的应用3R与QSAR统计分析工具数学模型构建第3页/共21页马斌：R在QSAR中的应用4QSAR中需要的数学方法R中的各种包提供几乎所有QSAR研究需要的数学方法如果将这些方法集合为一个QSAR专用的包，无论对R的发展还是QSAR的应用都是有利的第4页/共21页马斌：R在QSAR中的应用5QSAR过程与R1.结构参数的计算2.R的meta-analysis包1.主观筛选2.客观筛选1.多元分析2.神经网络1.交叉验证法 2.蒙特卡罗验证 1.输入未知参数第5页/共21页马斌：R在QSAR中的应用6R中包含QSAR数学方法的包逐步回归stepPLrstepwis
3、e前进法forwardlars遗传算法rfwdmvgafitgenalg多元分析lmtestbpcadepthheplotsICSLabdsvLeaps神经网络NeuralNeuralnet偏最小二乘法GplsPlsLsplsPlsPlsgenomiPplsgplsSpls第6页/共21页马斌：R在QSAR中的应用7R是QSAR的有力工具，但是需要的数学方法分散在大量不同的包中，就如无数的宝藏埋藏在世界各地。第7页/共21页马斌：R在QSAR中的应用8应用R进行QSAR的一个尝试以R进行多环芳烃（PAHs）在植物根际消解效果的QSAR分析为例子来说明R再QSAR中的应用。第8页/共21页马斌：
4、R在QSAR中的应用9为什么研究多环芳烃（PAHs）第9页/共21页马斌：R在QSAR中的应用10数据库的建立活性参数是PAHs在植物根际消解效果，采用meta-analysis中常用的效应值(effect sizes,d)表示。d=ln(E/C)分子结构参数通过Dragon 5(Talanet)计算得到965个结构参数，包括拓扑参数和理化参数。第10页/共21页马斌：R在QSAR中的应用11参数的筛选用主观选择对结构参数进行筛选，去掉所有值为恒量的参数，然后计算余下参数的相关系数矩阵，去掉相关系数大于等于0.95的两个参数中的一个。采用R极大提高筛选效率 第11页/共21页马斌：R在QSAR
5、中的应用12R中结构参数的主观筛选（去除恒量）mol.structuredim.molanfor(k in 1:(dim.mol2-n)if(mol.structure1,k=mol.structurea,k)+mol.structure-mol.structure,-k;+n-n+1;+kwrite.table(mol.structure,file=molstruc1.csv,sep=,)第12页/共21页马斌：R在QSAR中的应用13R中结构参数的主观筛选（去除高相关性参数）cor.matrixdim.cormfor(i in 1:(dim.cor1-1)for(j in(i+1):(di
6、m.cor2)if(cor.matrixi,j=0.95)+mol.structure-mol.structure,-i;+idim(mol.structure)write.table(cor.matrix,file=matrix.csv,sep=,)write.table(mol.structure,file=molstruc2.csv,sep=,)第13页/共21页马斌：R在QSAR中的应用14PLS构建模型#modeling QSAR by PLSlibrary(pls)lnRdataQSAR#cross-validation the QSAR modelcrossvalRMSEPR2p
7、lot(crossval)plot(RMSEP)plot(R2)plot(QSAR,ncomp=15)第15页/共21页马斌：R在QSAR中的应用16RMSE和R2R2RMSE第16页/共21页马斌：R在QSAR中的应用17优化预测模型第17页/共21页马斌：R在QSAR中的应用18用于预测和预报#predict unstudied compoundperylenenew.mol-predict(QSAR,type=response,newdata=perylene)输入还没有研究资料的perylene的结构参数，得到效应值为-0.86第18页/共21页马斌：R在QSAR中的应用19结论与展望
8、R在QSAR分析中的最大特点是快捷和简便。QSAR的模型构建、验证和应用过程中都有多种方法可以选择，而这些方法目前都分布在不同的包中收集和整理各种常用的QSAR用到的方法，编写QSAR常用过程的函数，并开发出不断更新的包就尤为重要本文为R的QSAR包作出了一个开端第19页/共21页马斌：R在QSAR中的应用20致谢感谢第一届R会议的所有贡献者本文的完成要感谢中国人民大学谢益辉的鼓励和杭州电子科技大学郑冰的大力帮助本文得到国家863计划（2007AA061101）、国家自然基金（40671092,20707020）和浙江省重点科技项目（2008C13024-3）的资助 第20页/共21页马斌：R在QSAR中的应用21感谢您的观看！第21页/共21页