The Method of Random Skewers

Ecological and evolutionary studies are often concerned with the properties of covariance matrices. The method of random skewers (RS method) has been used compare a matrix to an a priori vector or to compare two matrices. The method involves multiplying a matrix by many random vectors drawn from a uniform distribution over all possible vector directions. The comparisons are usually made using the average angle (or cosine) of the response vectors to an a priori vector or to the response vectors corresponding from another matrix. Angles are usually constrained to the interval 0°–90° because the distribution of response vectors is bipolar bimodal. The size of the average angle or cosine depends strongly on the relative sizes of the eigenvalues (especially the first). The distribution of angles between pairs of response vectors from two covariance matrices is more complicated because it depends on the differences in orientation of the eigenvectors and the relative sizes of the eigenvalues of the both matrices. The average absolute value of the angles between these pairs of response vectors depends on the relative sizes of the eigenvalues of the matrices making it difficult to interpret its meaning without knowledge of the eigenvalues and eigenvectors of the two matrices. Thus, it is simpler to just directly compare matrices in terms of these quantities.     

Predicting the QRS complex and detecting small changes using principal component analysis.

In this paper, a new method for QRS complex analysis and estimation based on principal component analysis (PCA) and polynomial fitting techniques is presented. Multi-channel ECG signals were recorded and QRS complexes were obtained from every channel and aligned perfectly in matrices. For every channel, the covariance matrix was calculated from the QRS complex data matrix of many heartbeats. Then the corresponding eigenvectors and eigenvalues were calculated and reconstruction parameter vectors were computed by expansion of every beat in terms of the principal eigenvectors. These parameter vectors show short-term fluctuations that have to be discriminated from abrupt changes or long-term trends that might indicate diseases. For this purpose, first-order poly-fit methods were applied to the elements of the reconstruction parameter vectors. In healthy volunteers, subsequent QRS complexes were estimated by calculating the corresponding reconstruction parameter vectors derived from these functions. The similarity, absolute error and RMS error between the original and predicted QRS complexes were measured. Based on this work, thresholds can be defined for changes in the parameter vectors that indicate diseases.     

Perils of parsimony: properties of reduced-rank estimates of genetic covariance matrices

Eigenvalues and eigenvectors of covariance matrices are important statistics for multivariate problems in many applications, including quantitative genetics. Estimates of these quantities are subject to different types of bias. This article reviews and extends the existing theory on these biases, considering a balanced one-way classification and restricted maximum-likelihood estimation. Biases are due to the spread of sample roots and arise from ignoring selected principal components when imposing constraints on the parameter space, to ensure positive semidefinite estimates or to estimate covariance matrices of chosen, reduced rank. In addition, it is shown that reduced-rank estimators that consider only the leading eigenvalues and -vectors of the "between-group" covariance matrix may be biased due to selecting the wrong subset of principal components. In a genetic context, with groups representing families, this bias is inverse proportional to the degree of genetic relationship among family members, but is independent of sample size. Theoretical results are supplemented by a simulation study, demonstrating close agreement between predicted and observed bias for large samples. It is emphasized that the rank of the genetic covariance matrix should be chosen sufficiently large to accommodate all important genetic principal components, even though, paradoxically, this may require including a number of components with negligible eigenvalues. A strategy for rank selection in practical analyses is outlined.     

人类群体遗传结构的协方差阵主成分分析方法

目的:探讨基因频率矩阵的中心化(或均值化)协方差阵主成分分析方法在人类群体遗传结构研究中的适用性和合理性。方法:从基因频率矩阵的结构特征入手,分析中心化、均值化协方差阵主成分分析与标准化相关阵主成分分析在特征根、特征向量以及降维效果等方面的差异,并通过实例比较不同方法在解释群体遗传结构特征上合理性。结果:中心化(或均值化)协方差阵的主成分不仅反映了基因变异程度的“方差信息量权”,而且反映了基因间相互影响程度的“相关信息量权”;标准化相关阵的主成分反映的仅是“相关信息量权”,不包括“方差信息量权”。通过比较中国26个汉族人群HLA-A基因座中心化协方差阵和标准化相关阵2种主成分分析结果,证实中心化协方差阵主成分分析方法在特征根与特征向量、保留主成分的个数和对主成分的群体遗传学解释的合理性等方面均优于标准化相关阵主成分分析方法。结论:在对群体遗传结构进行主成分分析时,应使用中心化(或均值化)变换消除基因频率矩阵中量级的影响,然后在用其协方差阵提取主成分。     

Measuring the structural similarity of two communities composed of the same species

