# 布尔网络

State space of a Boolean Network with N=4 nodes and K=1 links per node. Nodes can be either switched on (red) or off (blue). Thin (black) arrows symbolise the inputs of the Boolean function which is a simple "copy"-function for each node. The thick (grey) arrows show what a synchronous update does. Altogether there are 6 (orange) attractors, 4 of them are fixed points.

nodes and K=1 links per node. Nodes can be either switched on (red) or off (blue). Thin (black) arrows symbolise the inputs of the Boolean function which is a simple "copy"-function for each node. The thick (grey) arrows show what a synchronous update does. Altogether there are 6 (orange) attractors, 4 of them are fixed points.]]

A Boolean network consists of a discrete set of boolean variables each of which has a Boolean function (possibly different for each variable) assigned to it which takes inputs from a subset of those variables and output that determines the state of the variable it is assigned to. This set of functions in effect determines a topology (connectivity) on the set of variables, which then become nodes in a network. Usually, the dynamics of the system is taken as a discrete time series where the state of the entire network at time t+1 is determined by evaluating each variable's function on the state of the network at time t. This may be done synchronously or asynchronously.

A Boolean network consists of a discrete set of boolean variables each of which has a Boolean function (possibly different for each variable) assigned to it which takes inputs from a subset of those variables and output that determines the state of the variable it is assigned to. This set of functions in effect determines a topology (connectivity) on the set of variables, which then become nodes in a network. Usually, the dynamics of the system is taken as a discrete time series where the state of the entire network at time t+1 is determined by evaluating each variable's function on the state of the network at time t. This may be done synchronously or asynchronously.

Boolean networks have been used in biology to model regulatory networks. Although Boolean networks are a crude simplification of genetic reality where genes are not simple binary switches, there are several cases where they correctly capture the correct pattern of expressed and suppressed genes.

Boolean networks have been used in biology to model regulatory networks. Although Boolean networks are a crude simplification of genetic reality where genes are not simple binary switches, there are several cases where they correctly capture the correct pattern of expressed and suppressed genes.

The seemingly mathematical easy (synchronous) model was only fully understood in the mid 2000s.

## Classical model

A Boolean network is a particular kind of sequential dynamical system, where time and states are discrete, i.e. both the set of variables and the set of states in the time series each have a bijection onto an integer series. Such systems are like cellular automata on networks, except for the fact that when they are set up each node has a rule that is randomly chosen from all 2模板:Sup possible ones with K inputs. With K=2 class 2 behavior tends to dominate. But for K>2, the behavior one sees quickly approaches what is typical for a random mapping in which the network representing the evolution of the 2模板:Sup states of the N underlying nodes is itself connected essentially randomly.

A Boolean network is a particular kind of sequential dynamical system, where time and states are discrete, i.e. both the set of variables and the set of states in the time series each have a bijection onto an integer series. Such systems are like cellular automata on networks, except for the fact that when they are set up each node has a rule that is randomly chosen from all 2}} possible ones with K inputs. With K=2 class 2 behavior tends to dominate. But for K>2, the behavior one sees quickly approaches what is typical for a random mapping in which the network representing the evolution of the 2 states of the N underlying nodes is itself connected essentially randomly.

A random Boolean network (RBN) is one that is randomly selected from the set of all possible boolean networks of a particular size, N. One then can study statistically, how the expected properties of such networks depend on various statistical properties of the ensemble of all possible networks. For example, one may study how the RBN behavior changes as the average connectivity is changed.

A random Boolean network (RBN) is one that is randomly selected from the set of all possible boolean networks of a particular size, N. One then can study statistically, how the expected properties of such networks depend on various statistical properties of the ensemble of all possible networks. For example, one may study how the RBN behavior changes as the average connectivity is changed.

The first Boolean networks were proposed by Stuart A. Kauffman in 1969, as random models of genetic regulatory networks but their mathematical understanding only started in the 2000s.

The first Boolean networks were proposed by Stuart A. Kauffman in 1969, as random models of genetic regulatory networks but their mathematical understanding only started in the 2000s.

1969年，Stuart A. Kauffman提出了第一个布尔网络，作为遗传调控网络的随机模型，但其数学理解在2000年才开始。

## Attractors

