The Central Idea

A classifier uses a feature vector \(\mathbf{x}\) to estimate how plausible each class is, then applies a rule that selects a class.

\[ \boxed{ \mathbf{x} \longrightarrow \text{estimated class probabilities} \longrightarrow \text{predicted class} } \]

Probability estimation and class assignment are related, but they are not the same task. The distinction is essential because different decisions can use the same probability model with different thresholds.

Classification versus Regression

Feature Regression Classification
Response Quantitative value Category or class
Typical output Predicted numerical value Class probability and label
Example Predict a house price Predict whether a payment will be late
Common evaluation Residuals, MAE, RMSE Confusion matrix, log loss, ROC or PR measures

Binary Classification

A binary response has two classes. We often encode the event of interest as \(Y=1\) and the other outcome as \(Y=0\). Examples include

The labels \(0\) and \(1\) are codes. They do not imply that one class is intrinsically more valuable than the other.

The Positive Class

The class designated \(Y=1\) is commonly called the positive class. Positive does not necessarily mean desirable. Fraud, equipment failure, and disease may all be coded positive because they are the events an analysis is designed to detect.

State the positive class before reporting sensitivity, specificity, precision, or recall. These measures change meaning if the class labels are reversed.

Conditional Class Probability

For binary classification, the central statistical quantity is

\[ p(\mathbf{x}) = P(Y=1\mid \mathbf{X}=\mathbf{x}). \]

This is the probability of the positive class among cases having features \(\mathbf{x}\), under the model and the population represented by the data.

\[ 0\le p(\mathbf{x})\le 1. \]

Why Ordinary Linear Regression Is Problematic

If a \(0\)-\(1\) response is modeled by

\[ \hat{p}=b_0+b_1x, \]

the fitted values can fall below \(0\) or above \(1\). The error variance also changes with the probability because

\[ \operatorname{Var}(Y\mid \mathbf{x}) = p(\mathbf{x})[1-p(\mathbf{x})]. \]

A probability model should respect the probability scale and its changing variance.

Odds

A probability \(p\) can be expressed as odds:

\[ \boxed{ \text{odds} = \frac{p}{1-p} }. \]

If \(p=0.75\), the odds are

\[ \frac{0.75}{0.25}=3. \]

This means three expected positive outcomes for every one expected negative outcome among comparable cases. It does not mean a probability of \(3\).

From Odds Back to Probability

If the odds are \(o\), then

\[ \boxed{ p=\frac{o}{1+o} }. \]

Odds of \(4\) correspond to

\[ p=\frac{4}{1+4}=0.80. \]

Log-Odds

The logarithm of the odds is called the logit:

\[ \boxed{ \operatorname{logit}(p) = \log\left(\frac{p}{1-p}\right) }. \]

Probabilities between \(0\) and \(1\) map to log-odds anywhere on the real line. This makes it possible to model log-odds with a linear predictor.

Logistic Regression

Logistic regression assumes that the log-odds of the positive class are linear in the predictors:

\[ \boxed{ \log\left(\frac{p(\mathbf{x})}{1-p(\mathbf{x})}\right) = \beta_0+\beta_1x_1+\cdots+\beta_kx_k }. \]

The right side may also include indicator variables, interactions, polynomial terms, or spline terms. “Linear” refers to the coefficients, not necessarily to a straight relationship on the probability scale.

The Logistic Function

Let

\[ \eta = \beta_0+\beta_1x_1+\cdots+\beta_kx_k. \]

Solving the logit equation for \(p\) gives

\[ \boxed{ p(\mathbf{x}) = \frac{e^\eta}{1+e^\eta} = \frac{1}{1+e^{-\eta}} }. \]

The logistic function converts any real value of \(\eta\) into a valid probability.

The Logistic Curve

An S-shaped logistic probability curve The horizontal axis shows the linear predictor eta, and the vertical axis shows probability. The curve rises from near zero, crosses one half when eta is zero, and approaches one. η p 0.5 0
The logistic function approaches \(0\) and \(1\) without crossing those bounds. Its probability is \(0.5\) when the linear predictor is \(0\).

Interpreting a Logistic Coefficient

Holding all other predictors fixed, a one-unit increase in \(x_j\) changes the log-odds by \(\beta_j\) and multiplies the odds by

\[ \boxed{ e^{\beta_j} }. \]

