Modeling bacterial growth with differential equations
I hope this is the right place for this question. I am working on building a growth model for bacteria for a risk assessment, and would like to move the growth model past static temperature conditions. The primary growth model, using the modified Gompertz equation, is:
$$log_{10} R = A_g cdot exp(-exp{mucdot exp(1)cdot(l-t)/A_g +1}), $$
where $R$ is the relative population. There are secondary models for the three parameters $A_g$, $mu$ and $l$ (i.e., the max number of bacteria at the stationary phase; the growth rate; and the lag phase, respectively). For example, for mu the equation is:
$$sqrt{mu}=0.0421cdot(T-12.0570),$$ where $T$ is temperature.
All of the secondary models are dependent upon temperature, since bacteria obviously favor warmer conditions and grow faster, up to a certain point.
My trouble is that one paper specifying a different growth model (logistic) was able to model bacterial growth under ANY temperature regime you can think of by deriving a differential equation for their primary model. I have little experience in differential equations and sadly, I am actually just not even sure which part of my model would need to be integrated. Since I am interested in how the changing parameters (due to changing temperature) lead to growth in the end, I would think I would need to integrate with respect to $T$? But I am unsure if this means solving the differential equations of the second-order parameters (because they are the only ones that directly have $T$ in them) or if that means doing something else to the primary equation?
For reference, the paper I am using most often is:
Corradini, M. G., Amézquita, A., Normand, M. D., & Peleg, M. (2006). Modeling and predicting non-isothermal microbial growth using general purpose software. International journal of food microbiology, 106(2), 223-228.
(Open source PDF available through Google Scholar)
differential-equations applications
asked Jan 30 '14 at 20:14
HFBrowningHFBrowning
1136
|
show 1 more comment
I hope this is the right place for this question. I am working on building a growth model for bacteria for a risk assessment, and would like to move the growth model past static temperature conditions. The primary growth model, using the modified Gompertz equation, is:
$$log_{10} R = A_g cdot exp(-exp{mucdot exp(1)cdot(l-t)/A_g +1}), $$
where $R$ is the relative population. There are secondary models for the three parameters $A_g$, $mu$ and $l$ (i.e., the max number of bacteria at the stationary phase; the growth rate; and the lag phase, respectively). For example, for mu the equation is:
$$sqrt{mu}=0.0421cdot(T-12.0570),$$ where $T$ is temperature.
All of the secondary models are dependent upon temperature, since bacteria obviously favor warmer conditions and grow faster, up to a certain point.
My trouble is that one paper specifying a different growth model (logistic) was able to model bacterial growth under ANY temperature regime you can think of by deriving a differential equation for their primary model. I have little experience in differential equations and sadly, I am actually just not even sure which part of my model would need to be integrated. Since I am interested in how the changing parameters (due to changing temperature) lead to growth in the end, I would think I would need to integrate with respect to $T$? But I am unsure if this means solving the differential equations of the second-order parameters (because they are the only ones that directly have $T$ in them) or if that means doing something else to the primary equation?
For reference, the paper I am using most often is:
Corradini, M. G., Amézquita, A., Normand, M. D., & Peleg, M. (2006). Modeling and predicting non-isothermal microbial growth using general purpose software. International journal of food microbiology, 106(2), 223-228.
(Open source PDF available through Google Scholar)
differential-equations applications
asked Jan 30 '14 at 20:14
HFBrowningHFBrowning
1136
There is the Gompertz equation, and the Gompertz differential equation (I'm relying on wikipedia for this). The former is the solution of the latter. Perhaps you can take your formulation and from it create a "modified Gompertz differential equation"?
– Jeff Snider
Jan 30 '14 at 21:04
This might be a good route, if I understood more. Basically I think I am looking for a conceptual explanation of what needs to happen, rather than mechanistic? I've taken a fair amount of calculus but it was a while ago, and I'm getting stuck in the weeds on this one. For instance: I know that if I wanted to calculate overall growth for a simple varying temperature, I could just modify the second-order parameters, stick them into the primary model, and see much growth happened for each sequential time period. But I want to do that for something consistently decreasing, and not sure (cont)
– HFBrowning
Jan 30 '14 at 21:19
what I need to integrate/or differentiate with respect to. Plus the first option is an excessive amount of work. Here's a pic of what's "easy" to me, and then what I actually want to do. So I think it's a very easy question for math people probably. [IMG]i60.tinypic.com/24y4wsy.png[/IMG]
– HFBrowning
Jan 30 '14 at 21:22
In general when you have a complex differential equation and you're more interested in the result than the math, solving it is best done by a numerical package. There are a lot of different packages to choose from. The first thing you need is to formulate your problem as a differential equation, where your equation tells you the rate of change rather than the current population size.
– Jeff Snider
Jan 30 '14 at 21:26
Jeff - exactly! Your last sentence is what I am struggling with but maybe didn't articulate well.
– HFBrowning
Jan 30 '14 at 21:28
|
show 1 more comment
I hope this is the right place for this question. I am working on building a growth model for bacteria for a risk assessment, and would like to move the growth model past static temperature conditions. The primary growth model, using the modified Gompertz equation, is:
$$log_{10} R = A_g cdot exp(-exp{mucdot exp(1)cdot(l-t)/A_g +1}), $$
where $R$ is the relative population. There are secondary models for the three parameters $A_g$, $mu$ and $l$ (i.e., the max number of bacteria at the stationary phase; the growth rate; and the lag phase, respectively). For example, for mu the equation is:
$$sqrt{mu}=0.0421cdot(T-12.0570),$$ where $T$ is temperature.
All of the secondary models are dependent upon temperature, since bacteria obviously favor warmer conditions and grow faster, up to a certain point.
My trouble is that one paper specifying a different growth model (logistic) was able to model bacterial growth under ANY temperature regime you can think of by deriving a differential equation for their primary model. I have little experience in differential equations and sadly, I am actually just not even sure which part of my model would need to be integrated. Since I am interested in how the changing parameters (due to changing temperature) lead to growth in the end, I would think I would need to integrate with respect to $T$? But I am unsure if this means solving the differential equations of the second-order parameters (because they are the only ones that directly have $T$ in them) or if that means doing something else to the primary equation?
For reference, the paper I am using most often is:
Corradini, M. G., Amézquita, A., Normand, M. D., & Peleg, M. (2006). Modeling and predicting non-isothermal microbial growth using general purpose software. International journal of food microbiology, 106(2), 223-228.
(Open source PDF available through Google Scholar)
differential-equations applications
asked Jan 30 '14 at 20:14
HFBrowningHFBrowning
1136
I hope this is the right place for this question. I am working on building a growth model for bacteria for a risk assessment, and would like to move the growth model past static temperature conditions. The primary growth model, using the modified Gompertz equation, is:
$$log_{10} R = A_g cdot exp(-exp{mucdot exp(1)cdot(l-t)/A_g +1}), $$
where $R$ is the relative population. There are secondary models for the three parameters $A_g$, $mu$ and $l$ (i.e., the max number of bacteria at the stationary phase; the growth rate; and the lag phase, respectively). For example, for mu the equation is:
$$sqrt{mu}=0.0421cdot(T-12.0570),$$ where $T$ is temperature.
All of the secondary models are dependent upon temperature, since bacteria obviously favor warmer conditions and grow faster, up to a certain point.
My trouble is that one paper specifying a different growth model (logistic) was able to model bacterial growth under ANY temperature regime you can think of by deriving a differential equation for their primary model. I have little experience in differential equations and sadly, I am actually just not even sure which part of my model would need to be integrated. Since I am interested in how the changing parameters (due to changing temperature) lead to growth in the end, I would think I would need to integrate with respect to $T$? But I am unsure if this means solving the differential equations of the second-order parameters (because they are the only ones that directly have $T$ in them) or if that means doing something else to the primary equation?
For reference, the paper I am using most often is:
Corradini, M. G., Amézquita, A., Normand, M. D., & Peleg, M. (2006). Modeling and predicting non-isothermal microbial growth using general purpose software. International journal of food microbiology, 106(2), 223-228.
(Open source PDF available through Google Scholar)
differential-equations applications
differential-equations applications
asked Jan 30 '14 at 20:14
HFBrowningHFBrowning
1136
asked Jan 30 '14 at 20:14
HFBrowningHFBrowning
1136
edited Jan 30 '14 at 21:00
HFBrowning
asked Jan 30 '14 at 20:14
HFBrowningHFBrowning
1136
asked Jan 30 '14 at 20:14
HFBrowningHFBrowning
1136
asked Jan 30 '14 at 20:14
HFBrowningHFBrowning
1136
1136
There is the Gompertz equation, and the Gompertz differential equation (I'm relying on wikipedia for this). The former is the solution of the latter. Perhaps you can take your formulation and from it create a "modified Gompertz differential equation"?
– Jeff Snider
Jan 30 '14 at 21:04
This might be a good route, if I understood more. Basically I think I am looking for a conceptual explanation of what needs to happen, rather than mechanistic? I've taken a fair amount of calculus but it was a while ago, and I'm getting stuck in the weeds on this one. For instance: I know that if I wanted to calculate overall growth for a simple varying temperature, I could just modify the second-order parameters, stick them into the primary model, and see much growth happened for each sequential time period. But I want to do that for something consistently decreasing, and not sure (cont)
– HFBrowning
Jan 30 '14 at 21:19
what I need to integrate/or differentiate with respect to. Plus the first option is an excessive amount of work. Here's a pic of what's "easy" to me, and then what I actually want to do. So I think it's a very easy question for math people probably. [IMG]i60.tinypic.com/24y4wsy.png[/IMG]
– HFBrowning
Jan 30 '14 at 21:22
In general when you have a complex differential equation and you're more interested in the result than the math, solving it is best done by a numerical package. There are a lot of different packages to choose from. The first thing you need is to formulate your problem as a differential equation, where your equation tells you the rate of change rather than the current population size.
– Jeff Snider
Jan 30 '14 at 21:26
Jeff - exactly! Your last sentence is what I am struggling with but maybe didn't articulate well.
– HFBrowning
Jan 30 '14 at 21:28
|
show 1 more comment
There is the Gompertz equation, and the Gompertz differential equation (I'm relying on wikipedia for this). The former is the solution of the latter. Perhaps you can take your formulation and from it create a "modified Gompertz differential equation"?
– Jeff Snider
Jan 30 '14 at 21:04
This might be a good route, if I understood more. Basically I think I am looking for a conceptual explanation of what needs to happen, rather than mechanistic? I've taken a fair amount of calculus but it was a while ago, and I'm getting stuck in the weeds on this one. For instance: I know that if I wanted to calculate overall growth for a simple varying temperature, I could just modify the second-order parameters, stick them into the primary model, and see much growth happened for each sequential time period. But I want to do that for something consistently decreasing, and not sure (cont)
– HFBrowning
Jan 30 '14 at 21:19
what I need to integrate/or differentiate with respect to. Plus the first option is an excessive amount of work. Here's a pic of what's "easy" to me, and then what I actually want to do. So I think it's a very easy question for math people probably. [IMG]i60.tinypic.com/24y4wsy.png[/IMG]
– HFBrowning
Jan 30 '14 at 21:22
In general when you have a complex differential equation and you're more interested in the result than the math, solving it is best done by a numerical package. There are a lot of different packages to choose from. The first thing you need is to formulate your problem as a differential equation, where your equation tells you the rate of change rather than the current population size.
– Jeff Snider
Jan 30 '14 at 21:26
Jeff - exactly! Your last sentence is what I am struggling with but maybe didn't articulate well.
– HFBrowning
Jan 30 '14 at 21:28
There is the Gompertz equation, and the Gompertz differential equation (I'm relying on wikipedia for this). The former is the solution of the latter. Perhaps you can take your formulation and from it create a "modified Gompertz differential equation"?
– Jeff Snider
Jan 30 '14 at 21:04
There is the Gompertz equation, and the Gompertz differential equation (I'm relying on wikipedia for this). The former is the solution of the latter. Perhaps you can take your formulation and from it create a "modified Gompertz differential equation"?
– Jeff Snider
Jan 30 '14 at 21:04
This might be a good route, if I understood more. Basically I think I am looking for a conceptual explanation of what needs to happen, rather than mechanistic? I've taken a fair amount of calculus but it was a while ago, and I'm getting stuck in the weeds on this one. For instance: I know that if I wanted to calculate overall growth for a simple varying temperature, I could just modify the second-order parameters, stick them into the primary model, and see much growth happened for each sequential time period. But I want to do that for something consistently decreasing, and not sure (cont)
– HFBrowning
Jan 30 '14 at 21:19
This might be a good route, if I understood more. Basically I think I am looking for a conceptual explanation of what needs to happen, rather than mechanistic? I've taken a fair amount of calculus but it was a while ago, and I'm getting stuck in the weeds on this one. For instance: I know that if I wanted to calculate overall growth for a simple varying temperature, I could just modify the second-order parameters, stick them into the primary model, and see much growth happened for each sequential time period. But I want to do that for something consistently decreasing, and not sure (cont)
– HFBrowning
Jan 30 '14 at 21:19
what I need to integrate/or differentiate with respect to. Plus the first option is an excessive amount of work. Here's a pic of what's "easy" to me, and then what I actually want to do. So I think it's a very easy question for math people probably. [IMG]i60.tinypic.com/24y4wsy.png[/IMG]
– HFBrowning
Jan 30 '14 at 21:22
what I need to integrate/or differentiate with respect to. Plus the first option is an excessive amount of work. Here's a pic of what's "easy" to me, and then what I actually want to do. So I think it's a very easy question for math people probably. [IMG]i60.tinypic.com/24y4wsy.png[/IMG]
– HFBrowning
Jan 30 '14 at 21:22
In general when you have a complex differential equation and you're more interested in the result than the math, solving it is best done by a numerical package. There are a lot of different packages to choose from. The first thing you need is to formulate your problem as a differential equation, where your equation tells you the rate of change rather than the current population size.
– Jeff Snider
Jan 30 '14 at 21:26
In general when you have a complex differential equation and you're more interested in the result than the math, solving it is best done by a numerical package. There are a lot of different packages to choose from. The first thing you need is to formulate your problem as a differential equation, where your equation tells you the rate of change rather than the current population size.
– Jeff Snider
Jan 30 '14 at 21:26
Jeff - exactly! Your last sentence is what I am struggling with but maybe didn't articulate well.
– HFBrowning
Jan 30 '14 at 21:28
Jeff - exactly! Your last sentence is what I am struggling with but maybe didn't articulate well.
– HFBrowning
Jan 30 '14 at 21:28
|
show 1 more comment
1 Answer
1
active
oldest
votes
I have this modified Gompertz equation, how do I convert it back to a differential equation?
For some suitable constants $(a,b,c)$ determined by $(A_g,mu,l)$, one has $$lnln R=a-exp(b-ct)$$ hence $$frac{R'}{Rln R}=cexp(b-ct)=c(a-lnln R)$$ that is, $$R'=cRln R(a-lnln R)$$
answered Jan 22 '18 at 6:12
DidDid
246k23221455
add a comment |
Your Answer
StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f657643%2fmodeling-bacterial-growth-with-differential-equations%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
I have this modified Gompertz equation, how do I convert it back to a differential equation?
For some suitable constants $(a,b,c)$ determined by $(A_g,mu,l)$, one has $$lnln R=a-exp(b-ct)$$ hence $$frac{R'}{Rln R}=cexp(b-ct)=c(a-lnln R)$$ that is, $$R'=cRln R(a-lnln R)$$
answered Jan 22 '18 at 6:12
DidDid
246k23221455
add a comment |
I have this modified Gompertz equation, how do I convert it back to a differential equation?
For some suitable constants $(a,b,c)$ determined by $(A_g,mu,l)$, one has $$lnln R=a-exp(b-ct)$$ hence $$frac{R'}{Rln R}=cexp(b-ct)=c(a-lnln R)$$ that is, $$R'=cRln R(a-lnln R)$$
answered Jan 22 '18 at 6:12
DidDid
246k23221455
add a comment |
I have this modified Gompertz equation, how do I convert it back to a differential equation?
For some suitable constants $(a,b,c)$ determined by $(A_g,mu,l)$, one has $$lnln R=a-exp(b-ct)$$ hence $$frac{R'}{Rln R}=cexp(b-ct)=c(a-lnln R)$$ that is, $$R'=cRln R(a-lnln R)$$
answered Jan 22 '18 at 6:12
DidDid
246k23221455
I have this modified Gompertz equation, how do I convert it back to a differential equation?
For some suitable constants $(a,b,c)$ determined by $(A_g,mu,l)$, one has $$lnln R=a-exp(b-ct)$$ hence $$frac{R'}{Rln R}=cexp(b-ct)=c(a-lnln R)$$ that is, $$R'=cRln R(a-lnln R)$$
answered Jan 22 '18 at 6:12
DidDid
246k23221455
answered Jan 22 '18 at 6:12
DidDid
246k23221455
answered Jan 22 '18 at 6:12
DidDid
246k23221455
answered Jan 22 '18 at 6:12
DidDid
246k23221455
246k23221455
add a comment |
add a comment |
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f657643%2fmodeling-bacterial-growth-with-differential-equations%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
There is the Gompertz equation, and the Gompertz differential equation (I'm relying on wikipedia for this). The former is the solution of the latter. Perhaps you can take your formulation and from it create a "modified Gompertz differential equation"?
– Jeff Snider
Jan 30 '14 at 21:04
This might be a good route, if I understood more. Basically I think I am looking for a conceptual explanation of what needs to happen, rather than mechanistic? I've taken a fair amount of calculus but it was a while ago, and I'm getting stuck in the weeds on this one. For instance: I know that if I wanted to calculate overall growth for a simple varying temperature, I could just modify the second-order parameters, stick them into the primary model, and see much growth happened for each sequential time period. But I want to do that for something consistently decreasing, and not sure (cont)
– HFBrowning
Jan 30 '14 at 21:19
what I need to integrate/or differentiate with respect to. Plus the first option is an excessive amount of work. Here's a pic of what's "easy" to me, and then what I actually want to do. So I think it's a very easy question for math people probably. [IMG]i60.tinypic.com/24y4wsy.png[/IMG]
– HFBrowning
Jan 30 '14 at 21:22
In general when you have a complex differential equation and you're more interested in the result than the math, solving it is best done by a numerical package. There are a lot of different packages to choose from. The first thing you need is to formulate your problem as a differential equation, where your equation tells you the rate of change rather than the current population size.
– Jeff Snider
Jan 30 '14 at 21:26
Jeff - exactly! Your last sentence is what I am struggling with but maybe didn't articulate well.
– HFBrowning
Jan 30 '14 at 21:28