mle: Compute MLE based on R0 and R

In stephaneguerrier/CPreval: TO DO

Description

Proportion estimated using the MLE and confidence intervals based the asymptotic distribution of the estimator as well as the Clopper–Pearson approach.

Usage

mle(
  R0,
  R,
  pi0,
  n,
  alpha0 = 0,
  alpha = 0,
  beta0 = 0,
  beta = 0,
  gamma = 0.05,
  simulation = FALSE,
  ...
)

Arguments

R0	A numeric that provides the people of positive people in the sample that were known to be positive.
R	A numeric that provides the people of positive people in the sample.
n	A numeric that provides the sample size.
alpha0	A numeric that provides the False Negative (FN) rate for the sample R0. Default value is 0.
alpha	A numeric that provides the False Negative (FN) rate for the sample R. Default value is 0.
beta0	A numeric that provides the False Positive (FP) rate for the sample R0. Default value is 0.
beta	A numeric that provides the False Positive (FP) rate for the sample R. Default value is 0.
gamma	A numeric that used to compute a (1 - gamma) confidence region for the proportion. Default value is 0.05.
...	Additional arguments.

Value

A CPreval object with the structure:

• estimate: Estimated proportion

• sd: Estimated standard error of the estimator

• ci_asym: Asymptotic confidence interval

• ci_cp: Confidence interval based on the Clopper–Pearson approach

• gamma: Confidence level (i.e. 1 - gamma) for confidence interval

• method: Estimation method (in this case sample survey)

• measurement: A vector with (alpha0, alpha, beta0, beta)

• ...: Additional parameters

Author(s)

Stephane Guerrier

Examples

# Samples without measurement error
X = sim_Rs(p = 3/100, pi0 = 1/100, n = 1500, seed = 18)
mle(R0 = X$R0, R = X$R, pi0 = X$pi0, n = X$n)

# With measurement error
X = sim_Rs(p = 30/1000, pi0 = 10/1000, n = 1500, alpha0 = 0.01,
alpha = 0.01, beta0 = 0.05, beta = 0.05, seed = 18)
mle(R0 = X$R0, R = X$R, pi0 = X$pi0, n = X$n)
mle(R0 = X$R0, R = X$R, pi0 = X$pi0, n = X$n, alpha0 = 0.01,
alpha = 0.01, beta0 = 0.05, beta = 0.05)

stephaneguerrier/CPreval documentation built on June 6, 2020, 2:28 a.m.