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.     

黄河首曲玛曲县高寒湿地景观格局演变

高寒湿地是青藏高原地区最重要的生态水源涵养区之一, 也是局部气候的有效调节者, 其动态变化与成因亟待深入研究。该研究基于遥感图像分析、地理信息系统空间分析和景观生态指数分析结合的方法, 以黄河首曲玛曲县高寒湿地为研究对象, 对1995-2018年6期湿地的动态变化进行研究。结果表明, 研究区湿地在1995-2010年间不断退化, 1995-2010年湿地面积总共减少了18 680.31 hm²。在2010-2018年间黄河首曲高寒湿地面积有所增加, 但与20世纪90年代相比, 21世纪初开始湿地的面积普遍呈现下降趋势; 1995-2010年湿地斑块数不断增加, 斑块密度不断增大, 平均斑块面积下降, 景观的破碎度升高; 2010-2015年湿地斑块数和斑块密度减少, 2015-2018年湿地斑块数和斑块密度增加, 平均斑块面积先增大后减小, 景观的破碎度先降低后升高。1995-2010年研究区高寒湿地景观Shannon多样性指数和Shannon均匀度指数均呈现下降的趋势, 湿地的景观结构趋于简单, 景观类型分布更加集中。2010-2018年湿地景观Shannon多样性指数和Shannon均匀度指数均呈现上升趋势, 湿地的景观结构趋于复杂, 景观类型增加且分布更加分散。进一步的驱动力分析表明, 引起黄河首曲高寒湿地景观格局演变的主要因素是蒸发量和降水量, 其次是人口数量和大牲畜数量等人类活动影响。气候因子是影响黄河首曲高寒湿地面积变化的主要原因, 过度的人类经济活动在一定程度上加剧了湿地的变化。     

Ageing Effect Evaluation on HD-sEMG Signals Using CCA Approach

ObjectivesThe objective of the proposed study is to exploit the technology of high-density surface electromyography (HD-sEMG), in order to evaluate the muscle activation in young and elderly subjects during a daily life gesture, namely, Sit To Stand (STS), using wireless connected ambulatory equipment (TMSi©) and Blind Source Separation (BSS) approach with Canonical Correlation Analysis (CCA).Materials and methodsSixteen subjects participated (50% females) divided into two categories (‘H1’: young (30.62 yrs ±5.92, 23.95 kg/m² ±3.08), versus ‘H2’: old (61.87 yrs ±7.98, 23.4 kg/m² ±3.38)), in the recording of HD-sEMG signals, using 32-electrodes square grids (4×8), during Sit To Stand (STS) motion, three times at spontaneous speed. The studied muscle is the Rectus Femoris (RF) muscle. The recorded HD-sEMG signals were analyzed with CCA approach to extract correlation coefficient sets according to two age categories (young versus old), in order to evaluate its discriminating power with ageing. Statistical tests (t-test) were used to evaluate the discrimination for these two categories.ResultsThe calculation of CCA correlation coefficients showed a significant difference between young and old category concerning the mean CCA correlation coefficient (P<0.001***) and also the standard deviation of the CCA correlation coefficients (P<0.0001****).ConclusionThe obtained results are promising and indicate a clear difference between the obtained source variability using CCA method between the young and the old tested subjects during daily life motion. Furthermore, these estimated sources seem to be impacted by both anatomical and functional modifications with ageing.     

Cluster analysis of pressure pain threshold maps from the trapezius muscle

The aim of this study was to investigate and present a new mapping method to describe muscle pain sensitivity based on the assessment of pressure pain threshold (PPT) over the trapezius muscle. PPT data were recorded from 36 points in 20 healthy males using a standardised grid. Points were clustered using the K-means algorithm with a fixed initialisation procedure. The total number of clusters was determined on the basis of (1) R ² evaluation of the clustering outcome compared against a desired 95% reduction in variance criterion and (2) the number of empty clusters. A minimum of three clusters were found which fulfilled the criteria. The proposed method enables the identification of a relation between the muscle subdivisions and pressure pain sensitivity within the trapezius. Further, the cluster analysis will enable the study of differences in pain sensitivity distributions between patients and controls and quantify the effect of intervention (physical or pharmacological treatments).     

