get_tresholding_estimator_bias: Get bias of thresholding estimator

In timothy-barry/glmeiv: Implements the GLM-EIV model and method

Get bias of thresholding estimator

Description

Assumes without loss of generality that m_intercept and g_intercept are 0. Also, assumes the Bayes-optimal threshold has been chosen.

Usage

get_tresholding_estimator_bias(
  m_perturbation,
  g_perturbation,
  g_intercept,
  pi,
  c,
  return_bias = TRUE
)

Arguments

m_perturbation	perturbation coefficient for mRNA model
g_perturbation	perturbation coefficient for gRNA model
g_intercept	intercept term in gRNA model
pi	probability of perturbation
c	the threshold
return_bias	if TRUE, returns (absolute) bias, defined as ; if FALSE, returns attenuation fraction

Details

Assumes errors are N(0,1).

Value

the bias of the thresholding estimator

Examples

m_perturbation <- 1
g_perturbation <- 1
g_intercept <- 0.5
pi <- 0.5
c <- 1
get_tresholding_estimator_bias(m_perturbation, g_perturbation, g_intercept, pi, c)

timothy-barry/glmeiv documentation built on Jan. 30, 2024, 3:46 p.m.