simulate_pheno_ts: Compute phenotypic dynamic Time Series of trait + demographic...

In ashander/phenoecosim: Simulations of eco-evo trait and demographic dynamics under changing environments

Description

Compute phenotypic dynamic Time Series of trait + demographic change under stabilizing selection as function of environment (after Lande Chevin)

Usage

simulate_pheno_ts(T, X, params, env_args, growth_fun = "independent",
  poisson = FALSE, varying_g = FALSE)

Arguments

T	end time, assuming start time of 1
X	parameters (z, a, b, wbar, logN, theta)
params	a list with (gamma_sh, omegaz, A, B, R0, var_a, Vb, delta, sigma_xi, rho_tau, fractgen)
env_args	extra args for env.fn
growth_fun	density-dependence ("independent", "gompertz", "thetalogistic", "ceiling")
poisson	(logical) return N(t+1) = Poisson(f(N(t)). NOTE: only implemented for independent and thetalog!
varying_g	varying genetic variance using stochastic house of cards (SHC) approximation. See Details.

Details

NB - for now assume Tchange = 0 and demography after CL 2010 For the SHC after Burger and Lynch (1995) and Kopp and Matuszewski (2013) \[ \sigma_SHC() = \frac2 V_m V_s1 + \fracV_s\alpha ^ 2 N_e \] For $\alpha^2$ the variance of the effect of new mutations, $V_m$ the genomic mutation rate per generation, and $N_e \approx 2 R_0 N / (2 R_0 - 1)$. As population size increases, the function approaches a limit of $2 V_s V_m$ $V_s = \omega^2 + \sigma_e^2$ (we hard code the latter as 1) is the width of the fitness function.

Value

a long matrix with columns zbar, abar, bbar, Wbar, Npop, theta

ashander/phenoecosim documentation built on May 10, 2019, 1:52 p.m.