topic Re: Non normal distribution in regression in Statistical Procedures

Non normal distribution in regression

Babloo — Mon, 25 Jan 2016 10:45:45 GMT

May I request someone to shed some light on stattistical test to be conducted when erros in regression don't follow the normal distribution?

Re: Non normal distribution in regression

Rick_SAS — Mon, 25 Jan 2016 11:01:40 GMT

Look at the documentation for the GENMOD procedure, which includes sections about Goodness-of-Fit tests and related statistics. The doc for PROC GENMOD also explain estimates and contrasts.

If you provide more information about your model, more can be said.

Re: Non normal distribution in regression

Babloo — Mon, 25 Jan 2016 11:31:44 GMT

Thank you for your response Rick.

It is a general question which I came across , hence anticipating the simple answer in layman's term rather than bookish languague.

Re: Non normal distribution in regression

Rick_SAS — Mon, 25 Jan 2016 11:59:03 GMT

You might enjoy this graphical comparison of the assumptions for error distributions in linear and nonlinear models:

Re: Non normal distribution in regression

Babloo — Mon, 25 Jan 2016 12:14:54 GMT

Thanks again Rick.

So can I assume that answer to my question is 'proc genmod'?

I've also an another novice question - I know to find whether data is following normal distribution or not , but I don't know how to find whether error is following normal distribution or not. May I request you to guide me on this?

Re: Non normal distribution in regression

Rick_SAS — Mon, 25 Jan 2016 12:41:26 GMT

PROC GENMOD is a good place to start for fitting models of this type. There are alternatives, especially if you think the errors are correlated (as in a time series), but I don't want to overwhelm you with too many options.

I recommend that you do an internet search for "SAS" and "regression diagnostics" or "diagnostic plots". This is a deep area that is worth learning about.

A (very) short answer is that when the response variable is contnuous, you can examine the error distribution by fitting a model and then plotting the distribution of the raw residuals. Most SAS procedures, including GENMOD, have an OUTPUT statement that enables you to write the residual values to a data set. The simplest plot is a histogram of the residuals. Does the histogram look approximately "bell shaped"?

You can also plot the raw residuals versus each of the explanatory variables. If any of the plots look "fan shaped" (the size of the residuals depend on an X), that indicates that the model is not capturing the variation in the data. If so, many practitioners try to fit a more sophisticated model.

Re: Non normal distribution in regression

Reeza — Mon, 25 Jan 2016 15:25:06 GMT

Error is data...you should be able to isolate your error terms if you run a regression model.

@Babloo wrote:

I know to find whether data is following normal distribution or not , but I don't know how to find whether error is following normal distribution or not.

Re: Non normal distribution in regression

Ksharp — Tue, 26 Jan 2016 02:56:42 GMT

Just a thought.

You can use proc univariate to check X and Y to see whether they are all normal distribution or not . If they were all conform to normal then you can say residual term is normal distribution.

Re: Non normal distribution in regression

plf515 — Tue, 26 Jan 2016 20:55:15 GMT

If the dependent variable is continuous but the assumptions of OLS regression are not met regarding normality of residuals, then I suggest PROC ROBUSTREG and PROC QUANTREG both of which relax those assumptions. I've written papers on these for last years SGF.

Re: Non normal distribution in regression

Babloo — Wed, 27 Jan 2016 06:21:50 GMT

Are you saying that if my data follows a normal distribution then error in the data will also follow a normal distribution? If not, may I request you to write a simple SAS code to demonstarte normal distribution for errors?

Re: Non normal distribution in regression

Ksharp — Wed, 27 Jan 2016 07:34:19 GMT

Yes. According to Statistical Theory , any linear combination of normal variables is also normal distribution. Therefore,

Y-X= epsilon , if Y and X all conform to normal distribution then epsilon also conform normal.

Otherwise, you could use other Robust Regression Method as other suggest .

This is just my two cents.

Re: Non normal distribution in regression

plf515 — Wed, 27 Jan 2016 12:23:58 GMT

But the converse isn't true. That is, you can have Y be non-normal and still have normal residual

Re: Non normal distribution in regression

lvm — Wed, 27 Jan 2016 15:35:06 GMT

Look at it this way. If Y is dependent (conditional) on X, then it is irrelevant to test whether Y is normally distriubuted (independent of X). That is, using proc univariate to assess normality of Y is meaningless. You want to check the normality of the residuals, or better, the normality of the studentized residuals. This is automatically done in graphic form by several procedures.

Re: Non normal distribution in regression

Reeza — Wed, 27 Jan 2016 15:39:32 GMT

But OLS doesn't have an assumption that the X and Y are normally distributed, only the errors. More an assumption that they're random rather than systematic.

Re: Non normal distribution in regression

Reeza — Wed, 27 Jan 2016 15:43:23 GMT

Look at the Fit Diagnostics panel from Proc Reg. I think it's produced by default these days.

These charts help assess normality of the Residuals. You could also extract the residuals from proc reg and pass them to proc NPAR1WAY which has a bunch of tests for normality.

Re: Non normal distribution in regression

lvm — Wed, 27 Jan 2016 15:47:28 GMT

I didn't say anything about a distriubtion for X. OLS does not rely on normality of Y|X for parameter estimation, but for inference (SEs, etc.), it is assumed that the error term is normally distributed. This is the same thing as Y|X being normally distributed (Y given X, not marginal Y). One assess the normality of Y|X using the residuals.

Re: Non normal distribution in regression

Reeza — Wed, 27 Jan 2016 15:52:06 GMT

@lvm I was referring to @Ksharp response in my post about normality, I agree with your response.