Summary The primary purpose of this paper is to propose empirical measures of the structural differences between two communities of plants or animals composed of the same species. Structure is defined to consist of; 1) the species in the community, 2) the pattern of interactions as represented by the covariance or correlation matrix of successive observations on each species, and 3) the mean abundances of each species in each of the two communities. Statistical tests are proposed for testing whether the covariance matrices and the vectors of mean densities for each community are equal and empirical measures of the differences between the covariance matrices and mean vectors are proposed. Given unequal covariance or correlation matrices the factor analysis model is used to derive empirical measures of the degree to which each variable of the ecosystem is responsible for the observed defferences in the pattern of interactions in each community. These tests and measures were applied to data gathered byHunter (1966) on the abundances of six species ofDrosophila censused monthly over a period of approximately two and a half years in two adjacent, but different habitats near Bogota, colombia. The two covariance matrices were significantly different indicating different patterns of interactions in the twoDrosophila communities.     

Persistence of changes in the genetic covariance matrix after a bottleneck

Abstract.— Genetic variance, phenotypic variance, and the genetic covariance matrix ( G ) can change as a result of genetic drift. These changes will persist over time to some extent and will continue if population size remains relatively small. Nine populations founded by a single pair of Drosophila melanogaster were measured for a series of six morphological characteristics for a large number of parent-offspring families at both the third generation after the bottlenecks and after 20 generations. From these data, the phenotypic variance, additive genetic variance, and G were estimated for each line at each generation. Phenotypic and genetic variances were highly correlated over time, so that the measurements made at the third generation were predictive of the state of the population 17 generations later. Genetic covariances were also somewhat stable over time; however, the G matrices of some lines changed significantly over the intervening generations. This change did not return the populations toward their original state before the population bottlenecks. We conclude that the genetic covariance matrix can change as a result of mild genetic drift over a short span of time.     

Neighborhood Dependence in Bayesian Spatial Models

The conditional autoregressive model and the intrinsic autoregressive model are widely used as prior distribution for random spatial effects in Bayesian models. Several authors have pointed out impractical or counterintuitive consequences on the prior covariance matrix or the posterior covariance matrix of the spatial random effects. This article clarifies many of these puzzling results. We show that the neighborhood graph structure, synthesized in eigenvalues and eigenvectors structure of a matrix associated with the adjacency matrix, determines most of the apparently anomalous behavior. We illustrate our conclusions with regular and irregular lattices including lines, grids, and lattices based on real maps.     

Selection Response Decomposition (SRD): A New Tool for Dissecting Differences and Similarities Between Matrices

Genetic and phenotypic variance/covariance matrices are a fundamental measure of the amount of variation and the pattern of association among traits for current investigations in evolutionary biology. Still, few methods have been developed to accomplish the goal of pinpointing in which traits two matrices differ most, hampering further works on the field. We here described a novel method for dissecting matrix comparisons. This method is called Selection Response Decomposition and is an extension of the random skewers in the sense that evolutionary responses produced by known simulated selection vectors are unfolded and then compared in terms of the direct and indirect responses to selection for any trait. We also applied the method in diverse case studies, illustrating its potential. Both theoretical matrices and empirical biological data were used in the comparisons made. In the theoretical ones, the method was able to determine exactly which traits were responsible for the known a priori differences between the matrices, as well as where matrices remained similar to each other. Similar support could be observed in comparisons carried on between matrices produced from empirical biological data, since reasonable and detailed interpretations could be made regarding matrix comparisons. SRD represents an excellent tool for matrix comparisons and should provide quantitative evolutionary biology with a new method for analyzing and comparing variance/covariance patterns.     

Genetic evaluation methods for populations with dominance and inbreeding

Summary The effect of inbreeding on mean and genetic covariance matrix for a quantitative trait in a population with additive and dominance effects is shown. This genetic covariance matrix is a function of five relationship matrices and five genetic parameters describing the population. Elements of the relationship matrices are functions of Gillois (1964) identity coefficients for the four genes at a locus in two individuals. The equivalence of the path coefficient method (Jacquard 1966) and the tabular method (Smith and Mäki-Tanila 1990) to compute the covariance matrix of additive and dominance effects in a population with inbreeding is shown. The tabular method is modified to compute relationship matrices rather than the covariance matrix, which is trait dependent. Finally, approximate and exact Best Linear Unbiased Predictions (BLUP) of additive and dominance effects are compared using simulated data with inbreeding but no directional selection. The trait simulated was affected by 64 unlinked biallelic loci with equal effect and complete dominance. Simulated average inbreeding levels ranged from zero in generation one to 0.35 in generation five. The approximate method only accounted for the effect of inbreeding on mean and additive genetic covariance matrix, whereas the exact accounted for all of the changes in mean and genetic covariance matrix due to inbreeding. Approximate BLUP, which is computable for large populations where exact BLUP is not feasible, yielded unbiased predictions of additive and dominance effects in each generation with only slightly reduced accuracies relative to exact BLUP.     

