Summarize Attractors

Description

This function summarizes the posterior probability of possible attractors.

Usage

1
attractorSummary(tpost, post.prob.limit = 0.01, wildtype = TRUE)

Arguments

tpost

a ternaryPost object

post.prob.limit

the minimum posterior probability for an attractor to be listed

wildtype

if TRUE, the wildtype attractors are summarized; if FALSE, the perturbed attractors are summarized.

Value

The function returns a matrix of attractors and posterior probabilities for each perturbation.

Author(s)

Matthew N. McCall and Anthony Almudevar

See Also

Almudevar A, McCall MN, McMurray H, Land H (2011). Fitting Boolean Networks from Steady State Perturbation Data, Statistical Applications in Genetics and Molecular Biology, 10(1): Article 47.

Examples

1
2
3
4
5
6
ssObj <- matrix(c(1,1,1,0,1,1,0,0,1),nrow=3)
pObj <- matrix(c(1,0,0,0,1,0,0,0,1),nrow=3)
rownames(ssObj) <- rownames(pObj) <- colnames(ssObj) <- colnames(pObj) <- c("Gene1","Gene2","Gene3")
tnfitObj <- tnetfit(ssObj, pObj)
tnpostObj <- tnetpost(tnfitObj, mdelta=10, msample=10)
attractorSummary(tnpostObj)