Log–log plot

A log–log plot of y = x (blue), y = x2 (green), and y = x3 (red).
Note the logarithmic scale markings on each of the axes, and that the log x and log y axes (where the logarithms are 0) are where x and y themselves are 1.

In science and engineering, a log–log graph or log–log plot is a two-dimensional graph of numerical data that uses logarithmic scales on both the horizontal and vertical axes. Monomials – relationships of the form y=ax^k – appear as straight lines in a log–log graph, with the power and constant term corresponding to slope and intercept of the line, and thus these graphs are very useful for recognizing these relationships and estimating parameters. Any base can be used for the logarithm, though most common are 10, e, and 2.

Relation with monomials

Given a monomial equation y=ax^k, taking the logarithm of the equation (with any base) yields:

\log y = k \log x + \log a.

Setting X = \log x and Y = \log y, which corresponds to using a log–log graph, yields the equation:

Y = mX + b

where m = k is the slope of the line (gradient) and b = log a is the intercept on the (log y)-axis, meaning where log x = 0, so, reversing the logs, a is the y value corresponding to x = 1.[1]

Equations

The equation for a line on a log–log scale would be:

 \log_{10}F(x) = m \log_{10}x + b,
 F(x) = x^m\cdot10^b,

where m is the slope and b is the intercept point on the log plot.

Slope of a log–log plot

Finding the slope of a log–log plot using ratios

To find the slope of the plot, two points are selected on the x-axis, say x1 and x2. Using the above equation:

 \log[F (x_1)] = m \log (x_1) + b, \,

and

 \mathrm {log}[F (x_2)] = m \log (x_2) + b. \,

The slope m is found taking the difference:

 m = \frac { \mathrm {log} (F_2) - \mathrm {log} (F_1)} { \log(x_2) - \log(x_1) } = \frac {\log (F_2/F_1)}{\log(x_2/x_1)}, \,

where F1 is shorthand for F ( x1 ) and F2 is shorthand for F ( x2 ). The figure at right illustrates the formula. Notice that the slope in the example of the figure is negative. The formula also provides a negative slope, as can be seen from the following property of the logarithm:

 \log(x_1/x_2) = -\log(x_2/x_1). \,

Finding the function from the log–log plot

The above procedure now is reversed to find the form of the function F(x) using its (assumed) known log–log plot. To find the function F, pick some fixed point (x0, F0), where F0 is shorthand for F(x0), somewhere on the straight line in the above graph, and further some other arbitrary point (x1, F1) on the same graph. Then from the slope formula above:

 m = \frac {\log (F_1 / F_0)}{\log(x_1 / x_0)}

which leads to

 \log(F_1 / F_0) = m \log(x_1 / x_0) = \log[(x_1 / x_0)^m ]. \,

Notice that 10log10(F1) = F1. Therefore, the logs can be inverted to find:

 \frac{F_1}{F_0} = \left(\frac{x_1}{x_0}\right)^m

or

F_1 = \frac{F_0}{x_0^m} \,\, x^m, \,

which means that

 F(x) = \mathrm{constant}\cdot x^m.

In other words, F is proportional to x to the power of the slope of the straight line of its log–log graph. Specifically, a straight line on a log–log plot containing points (F0, x0) and (F1, x1) will have the function:

 F(x) = {F_0}\left(\frac{x}{x_0}\right)^\frac {\log (F_1/F_0)}{\log(x_1/x_0)},

Of course, the inverse is true too: any function of the form

 F(x) = \mathrm{constant} \cdot x^m

will have a straight line as its log–log graph representation, where the slope of the line is m.

Finding the area under a straight-line segment of log–log plot

To calculate the area under a continuous, straight-line segment of a log–log plot (or estimating an area of an almost-straight line), take the function defined previously

 F(x) = \mathrm{constant}\cdot x^m.

and integrate it. Since it is only operating on a definite integral (two defined endpoints), the area A under the plot takes the form

 A(x) = \int_{x_0}^{x_1} F(x) dx = \frac{\mathrm{constant}}{m+1}\cdot x^{m+1}: [x_0,x_1]

Rearranging the original equation and plugging in the fixed point values, it is found that

 \mathrm{constant} = \frac{F_0}{x_0^m}

