The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

Some hypothesis tests for the covariance matrix when the dimension is large compared to the sample size

This paper analyzes whether standard covariance matrix tests work whendimensionality is large, and in particular larger than sample size. Inthe latter case, the singularity of the sample covariance matrix makeslikelihood ratio tests degenerate, but other tests based on quadraticforms of sample covariance matrix eigenvalues remain well-defined. Westudy the consistency property and limiting distribution of these testsas dimensionality and sample size go to infinity together, with theirratio converging to a finite non-zero limit. We find that the existingtest for sphericity is robust against high dimensionality, but not thetest for equality of the covariance matrix to a given matrix. For thelatter test, we develop a new correction to the existing test statisticthat makes it robust against high dimensionality.

Honey, I shrunk the sample covariance matrix

The central message of this paper is that nobody should be using the samplecovariance matrix for the purpose of portfolio optimization. It containsestimation error of the kind most likely to perturb a mean-varianceoptimizer. In its place, we suggest using the matrix obtained from thesample covariance matrix through a transformation called shrinkage. Thistends to pull the most extreme coefficients towards more central values,thereby systematically reducing estimation error where it matters most.Statistically, the challenge is to know the optimal shrinkage intensity,and we give the formula for that. Without changing any other step in theportfolio optimization process, we show on actual stock market data thatshrinkage reduces tracking error relative to a benchmark index, andsubstantially increases the realized information ratio of the activeportfolio manager.

Breeding system and genetic variance in the monogamous, semi-social shrew, Crocidura russula

The population-genetic consequences of monogamy and male philopatry (a rare breeding system in mammals) were investigated using microsatellite markers in the semisocial and anthropophilic shrew Crocidura russula. A hierarchical sampling design over a 16-km geographical transect revealed a large genetic diversity (h = 0.813) with significant differentiation among subpopulations (F-ST = 5-6%), which suggests an exchange of 4.4 migrants per generation. Demic effective-size estimates were very high, due both to this limited gene inflow and to the inner structure of subpopulations. These were made of 13-20 smaller units (breeding groups), comprising an estimate of four breeding pairs each. Members of the same breeding groups displayed significant coancestries (F-LS = 9-10%), which was essentially due to strong male kinship: syntopic males were on average related at the half-sib level. Female dispersal among breeding groups was not complete (similar to 39%), and insufficient to prevent inbreeding. From our results, the breeding strategy of C. russula seems less efficient than classical mammalian systems (polygyny and male dispersal) in disentangling coancestry from inbreeding, but more so in retaining genetic variance.

Covariance matrix formula for Birnbaum-Saunders regression models

The Birnbaum-Saunders regression model is commonly used in reliability studies. We derive a simple matrix formula for second-order covariances of maximum-likelihood estimators in this class of models. The formula is quite suitable for computer implementation, since it involves only simple operations on matrices and vectors. Some simulation results show that the second-order covariances can be quite pronounced in small to moderate sample sizes. We also present empirical applications.

Estimation of genetic variance for macro- and micro-environmental sensitivity using double hierarchical generalized linear models

Background: Genetic variation for environmental sensitivity indicates that animals are genetically different in their response to environmental factors. Environmental factors are either identifiable (e.g. temperature) and called macro-environmental or unknown and called micro-environmental. The objectives of this study were to develop a statistical method to estimate genetic parameters for macro- and micro-environmental sensitivities simultaneously, to investigate bias and precision of resulting estimates of genetic parameters and to develop and evaluate use of Akaike’s information criterion using h-likelihood to select the best fitting model. Methods: We assumed that genetic variation in macro- and micro-environmental sensitivities is expressed as genetic variance in the slope of a linear reaction norm and environmental variance, respectively. A reaction norm model to estimate genetic variance for macro-environmental sensitivity was combined with a structural model for residual variance to estimate genetic variance for micro-environmental sensitivity using a double hierarchical generalized linear model in ASReml. Akaike’s information criterion was constructed as model selection criterion using approximated h-likelihood. Populations of sires with large half-sib offspring groups were simulated to investigate bias and precision of estimated genetic parameters. Results: Designs with 100 sires, each with at least 100 offspring, are required to have standard deviations of estimated variances lower than 50% of the true value. When the number of offspring increased, standard deviations of estimates across replicates decreased substantially, especially for genetic variances of macro- and micro-environmental sensitivities. Standard deviations of estimated genetic correlations across replicates were quite large (between 0.1 and 0.4), especially when sires had few offspring. Practically, no bias was observed for estimates of any of the parameters. Using Akaike’s information criterion the true genetic model was selected as the best statistical model in at least 90% of 100 replicates when the number of offspring per sire was 100. Application of the model to lactation milk yield in dairy cattle showed that genetic variance for micro- and macro-environmental sensitivities existed. Conclusion: The algorithm and model selection criterion presented here can contribute to better understand genetic control of macro- and micro-environmental sensitivities. Designs or datasets should have at least 100 sires each with 100 offspring.

