Set of possible priors for the number of changepoints (CPs) or incoming edges.

Description

Set of 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

data(priors)

Format

A matrix of 96 rows by 44 columns (kmax, alpha, beta and the probability for k=0 to 40 according to the chosen values of kmax, alpha and beta).

References

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

See Also

choosePriors,ARTIVAsubnet, ARTIVAnet

Examples

# See some prior probability 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/priors.Rd_%03d_medium.png", width=480, height=480)
> ### Name: priors
> ### Title: Set of possible priors for the number of changepoints (CPs) or
> ###   incoming edges.
> ### Aliases: priors
> ### Keywords: datasets
> 
> ### ** Examples
> 
> # See some prior probability density when the maximal number of
> # CPs/parents is equal to 5
> data(priors)
> choosePriors(kmax=5,priors)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>