chooseNetwork {BoolNet}

Extract a single Boolean network from a probabilistic Boolean network

Description

Creates a BooleanNetwork object with exactly one function per gene by extracting a specified set of transition functions from a ProbabilisticBooleanNetwork or BooleanNetworkCollection object.

Usage

chooseNetwork(probabilisticNetwork, 
              functionIndices, 
              dontCareValues=NULL, 
              readableFunctions=FALSE)

Arguments

probabilisticNetwork

A ProbabilisticBooleanNetwork or BooleanNetworkCollection object as returned by reconstructNetwork or loadNetwork

functionIndices

A vector of function indices with one entry for each gene

dontCareValues

If probabilisticNetwork is of class BooleanNetworkCollection, this specifies the values to fill in for "don't care" (*) values in the truth tables of the transition functions. This is a list containing one vector of Boolean values for each gene. The lengths of the vectors must coincide with the numbers of "don't care" values in the functions.

readableFunctions

If probabilisticNetwork is of class BooleanNetworkCollection, the string representations of the transition functions must be refreshed after filling in values for the "don't care" entries. This parameter specifies if readable DNF representations of the transition function truth tables are generated and displayed when the network is printed. If set to FALSE, the truth table result column is displayed. If set to "canonical", a canonical Disjunctive Normal Form is generated from each truth table. If set to "short", the canonical DNF is minimized by joining terms (which can be time-consuming for functions with many inputs). If set to TRUE, a short DNF is generated for functions with up to 12 inputs, and a canonical DNF is generated for functions with more than 12 inputs.

Value

Returns an object of class BooleanNetwork consisting of the transition functions whose indices were specified in functionIndices. The class BooleanNetwork is described in more detail in loadNetwork.

Constant genes are automatically fixed (e.g. knocked-out or over-expressed). This means that they are always set to the constant value, and states with the complementary value are not considered in transition tables etc. If you would like to change this behaviour, use fixGenes to reset the fixing.

See Also

reconstructNetwork, loadNetwork

Examples

  ## Not run: 
# load example data
data(examplePBN)

# extract a unique network
# - always use the first function
net <- chooseNetwork(examplePBN, rep(1, length(examplePBN$genes)))

# get attractors from this network
print(getAttractors(net))

## End(Not run)