Using Logistic Growth Models

The spread of a virus such as COVID-19 depends on how many people have the virus and how many people are left in the population to which the virus can spread. This behavior can be modeled using a logistic growth model. In this section, you will learn about the different components of a logistic growth model and what behaviors it is used to model.

Using Logistic Growth Models

Exponential growth cannot continue forever. Exponential models, while they may be useful in the short term, tend to fall apart the longer they continue. Consider an aspiring writer who writes a single line on day one and plans to double the number of lines she writes each day for a month. By the end of the month, she must write over 17 billion lines, or one-half-billion pages. It is impractical, if not impossible, for anyone to write that much in such a short period of time. Eventually, an exponential model must begin to approach some limiting value, and then the growth is forced to slow. For this reason, it is often better to use a model with an upper bound instead of an exponential growth model, though the exponential growth model is still useful over a short term, before approaching the limiting value.

The logistic growth model is approximately exponential at first, but it has a reduced rate of growth as the output approaches the model's upper bound, called the carrying capacity. For constants a, b, and  c, the logistic growth of a population over time x is represented by the model

f(t)=\dfrac{c}{1+ae^{−bx}}

The graph in Figure 6 shows how the growth rate changes over time. The graph increases from left to right, but the growth rate only increases until it reaches its point of maximum growth rate, at which point the rate of increase decreases.

Graph of f(t)=c/(1+ae^(-tx)). The carrying capacity is the asymptote at y=c. The initial value of population is (0, c/(1+a)).

Figure 6

Logistic Growth

The logistic growth model is

f(t)=\dfrac{c}{1+ae^{−bx}}

where
  • \frac{c}{1+a} is the initial value
  • c is the carrying capacity, or limiting value
  • b is a constant determined by the rate of growth.

Example 6
Using the Logistic-Growth Model
An influenza epidemic spreads through a population rapidly, at a rate that depends on two factors: The more people who have the flu, the more rapidly it spreads, and also the more uninfected people there are, the more rapidly it spreads. These two factors make the logistic model a good one to study the spread of communicable diseases. And, clearly, there is a maximum value for the number of people infected: the entire population.

For example, at time t=0 there is one person in a community of 1,000 people who has the flu. So, in that community, at most 1,000 people can have the flu. Researchers find that for this particular strain of the flu, the logistic growth constant is b=0.6030. Estimate the number of people in this community who will have had this flu after ten days. Predict how many people in this community will have had this flu after a long period of time has passed.

Solution
We substitute the given data into the logistic growth model

f(t)=\dfrac{c}{1+ae^{−bx}}

Because at most 1,000 people, the entire population of the community, can get the flu, we know the limiting value is c=1000.  To find a, we use the formula that the number of cases at time t=0 is \frac{c}{1+a}=1, from which it follows that a=999. This model predicts that, after ten days, the number of people who have had the flu is f(t)=\dfrac{1000}{1+999e^{−0.6030x}}≈293.8. Because the actual number must be a whole number (a person has either had the flu or not) we round to 294. In the long term, the number of people who will contract the flu is the limiting value, c=1000.

Analysis
Remember that, because we are dealing with a virus, we cannot predict with certainty the number of people infected. The model only approximates the number of people infected and will not give us exact or actual values.

The graph in Figure 7 gives a good picture of how this model fits the data.
Graph of f(x)=1000/(1+999e^(-0.5030x)) with the y-axis labeled as “Cases” and the x-axis labeled as “Days”. There was 1 case
Figure 7 The graph of f(t)=\dfrac{1000}{1+999e^{−0.6030x}}

Try It #5
Using the model in Example 6, estimate the number of cases of flu on day 15.


Choosing an Appropriate Model for Data

Now that we have discussed various mathematical models, we need to learn how to choose the appropriate model for the raw data we have. Many factors influence the choice of a mathematical model, among which are experience, scientific laws, and patterns in the data itself. Not all data can be described by elementary functions. Sometimes, a function is chosen that approximates the data over a given interval. For instance, suppose data were gathered on the number of homes bought in the United States from the years 1960 to 2013. After plotting these data in a scatter plot, we notice that the shape of the data from the years 2000 to 2013 follow a logarithmic curve. We could restrict the interval from 2000 to 2010, apply regression analysis using a logarithmic model, and use it to predict the number of home buyers for the year 2015.

