This function computes the expected log expression under the nonlinear genotype-by-treatment model used by map_qtl.

get_mu(g, t, b0, b1, b2, b3)

Arguments

g

A vector of the genotype coded as (0, 1, 2) or the imputation-based allelic dosage.

t

A vector of treatment indicators (0 or 1).

b0

A scalar specifying the intercept.

b1

A scalar specifying the genotype effect.

b2

A scalar specifying the treatment effect.

b3

A scalar specifying the genotype-by-treatment interaction effect.

Value

A numeric vector of expected log expression values.