-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathODE_MacrophagePolarizationProbFinalComplete.R
69 lines (65 loc) · 12.7 KB
/
ODE_MacrophagePolarizationProbFinalComplete.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
MacrophagePolarizationProbFinalCompletefunction(t, state, parameters) {
with(as.list(c(state, parameters)),{
# Input Nodes
w_STAT3 = ((((((((IL10)+(A2a)-(IL10)*(A2a))+(IL10e)-((IL10)+(A2a)-(IL10)*(A2a))*(IL10e))+((((IL6e)*(IL4e))*(1-IFNG)))-(((IL10)+(A2a)-(IL10)*(A2a))+(IL10e)-((IL10)+(A2a)-(IL10)*(A2a))*(IL10e))*((((IL6e)*(IL4e))*(1-IFNG))))+(MCSF)-((((IL10)+(A2a)-(IL10)*(A2a))+(IL10e)-((IL10)+(A2a)-(IL10)*(A2a))*(IL10e))+((((IL6e)*(IL4e))*(1-IFNG)))-(((IL10)+(A2a)-(IL10)*(A2a))+(IL10e)-((IL10)+(A2a)-(IL10)*(A2a))*(IL10e))*((((IL6e)*(IL4e))*(1-IFNG))))*(MCSF)))*(1-((((((SOCS3)+(IFNG)-(SOCS3)*(IFNG))+(IFNB)-((SOCS3)+(IFNG)-(SOCS3)*(IFNG))*(IFNB))+(NFKB)-(((SOCS3)+(IFNG)-(SOCS3)*(IFNG))+(IFNB)-((SOCS3)+(IFNG)-(SOCS3)*(IFNG))*(IFNB))*(NFKB))+(GMCSF)-((((SOCS3)+(IFNG)-(SOCS3)*(IFNG))+(IFNB)-((SOCS3)+(IFNG)-(SOCS3)*(IFNG))*(IFNB))+(NFKB)-(((SOCS3)+(IFNG)-(SOCS3)*(IFNG))+(IFNB)-((SOCS3)+(IFNG)-(SOCS3)*(IFNG))*(IFNB))*(NFKB))*(GMCSF))+(IgG)-(((((SOCS3)+(IFNG)-(SOCS3)*(IFNG))+(IFNB)-((SOCS3)+(IFNG)-(SOCS3)*(IFNG))*(IFNB))+(NFKB)-(((SOCS3)+(IFNG)-(SOCS3)*(IFNG))+(IFNB)-((SOCS3)+(IFNG)-(SOCS3)*(IFNG))*(IFNB))*(NFKB))+(GMCSF)-((((SOCS3)+(IFNG)-(SOCS3)*(IFNG))+(IFNB)-((SOCS3)+(IFNG)-(SOCS3)*(IFNG))*(IFNB))+(NFKB)-(((SOCS3)+(IFNG)-(SOCS3)*(IFNG))+(IFNB)-((SOCS3)+(IFNG)-(SOCS3)*(IFNG))*(IFNB))*(NFKB))*(GMCSF))*(IgG)))))
w_SMAD23 = (((((TGFBe)+(TGFB)-(TGFBe)*(TGFB)))*(1-((IFNG)+(SOCS3)-(IFNG)*(SOCS3)))))
w_STAT6 = ((((IL4e)+(SMAD23)-(IL4e)*(SMAD23))+(HIF1A)-((IL4e)+(SMAD23)-(IL4e)*(SMAD23))*(HIF1A)))*(1-(((STAT1)+(TNFA)-(STAT1)*(TNFA))+(TNFAe)-((STAT1)+(TNFA)-(STAT1)*(TNFA))*(TNFAe)))
w_SOCS1 = (((STAT6)+(IgG)-(STAT6)*(IgG)))
w_TLR4 = (((TLR4)*(1-((IgG)+(GCGCR)-(IgG)*(GCGCR)))))
w_NFKB = ((((((TLR4)+(GMCSF)-(TLR4)*(GMCSF))+(TNFA)-((TLR4)+(GMCSF)-(TLR4)*(GMCSF))*(TNFA))+(TNFAe)-(((TLR4)+(GMCSF)-(TLR4)*(GMCSF))+(TNFA)-((TLR4)+(GMCSF)-(TLR4)*(GMCSF))*(TNFA))*(TNFAe))+(IL1B)-((((TLR4)+(GMCSF)-(TLR4)*(GMCSF))+(TNFA)-((TLR4)+(GMCSF)-(TLR4)*(GMCSF))*(TNFA))+(TNFAe)-(((TLR4)+(GMCSF)-(TLR4)*(GMCSF))+(TNFA)-((TLR4)+(GMCSF)-(TLR4)*(GMCSF))*(TNFA))*(TNFAe))*(IL1B)))*(1-((((((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))+(IgG)-((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))*(IgG))+(GCGCR)-(((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))+(IgG)-((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))*(IgG))*(GCGCR))+(HIF1A)-((((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))+(IgG)-((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))*(IgG))+(GCGCR)-(((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))+(IgG)-((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))*(IgG))*(GCGCR))*(HIF1A))+(MCSF)-(((((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))+(IgG)-((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))*(IgG))+(GCGCR)-(((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))+(IgG)-((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))*(IgG))*(GCGCR))+(HIF1A)-((((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))+(IgG)-((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))*(IgG))+(GCGCR)-(((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))+(IgG)-((((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))+(STAT3)-(((IL10)+(IL10e)-(IL10)*(IL10e))+(STAT6)-((IL10)+(IL10e)-(IL10)*(IL10e))*(STAT6))*(STAT3))*(IgG))*(GCGCR))*(HIF1A))*(MCSF)))
w_STAT1 = ((((IFNB)+(IFNG)-(IFNB)*(IFNG))+(((IL6e)*(IFNG)))-((IFNB)+(IFNG)-(IFNB)*(IFNG))*(((IL6e)*(IFNG)))))*(1-((((STAT3)+(SOCS1)-(STAT3)*(SOCS1))+(STAT6)-((STAT3)+(SOCS1)-(STAT3)*(SOCS1))*(STAT6))+(IgG)-(((STAT3)+(SOCS1)-(STAT3)*(SOCS1))+(STAT6)-((STAT3)+(SOCS1)-(STAT3)*(SOCS1))*(STAT6))*(IgG)))
w_SOCS3 = (STAT1)
w_AP1 = ((((((((GMCSF)+(TNFAe)-(GMCSF)*(TNFAe))+(TNFA)-((GMCSF)+(TNFAe)-(GMCSF)*(TNFAe))*(TNFA))+(IL6e)-(((GMCSF)+(TNFAe)-(GMCSF)*(TNFAe))+(TNFA)-((GMCSF)+(TNFAe)-(GMCSF)*(TNFAe))*(TNFA))*(IL6e))+(IL1B)-((((GMCSF)+(TNFAe)-(GMCSF)*(TNFAe))+(TNFA)-((GMCSF)+(TNFAe)-(GMCSF)*(TNFAe))*(TNFA))+(IL6e)-(((GMCSF)+(TNFAe)-(GMCSF)*(TNFAe))+(TNFA)-((GMCSF)+(TNFAe)-(GMCSF)*(TNFAe))*(TNFA))*(IL6e))*(IL1B)))*(1-((GCGCR)+(STAT6)-(GCGCR)*(STAT6)))))
w_IL10 = (((((((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))+(((STAT3)*(IL6e)))-((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))*(((STAT3)*(IL6e))))+(((TLR4)*(A2a)))-(((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))+(((STAT3)*(IL6e)))-((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))*(((STAT3)*(IL6e))))*(((TLR4)*(A2a))))+(MCSF)-((((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))+(((STAT3)*(IL6e)))-((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))*(((STAT3)*(IL6e))))+(((TLR4)*(A2a)))-(((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))+(((STAT3)*(IL6e)))-((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))*(((STAT3)*(IL6e))))*(((TLR4)*(A2a))))*(MCSF))+(ERK)-(((((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))+(((STAT3)*(IL6e)))-((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))*(((STAT3)*(IL6e))))+(((TLR4)*(A2a)))-(((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))+(((STAT3)*(IL6e)))-((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))*(((STAT3)*(IL6e))))*(((TLR4)*(A2a))))+(MCSF)-((((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))+(((STAT3)*(IL6e)))-((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))*(((STAT3)*(IL6e))))+(((TLR4)*(A2a)))-(((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))+(((STAT3)*(IL6e)))-((IL10e)+(((IL4e)*(STAT6)))-(IL10e)*(((IL4e)*(STAT6))))*(((STAT3)*(IL6e))))*(((TLR4)*(A2a))))*(MCSF))*(ERK)))*(1-(((((((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))+(IFNG)-(((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))*(IFNG))+(IFNB)-((((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))+(IFNG)-(((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))*(IFNG))*(IFNB))+(GMCSF)-(((((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))+(IFNG)-(((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))*(IFNG))+(IFNB)-((((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))+(IFNG)-(((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))*(IFNG))*(IFNB))*(GMCSF))+(IgG)-((((((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))+(IFNG)-(((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))*(IFNG))+(IFNB)-((((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))+(IFNG)-(((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))*(IFNG))*(IFNB))+(GMCSF)-(((((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))+(IFNG)-(((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))*(IFNG))+(IFNB)-((((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))+(IFNG)-(((STAT1)+(NFKB)-(STAT1)*(NFKB))+(SOCS3)-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(SOCS3))*(IFNG))*(IFNB))*(GMCSF))*(IgG)))
w_TNFA = ((((((TNFAe)+(AP1)-(TNFAe)*(AP1))+(((IL6e)*(AP1)))-((TNFAe)+(AP1)-(TNFAe)*(AP1))*(((IL6e)*(AP1))))+(((Fra1)*(AP1)))-(((TNFAe)+(AP1)-(TNFAe)*(AP1))+(((IL6e)*(AP1)))-((TNFAe)+(AP1)-(TNFAe)*(AP1))*(((IL6e)*(AP1))))*(((Fra1)*(AP1))))*(1-(((STAT1)+(NFKB)-(STAT1)*(NFKB))+(((TLR4)*(A2a)))-((STAT1)+(NFKB)-(STAT1)*(NFKB))*(((TLR4)*(A2a)))))))
w_TGFB = ((((((TGFBe)+(((IL4e)*(STAT6)))-(TGFBe)*(((IL4e)*(STAT6))))+(((IL10e)*(STAT3)))-((TGFBe)+(((IL4e)*(STAT6)))-(TGFBe)*(((IL4e)*(STAT6))))*(((IL10e)*(STAT3))))+(((IL10)*(STAT3)))-(((TGFBe)+(((IL4e)*(STAT6)))-(TGFBe)*(((IL4e)*(STAT6))))+(((IL10e)*(STAT3)))-((TGFBe)+(((IL4e)*(STAT6)))-(TGFBe)*(((IL4e)*(STAT6))))*(((IL10e)*(STAT3))))*(((IL10)*(STAT3))))+(A2a)-((((TGFBe)+(((IL4e)*(STAT6)))-(TGFBe)*(((IL4e)*(STAT6))))+(((IL10e)*(STAT3)))-((TGFBe)+(((IL4e)*(STAT6)))-(TGFBe)*(((IL4e)*(STAT6))))*(((IL10e)*(STAT3))))+(((IL10)*(STAT3)))-(((TGFBe)+(((IL4e)*(STAT6)))-(TGFBe)*(((IL4e)*(STAT6))))+(((IL10e)*(STAT3)))-((TGFBe)+(((IL4e)*(STAT6)))-(TGFBe)*(((IL4e)*(STAT6))))*(((IL10e)*(STAT3))))*(((IL10)*(STAT3))))*(A2a)))*(1-((NFKB)+(STAT1)-(NFKB)*(STAT1)))
w_Fra1 = (Fra1)+(TLR4)-(Fra1)*(TLR4)
w_HIF1A = (((Hipoxia)+(TGFB)-(Hipoxia)*(TGFB))+(TGFBe)-((Hipoxia)+(TGFB)-(Hipoxia)*(TGFB))*(TGFBe))+(((TLR4)*(A2a)))-(((Hipoxia)+(TGFB)-(Hipoxia)*(TGFB))+(TGFBe)-((Hipoxia)+(TGFB)-(Hipoxia)*(TGFB))*(TGFBe))*(((TLR4)*(A2a)))
w_ERK = (IgG)
w_IL4e = (IL4e)
w_IL10e = ((IL10e)+(STAT3)-(IL10e)*(STAT3))+(STAT6)-((IL10e)+(STAT3)-(IL10e)*(STAT3))*(STAT6)
w_TNFAe = ((TNFAe)+(AP1)-(TNFAe)*(AP1))+(NFKB)-((TNFAe)+(AP1)-(TNFAe)*(AP1))*(NFKB)
w_TGFBe = ((TGFBe)+(STAT3)-(TGFBe)*(STAT3))+(STAT6)-((TGFBe)+(STAT3)-(TGFBe)*(STAT3))*(STAT6)
w_IgG = (IgG)
w_A2a = (A2a)
w_IFNG = (IFNG)+(STAT1)-(IFNG)*(STAT1)
w_IFNB = (IFNB)+(STAT1)-(IFNB)*(STAT1)
w_IL1B = ((IL1B)+(AP1)-(IL1B)*(AP1))+(NFKB)-((IL1B)+(AP1)-(IL1B)*(AP1))*(NFKB)
w_IL6e = ((IL6e)+(Fra1)-(IL6e)*(Fra1))+(AP1)-((IL6e)+(Fra1)-(IL6e)*(Fra1))*(AP1)
w_MCSF = (MCSF)
w_GMCSF = (GMCSF)
w_Hipoxia = (Hipoxia)
w_GCGCR = (GCGCR)
# Rates of Change
dSTAT3 = 1/(1+(exp(-h*(w_STAT3-b)))) - (alphaSTAT3*STAT3)
dSMAD23 = 1/(1+(exp(-h*(w_SMAD23-b)))) - (alphaSMAD23*SMAD23)
dSTAT6 = 1/(1+(exp(-h*(w_STAT6-b)))) - (alphaSTAT6*STAT6)
dSOCS1 = 1/(1+(exp(-h*(w_SOCS1-b)))) - (alphaSOCS1*SOCS1)
dTLR4 = 1/(1+(exp(-h*(w_TLR4-b)))) - (alphaTLR4*TLR4)
dNFKB = 1/(1+(exp(-h*(w_NFKB-b)))) - (alphaNFKB*NFKB)
dSTAT1 = 1/(1+(exp(-h*(w_STAT1-b)))) - (alphaSTAT1*STAT1)
dSOCS3 = 1/(1+(exp(-h*(w_SOCS3-b)))) - (alphaSOCS3*SOCS3)
dAP1 = 1/(1+(exp(-h*(w_AP1-b)))) - (alphaAP1*AP1)
dIL10 = 1/(1+(exp(-h*(w_IL10-b)))) - (alphaIL10*IL10)
dTNFA = 1/(1+(exp(-h*(w_TNFA-b)))) - (alphaTNFA*TNFA)
dTGFB = 1/(1+(exp(-h*(w_TGFB-b)))) - (alphaTGFB*TGFB)
dFra1 = 1/(1+(exp(-h*(w_Fra1-b)))) - (alphaFra1*Fra1)
dHIF1A = 1/(1+(exp(-h*(w_HIF1A-b)))) - (alphaHIF1A*HIF1A)
dERK = 1/(1+(exp(-h*(w_ERK-b)))) - (alphaERK*ERK)
dIL4e = 1/(1+(exp(-h*(w_IL4e-b)))) - (alphaIL4e*IL4e)
dIL10e = 1/(1+(exp(-h*(w_IL10e-b)))) - (alphaIL10e*IL10e)
dTNFAe = 1/(1+(exp(-h*(w_TNFAe-b)))) - (alphaTNFAe*TNFAe)
dTGFBe = 1/(1+(exp(-h*(w_TGFBe-b)))) - (alphaTGFBe*TGFBe)
dIgG = 1/(1+(exp(-h*(w_IgG-b)))) - (alphaIgG*IgG)
dA2a = 1/(1+(exp(-h*(w_A2a-b)))) - (alphaA2a*A2a)
dIFNG = 1/(1+(exp(-h*(w_IFNG-b)))) - (alphaIFNG*IFNG)
dIFNB = 1/(1+(exp(-h*(w_IFNB-b)))) - (alphaIFNB*IFNB)
dIL1B = 1/(1+(exp(-h*(w_IL1B-b)))) - (alphaIL1B*IL1B)
dIL6e = 1/(1+(exp(-h*(w_IL6e-b)))) - (alphaIL6e*IL6e)
dMCSF = 1/(1+(exp(-h*(w_MCSF-b)))) - (alphaMCSF*MCSF)
dGMCSF = 1/(1+(exp(-h*(w_GMCSF-b)))) - (alphaGMCSF*GMCSF)
dHipoxia = 1/(1+(exp(-h*(w_Hipoxia-b)))) - (alphaHipoxia*Hipoxia)
dGCGCR = 1/(1+(exp(-h*(w_GCGCR-b)))) - (alphaGCGCR*GCGCR)
list(c(dSTAT3,dSMAD23,dSTAT6,dSOCS1,dTLR4,dNFKB,dSTAT1,dSOCS3,dAP1,dIL10,dTNFA,dTGFB,dFra1,dHIF1A,dERK,dIL4e,dIL10e,dTNFAe,dTGFBe,dIgG,dA2a,dIFNG,dIFNB,dIL1B,dIL6e,dMCSF,dGMCSF,dHipoxia,dGCGCR))
})
}