基于Progressive多序列比对方法的求解多序列比对的启发式算法

在生物信息学研究中,生物序列比对问题占有重要的地位。多序列比对问题是一个NPC问题,由于时间和空间的限制不能够求出精确解。文中简要介绍了Feng和Doolittle提出的多序列比对算法的基本思想,并改进了该算法使之具有更好的比对精度。实验结果表明,新算法对解决一般的progressive多序列比对方法中遇到的局部最优问题有较好的效果。     

Use of the EM algorithm to detect QTL affecting multiple-traits in an across half-sib family analysis

QTL detection experiments in livestock species commonly use the half-sib design. Each male is mated to a number of females, each female producing a limited number of progeny. Analysis consists of attempting to detect associations between phenotype and genotype measured on the progeny. When family sizes are limiting experimenters may wish to incorporate as much information as possible into a single analysis. However, combining information across sires is problematic because of incomplete linkage disequilibrium between the markers and the QTL in the population. This study describes formulæ for obtaining MLEs via the expectation maximization (EM) algorithm for use in a multiple-trait, multiple-family analysis. A model specifying a QTL with only two alleles, and a common within sire error variance is assumed. Compared to single-family analyses, power can be improved up to fourfold with multi-family analyses. The accuracy and precision of QTL location estimates are also substantially improved. With small family sizes, the multi-family, multi-trait analyses reduce substantially, but not totally remove, biases in QTL effect estimates. In situations where multiple QTL alleles are segregating the multi-family analysis will average out the effects of the different QTL alleles.     

基于遗传算法和支持向量机的乳杆菌发酵过程建模

由于生化反应过程的复杂性和高度非线性,多数简单的数学模型不能准确描述。该文基于Matlab软件,利用改进的支持向量机(υ-SVR)对植物乳酸杆菌发酵这一典型生化过程进行研究,应用遗传算法估计模型最优参数,建立植物乳杆菌的菌体密度预测模型。同时建立传统的logistic动力学模型以进行比较。结果表明,采用结合遗传算法的υ-SVR预测模型拟合误差小,皮尔森相关系数(R)更高,可以较好地预测乳酸杆菌的发酵过程,为其优化控制及放大提供依据。     

Nonparametric estimation of ROC curves in the absence of a gold standard

In the evaluation of diagnostic accuracy of tests, a gold standard on the disease status is required. However, in many complex diseases, it is impossible or unethical to obtain such a gold standard. If an imperfect standard is used, the estimated accuracy of the tests would be biased. This type of bias is called imperfect gold standard bias. In this article we develop a nonparametric maximum likelihood method for estimating ROC curves and their areas of ordinal-scale tests in the absence of a gold standard. Our simulation study shows that the proposed estimators for the ROC curve areas have good finite-sample properties in terms of bias and mean squared error. Further simulation studies show that our nonparametric approach is comparable to the binormal parametric method, and is easier to implement. Finally, we illustrate the application of the proposed method in a real clinical study on assessing the accuracy of seven specific pathologists in detecting carcinoma in situ of the uterine cervix.     

On computation of semiparametric maximum likelihood estimators with shape constraints

Large sample theory of semiparametric models based on maximum likelihood estimation (MLE) with shape constraint on the nonparametric component is well studied. Relatively less attention has been paid to the computational aspect of semiparametric MLE. The computation of semiparametric MLE based on existing approaches such as the expectation‐maximization (EM) algorithm can be computationally prohibitive when the missing rate is high. In this paper, we propose a computational framework for semiparametric MLE based on an inexact block coordinate ascent (BCA) algorithm. We show theoretically that the proposed algorithm converges. This computational framework can be applied to a wide range of data with different structures, such as panel count data, interval‐censored data, and degradation data, among others. Simulation studies demonstrate favorable performance compared with existing algorithms in terms of accuracy and speed. Two data sets are used to illustrate the proposed computational method. We further implement the proposed computational method in R package BCA1SG , available at CRAN.     

A neural network model for the prediction of membrane-spanning amino acid sequences

The architecture and weights of an artificial neural network model that predicts putative transmembrane sequences have been developed and optimized by the algorithm of structure evolution. The resulting filter is able to classify membrane/nonmembrane transition regions in sequences of integral human membrane proteins with high accuracy. Similar results have been obtained for both training and test set data, indicating that the network has focused on general features of transmembrane sequences rather than specializing on the training data. Seven physicochemical amino acid properties have been used for sequence encoding. The predictions are compared to hydrophobicity plots.     

气候变化下青藏高原全缘叶绿绒蒿的潜在分布变迁

青藏高原物种丰富且属于气候变化敏感区,研究气候变化对青藏高原物种的潜在分布影响,对于该区域物种多样性保护具有重要意义。该研究以一级濒危藏药植物全缘叶绿绒蒿为研究对象,利用加权平均算法(weighted average algorithm, WAA)构建随机森林(RF)、灵活判别分析(FDA)及人工神经网络(ANN)的集成模型,同时对比分析了WAA模型和不同生态位模型的预测精度。最后利用WAA模型预测了全缘叶绿绒蒿在当前(1970~2000年平均)和未来(2041~2060年平均)气候情景下的潜在分布,其中未来气候考虑了2种“共享社会经济路径”(SSP2-45和SSP5-85)。结果显示:(1) WAA模型的预测表明,基于RF、FDA和ANN的集成模型的AUC值为0.926,在AUC值最高RF模型的基础上提高了3%,在FDA和ANN模型的AUC值的基础上均提高了5%。(2) WAA模型确定,全缘叶绿绒蒿的潜在分布对年降水量和最暖季降水量最为敏感,其次是最热月份最高气温,同时对最湿月份降水量以及等温性表现出较低的敏感性。(3)当前全缘叶绿绒蒿潜在分布区主要分布在甘肃西南部、青海东部至南部、四川西部和西北部、云南西北部和东北部、西藏东部。(4)未来气候变化下青藏高原全缘叶绿绒蒿潜在分布预测表明,在2050年SSP2-45情景下,全缘叶绿绒蒿的潜在分布区大小与当前潜在分布区大小基本相同,但整体向西北方向高海拔高纬度地区迁移;在SSP5-85情景下,全缘叶绿绒蒿的潜在分布区明显收缩,且向西北高纬度高海拔地区延伸的趋势更加明显。