If \(\beta_j=0.40\), the odds ratio is

\[ e^{0.40}\approx 1.492. \]

The fitted odds are about \(49.2\%\) higher for a one-unit increase in \(x_j\), conditional on the other predictors. This is not a constant \(49.2\)-percentage-point increase in probability.

Probability Changes Are Not Constant

In a one-predictor logistic model, the fitted probability slope is

\[ \frac{dp}{dx} = \beta_1p(1-p). \]

The same coefficient can produce a small probability change near \(p=0\) or \(p=1\) and a larger change near \(p=0.5\). Report probability changes at relevant predictor values when stakeholders think in probabilities.

Indicator Predictors

Suppose \(G=1\) identifies a group and \(G=0\) is the reference group:

\[ \log\left(\frac{p}{1-p}\right) = \beta_0+\beta_1G. \]

Then \(e^{\beta_1}\) is the fitted odds ratio comparing \(G=1\) with \(G=0\). With additional predictors, the comparison holds those predictors fixed.

Interactions and Nonlinear Terms

Classification models can represent changing effects:

\[ \operatorname{logit}(p) = \beta_0+\beta_1x+\beta_2G+\beta_3xG+\beta_4x^2. \]

The interaction allows the association between \(x\) and the log-odds to differ by group. The quadratic term permits curvature. Preserve lower-order terms when including interactions or powers unless a strong structural reason supports a restriction.

Estimating the Coefficients

Logistic regression is ordinarily fitted by maximum likelihood. For observations \(i=1,\ldots,n\), with \(y_i\in\{0,1\}\), the likelihood is

\[ L(\boldsymbol{\beta}) = \prod_{i=1}^{n} p_i^{y_i}(1-p_i)^{1-y_i}. \]

The fitted coefficients make the observed class outcomes as likely as possible under the model. Numerical algorithms perform the optimization.

Log Loss

The average negative log-likelihood is called log loss or cross-entropy:

\[ \boxed{ -\frac{1}{n} \sum_{i=1}^{n} \left[ y_i\log(\hat{p}_i) +(1-y_i)\log(1-\hat{p}_i) \right] }. \]

Smaller values are better. Log loss rewards accurate probabilities and penalizes confident errors severely.

A Worked Probability Model

Suppose a fitted model for late payment uses debt-to-income percentage \(x\):

\[ \operatorname{logit}(\hat{p}) = -4+0.10x. \]

For an applicant with \(x=30\),

\[ \hat{\eta} = -4+0.10(30) = -1, \] \[ \hat{p} = \frac{1}{1+e^1} \approx 0.269. \]

The model estimates a \(26.9\%\) probability of late payment for a case with this predictor value.

Interpret the Worked Coefficient

A one-percentage-point increase in debt-to-income percentage multiplies the fitted odds of late payment by

\[ e^{0.10}\approx 1.105. \]

Thus, the fitted odds increase by about \(10.5\%\) per percentage point, within the model and holding any other included predictors fixed.

The \(50\%\) Probability Point

The logistic probability is \(0.5\) when the linear predictor is \(0\). In the worked model,

\[ -4+0.10x=0, \]

so

\[ x=40. \]

This point describes the model. It is not automatically the best decision threshold or an operational boundary.

From Probability to a Class

A threshold \(t\) converts the fitted probability into a label:

\[ \hat{Y} = \begin{cases} 1, & \hat{p}\ge t,\\ 0, & \hat{p}<t. \end{cases} \]

The conventional threshold \(t=0.5\) is only one option. The appropriate threshold depends on the purpose of the decision and the consequences of each error.

Changing the Threshold

Threshold selection is a decision-policy choice, not a hidden property of the fitted model.

The Confusion Matrix

For a fixed threshold, compare predicted and actual classes:

Prediction Actual positive Actual negative
Predicted positive True positive (TP) False positive (FP)
Predicted negative False negative (FN) True negative (TN)

A Worked Confusion Matrix

Suppose a screening classifier is evaluated on \(1{,}000\) cases:

Prediction Condition present Condition absent Row total
Positive TP \(=72\) FP \(=108\) \(180\)
Negative FN \(=18\) TN \(=802\) \(820\)
Column total \(90\) \(910\) \(1{,}000\)

Sensitivity or Recall

Sensitivity is the proportion of actual positives correctly identified:

