mle_pe: Find MLEs for microbiome mixture

In kbroman/mbmixture: Microbiome Mixture Analysis

Find MLEs for microbiome mixture

Description

Find joint MLEs of p and e for microbiome mixture model

Usage

mle_pe(
  tab,
  interval = c(0, 1),
  tol = 0.000001,
  check_boundary = FALSE,
  SE = FALSE
)

Arguments

tab	Dataset of read counts as 3d array of size 3x3x2, genotype in first sample x genotype in second sample x allele in read.
interval	Interval to which each parameter should be constrained
tol	Tolerance for convergence
check_boundary	If TRUE, explicitly check the boundaries of interval.
SE	If TRUE, get estimated standard errors.

Value

A vector containing the estimates of p and e along with the evaluated log likelihood and likelihood ratio test statistics for the hypotheses p=0 and p=1.

Examples

data(mbmixdata)
mle_pe(mbmixdata)

kbroman/mbmixture documentation built on Nov. 28, 2024, 6:38 p.m.