|
Matlab function expressing the boolean network. |
||
| function new = rules toy model(setnow, lnow, nodes) |
||
| new = zeros(lnow,1); |
||
| Lig1 = 1; |
Lig2 = 2; |
A1 = 3; |
| A2 = 4; |
Inputk3r = 5; |
K1 = 6; |
| K2 = 7; |
K3p1 = 8; |
K3p2 = 9; |
| K3r = 10; |
Ph1 = 11; |
Ph2 = 12; |
| R1 = 13; |
R2 = 14; |
Geneexp = 15; |
| Tf1 = 16; |
||
| for i = 1:lnow |
||
| inode = setnow(i) ; |
||
| switch inode |
||
| case A1 |
new(i) = nodes(R1); |
|
| case A2 |
new(i) = and((nodes(R1) == 0), nodes(R2)); |
|
| case K1 |
new(i) = and((nodes(Ph1) == 0), nodes(A1)); |
|
| case K2 |
new(i) = or(and(nodes(A2), nodes(Ph2)), and(nodes(K1), nodes(A2))); |
|
| case K3p1 |
new(i) = and(nodes(K1), nodes(K3r)); |
|
| case K3p2 |
new(i) = and(nodes(K3r), nodes(K2)); |
|
| case K3r |
new(i) = nodes(Inputk3r); |
|
| case Ph1 |
new(i) = nodes(K3p1); |
|
| case Ph2 |
new(i) = nodes(A1); |
|
| case R1 |
new(i) = nodes(Lig1); |
|
| case R2 |
new(i) = nodes(Lig2); |
|
| case Geneexp |
new(i) = nodes(Tf1); |
|
| case Tf1 |
new(i) = or(or(and(nodes(K3p1), nodes(K3p2)), nodes(K3p2)), nodes(K3p1)); |
|
| end |
||
| end |
||
|
Alternative export of the toy model for simulations in Matlab, e.g. as used in [2]. | ||
Saez-Rodriguez et al. BMC Bioinformatics 2006 7:506 doi:10.1186/1471-2105-7-506 |
||