PARma: identification of microRNA target sites in AGO-PAR-CLIP data

PARma is a complete data analysis software for AGO-PAR-CLIP experiments to identify target sites of microRNAs as well as the microRNA binding to these sites. It integrates specific characteristics of the experiments into a generative model. The model and a novel pattern discovery tool are iteratively applied to data to estimate seed activity probabilities, cluster confidence scores and to assign the most probable microRNA. Based on differential PAR-CLIP analysis and comparison to RIP-Chip data, we show that PARma is more accurate than existing approaches. PARma is available from http://www.bio.ifi.lmu.de/PARma     

Prediction of methionine oxidation risk in monoclonal antibodies using a machine learning method

ABSTRACT

Monoclonal antibodies (mAbs) have become a major class of protein therapeutics that target a spectrum of diseases ranging from cancers to infectious diseases. Similar to any protein molecule, mAbs are susceptible to chemical modifications during the manufacturing process, long-term storage, and in vivo circulation that can impair their potency. One such modification is the oxidation of methionine residues. Chemical modifications that occur in the complementarity-determining regions (CDRs) of mAbs can lead to the abrogation of antigen binding and reduce the drug’s potency and efficacy. Thus, it is highly desirable to identify and eliminate any chemically unstable residues in the CDRs during the therapeutic antibody discovery process. To provide increased throughput over experimental methods, we extracted features from the mAbs’ sequences, structures, and dynamics, used random forests to identify important features and develop a quantitative and highly predictive in silico methionine oxidation model.     

Apoptosis detection via automated algorithms to analyze biomarker translocation in reporter cells

Rapid, efficient, and robust quantitative analyses of dynamic apoptotic events are essential in a high-throughput screening workflow. Currently used methods have several bottlenecks, specifically, limitations in available fluorophores for downstream assays and misinterpretation of statistical image data analysis. In this study, we developed cytochrome-C (Cyt-C) and caspase-3/-8 reporter cell lines using lung (PC9) and breast (T47D) cancer cells, and characterized them from the response to apoptotic stimuli. In these two reporter cell lines, the spatial fluorescent signal translocation patterns served as reporters of activations of apoptotic events, such as Cyt-C release and caspase-3/-8 activation. We also developed a vision-based, tunable, automated algorithm in MATLAB to implement the robust and accurate analysis of signal translocation in single or multiple cells. Construction of the reporter cell lines allows live monitoring of apoptotic events without the need for any other dyes or fixatives. Our algorithmic implementation forgoes the use of simple image statistics for more robust analytics. Our optimized algorithm can achieve a precision greater than 90% and a sensitivity higher than 85%. Combining our automated algorithm with reporter cells bearing a single-color dye/fluorophore, we expect our approach to become an integral component in the high-throughput drug screening workflow.     

Comparing methods for metabolic network analysis and an application to metabolic engineering

Bioinformatics tools have facilitated the reconstruction and analysis of cellular metabolism of various organisms based on information encoded in their genomes. Characterization of cellular metabolism is useful to understand the phenotypic capabilities of these organisms. It has been done quantitatively through the analysis of pathway operations. There are several in silico approaches for analyzing metabolic networks, including structural and stoichiometric analysis, metabolic flux analysis, metabolic control analysis, and several kinetic modeling based analyses. They can serve as a virtual laboratory to give insights into basic principles of cellular functions. This article summarizes the progress and advances in software and algorithm development for metabolic network analysis, along with their applications relevant to cellular physiology, and metabolic engineering with an emphasis on microbial strain optimization. Moreover, it provides a detailed comparative analysis of existing approaches under different categories.