Porisma

Worked solution

Compare pollution models and audit a damaged regression fragment

(a) Traffic engineers compare traffic flow on a busy main road with air pollution at a nearby monitoring station. Each hour, sensors report the average traffic flow xx (vehicles per hour) for the preceding hour; the station supplies a pollution reading yy in a suitable unit, with higher readings indicating more pollution. A random sample of eight hours gives:

xxyy
17961.0
19182.2
21203.5
23154.2
23684.3
24204.5
25884.9
28505.5

(i) Draw the scatter diagram, labelling the axes.

Two proposed models are y=a+bxy=a+bx and y2=c+dxy^2=c+dx.

(ii) Find the correlation coefficient between (a) xx and yy, and (b) xx and y2y^2.

(iii) Use (i) and (ii) to explain which model is better.

(iv) Find a suitable regression equation and use it to estimate yy when x=2000x=2000.

(v) Comment on the validity of the headline “Heavy traffic causes air pollution” in light of these results.

(b) A partly destroyed research paper retains these statements: the mean of xx is 4.4; the regression of yy on xx is y=2.5x+3.8y=2.5x+3.8; the regression of xx on yy begins x=1.5y−…x=1.5y-\ldots. Calculate the missing constant at the end of the second equation.

[2017/JJC/II/10]

CJC · Discussion · Question 6

(a)(i) Preserve the hourly paired observations

Answer

The eight observations rise with decreasing apparent slope.

Plot average traffic flow horizontally and the corresponding hourly pollution reading vertically. The original response shows a curved increasing pattern; the independently drawn second panel shows the squared-response transformation considered below. The sampling unit is an hour, not eight different stations or eight vehicles.

How did this part go?

(a)(ii)(a) Correlate the original variables

Answer

rx,y=0.959r_{x,y}=0.959 (3 s.f.).

Using all eight pairs gives 0.9588684475…0.9588684475\ldots, matching the printed 0.9590.959. Despite the large positive value, inspect the curved scatter before selecting a straight response model.

How did this part go?

(a)(ii)(b) Correlate flow with squared pollution

Answer

rx,y2=0.995r_{x,y^2}=0.995 (3 s.f.).

Square each pollution reading while retaining its matching traffic flow. The coefficient is 0.9951588471…0.9951588471\ldots, matching the printed 0.9950.995. This transformation is of the response, not of flow.

How did this part go?

(a)(iii) Compare the two proposed forms

Answer

y2=c+dxy^2=c+dx is the better of these two sample models.

The original scatter curves and flattens, while the squared-response scatter is close to a straight rising line and has a correlation closer to 1. Together these support the transformed candidate. This compares fit to this sample; it does not establish a causal mechanism or a universally valid law.

How did this part go?

(a)(iv) Fit squared pollution on flow

Answer

Predicted reading y=2.8y=2.8 (2 s.f.).

The full independently calculated regression is y2=c+dxy^2=c+dx, with

c=−47.98514143459646,d=0.027897204434110024.\begin{aligned}c&=-47.98514143459646,\\d&=0.027897204434110024.\end{aligned}

matching the printed rounded form y2=0.0279x−48.0y^2=0.0279x-48.0. With the unrounded coefficients,

y=c+2000d=2.794506653…≈2.8.\begin{aligned}&y\\ & =\sqrt{c+2000d}\\ & =2.794506653\ldots\approx2.8.\end{aligned}

Take the positive root, consistent with the supplied positive pollution scale and observations. 2000 vehicles per hour is within 1796–2850; this is interpolation. The source prints 2.8 without a precision instruction, so that displayed result is preserved; 2 s.f. states its precision rather than inventing a printed requirement.

How did this part go?

(a)(v) Evaluate the headline

Answer

The results show association; they do not by themselves establish the causal headline.

The sampled hourly observations do not control possible factors such as weather or other pollution sources. Those are possible unmeasured influences, not data supplied by the question. A transformed coefficient close to 1 does not establish that changing traffic alone causes the observed pollution change.

How did this part go?

(b) Compute the formal constant and check consistency

Answer

The missing number is 17.817.8, giving formally x=1.5y−17.8x=1.5y-17.8; the supplied slopes are inconsistent.

Both genuine regression lines pass through the mean point. From the first equation,

yˉ=2.5(4.4)+3.8=14.8.\begin{aligned}&\bar y\\ & =2.5(4.4)+3.8\\ & =14.8.\end{aligned}

Write the fragment as x=1.5y−kx=1.5y-k; then

4.4=1.5(14.8)−k⇒k=17.8,4.4=1.5(14.8)-k\quad\Rightarrow\quad k=17.8,

which reproduces the printed missing-number answer. However, for a nondegenerate paired dataset, the product of the two regression slopes equals r2≤1r^2\le1. Here 2.5(1.5)=3.75>12.5(1.5)=3.75>1, so these cannot both be the least-squares regressions of the same data. The formal mean-point completion is retained, but the source fragment is mathematically inconsistent as stated. No corrected slope or replacement dataset is invented.

How did this part go?

Compare your sketch

Expand diagram

Use the scrollbars or arrow keys to inspect details.

All eight traffic/pollution pairs 1796/1,1918/2.2,2120/3.5,2315/4.2,2368/4.3,2420/4.5,2588/4.9,2850/5.5. Second panel retains the same pairs with response squared and full-precision fit y²−47.9851414346+0.0278972044341x. Positive y branch used for the prediction; no extra measurements invented.
All eight traffic/pollution pairs 1796/1,1918/2.2,2120/3.5,2315/4.2,2368/4.3,2420/4.5,2588/4.9,2850/5.5. Second panel retains the same pairs with response squared and full-precision fit y²−47.9851414346+0.0278972044341x. Positive y branch used for the prediction; no extra measurements invented.
All eight traffic/pollution pairs 1796/1,1918/2.2,2120/3.5,2315/4.2,2368/4.3,2420/4.5,2588/4.9,2850/5.5. Second panel retains the same pairs with response squared and full-precision fit y²−47.9851414346+0.0278972044341x. Positive y branch used for the prediction; no extra measurements invented.