# NOT RUN {
## Includes extra variables
vars = list(Intercept = -10.6,
distance = c(0, 15, 50),
volume = c(25, 50),
biomass = 100,
alive = 1,
tech_rep = 1:10,
rep = 1:3, Cq = 1)
## Intercept only
ans = sim_eDNA_lm(Cq ~ 1, vars,
betas = c(intercept = -15),
sigma_ln_eDNA = 1e-5,
std_curve_alpha = 21.2, std_curve_beta = -1.5)
print(ans)
ans = sim_eDNA_lm(Cq ~ distance + volume, vars,
betas = c(intercept = -10.6, distance = -0.05, volume = 0.1),
sigma_ln_eDNA = 1, std_curve_alpha = 21.2, std_curve_beta = -1.5)
# }
Run the code above in your browser using DataLab