gaga (version 2.18.0)

getpar: Extract hyper-parameter estimates from a gagafit or nnfit object

Description

Extracts the hyper-parameter estimates from a gagafit or nnfit object and puts them in a list.

Usage

getpar(fit)

Arguments

fit

Object of class gagafit or nnfit.

Value

For gagafit objects, a list with components:

a0

Estimated value of hyper-parameter a0

nu

Estimated value of hyper-parameter nu

balpha

Estimated value of hyper-parameter balpha

nualpha

Estimated value of hyper-parameter nualpha

probclus

Estimated cluster probabilities

probpat

Estimated prior probability of each expression pattern

For nnfit objects, a vector with elements mu0, tau02, v0, sigma02, probclus and probpat. These are the hierarchical N(mu0,tau0^2) * IG(.5*v0,.5*v0*sigma0^2) prior parameter estimates.

Details

This routine simply evaluates the component parest from a gagafit or nnfit object, which causes an error if this component is not available. This routine is used internally by a number of other routines.

References

Rossell D. (2009) GaGa: a Parsimonious and Flexible Model for Differential Expression Analysis. Annals of Applied Statistics, 3, 1035-1051.

Yuan, M. and Kendziorski, C. (2006). A unified approach for simultaneous gene clustering and differential expression identification. Biometrics 62(4): 1089-1098.

See Also

fitGG, fitNN, parest