Three kinds of functions that are often useful in mathematical models are linear functions, exponential functions, and logarithmic functions. If the data lies on a straight line, or seems to lie approximately along a straight line, a linear model may be best. If the data is non-linear, we often consider an exponential or logarithmic model, though other models, such as quadratic models, may also be considered.

In choosing between an exponential model and a logarithmic model, we look at the way the data curves. This is called the concavity. If we draw a line between two data points, and all (or most) of the data between those two points lies above that line, we say the curve is concave down. We can think of it as a bowl that bends downward and therefore cannot hold water. If all (or most) of the data between those two points lies below the line, we say the curve is concave up. In this case, we can think of a bowl that bends upward and can therefore hold water. An exponential curve, whether rising or falling, whether representing growth or decay, is always concave up away from its horizontal asymptote. A logarithmic curve is always concave away from its vertical asymptote. In the case of positive data, which is the most common case, an exponential curve is always concave up, and a logarithmic curve always concave down.

A logistic curve changes concavity. It starts out concave up and then changes to concave down beyond a certain point, called a point of inflection.

After using the graph to help us choose a type of function to use as a model, we substitute points, and solve to find the parameters. We reduce round-off error by choosing points as far apart as possible.


Example 7
Choosing a Mathematical Model

Does a linear, exponential, logarithmic, or logistic model best fit the values listed in Table 1? Find the model, and use a graph to check your choice.

x 1 2 3 4 5 6 7 8 9
y 0 1.386 2.197 2.773 3.219 3.584 3.892 4.159 4.394

Table 1

Solution
First, plot the data on a graph as in Figure 8. For the purpose of graphing, round the data to two decimal places.

Graph of the previous table's values.
Figure 8

Clearly, the points do not lie on a straight line, so we reject a linear model. If we draw a line between any two of the points, most or all of the points between those two points lie above the line, so the graph is concave down, suggesting a logarithmic model. We can try y=aln(bx). Plugging in the first point, (1,0), gives 0=alnb. We reject the case that a=0 (if it were, all outputs would be 0), so we know ln(b)=0. Thus b=1 and y=aln(x). Next we can use the point (9,4.394) to solve for a:

y=aln(x)
4.394=aln(9)
a=\dfrac{4.394}{ln(9)}

Because a=\dfrac{4.394}{ln(9)}≈2, an appropriate model for the data is y=2ln(x).

To check the accuracy of the model, we graph the function together with the given points as in Figure 9.

Graph of previous table's values showing that it fits the function y=2ln(x) with an asymptote at x=0.

Figure 9 The graph of y=2lnx.

We can conclude that the model is a good fit to the data.
Compare Figure 9 to the graph of y=ln(x^2) shown in Figure 10.

Graph of previous table's values showing that it fits the function y=2ln(x) with an asymptote at x=0.

Figure 10 The graph of y=ln(x^2)

The graphs appear to be identical when x > 0.  A quick check confirms this conclusion: y=ln(x^2)=2ln(x) for x > 0.

However, if x < 0, the graph of y=ln(x^2) includes a "extra" branch, as shown in Figure 11. This occurs because, while y=2ln(x)  cannot have negative values in the domain (as such values would force the argument to be negative), the function y=ln(x^2)  can have negative domain values.


Graph of y=ln(x^2).
Figure 11

Try It #6
Does a linear, exponential, or logarithmic model best fit the data in Table 2? Find the model.

  
1 2 3 4 5 6 7 8 9
y   
3.297 5.437 8.963 14.778 24.365 40.172 66.231 109.196 180.034

Table 2


Source: Rice University, https://openstax.org/books/college-algebra/pages/6-7-exponential-and-logarithmic-models
Creative Commons License This work is licensed under a Creative Commons Attribution 4.0 License.