Since a Boolean network has only 2N possible states, a trajectory will sooner or later reach a previously visited state, and thus, since the dynamics are deterministic, the trajectory will fall into a steady state or cycle called an attractor (though in the broader field of dynamical systems a cycle is only an attractor if perturbations from it lead back to it). If the attractor has only a single state it is called a point attractor, and if the attractor consists of more than one state it is called a cycle attractor. The set of states that lead to an attractor is called the basin of the attractor. States which occur only at the beginning of trajectories (no trajectories lead to them), are called garden-of-Eden states and the dynamics of the network flow from these states towards attractors. The time it takes to reach an attractor is called transient time.

Since a Boolean network has only 2N possible states, a trajectory will sooner or later reach a previously visited state, and thus, since the dynamics are deterministic, the trajectory will fall into a steady state or cycle called an attractor (though in the broader field of dynamical systems a cycle is only an attractor if perturbations from it lead back to it). If the attractor has only a single state it is called a point attractor, and if the attractor consists of more than one state it is called a cycle attractor. The set of states that lead to an attractor is called the basin of the attractor. States which occur only at the beginning of trajectories (no trajectories lead to them), are called garden-of-Eden states and the dynamics of the network flow from these states towards attractors. The time it takes to reach an attractor is called transient time.

garden-of-Eden states 网络的动态从这些状态流向吸引子。到达吸引子所需的时间称为瞬时 transient time 。


With growing computer power and increasing understanding of the seemingly simple model, different authors gave different estimates for the mean number and length of the attractors, here a brief summary of key publications.

With growing computer power and increasing understanding of the seemingly simple model, different authors gave different estimates for the mean number and length of the attractors, here a brief summary of key publications.

{ | class = “ wikitable sortable”
Author Author 作者 Year Year 年份 Mean attractor length Mean attractor length 平均吸引长度 Mean attractor number Mean attractor number 平均吸引子数 comment comment 评论
Kauffmann  Kauffmann

1969 1969 1969 $\displaystyle{ \langle A\rangle\sim \sqrt{N} }$ $\displaystyle{ \langle A\rangle\sim \sqrt{N} }$ < math > langle a rangle sim sqrt { n } </math > $\displaystyle{ \langle\nu\rangle\sim \sqrt{N} }$ $\displaystyle{ \langle\nu\rangle\sim \sqrt{N} }$ < math > langle nu rangle sim sqrt { n } </math >
Bastolla/ Parisi Bastolla/ Parisi Bastolla/ Parisi 1998 1998 1998 faster than a power law, $\displaystyle{ \langle A\rangle \gt N^x \forall x }$ faster than a power law, $\displaystyle{ \langle A\rangle \gt N^x \forall x }$ 比幂定律快，< math > langle a rangle > n ^ x for all x </math > faster than a power law, $\displaystyle{ \langle\nu\rangle \gt N^x \forall x }$ faster than a power law, $\displaystyle{ \langle\nu\rangle \gt N^x \forall x }$ 比幂定律快，< math > langle nu rangle > n ^ x for all x </math > first numerical evidences first numerical evidences

Bilke/ Sjunnesson Bilke/ Sjunnesson Bilke/Sjunnesson 2002 2002 2002 linear with system size, $\displaystyle{ \langle\nu\rangle \sim N }$ linear with system size, $\displaystyle{ \langle\nu\rangle \sim N }$ 与系统大小成线性关系
Socolar/Kauffman Socolar/Kauffman Socolar/Kauffman 2003 2003 2003 faster than linear, $\displaystyle{ \langle\nu\rangle \gt N^x }$ with $\displaystyle{ x \gt 1 }$ faster than linear, $\displaystyle{ \langle\nu\rangle \gt N^x }$ with $\displaystyle{ x \gt 1 }$

Samuelsson/Troein Samuelsson/Troein Samuelsson/Troein 2003 2003 2003 superpolynomial growth, $\displaystyle{ \langle\nu\rangle \gt N^x \forall x }$ superpolynomial growth, $\displaystyle{ \langle\nu\rangle \gt N^x \forall x }$ 超多项式生长，< math > langle nu rangle > n ^ x for all x </math > mathematical proof mathematical proof

Mihaljev/Drossel Mihaljev/Drossel Mihaljev/Drossel 2005 2005 2005 faster than a power law, $\displaystyle{ \langle A\rangle \gt N^x \forall x }$ faster than a power law, $\displaystyle{ \langle A\rangle \gt N^x \forall x }$ 比幂定律快，< math > langle a rangle > n ^ x for all x </math > faster than a power law, $\displaystyle{ \langle\nu\rangle \gt N^x \forall x }$ faster than a power law, $\displaystyle{ \langle\nu\rangle \gt N^x \forall x }$ 比幂定律快，< math > langle nu rangle > n ^ x for all x </math >

|}

