MLE {PoDBAY}

Maximum Likelihood estimation

Description

Function calculates the log likelihood value which is used after the initial guesses of the parameters are set in the PoDMLE function.

Usage

MLE(params,
    nondiseasedTiters,
    diseasedTiters,
    adjustTiters = FALSE,
    adjustFrom = log2(10),
    adjustTo = log2(5))

Arguments

params

named numeric vector: PoD curve parameters (et50, slope, pmax)

nondiseasedTiters

numeric vector: non-diseased subjects titers

diseasedTiters

numeric vector: diseased subjects titers

adjustTiters

boolean: set to TRUE if titer values should be adjusted, for details see PoD function

adjustFrom

numeric: value specifying the detection limit, all values below the detection limit will be adjusted to adjustTo value

adjustTo

numeric: value to which titers below the detection limit will be adjusted

Details

MLE function is used inside of PoDMLE function and esimates the PoD curve parameters.

Based on the provided titers for diseased and non-diseased subjects the PoD curve parameters which maximize the log likelihood are chosen as optimal estimates of parameters.

Value

log likelihood, numeric value

Examples

# Data preparation
data(diseased)
data(nondiseased)
data(PoDParams)

# MLE calculation
MLE(PoDParams, nondiseased$titers, diseased$titers)