X | Y = y ~ N(mu_x + rho sd_x (y - mu_y) / sd_y, sd_x sqrt(1 - rho^2)) and symmetrically for Y | X (Hedderich, Sachs & Reynarowych 2023, p. 324).
Usage
bvncnd(value, mu_x, mu_y, sd_x, sd_y, rho, given = c("y", "x"))X | Y = y ~ N(mu_x + rho sd_x (y - mu_y) / sd_y, sd_x sqrt(1 - rho^2)) and symmetrically for Y | X (Hedderich, Sachs & Reynarowych 2023, p. 324).
bvncnd(value, mu_x, mu_y, sd_x, sd_y, rho, given = c("y", "x"))