Fitting Linear Models to Data

Drawing and Interpreting Scatter Plots

A scatter plot is a graph of plotted points that may show a relationship between two sets of data. If the relationship is from a linear model, or a model that is nearly linear, the professor can draw conclusions using his knowledge of linear functions. Figure 1 shows a sample scatter plot.

Figure 1 A scatter plot of age and final exam score variables

Notice this scatter plot does not indicate a linear relationship. The points do not appear to follow a trend. In other words, there does not appear to be a relationship between the age of the student and the score on the final exam.

Example 1

Using a Scatter Plot to Investigate Cricket Chirps

Table 1 shows the number of cricket chirps in 15 seconds, for several different air temperatures, in degrees Fahrenheit5. Plot this data, and determine whether the data appears to be linearly related.

5

Selected data from http://classic.globe.gov/fsl/scientistsblog/2007/10/. Retrieved Aug 3, 2010

Chirps	44	35	20.4	33	31	35	18.5	37	26
Temperature	80.5	70.5	57	66	68	72	52	73.5	53

Table 1 Cricket Chirps vs Air Temperature

Show solution

Solution

Plotting this data, as depicted in Figure 2 suggests that there may be a trend. We can see from the trend in the data that the number of chirps increases as the temperature increases. The trend appears to be roughly linear, though certainly not perfectly so.

Figure 2

Finding the Line of Best Fit

Once we recognize a need for a linear function to model that data, the natural follow-up question is “what is that linear function?” One way to approximate our linear function is to sketch the line that seems to best fit the data. Then we can extend the line until we can verify the y-intercept. We can approximate the slope of the line by extending it until we can estimate the $\frac{\text{rise}}{\text{run}}.$

Example 2

Finding a Line of Best Fit

Find a linear function that fits the data in Table 1 by “eyeballing” a line that seems to fit.

Show solution

Solution

On a graph, we could try sketching a line. Using the starting and ending points of our hand drawn line, points (0, 30) and (50, 90), this graph has a slope of

$$\begin{array}{ccc}\hfill m& =& \frac{60}{50}\hfill \\ \hfill & =& 1.2\hfill \end{array}$$

and a y-intercept at 30. This gives an equation of

$$T(c)=1.2c+30$$

where $c$ is the number of chirps in 15 seconds, and $T\left(c\right)$ is the temperature in degrees Fahrenheit. The resulting equation is represented in Figure 3.

Figure 3

Analysis

This linear equation can then be used to approximate answers to various questions we might ask about the trend.

Recognizing Interpolation or Extrapolation

While the data for most examples does not fall perfectly on the line, the equation is our best guess as to how the relationship will behave outside of the values for which we have data. We use a process known as interpolation when we predict a value inside the domain and range of the data. The process of extrapolation is used when we predict a value outside the domain and range of the data.

Figure 4 compares the two processes for the cricket-chirp data addressed in Example 2. We can see that interpolation would occur if we used our model to predict temperature when the values for chirps are between 18.5 and 44. Extrapolation would occur if we used our model to predict temperature when the values for chirps are less than 18.5 or greater than 44.

There is a difference between making predictions inside the domain and range of values for which we have data and outside that domain and range. Predicting a value outside of the domain and range has its limitations. When our model no longer applies after a certain point, it is sometimes called model breakdown. For example, predicting a cost function for a period of two years may involve examining the data where the input is the time in years and the output is the cost. But if we try to extrapolate a cost when $x=50,$ that is in 50 years, the model would not apply because we could not account for factors fifty years in the future.

Figure 4 Interpolation occurs within the domain and range of the provided data whereas extrapolation occurs outside.

Interpolation and Extrapolation

Different methods of making predictions are used to analyze data.

The method of interpolation involves predicting a value inside the domain and/or range of the data.
The method of extrapolation involves predicting a value outside the domain and/or range of the data.
Model breakdown occurs at the point when the model no longer applies.

Example 3

Understanding Interpolation and Extrapolation

Use the cricket data from Table 1 to answer the following questions:

1. ⓐWould predicting the temperature when crickets are chirping 30 times in 15 seconds be interpolation or extrapolation? Make the prediction, and discuss whether it is reasonable.
2. ⓑWould predicting the number of chirps crickets will make at 40 degrees be interpolation or extrapolation? Make the prediction, and discuss whether it is reasonable.

Show solution

Solution

1. ⓐThe number of chirps in the data provided varied from 18.5 to 44. A prediction at 30 chirps per 15 seconds is inside the domain of our data, so would be interpolation. Using our model:
   $$\begin{array}{ccc}\hfill \text{T (30)}& =& \text{30 + 1.2(30)}\hfill \\ & =& \text{66 degrees}\hfill \end{array}$$

   Based on the data we have, this value seems reasonable.

2. ⓑThe temperature values varied from 52 to 80.5. Predicting the number of chirps at 40 degrees is extrapolation because 40 is outside the range of our data. Using our model:
   $$\begin{array}{l}40=30+1.2c\hfill \\ 10=1.2c\hfill \\ c\approx 8.33\hfill \end{array}$$