Proportionality of Diagonal Blocks with Blockwise Independence

In this paper, the tests of similarities among group covariance matrices and the differences among block covariance matrices within a group under repeated measurement model are studied. There are nine hierarchical nested structures of covariance matrices which have been tested. The likelihood ratio tests have been derived for these nine hierarchically structured models. An algorithm for determining the numerical solution of the corresponding maximum likelihood equations is also given.     

QUANTITATIVE GENETICS OF MORPHOLOGICAL DIFFERENTIATION IN PEROMYSCUS. I. TESTS OF THE HOMOGENEITY OF GENETIC COVARIANCE STRUCTURE AMONG SPECIES AND SUBSPECIES

Phenotypic and additive genetic covariance matrices were estimated for 15 morphometric characters in three species and subspecies of Peromyscus. Univariate and multivariate ANOVAs indicate these groups are highly diverged in all characters, P. leucopus having the largest body size, P. maniculatus bairdii the smallest, and P. maniculatus nebrascensis being intermediate. Comparing the structure of P and G within each taxon revealed significant similarities in all three cases. This proportionality was strong enough to justify using P in the place of G to analyze evolutionary processes using quantitative genetic models when G can not be estimated, as in fossil material. However, the similarity between genetic and phenotypic covariance structures is sufficiently low that estimates of the genetic parameters should be used when possible. The additive genetic covariance matrices were compared to examine the assumption that they remain constant during evolution, an assumption which underlies many applications of quantitative-genetic models. While matrix permutation tests indicated statistically significant proportionality between the genetic covariance structures of the two P. maniculatus subspecies, there is no evidence of significant genetic structural similarity between species. This result suggests that the assumption of constant genetic covariance structure may be valid only within species. (It does not, however, necessarily imply a causal relationship between speciation and heterogeneity of genetic covariance structures.) The low matrix correlation for the two P. maniculatus subspecies' genetic covariance matrices indicates G may not be functionally constant, even within species. The lack of similarity observed here may be due partly to sampling variation.     

Direct product-matrix method treatment of macromolecular binding

A major limitation of the well-known matrix method treatment of linear polymer systems is the large size of the transfer matrix required to account simultaneously for such factors as polymer conformation, polymer sequence, ligand binding, etc. This paper shows how to combine such factors by using a direct product of small matrices, one for each factor, to generate the large transfer matrix. The properties of the direct product operation can be used to simplify considerably the calculation of eigenvalues and eigenvectors of the transfer matrix. This permits one to develop fairly simply the statistical thermodynamics of relatively complicated macromolecular binding systems involving a large number of states for the individual polymer sites. Applications include treatments of sequence-specific binding of ligands, sequence-specific conformational changes, and the combination of the two processes.     

Constraints and missing reactions in the urea cycle.

The stoichiometric relations in a series of biochemical reactions are summarized by a stoichiometric number matrix (with a column for each reaction) and a conservation matrix (with a row for each constraint). These two matrices for a series or cycle of biochemical reactions are related because the columns of the stoichiometric number matrix are in the null space of the conservation matrix, and the rows of the transpose of the conservation matrix are in the null space of the transpose of the stoichiometric number matrix. The conservation matrix for a system of biochemical reactions is of interest because it shows the nature of the constraints in addition to the conservation of atoms and groups. Constraints beyond those for the conservation of atoms and groups indicate "missing reactions" that do not occur because the enzymes involved couple reactions that could occur and still conserve atoms and groups. The interpretation of conservation matrices and stoichiometric matrices for a reaction system is complicated by the fact that they are not unique. However, their row-reduced forms are unique, as are their dimensions, which represent the number of reactants and number of independent reactions. Two matrices that look different contain the same information if they have the same row-reduced form. The urea cycle, which involves five enzyme-catalyzed reactions, and its net reaction are discussed in terms of the linear constraints produced by enzyme catalysis. A procedure to obtain a set of conservation equations that will yield the correct net reaction is described.     

Testing substitution models within a phylogenetic tree

Phylogenetic tree reconstruction frequently assumes the homogeneity of the substitution process over the whole tree. To test this assumption statistically, we propose a test based on the sample covariance matrix of the set of substitution rate matrices estimated from pairwise sequence comparison. The sample covariance matrix is condensed into a one-dimensional test statistic Delta = sum ln(1 + delta(i)), where delta(i) are the eigenvalues of the sample covariance matrix. The test does not assume a specific mutational model. It analyses the variation in the estimated rate matrices. The distribution of this test statistic is determined by simulations based on the phylogeny estimated from the data. We study the power of the test under various scenarios and apply the test to X chromosome and mtDNA primate sequence data. Finally, we demonstrate how to include rate variation in the test.     