Substituting back into the integral, you find that for A over x0 to x1

 A = \frac{F_0/x_0^m}{m+1}\cdot (x_1^{m+1}-x_0^{m+1})
 \log A = \log \left[\frac{F_0 / x_0^m}{m+1}\cdot (x_1^{m+1}-x_0^{m+1})\right] = \log \frac{F_0}{m+1} - \log \frac{1}{x_0^m} + \log (x_1^{m+1}-x_0^{m+1})
 \log A = \log \frac{F_0}{m+1} + \log \left(\frac{x_1^{m+1}-x_0^{m+1}}{x_0^m}\right) = \log \frac{F_0}{m+1} + \log \left(\frac{x_1^m}{x_0^m}\cdot x_1 - \frac{x_0^{m+1}}{x_0^m}\right)

Therefore:  A = \frac{F_0}{m+1} \cdot \left[x_1 \cdot \left(\frac {x_1}{x_0}\right)^m - x_0\right]

For m=-1, the integral becomes  A_{(m=-1)} = \int_{x_0}^{x_1} F(x) = \int_{x_0}^{x_1} \frac {\mathrm{constant}}{x} = \frac{F_0}{x_0^{-1}} \int_{x_0}^{x_1} \frac {1}{x} = F_0 \cdot x_0 \cdot \ln x: [x_0,x_1]

 A_{(m=-1)} = F_0 \cdot x_0 \cdot \ln \frac{x_1}{x_0}

Applications

These graphs are useful when the parameters a and b need to be estimated from numerical data. Specifications such as this are used frequently in economics.

One example is the estimation of money demand functions based on inventory theory, in which it can be assumed that money demand at time t is given by

M_t = AR_t^bY_t^cU_t,

where M is the real quantity of money held by the public, R is the rate of return on an alternative, higher yielding asset in excess of that on money, Y is the public's real income, U is an error term assumed to be lognormally distributed, A is a scale parameter to be estimated, and b and c are elasticity parameters to be estimated. Taking logs yields

m_t = a + br_t + cy_t + u_t,

where m = log M, a = log A, r = log R, y = log Y, and u = log U with u being normally distributed. This equation can be estimated using ordinary least squares.

Another economic example is the estimation of a firm's Cobb–Douglas production function, which is the right side of the equation

Q_t=AN_t^{\alpha}K_t^{\beta}U_t,

in which Q is the quantity of output that can be produced per month, N is the number of hours of labor employed in production per month, K is the number of hours of physical capital utilized per month, U is an error term assumed to be lognormally distributed, and A, \alpha, and \beta are parameters to be estimated. Taking logs gives the linear regression equation

q_t = a + \alpha n_t + \beta k_t + u_t

where q = log Q, a=log A, n=log N, k=log K, and u=log U.

Log–log regression can also be used to estimate the fractal dimension of a naturally occurring fractal.

However, going in the other direction – observing that data appears as an approximate line on a log–log scale and concluding that the data follows a power law – is invalid.[2]

In fact, many other functional forms appear approximately linear on the log–log scale, and simply evaluating the goodness of fit of a linear regression on logged data using the coefficient of determination (R2) may be invalid, as the assumptions of the linear regression model, such as Gaussian error, may not be satisfied; in addition, tests of fit of the log–log form may exhibit low statistical power, as these tests may have low likelihood of rejecting power laws in the presence of other true functional forms. While simple log–log plots may be instructive in detecting possible power laws, and have been used dating back to Pareto in the 1890s, validation as a power laws requires more sophisticated statistics.[2]

These graphs are also extremely useful when data are gathered by varying the control variable along an exponential function, in which case the control variable x is more naturally represented on a log scale, so that the data points are evenly spaced, rather than compressed at the low end. The output variable y can either be represented linearly, yielding a lin–log graph (log x, y), or its logarithm can also be taken, yielding the log–log graph (log x, log y).

Bode plot (a graph of the frequency response of a system) is also log-log plot.

See also

References

  1. M. Bourne Graphs on Logarithmic and Semi-Logarithmic Paper (www.intmath.com)
  2. 1 2 Clauset, A.; Shalizi, C. R.; Newman, M. E. J. (2009). "Power-Law Distributions in Empirical Data". SIAM Review 51 (4): 661–703. arXiv:0706.1062. Bibcode:2009SIAMR..51..661C. doi:10.1137/070710111.
This article is issued from Wikipedia - version of the Thursday, April 28, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.