## Stability

In dynamical systems theory, the structure and length of the attractors of a network corresponds to the dynamic phase of the network. The stability of Boolean networks depends on the connections of their nodes. A Boolean network can exhibit stable, critical or chaotic behavior. This phenomenon is governed by a critical value of the average number of connections of nodes ($\displaystyle{ K_{c} }$), and can be characterized by the Hamming distance as distance measure. In the unstable regime, the distance between two initially close states on average grows exponentially in time, while in the stable regime it decreases exponentially. In this, with "initially close states" one means that the Hamming distance is small compared with the number of nodes ($\displaystyle{ N }$) in the network.

In dynamical systems theory, the structure and length of the attractors of a network corresponds to the dynamic phase of the network. The stability of Boolean networks depends on the connections of their nodes. A Boolean network can exhibit stable, critical or chaotic behavior. This phenomenon is governed by a critical value of the average number of connections of nodes ($\displaystyle{ K_{c} }$), and can be characterized by the Hamming distance as distance measure. In the unstable regime, the distance between two initially close states on average grows exponentially in time, while in the stable regime it decreases exponentially. In this, with "initially close states" one means that the Hamming distance is small compared with the number of nodes ($\displaystyle{ N }$) in the network.

For N-K-model the network is stable if $\displaystyle{ K\lt K_{c} }$, critical if $\displaystyle{ K=K_{c} }$, and unstable if $\displaystyle{ K\gt K_{c} }$.

For N-K-model the network is stable if $\displaystyle{ K\lt K_{c} }$, critical if $\displaystyle{ K=K_{c} }$, and unstable if $\displaystyle{ K\gt K_{c} }$.

The state of a given node $\displaystyle{ n_{i} }$ is updated according to its truth table, whose outputs are randomly populated. $\displaystyle{ p_{i} }$ denotes the probability of assigning an off output to a given series of input signals.

The state of a given node $\displaystyle{ n_{i} }$ is updated according to its truth table, whose outputs are randomly populated. $\displaystyle{ p_{i} }$ denotes the probability of assigning an off output to a given series of input signals.

If $\displaystyle{ p_{i}=p=const. }$ for every node, the transition between the stable and chaotic range depends on $\displaystyle{ p }$. According to Bernard Derrida and Yves Pomeau

If $\displaystyle{ p_{i}=p=const. }$ for every node, the transition between the stable and chaotic range depends on $\displaystyle{ p }$. According to Bernard Derrida and Yves Pomeau

, the critical value of the average number of connections is $\displaystyle{ K_{c}=1/[2p(1-p)] }$.

， 平均连接数的临界值为 $\displaystyle{ K_{c}=1/[2p(1-p)] }$

If $\displaystyle{ K }$ is not constant, and there is no correlation between the in-degrees and out-degrees, the conditions of stability is determined by $\displaystyle{ \langle K^{in}\rangle }$ The network is stable if $\displaystyle{ \langle K^{in}\rangle \lt K_{c} }$, critical if $\displaystyle{ \langle K^{in}\rangle =K_{c} }$, and unstable if $\displaystyle{ \langle K^{in}\rangle \gt K_{c} }$.

If $\displaystyle{ K }$ is not constant, and there is no correlation between the in-degrees and out-degrees, the conditions of stability is determined by $\displaystyle{ \langle K^{in}\rangle }$ The network is stable if $\displaystyle{ \langle K^{in}\rangle \lt K_{c} }$, critical if $\displaystyle{ \langle K^{in}\rangle =K_{c} }$, and unstable if $\displaystyle{ \langle K^{in}\rangle \gt K_{c} }$.