\[ \boxed{ \text{sensitivity} = \text{recall} = \frac{\mathrm{TP}}{\mathrm{TP}+\mathrm{FN}} }. \]

In the example,

\[ \frac{72}{72+18} = \frac{72}{90} = 0.80. \]

The classifier detects \(80\%\) of the cases with the condition.

Specificity

Specificity is the proportion of actual negatives correctly identified:

\[ \boxed{ \text{specificity} = \frac{\mathrm{TN}}{\mathrm{TN}+\mathrm{FP}} }. \]

In the example,

\[ \frac{802}{802+108} = \frac{802}{910} \approx 0.881. \]

Precision or Positive Predictive Value

Precision is the proportion of predicted positives that are actually positive:

\[ \boxed{ \text{precision} = \frac{\mathrm{TP}}{\mathrm{TP}+\mathrm{FP}} }. \]

In the example,

\[ \frac{72}{72+108} = \frac{72}{180} = 0.40. \]

Forty percent of positive classifications are true positives in this evaluation population.

Negative Predictive Value

Negative predictive value is the proportion of predicted negatives that are actually negative:

\[ \boxed{ \text{NPV} = \frac{\mathrm{TN}}{\mathrm{TN}+\mathrm{FN}} }. \]

In the example,

\[ \frac{802}{802+18} = \frac{802}{820} \approx 0.978. \]

Accuracy

Accuracy is the overall proportion classified correctly:

\[ \boxed{ \text{accuracy} = \frac{\mathrm{TP}+\mathrm{TN}} {\mathrm{TP}+\mathrm{FP}+\mathrm{FN}+\mathrm{TN}} }. \]

In the example,

\[ \frac{72+802}{1000} = 0.874. \]

Why Accuracy Can Mislead

Only \(90\) of the \(1{,}000\) cases in the example are positive, so the positive prevalence is \(9\%\). A rule that predicts every case negative would have

\[ \frac{910}{1000}=91\% \]

accuracy while detecting none of the positive cases. Always compare a classifier with simple baselines and examine class-specific performance.

False-Positive and False-Negative Rates

\[ \text{false-positive rate} = \frac{\mathrm{FP}}{\mathrm{FP}+\mathrm{TN}} = 1-\text{specificity}, \] \[ \text{false-negative rate} = \frac{\mathrm{FN}}{\mathrm{FN}+\mathrm{TP}} = 1-\text{sensitivity}. \]

In the worked example, the false-positive rate is about \(0.119\), and the false-negative rate is \(0.20\).

Balanced Accuracy

Balanced accuracy gives equal weight to performance in the two actual classes:

\[ \boxed{ \text{balanced accuracy} = \frac{\text{sensitivity}+\text{specificity}}{2} }. \]

For the worked example,

\[ \frac{0.800+0.8813}{2} \approx 0.841. \]

The \(F_1\) Score

The \(F_1\) score is the harmonic mean of precision and recall:

\[ \boxed{ F_1 = 2\frac{(\text{precision})(\text{recall})} {\text{precision}+\text{recall}} = \frac{2\mathrm{TP}} {2\mathrm{TP}+\mathrm{FP}+\mathrm{FN}} }. \]

For the worked example,

\[ F_1 = \frac{144}{144+108+18} = \frac{144}{270} \approx 0.533. \]

\(F_1\) omits true negatives, so it should not be treated as a universal performance measure.

Conditioning Directions Matter

Sensitivity and precision answer different questions:

\[ \text{sensitivity} = P(\widehat{Y}=1\mid Y=1), \] \[ \text{precision} = P(Y=1\mid \widehat{Y}=1). \]

Reversing the condition reverses the question. High sensitivity does not guarantee high precision.

Prevalence Affects Predictive Values

Sensitivity and specificity condition on the actual class. Precision and negative predictive value condition on the predicted class and therefore change when class prevalence changes.

\[ P(Y=1\mid \widehat{Y}=1) = \frac{ P(\widehat{Y}=1\mid Y=1)P(Y=1) }{ P(\widehat{Y}=1) }. \]

A screening rule can have the same sensitivity and specificity in two populations but a lower precision where the condition is rarer.

Choose Metrics from the Decision

