ARTIVA (version 1.2.3)

choosePriors: Function to plot an overview of possible priors for the number of changepoints and/or edges.

Description

Plots an overview of some possible priors used in the ARTIVAsubnet function for the number of changepoints (resp. incoming edges) according to a given number of maximum changepoints maxCP (resp. incoming Edges maxPred) when parameters (alphaCP, betaCP for the CPs or alphaEdges, betaEdges for the edges) in function ARTIVAsubnet are set to default (alpha=1, beta=0.5). In the ARTIVAsubnet procedure, the number of CPs (respectively the number of incoming edges) is sampled from a truncated Poisson with mean lambda, where lambda is drawn from an Inverse Gamma distribution (alpha, beta), see Lebre et al. (2010) for more details.

Usage

choosePriors(kmax,priors)

Arguments

kmax

Maximum number of changepoints or incomins Edges (parents)

priors

Table describing the priors which can be loaded with data(priors)

Value

NULL, a graph is plotted.

References

S. Lebre, J. Becq, F. Devaux, M. P. H. Stumpf, G. Lelandais (2010) Statistical inference of the time-varying structure of gene-regulation networks, BMC Systems Biology, 4:130.

See Also

ARTIVAsubnet, ARTIVAnet

Examples

# See some prior probabbility density when the maximal number of
# CPs/parents is equal to 5
data(priors)
choosePriors(kmax=5,priors)