The conditions of stability are the same in the case of networks with scale-free topology where the in-and out-degree distribution is a power-law distribution: $\displaystyle{ P(K) \propto K^{-\gamma} }$, and $\displaystyle{ \langle K^{in} \rangle=\langle K^{out} \rangle }$, since every out-link from a node is an in-link to another.

The conditions of stability are the same in the case of networks with scale-free topology where the in-and out-degree distribution is a power-law distribution: $\displaystyle{ P(K) \propto K^{-\gamma} }$, and $\displaystyle{ \langle K^{in} \rangle=\langle K^{out} \rangle }$, since every out-link from a node is an in-link to another.

Sensitivity shows the probability that the output of the Boolean function of a given node changes if its input changes. For random Boolean networks,

Sensitivity shows the probability that the output of the Boolean function of a given node changes if its input changes. For random Boolean networks,

$\displaystyle{ q_{i}=2p_{i}(1-p_{i}) }$. In the general case, stability of the network is governed by the largest eigenvalue $\displaystyle{ \lambda_{Q} }$ of matrix $\displaystyle{ Q }$, where $\displaystyle{ Q_{ij}=q_{i}A_{ij} }$, and $\displaystyle{ A }$ is the adjacency matrix of the network. The network is stable if $\displaystyle{ \lambda_{Q}\lt 1 }$, critical if $\displaystyle{ \lambda_{Q}=1 }$, unstable if $\displaystyle{ \lambda_{Q}\gt 1 }$.

$\displaystyle{ q_{i}=2p_{i}(1-p_{i}) }$。在一般情况下，网络的稳定性由最大的特征值$\displaystyle{ \lambda_{Q} }$来控制。的矩阵 $\displaystyle{ Q }$，其中$\displaystyle{ Q_{ij}=q_{i}A_{ij} }$$\displaystyle{ A }$ 是网络的邻接矩阵。如果 $\displaystyle{ \lambda_{Q}\lt 1 }$，网络是稳定的；如果 $\displaystyle{ \lambda_{Q}=1 }$，网络是临界的；如果 $\displaystyle{ \lambda_{Q}\gt 1 }$，网络是不稳定的。

## Other topologies

One theme is to study different underlying graph topologies.

One theme is to study different underlying graph topologies.

• The homogeneous case simply refers to a grid which is simply the reduction to the famous Ising model.
• 同质情况 Homogeneous Case 只是指网格，它只是对著名的伊辛模型 lsing model 的还原。
• Scale-free topologies may be chosen for Boolean networks. One can distinguish the case where only in-degree distribution in power-law distributed, or only the out-degree-distribution or both.

## Other updating schemes

Classical Boolean networks (sometimes called CRBN, i.e. Classic Random Boolean Network) are synchronously updated. Motivated by the fact that genes don't usually change their state simultaneously, different alternatives have been introduced. A common classification is the following:

Classical Boolean networks (sometimes called CRBN, i.e. Classic Random Boolean Network) are synchronously updated. Motivated by the fact that genes don't usually change their state simultaneously, different alternatives have been introduced. A common classification is the following:

• Deterministic asynchronous updated Boolean networks (DRBNs) are not synchronously updated but a deterministic solution still exists. A node i will be updated when t ≡ Qi (mod Pi) where t is the time step.

t≡Qi(modPi) 其中 t 是时间步长时，i 节点将被更新。确定性异步更新布尔网络(DRBNs)不是同步更新，但确定性解仍然存在。当 t≡Qi(modPi) 时，i 节点将被更新，其中 t 是时间步长。

• The most general case is full stochastic updating (GARBN, general asynchronous random boolean networks). Here, one (or more) node(s) are selected at each computational step to be updated.
• 最一般的情况是完全随机更新（GARBN，一般异步随机布尔网络）。在这里，在每个计算步骤中选择一个（或多个）节点进行更新。
• The Partially-Observed Boolean Dynamical System (POBDS) signal model differs from all previous deterministic and stochastic Boolean network models by removing the assumption of direct observability of the Boolean state vector and allowing uncertainty in the observation process, addressing the scenario encountered in practice.

## Classification

• The Scalable Optimal Bayesian Classification developed an optimal classification of trajectories accounting for potential model uncertainty and also proposed a particle-based trajectory classification that is highly scalable for large networks with much lower complexity than the optimal solution.