DnaSAM: Software to perform neutrality testing for large datasets with complex null models

Patterns of DNA sequence polymorphisms can be used to understand the processes of demography and adaptation within natural populations. High-throughput generation of DNA sequence data has historically been the bottleneck with respect to data processing and experimental inference. Advances in marker technologies have largely solved this problem. Currently, the limiting step is computational, with most molecular population genetic software allowing a gene-by-gene analysis through a graphical user interface. An easy-to-use analysis program that allows both high-throughput processing of multiple sequence alignments along with the flexibility to simulate data under complex demographic scenarios is currently lacking. We introduce a new program, named DnaSAM, which allows high-throughput estimation of DNA sequence diversity and neutrality statistics from experimental data along with the ability to test those statistics via Monte Carlo coalescent simulations. These simulations are conducted using the ms program, which is able to incorporate several genetic parameters (e.g. recombination) and demographic scenarios (e.g. population bottlenecks). The output is a set of diversity and neutrality statistics with associated probability values under a user-specified null model that are stored in easy to manipulate text file.     

中国野生假俭草的研究现状

对近年来我国野生假俭草的研究现状进行了综述,并提出了我国野生假俭草发展前景及今后科研的方向。     

FLAGS: A Flexible and Adaptive Association Test for Gene Sets Using Summary Statistics

Genome-wide association studies (GWAS) have been widely used for identifying common variants associated with complex diseases. Despite remarkable success in uncovering many risk variants and providing novel insights into disease biology, genetic variants identified to date fail to explain the vast majority of the heritability for most complex diseases. One explanation is that there are still a large number of common variants that remain to be discovered, but their effect sizes are generally too small to be detected individually. Accordingly, gene set analysis of GWAS, which examines a group of functionally related genes, has been proposed as a complementary approach to single-marker analysis. Here, we propose a flexible and adaptive test for gene sets (FLAGS), using summary statistics. Extensive simulations showed that this method has an appropriate type I error rate and outperforms existing methods with increased power. As a proof of principle, through real data analyses of Crohn’s disease GWAS data and bipolar disorder GWAS meta-analysis results, we demonstrated the superior performance of FLAGS over several state-of-the-art association tests for gene sets. Our method allows for the more powerful application of gene set analysis to complex diseases, which will have broad use given that GWAS summary results are increasingly publicly available.     

Estimation in Longitudinal or Panel Data Models with Random‐Effect‐Based Missing Responses

Summary In studies with longitudinal or panel data, missing responses often depend on values of responses through a subject‐level unobserved random effect. Besides the likelihood approach based on parametric models, there exists a semiparametric method, the approximate conditional model (ACM) approach, which relies on the availability of a summary statistic and a linear or polynomial approximation to some random effects. However, two important issues must be addressed in applying ACM. The first is how to find a summary statistic and the second is how to estimate the parameters in the original model using estimates of parameters in ACM. Our study is to address these two issues. For the first issue, we derive summary statistics under various situations. For the second issue, we propose to use a grouping method, instead of linear or polynomial approximation to random effects. Because the grouping method is a moment‐based approach, the conditions we assumed in deriving summary statistics are weaker than the existing ones in the literature. When the derived summary statistic is continuous, we propose to use a classification tree method to obtain an approximate summary statistic for grouping. Some simulation results are presented to study the finite sample performance of the proposed method. An application is illustrated using data from the study of Modification of Diet in Renal Disease.     

生态学中的点格局研究概况及其在国内的应用

点格局分析是研究生态学格局的工具之一,近年来在生态学中的应用越来越多。为深入了解点格局分析方法在国内的研究与应用情况,以所总结的研究进展、一般步骤和基本要点为背景,分析评述了1996—2015年期间以点格局为主题的国内中文核心期刊文献。结果表明,在国内生态学格局研究中,应用研究占据主导地位,研究对象广泛,包括以树木为主的乔、灌、草等不同生活型的植物,甚至包括景观;基础研究,包括概括性统计量、零模型与点过程模型等方面,以及专用软件工具包的开发等研究薄弱。在应用中存在一定问题,主要表现为:概括性统计量使用单一,且以Ripley的K-函数及其变形为主;零模型(或点过程模型)是科学问题的统计表达,但是有一半以上的研究未明确给出零模型。建议在未来应用研究中重视多种统计量的组合使用和原假设的建立,在探讨热带、亚热带森林等具有复杂空间结构系统的多样性格局时,考虑对象的不同世代和系统的不同垂直层次,并加强多变量或三维概括性统计量的开发、点格局分析方法与动态过程模型的结合研究等工作。     

美国 FDA 新颁发政策与程序手册 ——“对基于问题审评的申报资料的药学审评”的分析解读

对美国 FDA 于 2014 年 11 月 18 日发布的政策与程序手册“对基于问题审评的申报资料的药学审评”进行介绍与讨论,包括发 布该手册的目的、背景、采用基于问题的药学审评的优势、相关政策、各自职责与程序等,以及原文附件中与原料药及制剂申报相关的 技术问题,以便于药学工作者了解与药学审评相关的问题的全貌。     

Estimating risk difference in multicenter studies under baseline-risk heterogeneity

In this paper, we consider the case of efficient estimation of the risk difference in a multicenter study allowing for baseline heterogeneity. We consider the optimally weighted estimator for the common risk difference and show that this estimator has considerable bias when the true weights (which are inversely proportional to the variances of the center-specific risk difference estimates) are replaced by their sample estimates. In addition, we propose a new estimator for this situation of the Mantel-Haenszel type that is unbiased and, in addition, has a smaller variance for small sample sizes within the study centers. Simulations illustrate these findings.