We can compare the regions of interpolation and extrapolation using Figure 5.

Figure 5

Analysis

Our model predicts the crickets would chirp 8.33 times in 15 seconds. While this might be possible, we have no reason to believe our model is valid outside the domain and range. In fact, generally crickets stop chirping altogether below around 50 degrees.

Try It #1

According to the data from Table 1, what temperature can we predict it is if we counted 20 chirps in 15 seconds?

Finding the Line of Best Fit Using a Graphing Utility

While eyeballing a line works reasonably well, there are statistical techniques for fitting a line to data that minimize the differences between the line and data values6. One such technique is called least squares regression and can be computed by many graphing calculators, spreadsheet software, statistical software, and many web-based calculators7. Least squares regression is one means to determine the line that best fits the data, and here we will refer to this method as linear regression.

6

Technically, the method minimizes the sum of the squared differences in the vertical direction between the line and the data values.

7

For example, http://www.shodor.org/unchem/math/lls/leastsq.html

How To

Given data of input and corresponding outputs from a linear function, find the best fit line using linear regression.

1. Enter the input in List 1 (L1).
2. Enter the output in List 2 (L2).
3. On a graphing utility, select Linear Regression (LinReg).

Example 4

Finding a Least Squares Regression Line

Find the least squares regression line using the cricket-chirp data in Table 2.

Show solution

Solution

1. Enter the input (chirps) in List 1 (L1).
2. Enter the output (temperature) in List 2 (L2). See Table 2.

   L1	44	35	20.4	33	31	35	18.5	37	26
   L2	80.5	70.5	57	66	68	72	52	73.5	53

   Table 2
3. On a graphing utility, select Linear Regression (LinReg). Using the cricket chirp data from earlier, with technology we obtain the equation:

$$T(c)=30.281+1.143c$$

Analysis

Notice that this line is quite similar to the equation we “eyeballed” but should fit the data better. Notice also that using this equation would change our prediction for the temperature when hearing 30 chirps in 15 seconds from 66 degrees to:

$$\begin{array}{l}T(30)=30.281+1.143(30)\hfill \\ =64.571\hfill \\ \approx 64.6\text{degrees}\hfill \end{array}$$

The graph of the scatter plot with the least squares regression line is shown in Figure 6.

Figure 6

Q&A

Will there ever be a case where two different lines will serve as the best fit for the data?

No. There is only one best fit line.

Distinguishing Between Linear and Nonlinear Models

As we saw above with the cricket-chirp model, some data exhibit strong linear trends, but other data, like the final exam scores plotted by age, are clearly nonlinear. Most calculators and computer software can also provide us with the correlation coefficient, which is a measure of how closely the line fits the data. Many graphing calculators require the user to turn a "diagnostic on" selection to find the correlation coefficient, which mathematicians label as $r$ The correlation coefficient provides an easy way to get an idea of how close to a line the data falls.

We should compute the correlation coefficient only for data that follows a linear pattern or to determine the degree to which a data set is linear. If the data exhibits a nonlinear pattern, the correlation coefficient for a linear regression is meaningless. To get a sense for the relationship between the value of $r$ and the graph of the data, Figure 7 shows some large data sets with their correlation coefficients. Remember, for all plots, the horizontal axis shows the input and the vertical axis shows the output.

Figure 7 Plotted data and related correlation coefficients. (credit: “DenisBoigelot,” Wikimedia Commons)

Fitting a Regression Line to a Set of Data

Once we determine that a set of data is linear using the correlation coefficient, we can use the regression line to make predictions. As we learned above, a regression line is a line that is closest to the data in the scatter plot, which means that only one such line is a best fit for the data.

Example 6

Using a Regression Line to Make Predictions

Gasoline consumption in the United States has been steadily increasing. Consumption data from 1994 to 2004 is shown in Table 3.8 Determine whether the trend is linear, and if so, find a model for the data. Use the model to predict the consumption in 2008.

8

http://www.bts.gov/publications/national_transportation_statistics/2005/html/table_04_10.html

Year	'94	'95	'96	'97	'98	'99	'00	'01	'02	'03	'04
Consumption (billions of gallons)	113	116	118	119	123	125	126	128	131	133	136

Table 3

The scatter plot of the data, including the least squares regression line, is shown in Figure 8.

Figure 8

Show solution

Solution

We can introduce a new input variable, $t,$ representing years since 1994.

The least squares regression equation is:

$$C(t)=113.318+2.209t$$

Using technology, the correlation coefficient was calculated to be 0.9965, suggesting a very strong increasing linear trend.

Using this to predict consumption in 2008 $(t=14),$

$$\begin{array}{l}C(14)=113.318+2.209(14)\hfill \\ =144.244\hfill \end{array}$$

The model predicts 144.244 billion gallons of gasoline consumption in 2008.