Skip to content

Commit aa3c648

Browse files
authored
Bertozzi_PNAS2020: update to match Table 1 of paper (#222)
* refactor into paper parameters; change infection reaction to avoid copasi warning * add optimized parameters to parameters table * checklist; use `petab.v1` * delete notes; info is outdated or in #99 * change tests to python 3.11
1 parent f08442c commit aa3c648

22 files changed

Lines changed: 190 additions & 1263 deletions

.github/workflows/tests.yml

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -8,7 +8,7 @@ jobs:
88
runs-on: ubuntu-latest
99
strategy:
1010
matrix:
11-
python-version: [3.9]
11+
python-version: [3.11]
1212

1313
steps:
1414
- name: Check out repository

Benchmark-Models/Bertozzi_PNAS2020/Bertozzi_PNAS2020.yaml

Lines changed: 0 additions & 13 deletions
This file was deleted.

Benchmark-Models/Bertozzi_PNAS2020/Bertozzi_PNAS2020_notes_on_model.txt

Lines changed: 0 additions & 5 deletions
This file was deleted.
Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,3 @@
1+
conditionId N_ I0_ R0_ gamma_
2+
u_CA 39560000 I0_CA R0_CA gamma_CA
3+
u_NY 19540000 I0_NY R0_NY gamma_NY

Benchmark-Models/Bertozzi_PNAS2020/experimentalCondition_Bertozzi_PNAS2020.tsv

Lines changed: 0 additions & 5 deletions
This file was deleted.

Benchmark-Models/Bertozzi_PNAS2020/measurementData_Bertozzi_PNAS2020.tsv

Lines changed: 0 additions & 139 deletions
This file was deleted.
Lines changed: 23 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,23 @@
1+
observableId simulationConditionId measurement time datasetId realDate
2+
y_I_CA u_CA 225 0 data_CA 2020-03-20
3+
y_I_CA u_CA 187 1 data_CA 2020-03-21
4+
y_I_CA u_CA 282 2 data_CA 2020-03-22
5+
y_I_CA u_CA 462 3 data_CA 2020-03-23
6+
y_I_CA u_CA 430 4 data_CA 2020-03-24
7+
y_I_CA u_CA 460 5 data_CA 2020-03-25
8+
y_I_CA u_CA 901 6 data_CA 2020-03-26
9+
y_I_CA u_CA 758 7 data_CA 2020-03-27
10+
y_I_CA u_CA 437 8 data_CA 2020-03-28
11+
y_I_CA u_CA 757 9 data_CA 2020-03-29
12+
y_I_CA u_CA 1286 10 data_CA 2020-03-30
13+
y_I_NY u_NY 2945 0 data_NY 2020-03-20
14+
y_I_NY u_NY 3400 1 data_NY 2020-03-21
15+
y_I_NY u_NY 4090 2 data_NY 2020-03-22
16+
y_I_NY u_NY 5084 3 data_NY 2020-03-23
17+
y_I_NY u_NY 4797 4 data_NY 2020-03-24
18+
y_I_NY u_NY 5160 5 data_NY 2020-03-25
19+
y_I_NY u_NY 7036 6 data_NY 2020-03-26
20+
y_I_NY u_NY 6999 7 data_NY 2020-03-27
21+
y_I_NY u_NY 7534 8 data_NY 2020-03-28
22+
y_I_NY u_NY 7237 9 data_NY 2020-03-29
23+
y_I_NY u_NY 7015 10 data_NY 2020-03-30
Lines changed: 101 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,101 @@
1+
<?xml version='1.0' encoding='UTF-8'?>
2+
<sbml xmlns="http://www.sbml.org/sbml/level3/version2/core" level="3" version="2">
3+
<model id="Bertozzi_PNAS2020" name="Bertozzi_PNAS2020">
4+
<annotation>
5+
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:bqbiol="http://biomodels.net/biology-qualifiers/" xmlns:bqmodel="http://biomodels.net/model-qualifiers/">
6+
<rdf:Description rdf:about="#benchmark_models_petab">
7+
<bqbiol:isDescribedBy>
8+
<rdf:Bag>
9+
<rdf:li rdf:resource="http://identifiers.org/doi/10.1073/pnas.2006520117"/>
10+
</rdf:Bag>
11+
</bqbiol:isDescribedBy>
12+
</rdf:Description>
13+
</rdf:RDF>
14+
</annotation>
15+
<listOfCompartments>
16+
<compartment id="compartment" size="1" constant="true"/>
17+
</listOfCompartments>
18+
<listOfSpecies>
19+
<species id="I_" compartment="compartment" initialConcentration="0.1" boundaryCondition="false" constant="false" hasOnlySubstanceUnits="true"/>
20+
<species id="R_" compartment="compartment" initialConcentration="0" boundaryCondition="false" constant="false" hasOnlySubstanceUnits="true"/>
21+
<species id="S_" compartment="compartment" initialConcentration="0.9" boundaryCondition="false" constant="false" hasOnlySubstanceUnits="true"/>
22+
</listOfSpecies>
23+
<listOfParameters>
24+
<parameter id="N_" value="1" constant="true"/>
25+
<parameter id="beta_N" value="0" constant="true"/>
26+
<parameter id="gamma_" value="0.1" constant="true"/>
27+
<parameter id="I0_" value="0.1" constant="true"/>
28+
<parameter id="R0_" value="0.1" constant="true"/>
29+
</listOfParameters>
30+
<listOfInitialAssignments>
31+
<initialAssignment symbol="I_">
32+
<math xmlns="http://www.w3.org/1998/Math/MathML">
33+
<apply>
34+
<times/>
35+
<ci> I0_ </ci>
36+
</apply>
37+
</math>
38+
</initialAssignment>
39+
<initialAssignment symbol="S_">
40+
<math xmlns="http://www.w3.org/1998/Math/MathML">
41+
<apply>
42+
<minus/>
43+
<ci> N_ </ci>
44+
<ci> I0_ </ci>
45+
</apply>
46+
</math>
47+
</initialAssignment>
48+
<initialAssignment symbol="beta_N">
49+
<math xmlns="http://www.w3.org/1998/Math/MathML">
50+
<apply>
51+
<divide/>
52+
<apply>
53+
<times/>
54+
<ci> R0_ </ci>
55+
<ci> gamma_ </ci>
56+
</apply>
57+
<ci> N_ </ci>
58+
</apply>
59+
</math>
60+
</initialAssignment>
61+
</listOfInitialAssignments>
62+
<listOfReactions>
63+
<reaction id="infection" reversible="false">
64+
<listOfReactants>
65+
<speciesReference species="S_" stoichiometry="1" constant="true"/>
66+
<speciesReference species="I_" stoichiometry="1" constant="true"/>
67+
</listOfReactants>
68+
<listOfProducts>
69+
<speciesReference species="I_" stoichiometry="2" constant="true"/>
70+
</listOfProducts>
71+
<kineticLaw>
72+
<math xmlns="http://www.w3.org/1998/Math/MathML">
73+
<apply>
74+
<times/>
75+
<ci> beta_N </ci>
76+
<ci> I_ </ci>
77+
<ci> S_ </ci>
78+
</apply>
79+
</math>
80+
</kineticLaw>
81+
</reaction>
82+
<reaction id="recovery" reversible="false">
83+
<listOfReactants>
84+
<speciesReference species="I_" stoichiometry="1" constant="true"/>
85+
</listOfReactants>
86+
<listOfProducts>
87+
<speciesReference species="R_" stoichiometry="1" constant="true"/>
88+
</listOfProducts>
89+
<kineticLaw>
90+
<math xmlns="http://www.w3.org/1998/Math/MathML">
91+
<apply>
92+
<times/>
93+
<ci> gamma_ </ci>
94+
<ci> I_ </ci>
95+
</apply>
96+
</math>
97+
</kineticLaw>
98+
</reaction>
99+
</listOfReactions>
100+
</model>
101+
</sbml>

0 commit comments

Comments
 (0)