Operational priority Metric to emphasize
Find most actual positives Sensitivity or recall
Limit false alarms Specificity or false-positive rate
Make positive flags trustworthy Precision
Evaluate probability quality Log loss, Brier score, and calibration
Compare ranking across thresholds ROC or precision–recall analysis

Decision Costs

Let \(C_{\mathrm{FP}}\) be the cost of a false positive and \(C_{\mathrm{FN}}\) the cost of a false negative. A decision policy can be evaluated by expected cost:

\[ \text{expected cost} = C_{\mathrm{FP}}P(\mathrm{FP}) + C_{\mathrm{FN}}P(\mathrm{FN}), \]

with additional terms when correct decisions also have material benefits or costs. Cost assumptions should be documented and sensitivity-tested.

A Cost-Based Threshold

In a simplified setting with correct calibrated probabilities and zero cost for correct decisions, predict positive when

\[ C_{\mathrm{FP}}(1-p) \le C_{\mathrm{FN}}p. \]

Solving gives

\[ \boxed{ p\ge \frac{C_{\mathrm{FP}}} {C_{\mathrm{FP}}+C_{\mathrm{FN}}} }. \]

If a false negative costs four times as much as a false positive, the simplified threshold is \(1/(1+4)=0.20\), not \(0.50\).

The ROC Curve

A receiver operating characteristic curve plots

\[ \text{true-positive rate} \quad\text{against}\quad \text{false-positive rate} \]

as the classification threshold varies. Each point represents a different sensitivity–specificity tradeoff. The curve evaluates ranking behavior across thresholds, not the quality of one selected operating threshold.

Area under the ROC Curve

The ROC area, or AUC, can be interpreted as the probability that a randomly selected positive case receives a higher score than a randomly selected negative case, with ties handled appropriately.

Precision–Recall Curves

A precision–recall curve plots precision against recall as the threshold changes. It focuses attention on positive predictions and can be especially informative when the positive class is rare.

Unlike ROC analysis, the baseline precision depends on positive-class prevalence. Comparisons across populations therefore require care.

Ranking Is Not Calibration

A model has strong discrimination if it tends to rank positive cases above negative cases. A model is calibrated if predicted probabilities agree with observed frequencies.

A model can rank cases perfectly while reporting probabilities that are systematically too high or too low. Evaluate both capabilities when the probabilities themselves guide decisions.

Calibration

Informally, a calibrated model should satisfy

\[ P(Y=1\mid \hat{p}\approx q) \approx q. \]

Among cases assigned probabilities near \(0.70\), about \(70\%\) should be positive over repeated comparable cases. Calibration plots compare observed frequencies with predicted probabilities across ranges of risk.

The Brier Score

For binary outcomes, the Brier score is the mean squared probability error:

\[ \boxed{ \text{Brier score} = \frac{1}{n} \sum_{i=1}^{n} (\hat{p}_i-y_i)^2 }. \]

Smaller values are better. Unlike label accuracy, the Brier score uses the full probability prediction and reflects both calibration and resolution.

Training and Test Performance

Performance calculated on the same data used to fit the model is usually optimistic. A credible evaluation uses observations not used to estimate the model:

\[ \boxed{ \text{fit on training data} \quad\longrightarrow\quad \text{evaluate on unseen data} }. \]

All feature selection, transformations, imputation, scaling, and threshold tuning must be learned without using the final test outcomes.

Data Leakage

Leakage occurs when information unavailable at prediction time, or information from the evaluation outcomes, enters model training. Examples include

Leakage can produce excellent reported metrics and poor real-world performance.

Respect the Data Structure

Random row splitting is not always appropriate. Use time-based splits for future prediction, group-based splits when observations from the same person or organization are related, and location-based evaluation when geographic transfer matters.

The evaluation design should reproduce the intended deployment boundary.

Class Imbalance

When one class is rare, a model may achieve high accuracy by favoring the common class. Responses include

Artificially changing class proportions can distort probability calibration unless the shift is explicitly corrected.

Complete and Quasi-Complete Separation

In logistic regression, separation occurs when a predictor or combination of predictors perfectly distinguishes the observed classes. Ordinary maximum-likelihood coefficient estimates may then grow without a finite limit.

Diagnose separation rather than interpreting enormous coefficients literally. Penalized estimation, revised predictors, or additional information may be appropriate.

Other Classification Methods

