emfit {EBarrays}R Documentation

Implements EM algorithm for gene expression mixture model

Description

Implements the EM algorithm for gene expression mixture model

Usage

emfit(data,
      family,
      hypotheses,
      theta.init = NULL, p.init = NULL,
      num.iter = 20,
      verbose = getOption("verbose"),
      trace = TRUE,
      optim.control = list())

Arguments

data a matrix or an an object of class "exprSet" or "ExpressionSet"
family an object of class ``ebarraysFamily'' or a characterstring which can be coerced to one. Currently, only the characters "GG" and "LNN" are valid. Other families can be supplied by constructing them explicitly.
hypotheses an object of class ``ebarraysPatterns'' representing the hypotheses of interest. Such patterns can be generated by the function ebPatterns
theta.init initial values of theta. Should be appropriate for the corresponding family. If omitted, defaults are calculated based on the family.
p.init initial values of mixture probabilities. If omitted, all hypotheses are initially taken as equally likely.
num.iter number of EM iterations
verbose logical or numeric (0,1,2) indicating desired level of information printed for the user
trace whether the result should have slots thetaTrace and probTrace, matrices giving values for all iterations
optim.control list passed unchanged to optim for finer control

Value

an object of formal class ``ebarraysEmFit'' with the following slots

family the family passed to emfit
hypotheses the hypotheses passed to emfit
thetaEst numeric vector holding the parameter estimates
probEst numeric vector holding estimates of the mixture probabilities
thetaTrace optional matrix holding intermediate values of thetaEst
probTrace optional matrix holding intermediate values of probEst


The show() method can be used to summarize the object. It can also be used to generate posterior probabilities using postprob

Note

emfit is a generic function dispatched on the arguments data, family and hypotheses.

Author(s)

Christina Kendziorski, Michael Newton and Deepayan Sarkar

References

UW Madison Biostatistics and Medical Informatics Technical Report 166, ``On parametric empirical Bayes methods for comparing multiple groups using replicated gene expression profiles,'' by C.M. Kendziorski, M.A. Newton, H. Lan, and M.N. Gould. See http://www.biostat.wisc.edu/~kendzior/

See Also

ebPatterns, ebarraysFamily-class

Examples

data(sample.exprSet.1) ## from Biobase
patterns <- ebPatterns(c("1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1",
                         "1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2"))
gg.fit <- emfit(data = sample.exprSet.1, family = "GG", hypotheses = patterns,
                verbose = TRUE)
##                theta.init = c(10, 1, 1), p.init = c(.99, .01))
show(gg.fit)

[Package EBarrays version 1.6.0 Index]