From Scientific Expectations to Regression Forms
A scientific description often tells us more than whether a response should increase or decrease. It may suggest that an effect begins only after a threshold, that each additional unit has a smaller effect than the previous one, or that the response decreases and then increases. These expectations can guide the construction of a small, interpretable set of candidate regression curves before any coefficients are estimated.
A candidate curve is a proposed form for the mean response as a function of a predictor. It is not a claim that every observation will lie on that curve, nor does choosing a form prove a scientific mechanism. The aim here is to translate a stated expectation into a mathematical shape, specify the domain where that shape makes sense, and interpret what its terms imply.
For the examples, suppose observations follow \(Y_i=m(x_i)+\varepsilon_i\), where \(x_i\) is the predictor value for observation \(i\), and \(\varepsilon_i\) is its error around the mean function. Assume the errors are independent, have conditional mean zero, and have a common finite variance \(\sigma^2\): \(E(\varepsilon_i\mid x_i)=0\) and \(\operatorname{Var}(\varepsilon_i\mid x_i)=\sigma^2\). These assumptions define a familiar regression setting; they do not determine the shape of \(m(x)\). The scientific expectation guides the candidate mean function, while the error assumptions describe variation around it. The examples below are hypothetical illustrations, not results from real studies.
Recognizing the Shape Implied by Subject Knowledge
Begin by describing the expected response in ordinary language: Does the response change gradually, change rate at a particular point, level off, or reverse direction? Then choose a function whose behavior matches that description over the predictor range that matters. The coefficients give the form flexibility, but the form itself imposes important restrictions.
| Scientific expectation | Possible candidate form | Important implication |
|---|---|---|
| A change in rate at a threshold \(c\) | \(m(x)=\beta_0+\beta_1x+\beta_2(x-c)_+\) | Continuous at \(c\); the slope changes there |
| Increasing response with diminishing increments | \(m(x)=\beta_0+\beta_1\log(x/x_0)\) | Requires \(x>0\); for \(\beta_1>0\), increases at a decreasing rate |
| A response that decreases and then increases | \(m(x)=\beta_0+\beta_1x+\beta_2x^2\), with \(\beta_2>0\) | Has a minimum at its turning point, if that point is in the relevant range |
Here \(\beta_0,\beta_1,\beta_2\) are coefficients and \(x_0\) is a positive reference value with the same units as \(x\). The hinge term is defined by \((x-c)_+=\max(x-c,0)\). A logarithm in a regression expression should have a dimensionless argument; the ratio \(x/x_0\) provides one. These forms are not interchangeable descriptions of “nonlinearity”: each encodes a different pattern of change.
Worked Example: A Response That Changes Rate at a Threshold
Suppose a hypothetical production process is expected to improve gradually with operating time, but an additional change in rate may begin after 6 hours. Let \(x\) be operating time in hours and let the mean response \(m(x)\) be output units per batch. Consider the candidate
For \(x<6\), the positive-part term is zero, so \(m(x)=8+0.4x\). The slope below the threshold is \(0.4\) output units per batch per hour. For \(x>6\), the hinge term equals \(x-6\), giving \(m(x)=8+0.4x+1.1(x-6)\); the slope is then \(0.4+1.1=1.5\) output units per batch per hour.
At \(x=4\), \(m(4)=8+0.4(4)+1.1(0)=9.6\) output units per batch. At the threshold, \(m(6)=8+0.4(6)+1.1(0)=10.4\). At \(x=8\), \(m(8)=8+0.4(8)+1.1(2)=13.4\). The curve is continuous at 6 hours, but it is not differentiable there: its left-hand slope is \(0.4\) and its right-hand slope is \(1.5\). Thus neither side’s slope should be described as the slope exactly at the threshold.
The hinge construction separates a baseline trend from a change in trend. The coefficient \(\beta_2\) is the change in slope at \(c\), not the slope after \(c\) by itself. If the expected response instead starts only after a threshold, subject knowledge might support a different form, such as a hinge term without the baseline linear trend. The choice should reflect the stated expectation rather than convenience alone.
Worked Example: Diminishing Increases on a Positive Predictor Scale
Suppose a hypothetical process uses \(x\) litres of an input, and a measured response \(m(x)\) is in response units. Subject knowledge suggests that increasing input raises the mean response, but with smaller gains at higher input levels. A candidate is
where \(\ln\) denotes the natural logarithm and 10 litres is the reference amount \(x_0\). At 10 litres, \(\ln(10/10)=0\), so \(m(10)=4\) response units. At 40 litres,
The mean response at 40 litres is therefore about 6.495 response units under this candidate. To check its shape, differentiate with respect to \(x\): \(m'(x)=1.8/x\). This derivative is positive for \(x>0\) and gets smaller as \(x\) increases. For example, it is \(0.18\) response units per litre at 10 litres and \(0.045\) at 40 litres. The curve increases while its local rate of increase diminishes. It is defined only for positive input; it cannot be used at zero litres.
The logarithmic form is one way to express diminishing increases, not a universal model of saturation. It keeps increasing for every positive \(x\), rather than approaching a finite upper limit. If scientific knowledge requires a strict upper bound or a plateau, that expectation calls for a candidate with a corresponding limiting behavior. The form should not be assigned a stronger interpretation than its mathematics supports.
Encoding a U-Shaped Expectation
A quadratic mean function can express a response that decreases and then increases, provided its curvature and turning point are appropriate. For \(m(x)=\beta_0+\beta_1x+\beta_2x^2\), the coefficient of \(x^2\) controls the direction of curvature. If \(\beta_2>0\), the curve bends upward. Completing the square shows where its minimum occurs:
Because a square is nonnegative and \(\beta_2>0\), the first term is smallest when \(x=-\beta_1/(2\beta_2)\); that is the quadratic’s minimum. A U-shaped curve over the predictor values of interest requires more than a positive quadratic coefficient: the minimum must also fall within that range. Otherwise the curve may appear only increasing or only decreasing over the observed region.
Worked Example: A Quadratic Candidate for a U-Shaped Response
Imagine a hypothetical setting in which a response is expected to be lowest near a moderate value of a predictor, and higher at values both below and above it. Let \(x\) be a centered predictor in index units and let \(m(x)\) be the mean response in score units. Consider
Here \(\beta_2=0.5>0\), so the quadratic bends upward. Its minimum is at
The corresponding mean is \(m(2)=20-2(2)+0.5(2^2)=18\) score units. At \(x=0\), \(m(0)=20\); at \(x=4\), \(m(4)=20-8+8=20\). The equal values on either side of 2 illustrate the symmetry of this quadratic, and both exceed its minimum. This candidate expresses a U shape over a range containing 2; it does not imply a U shape over every possible predictor range.
A Practical Construction Process
The following sequence keeps the curve tied to subject knowledge and makes its implications explicit before coefficients are estimated.
Describe how the mean response is expected to change, including any threshold, diminishing increments, or reversal in direction.
Use a hinge term for a change in slope, a logarithmic term for one kind of diminishing increase, or a quadratic with positive curvature for a U-shaped expectation.
Confirm that the expression is defined and that its implied behavior is meaningful over the predictor values where it will be used.
Identify slopes, slope changes, or turning points, and state their units where appropriate. Do not treat a coefficient as having the same interpretation in every form.
It can be useful to retain two or three scientifically plausible forms when the verbal expectation does not identify a unique function. For instance, “diminishing gains” does not specify exactly how quickly gains decline. Candidate forms should differ for a reason that can be explained, rather than being collected simply because they are mathematically available. At this stage, the task is to formulate the possibilities—not to declare which curve fits best.
Common Mistakes and Why Formulation Matters
- Using a threshold coefficient as the post-threshold slope. In the hinge model, the post-threshold slope is the baseline slope plus the hinge coefficient. At the threshold itself, the curve is continuous but may have no single derivative.
- Calling every increasing concave curve “saturation.” The logarithmic form has diminishing increases but does not approach a finite ceiling. Its actual limiting behavior matters.
- Assuming a positive quadratic coefficient guarantees a U shape in the data range. It guarantees upward curvature; the minimum must also be within the relevant predictor range for the curve to decrease and then increase there.
- Ignoring domains and units. A logarithmic predictor must be positive, and a log argument should be dimensionless. A slope’s units are response units per predictor unit.
- Reading a candidate as a causal claim. A mean-function shape describes an association under a regression model. By itself, it does not establish that changing the predictor causes the response to change.
Check Your Understanding
For each question, explain the shape or interpretation in context rather than naming a formula alone.
- In \(m(x)=a+bx+d(x-c)_+\), what are the slopes below and above \(c\), and what can happen to the derivative at \(c\)?
- For \(m(x)=a+b\ln(x/x_0)\) with \(b>0\), state the domain and explain why increases diminish as \(x\) grows.
- A quadratic has positive coefficient on \(x^2\), but its minimum lies below all predictor values of interest. Does it show a U shape over that range? Explain.
- Why is the ratio \(x/x_0\) preferable to taking the logarithm of a quantity stated only in physical units?
- Give one reason that a scientifically plausible candidate curve does not, by itself, establish a causal effect.