Abstract:Random effects are implemented for aster models using two approximations taken from Breslow and Clayton [J. Amer. Statist. Assoc. 88 (1993) 9-25]. Random effects are analytically integrated out of the Laplace approximation to the complete data log likelihood, giving a closed-form expression for an approximate missing data log likelihood. Third and higher derivatives of the complete data log likelihood with respect to the random effects are ignored, giving a closed-form expression for second derivatives of the approximate missing data log likelihood, hence approximate observed Fisher information. This method is applicable to any exponential family random effects model. It is implemented in the CRAN package aster (R Core Team [R: A Language and Environment for Statistical Computing (2012) R Foundation for Statistical Computing], Geyer [R package aster (2012) <a class="link-external link-http" href="http://cran.r-project.org/package=aster" rel="external noopener nofollow">this http URL</a>]). Applications are analyses of local adaptation in the invasive California wild radish (Raphanus sativus) and the slender wild oat (Avena barbata) and of additive genetic variance for fitness in the partridge pea (Chamaecrista fasciculata).

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to expand the Aster model by introducing random effects in order to better analyze biological adaptability and the influence of genetics on adaptability. Specifically, the paper focuses on: 1. **Analysis of local adaptability**: Study the adaptability of different populations in different environments, especially in the invasive species California wild radish (*Raphanus sativus*) and slender wild oat (*Avena barbata*). These analyses usually involve nested categorical variables (such as populations nested in regions and plots nested in sites), so random effects need to be dealt with. 2. **Analysis of additive genetic variation**: Study the influence of additive genetic variation of sunflower pea (*Chamaecrista fasciculata*) on adaptability. This involves using a random - effect model to estimate genetic variance components, thereby evaluating the influence of genetic differences between individuals on adaptability. ### Solutions To achieve the above goals, the paper adopts the following methods: - **Approximate Integrated Likelihood (AIL)**: Use the Laplace approximation to deal with the complex likelihood function of the exponential family model containing random effects. This method can analytically integrate out the random effects, thereby obtaining a closed - form approximate log - likelihood function for missing data. - **Ignore high - order derivatives**: To simplify the calculation, the paper assumes that the third - and higher - order derivatives of the complete - data log - likelihood function with respect to the random effects are zero, thereby obtaining approximate second - order derivatives and the observed Fisher information matrix. - **Parameter transformation**: Introduce the square root of the variance component to avoid the problem of zero variance components and make the optimization process more stable. ### Application examples The paper shows the application of these methods through three specific examples: 1. **California wild radish (*Raphanus sativus*)**: Studied the adaptability of different populations in different environments, especially the performance of the northern coastal and southern inland populations at two experimental sites. The results show that each population shows higher adaptability in its native environment, supporting the conclusion of local adaptability. 2. **Slender wild oat (*Avena barbata*)**: Studied the performance of two ecotypes (xerophytic and hygrophytic) in different years and experimental sites, and evaluated the changes in survival rate and inflorescence number. 3. **Sunflower pea (*Chamaecrista fasciculata*)**: Analyzed the influence of additive genetic variation on adaptability, especially by estimating genetic variance components to evaluate the genetic differences between individuals. ### Conclusions The paper expands the Aster model by introducing random effects, providing an effective method to analyze biological adaptability and the influence of genetics on adaptability. These methods perform well in dealing with complex data structures (such as nested categorical variables) and estimating genetic variance components, providing a powerful tool for ecological and evolutionary biology research.

Local adaptation and genetic effects on fitness: Calculations for exponential family models with random effects

Linear mixed effects models for non‐Gaussian continuous repeated measurement data

Combining envelope methodology and aster models for variance reduction in life history analyses

Application of a genomic model for high-dimensional chemometric analysis

Adaptive Fitting of Linear Mixed-Effects Models with Correlated Random-effects

The accuracy of predicting maladaptation to new environments with genomic data

Epistasis and the structure of fitness landscapes: are experimental fitness landscapes compatible with Fisher's Geometric model?

Metabolism, Growth, and the Energetic Definition of Fitness: A Quantitative Genetic Study in the Land Snail Cornu aspersum

On the plains and prairies of Minnesota: The role of mathematical statistics in biological explanation

Generalized Single Index Models and Jensen Effects on Reproduction and Survival

The limits of predicting maladaptation to future environments with genomic data

Improving a nonlinear Gompertz growth model using bird-specific random coefficients in two heritage chicken lines

Quantitative analyses of empirical fitness landscapes

Inference of the distribution of fitness effects of mutations is affected by single nucleotide polymorphism filtering methods, sample size and population structure

Adaptation to a changing environment: what me Normal?

Exponential family measurement error models for single-cell CRISPR screens

lme4breeding: enabling genetic evaluation in the era of genomic data

Random effects compound Poisson model to represent data with extra zeros

Likelihood-based Inference for Skewed Responses in a Crossover Trial Setup

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

Bayesian hierarchical modelling for inferring genetic interactions in yeast