THE COVARIANCE STRUCTURE OF LIFE-HISTORY CHARACTERS IN DAPHNIA PULEX

The genetic covariance structure for life-history characters in two populations of cyclically parthenogenetic Daphnia pulex indicates considerable positive correlation among important fitness components, apparently at odds with the expectation if antagonistic pleiotropy is the dominant cause of the maintanence of genetic variation. Although there is no genetic correlation between offspring size and offspring number, present growth and present reproduction are both strongly positively correlated genetically with future reproduction, and early maturity is genetically correlated with larger clutch size. Although the ubiquity of antagonistic pleiotropy has been recently questioned, there are peculiarities of cyclical parthenogenesis that could lead to positive life-history covariance even when negative covariance would be expected in a similar sexual species. These include the influence of nonadditive gene action on evolution in clonally reproducing organisms, and the periodic release of hidden genetic variance within populations of cyclical parthenogens. Examination of matrix similarity, using the bootstrap for distribution-free hypothesis testing, reveals no evidence to suggest that the genetic covariance matrices differ between the populations. However, there is considerable evidence that the phenotypic and environmental covariance matrices differ between populations. These results indicate approximate stability of the genetic covariance matrix within species, an important assumption of many phenotypic evolution models, but should caution against the use of phenotypic in place of genetic covariance matrices.     

Weighted Concordance Analysis

Using the variance stabilizing technique, a product multinomial model is introduced to generate a new statistic to test observers' uncertainty in a weighted concordance analysis. Distance matrices which follow some specific rules are obtained by linear combinations of hierarchical distance matrices whose elements are equal to 0 or 1 and unit diagonal. The new statistic is compared with the kappa statistic interpreted by considering the covariance matrix generated by the data. By rewriting the test statistic in a barycentric form, one demonstrates how to modify the barycentric coefficients to derive an adequate measure of the interobserver agreement. The methods are illustrated using two examples.     

Local feature analysis: a statistical theory for reproducible essential dynamics of large macromolecules

Multivariate statistical methods are widely used to extract functional collective motions from macromolecular molecular dynamics (MD) simulations. In principal component analysis (PCA), a covariance matrix of positional fluctuations is diagonalized to obtain orthogonal eigenvectors and corresponding eigenvalues. The first few eigenvectors usually correspond to collective modes that approximate the functional motions in the protein. However, PCA representations are globally coherent by definition and, for a large biomolecular system, do not converge on the time scales accessible to MD. Also, the forced orthogonalization of modes leads to complex dependencies that are not necessarily consistent with the symmetry of biological macromolecules and assemblies. Here, we describe for the first time the application of local feature analysis (LFA) to construct a topographic representation of functional dynamics in terms of local features. The LFA representations are low dimensional, and like PCA provide a reduced basis set for collective motions, but they are sparsely distributed and spatially localized. This yields a more reliable assignment of essential dynamics modes across different MD time windows. Also, the intrinsic dynamics of local domains is more extensively sampled than that of globally coherent PCA modes.     

Patterns of phenotypic covariation and correlation in modern humans as viewed from morphological integration

Proportionality of phenotypic and genetic distance is of crucial importance to adequately focus on population history and structure, and it depends on the proportionality of genetic and phenotypic covariance. Constancy of phenotypic covariances is unlikely without constancy of genetic covariation if the latter is a substantial component of the former. If phenotypic patterns are found to be relatively stable, the most probable explanation is that genetic covariance matrices are also stable. Factors like morphological integration account for such stability. Morphological integration can be studied by analyzing the relationships among morphological traits. We present here a comparison of phenotypic correlation and covariance structure among worldwide human populations. Correlation and covariance matrices between 47 cranial traits were obtained for 28 populations, and compared with design matrices representing functional and developmental constraints. Among-population differences in patterns of correlation and covariation were tested for association with matrices of genetic distances (obtained after an examination of 10 Alu-insertions) and with Mahalanobis distances (computed after craniometrical traits). All matrix correlations were estimated by means of Mantel tests. Results indicate that correlation and covariance structure in our species is stable, and that among-group correlation/covariance similarity is not related to genetic or phenotypic distance. Conversely, genetic and morphological distance matrices were highly correlated. Correlation and covariation patterns were largely associated with functional and developmental factors, which probably account for the stability of covariance patterns.