Method Core idea Notable tradeoff
Logistic regression Model class log-odds Interpretable but may need nonlinear features
Nearest neighbors Use nearby labeled cases Flexible but scale- and dimension-sensitive
Decision tree Split feature space into regions Readable but individual trees can be unstable
Random forest Average many randomized trees Strong flexibility with less direct interpretation
Support vector machine Find a separating margin Effective boundaries but probabilities need extra work

Classification Boundaries

A decision boundary is the set of feature values where the selected class changes. For a logistic model with threshold \(0.5\),

\[ \beta_0+\beta_1x_1+\beta_2x_2=0 \]

is a straight boundary in the \(x_1,x_2\) plane. Interactions, polynomial terms, splines, kernels, and trees can create curved or segmented boundaries.

Multiclass Classification

A multiclass response has \(K>2\) possible classes. The model estimates

\[ p_k(\mathbf{x}) = P(Y=k\mid \mathbf{x}), \qquad k=1,\ldots,K, \]

subject to

\[ p_k(\mathbf{x})\ge 0, \qquad \sum_{k=1}^{K}p_k(\mathbf{x})=1. \]

The highest-probability class may be selected, but costs or abstention rules can justify another policy.

Softmax Probabilities

Multinomial logistic regression assigns each class a score \(\eta_k\) and converts the scores to probabilities using

\[ \boxed{ p_k = \frac{e^{\eta_k}} {\sum_{j=1}^{K}e^{\eta_j}} }. \]

The probabilities are nonnegative and sum to \(1\). One class is commonly treated as a reference for identifiable coefficient interpretation.

Multiclass Evaluation

A \(K\times K\) confusion matrix records actual and predicted classes. Class-specific precision and recall can be combined in several ways:

State the averaging convention because the same model can receive materially different summaries.

Ordinal Classification

Some categories have a natural order, such as low, medium, and high. Treating them as unrelated classes discards that order; treating them as equally spaced numbers adds an assumption that may be false.

Ordinal regression models cumulative probabilities while respecting the ranking of the categories.

Abstention and Human Review

A system need not force a label for every case. It can defer uncertain or high-impact cases:

\[ \hat{Y} = \begin{cases} 0, & \hat{p}<t_L,\\ \text{review}, & t_L\le\hat{p}<t_U,\\ 1, & \hat{p}\ge t_U. \end{cases} \]

A review option can improve risk control, but its workload, delay, and reviewer consistency must be measured.

Fairness Requires More Than One Metric

Performance may differ across groups because of sample composition, measurement quality, historical decisions, prevalence, or model behavior. Evaluate relevant error rates, calibration, coverage, and downstream effects by group.

Fairness criteria can conflict mathematically when group prevalences differ. The appropriate governance standard depends on the application, rights at stake, and legal and ethical requirements.

Association Is Not Causation

A classifier predicts an observed label from available features. A coefficient, feature importance, or accurate prediction does not establish that changing a predictor would change the outcome.

Causal decisions require a defined intervention, suitable study design, and defensible assumptions beyond predictive accuracy.

Dataset Shift

A model can degrade when deployment differs from development:

Predeployment validation cannot eliminate the need for postdeployment monitoring.

Deployment Monitoring

A production classification system should track

Define alert thresholds, ownership, escalation paths, and retraining criteria before deployment.

A Practical Classification Workflow

  1. Define the outcome, positive class, prediction time, and unit of analysis.
  2. Specify the decision, stakeholders, error costs, and constraints.
  3. Audit label quality, feature availability, missingness, and class balance.
  4. Create training, validation, and test partitions that match deployment.
  5. Build a transparent baseline before adding complexity.
  6. Fit models without leaking evaluation information.
  7. Evaluate discrimination, calibration, class-specific errors, and uncertainty.
  8. Select and document the threshold or abstention policy.
  9. Test robustness across groups, time periods, and plausible shifts.
  10. Monitor the deployed system and maintain a rollback path.

Reporting Checklist

Common Error: Treating \(0.5\) as Universal

A \(0.5\) threshold assumes a particular decision structure. It can be unsuitable when classes are imbalanced, error costs differ, resources are constrained, or review is available.

Common Error: Reporting Only Accuracy

Accuracy hides which class is being missed and can reward trivial majority-class predictions. Report metrics aligned with the operational decision and include the confusion matrix at the selected threshold.