A class of improved heteroskedasticity-consistent covariance matrix estimators

The heteroskedasticity-consistent covariance matrix estimator proposed by White (1980), also known as HC0, is commonly used in practical applications and is implemented into a number of statistical software. Cribari–Neto, Ferrari & Cordeiro (2000) have developed a bias-adjustment scheme that delivers bias-corrected White estimators. There are several variants of the original White estimator that also commonly used by practitioners. These include the HC1, HC2 and HC3 estimators, which have proven to have superior small-sample behavior relative to White’s estimator. This paper defines a general bias-correction mechamism that can be applied not only to White’s estimator, but to variants of this estimator as well, such as HC1, HC2 and HC3. Numerical evidence on the usefulness of the proposed corrections is also presented. Overall, the results favor the sequence of improved HC2 estimators.

Control charts for monitoring the mean vector and the covariance matrix of bivariate processes

In this article, we propose new control charts for monitoring the mean vector and the covariance matrix of bivariate processes. The traditional tools used for this purpose are the T (2) and the |S| charts. However, these charts have two drawbacks: (1) the T (2) and the |S| statistics are not easy to compute, and (2) after a signal, they do not distinguish the variable affected by the assignable cause. As an alternative to (1), we propose the MVMAX chart, which only requires the computation of sample means and sample variances. As an alternative to (2), we propose the joint use of two charts based on the non-central chi-square statistic (NCS statistic), named as the NCS charts. Once the NCS charts signal, the user can immediately identify the out-of-control variable. In general, the synthetic MVMAX chart is faster than the NCS charts and the joint T (2) and |S| charts in signaling processes disturbances.

Monitoring the mean vector and the covariance ­matrix of multivariate processes with sample means and sample ranges

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A single chart with supplementary runs rules for monitoring the mean vector and the covariance matrix of multivariate processes

The MRMAX chart is a single chart based on the standardized sample means and sample ranges for monitoring the mean vector and the covariance matrix of multivariate processes. User's familiarity with the computation of these statistics is a point in favor of the MRMAX chart. As a single chart, the recently proposed MRMAX chart is very appropriate for supplementary runs rules. In this article, we compare the supplemented MRMAX chart and the synthetic MRMAX chart with the standard MRMAX chart. The supplementary and the synthetic runs rules enhance the performance of the MRMAX chart. © 2013 Elsevier Ltd.

Genetic Variants in Matrix Metalloproteinase-9 Gene Modify Metalloproteinase-9 Levels in Black Subjects

Altered matrix metalloproteinases (MMPs) levels are involved in cardiovascular diseases and increased MMP-9 levels enhance the cardiovascular risk in apparently healthy subjects. We investigated the effects of MMP-9 gene polymorphisms and haplotypes on the circulating MMP-9 levels in healthy black subjects and the effects of an MMP-2 polymorphism on the plasma MMP-2 concentrations. We studied 190 healthy subjects, nonsmokers, self-reported as blacks (18-63 years). Genotypes for the MMP-2 C-1306T polymorphism and the MMP-9 C-1562T, 90(CA)(14-24) and Q279R polymorphisms (rs243865, rs3918242, rs2234681, and rs17576, respectively) were determined by TaqMan (R) Allele Discrimination assay and real-time polymerase chain reaction or restriction fragment length polymorphism. Alleles for the 90(CA)(14-24) polymorphism were grouped as low (L) when there were < 21 and high (H) when there were >= 21 CA repeats. The plasma levels of MMP-2 and MMP-9 were determined by gelatin zymography. The software PHASE 2.1 was used to estimate the haplotypes frequencies. Although we found no effects of the MMP-9 C-1562T or the Q279R polymorphisms on MMP-9 levels, higher MMP-9 levels were associated with the HH genotype for the -90(CA)(14-24) polymorphism compared with the HL or LL genotypes. Lower MMP-9 levels were found in carriers of the CRL haplotype (combining the C, R, and L alleles for the MMP-9 polymorphisms) compared with the CRH haplotype. Consistent with this finding, the CRL haplotype was more commonly found in subjects with low MMP-9 levels. The MMP-2 C-1306T polymorphism had no effects on the plasma MMP-2 levels. Our results show that MMP-9 genetic variations modify MMP-9 levels in black subjects and may offer biochemical evidence implicating MMP-9 in the pathogenesis of cardiovascular diseases in blacks.