logistic.mles {genpwr}

Function to calculate MLE's for logistic models

Description

Finds the maximum likelihood estimates for a given 2x3 table under the specified genetic model.

Usage

logistic.mles(t, model)

Arguments

t

A 2x3 table of the joint probabilities of disease and genotype. Rows are case vs. control and columns are genotypes.

model

The assumed genetic model(s) used in testing: 'Dominant', 'Additive', 'Recessive'

Value

A vector of logistic regression model coefficients.

Examples

logistic.mles(rbind(c(0.2967437, 0.1806723, 0.02258404), 
 c(0.3432563, 0.1393277, 0.01741596)), "Dominant")