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

Value

NULL, a graph is plotted.

Author(s)

S. Lebre and G. Lelandais.

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)

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(ARTIVA)
Loading required package: MASS
Loading required package: igraph

Attaching package: 'igraph'

The following objects are masked from 'package:stats':

    decompose, spectrum

The following object is masked from 'package:base':

    union

Loading required package: gplots

Attaching package: 'gplots'

The following object is masked from 'package:stats':

    lowess

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ARTIVA/choosePriors.Rd_%03d_medium.png", width=480, height=480)
> ### Name: choosePriors
> ### Title: Function to plot an overview of possible priors for the number
> ###   of changepoints and/or edges.
> ### Aliases: choosePriors
> ### Keywords: util
> 
> ### ** Examples
> 
> # See some prior probabbility density when the maximal number of
> # CPs/parents is equal to 5
> data(priors)
> choosePriors(kmax=5,priors)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>