Common Error: Calling Scores Probabilities

Some algorithms produce ranking scores that are not calibrated probabilities. A number between \(0\) and \(1\) is not automatically a trustworthy risk estimate. Calibration must be evaluated on relevant unseen data.

Common Error: Tuning on the Test Set

Repeatedly selecting models or thresholds based on test performance turns the test set into training information. Reserve the final test set for one honest evaluation after modeling decisions are fixed.

Practice 1: Convert Probability and Odds

A model assigns probability \(p=0.20\) to the positive class. Find the odds. Then convert odds of \(3\) back to a probability.

Practice 2: Calculate a Logistic Probability

For

\[ \operatorname{logit}(\hat{p})=-3+0.5x, \]

find \(\hat{p}\) when \(x=4\).

Practice 3: Interpret an Odds Ratio

A logistic coefficient is \(\beta=0.25\). Interpret \(e^\beta\).

Practice 4: Change the Threshold

What generally happens to sensitivity and the false-positive rate when the decision threshold is lowered?

Practice 5: Calculate Metrics

A classifier produces \(\mathrm{TP}=45\), \(\mathrm{FP}=15\), \(\mathrm{FN}=5\), and \(\mathrm{TN}=135\). Find sensitivity, specificity, precision, and accuracy.

Practice 6: Diagnose Accuracy

In a dataset, \(99\%\) of cases are negative. A classifier predicts every case negative. Explain why its \(99\%\) accuracy is inadequate evidence of usefulness.

Practice 7: Select a Cost Threshold

Under the simplified cost model, a false positive costs \(2\) units and a false negative costs \(8\) units. Find the probability threshold.

Practice 8: Distinguish Ranking and Calibration

A model ranks nearly every positive case above every negative case, but cases assigned probability \(0.80\) are positive only \(50\%\) of the time. What is strong, and what is weak?

Practice 9: Identify Leakage

An analyst imputes missing values using the full dataset and then creates training and test sets. What is the problem?

Practice 10: Evaluate a Causal Claim

A classifier accurately predicts employee turnover from commute length. A manager concludes that shortening commutes will reduce turnover by the amount implied by the model. What is wrong with this conclusion?

Practice Answers

Practice 1

\[ \text{odds} = \frac{0.20}{0.80} = 0.25. \]

Odds of \(3\) correspond to

\[ p=\frac{3}{1+3}=0.75. \]

Practice 2

\[ \hat{\eta} = -3+0.5(4) = -1, \] \[ \hat{p} = \frac{1}{1+e^1} \approx 0.269. \]

Practice 3

\[ e^{0.25}\approx 1.284. \]

Holding the other included predictors fixed, a one-unit increase in the predictor multiplies the fitted odds by about \(1.284\), an increase of about \(28.4\%\). It is not a constant probability increase.

Practice 4

More cases are labeled positive. Sensitivity generally rises, and the false-positive rate generally also rises.

Practice 5

\[ \text{sensitivity} = \frac{45}{45+5} = 0.90, \] \[ \text{specificity} = \frac{135}{135+15} = 0.90, \] \[ \text{precision} = \frac{45}{45+15} = 0.75, \] \[ \text{accuracy} = \frac{45+135}{200} = 0.90. \]

Practice 6

The classifier has zero sensitivity and detects none of the positive cases. Its accuracy merely reflects the class imbalance and matches the trivial majority-class baseline.

Practice 7

\[ t = \frac{C_{\mathrm{FP}}} {C_{\mathrm{FP}}+C_{\mathrm{FN}}} = \frac{2}{2+8} = 0.20. \]

Practice 8

Discrimination or ranking is strong, but calibration is poor because the stated probabilities exceed the observed event frequency.

Practice 9

The imputation values use information from the test data. Imputation should be fitted on the training data and then applied unchanged to the test data, ordinarily within a complete modeling pipeline.

Practice 10

Predictive association does not identify the causal effect of changing commute length. Confounding, selection, reverse causation, and the feasibility of the intervention require separate causal analysis.

What to Remember

Where We Go Next

Classification emphasizes categorical outcomes and probability-based decisions. For quantitative outcomes, a point prediction alone does not communicate how uncertain a future response may be.

The next lesson, Prediction Intervals, develops intervals designed to contain future individual outcomes and distinguishes them from confidence intervals for a mean response.