Simulated multi-observation-per-species CTmax data

A simulated dataset mimicking a thermal tolerance study where each species has multiple measurements of critical thermal maximum (CTmax) taken at different acclimation temperatures. The data-generating process is

Usage

ctmax_sim

Format

A data frame with 1,464 rows and 3 variables:

species

Character. Species name matching tree300 tips.

acclim_temp

Numeric. Acclimation temperature (degrees C). This is an observation-level covariate that varies within species.

CTmax

Numeric. Critical thermal maximum (degrees C). Contains NA values for unobserved species and MCAR missingness.

Source

Simulated. See data-raw/make_ctmax_sim.R for the generation script.

Details

$$\begin{aligned} &CTmax_{ij} = 38 + phylo_i \\ &\quad + 0.50 \times acclim\_temp_{ij} \\ &\quad + \epsilon_{ij} \end{aligned}$$

where \(phylo_i\) follows Brownian motion on tree300, the within-species acclimation response ratio is 0.50, and \(\epsilon_{ij} \sim N(0, 1.5^2)\), where the residual standard deviation is 1.5 degrees Celsius. Thirty percent of species are entirely unobserved (all CTmax values are NA), and an additional 15% of remaining observations are missing at random.

This dataset demonstrates pigauto's multi-observation and observation-level covariate support. Use with tree300 and set species_col = "species" in impute.

See also

tree300, impute