topic Re: PROC REG - RESTRICTIONS in Statistical Procedures

PROC REG - RESTRICTIONS

AlexeyS — Sun, 23 Jun 2019 04:42:37 GMT

I want to run linear regression, but I have to add the next restriction : yhat>=0, i.e predicted variable have to be greater or equal to zero. I use proc reg. if I cannot use it, maybe SAS has another procedure that can do it.

Thank you

Re: PROC REG - RESTRICTIONS

Rick_SAS — Sun, 23 Jun 2019 10:14:31 GMT

Please tell us more about your problem. In particular:

1. Is the response variable an inherently positive quantity such as a count, a volume, or a weight? If it is a count, you might consider using PROC GENMOD and Poisson regression.

2. Are all of the observed responses positive? If so, the predicted values will also be positive on the convex hull of the observations. For an example, see the article "Truncate response surfaces."

Re: PROC REG - RESTRICTIONS

plf515 — Sun, 23 Jun 2019 10:48:27 GMT

Why do you want to restrict the predicted values to positive numbers?

What is the nature of your dependent variable?

If the DV is a count, then there is PROC COUNTREG which lets you do Poisson regression or negative binomial regression (I have never had a case where the assumptions of Poisson are valid - they may exist, but it may be as unlikely as a cavalry officer getting kicked to death by his horse. COUNTREG offers other options too.

If the DV is bounded on both ends (not just the low end at 0) then beta regression could work (you may have to divide the DV by its max value). You can do this in PROC GLIMMIX using the DIST = beta option on the model statement.

If neither of those are the case, you could try taking the log of the DV and then back transforming after the regression - but I don't necessarily recommend this.

Re: PROC REG - RESTRICTIONS

AlexeyS — Sun, 23 Jun 2019 11:00:20 GMT

My predicted dependent variable should be greater or equal to zero. It's like volume. I use proc reg, and didn't find an option to restrict yhat be ge zero

Re: PROC REG - RESTRICTIONS

AlexeyS — Sun, 23 Jun 2019 11:02:16 GMT

My dependent variable is greater or equal to zero. I want that predicted DV will be the same. My variable looks like volume. I use proc reg but I saw that this procedure don't have an option to restrict predicted value of DV.

Re: PROC REG - RESTRICTIONS

PaigeMiller — Sun, 23 Jun 2019 11:15:42 GMT

@AlexeyS wrote:

My predicted dependent variable should be greater or equal to zero. It's like volume. I use proc reg, and didn't find an option to restrict yhat be ge zero

It seems to me that @Rick_SAS has already provided the solution. If it doesn't work in your case, please give specific and detailed information about why it doesn't work.

Re: PROC REG - RESTRICTIONS

Rick_SAS — Sun, 23 Jun 2019 12:11:46 GMT

There is no way to restrict the predicted value in PROC REG. The predictions are determined by least squares regression.

If you are getting negative predicted values, then you need to examine other options. For example, the regression fit might be influenced by one or more outliers or high-leverage points that are pulling the regression line up or down. If so, you might try PROC ROBUSTREG for robust regression.

I also like Peter's suggestion. However, instead of performing an OLS fit of LOG(Y) and then taking EXP of the predicted values, you could use PROC GENMOD to fit a log-link model. For a discussion of the two models and why they are different, see "Error distributions and exponential regression models."

Re: PROC REG - RESTRICTIONS

Ksharp — Sun, 23 Jun 2019 12:37:59 GMT

I am not sure, you could try other distribution and LINKED function.

proc genmod data=sashelp.class;
model weight=height/dist=gamma link=log;
run;