- 96 Views
- Uploaded on
- Presentation posted in: General

Predictions from the regression model

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Predictions from the regression model

(Session 09)

At the end of this session, you will be able to

- assess the regression model in terms of its appropriateness as a prediction equation
- distinguish between two types of possible predictions and why their precision differs
- use software to produce predictors and corresponding confidence intervals
- understand dangers of predicting beyond values used in developing the predictor

Consider a health survey of a random sample of 20 districts in which one variable of interest was the perinatal mortality rate, i.e. no. of deaths in first 30 days after birth per 1000 live births;

In all districts, the number of health centres per 1000 households is known.

Suppose a health official wants to predict perinatal mortality in districts not covered by the survey.

-----------------------------------------

mort | Coef. Std.Err. t P>|t|

----------+------------------------------

hcentres | -1.9139 .4071 -4.70 0.000

const. | 38.111 4.703 8.10 0.000

-----------------------------------------

Adjusted R2 = 52.6%

The equation of the line is:

mort = 38.11 - 1.914 (hcentres)

- There are two types of predictions possible:
- What would be the average mortality rate in districts where there are only 9 health centres per 1000 households?
- A particular district has 9 health centres per 1000 hhs. What would be the predicted mortality rate in this district?
- The answer to both predictions are:
- mortality = 38.11 - 1.914(9)= 20.88
- But would their standard errors be the same?

The standard error for (1), call it “sepred”, can be obtained with statistical software.

The standard error for (2) is the square root of:

(sepred)2 + (residual M.S. in anova)

The addition of the RMS is to account for variation of individual values about their mean.

Note: For the case where there is only 1 regressor variable, a simple formula does exist (see final slide).

The standard error can be used to determine (say) 95% confidence intervals for each type of prediction.

Because individual predictions have a larger standard error, the corresponding 95% C.I. will be wider for such predictions.

The graphs below show the difference in precision (Stata output)

It is clear from above graphs that the confidence bands widen as predictions are made away from the mean of the x values.

This is also true of predictions made from results of a multiple linear regression model.

Hence extrapolating to make predictions beyond x’s used in the model fitting can be dangerous since they would be very imprecise!

One aim of the Adult Morbidity and Mortality Project (AMMP) in Tanzania was to investigate how health outcomes varied across different levels of poverty.

All households in several sentinel sites were monitored over many years…

But, initially there was no measure of income poverty to address above objective.

Question: Can a prediction model be developed to give a poverty measure using socio-economic variables?

------------------------------------------

lnexpdf | Coef. Std.Err. t P>|t|

-------------+----------------------------

hhsize | -.22667 .01213 -18.69 0.000

hhsize2 | .00956 .00085 11.24 0.000

water | .11726 .03007 3.90 0.000

qmeat | .11932 .0100 11.93 0.000

qmilk | .03158 .00601 5.26 0.000

iron | .20931 .03164 6.62 0.000

table | .18927 .03559 5.32 0.000

wheatf | .29468 .03232 9.12 0.000

seeds | .18994 .05095 3.73 0.000

num_meal | .15502 .03265 4.75 0.000

const. | 9.1540 .0930 98.41 0.000

------------------------------------------

lnexpdf (y) = 9.154 – 0.227(hhsize)

+ 0.0096(hhsize)2 + 0.117(water)

+ 0.119(qmeat) + 0.032(qmilk)

+ 0.209(iron) +0.189(table)

+ 0.295(wheatf) + 0.190(seeds)

+0.155(num_meal)

R2=58%. A plot of predicted values vs actual also indicates closeness of prediction.

terciles | terciles of prediction|of lnexpdf| Lowest Middle Highest | Total

----------+-------------------------+------

Lowest | 248 80 13 | 341

Middle | 95 182 74 | 351

Highest | 9 86 262 | 357

----------+-------------------------+------

Total | 352 348 349 | 1049

----------+-------------------------+------

Percent correctly classified

= [(248 + 182 + 262)/1049]*100 = 66%

---------------+-------------------+------

If really below| Below poverty line|

basic needs | on predictions |

poverty line? | No Yes | Total

---------------+-------------------+------

No | 838 42 | 880

Yes | 99 70 | 169

---------------+-------------------+------

Total | 937 112 | 1,049

---------------+-------------------+------

Percent correctly classified as below line

= [(838 + 70)/1049]*100 = 86.6%

The answer is “yes”.

Although at a household level, the predictions were not so good, health outcomes like mortality rates, percent immunised, etc., were measured in AMMP at community level.

The predictions, averaged over HHs in a community were satisfactory, e.g. further analysis showed that the average prediction percent error was under 6%.

- Select carefully a set of potential explanatory variables
- Decide the direction of influence expected (+ or – for regression coefficients)
- Use a suitable model selection approach to identify the subset of variables contributing significantly to explaining the variation in the key response(y)
- Conduct a residual analysis and take remedial action if problems occur

- Check whether the final model makes sense in terms of the explanatory variables included and the signs of their coefficients
- Consider whether there are other latent variables which may be non-measurable but influential! Conclusions should mention these!
- Examine how well the prediction performs, studying also the precision of the prediction
- Give the form of the model equation and corresponding conclusions

If s2 represents the Residual Mean Square in the anova, the std. error for the predicted mean =

Std. error for a predicted individual value =

Practical work follows to ensure learning objectives are achieved…