em.draw {metaArray}R Documentation

Plot of transformed expression produced by EM algorithm

Description

Given a numeric vector, a plot of four panels is drawn: 1) fitted mixture distribution 2) transformed expression against original expression 3) histogram of original expression 4) progression of log-likelihood during the fit

Usage

  em.draw(vec, cl, threshold=0.0001)
  #em.draw(vec, cl=1-metastasis, threshold=0.0001)

Arguments

vec A numeric vector, especially a particular row of expression matrix
cl A vector of 0s and 1s. Use 1 for normal phenotype and 0 for non-normal phenotype. If left blank, all samples will be labeled as normal phenotype. Normal component of mixture is estimated using samples with normal phenotype only. POE for samples with non-normal phenotype will be calculated after EM algorithm finishes ML estimation.
threshold Crieterion for convergence in likelihood.

Value

A plot of four panels will appear upon the call.

Author(s)

Debashis Ghosh <ghoshd@umich.edu> , Hyungwon Choi <hwchoi@umich.edu>

Examples



[Package metaArray version 1.12.0 Index]