Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Forecasting
- /
- How should I do the regression with partially miss...

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

09-05-2013 05:02 PM

Suppose I have three variables: Y, X1 and X2, while both Y and X1 have 100 observations, but X2 only has, say, 30 observations.

I want to estimate an equatino as Y=X1*b1+X2*b2, while utilizing all the information I have, i.e., I do not want to discard the 70 observations with missing X2s. How am I supposed to write the code?

Can I write it in this way:

prco model data=yx1x2;

parameters b1 b2;

if x2=. then

eq1=y-x1*b1;

else

eq1=y-x1*b1-x2*b2;

fit eq1;

run;

Behind the scene, how does SAS process this equation, I mean, what is the algorithm ?

Thank you very much.

Accepted Solutions

Solution

09-06-2013
03:59 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

09-06-2013 03:59 PM

All Replies

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

09-05-2013 05:06 PM

1. Can you impute your missing data? SAS has procedures for that.

2. Is your data missing at random or systematic and continuous or categorical? If categorical, can you include "Missing" as a category?

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

09-06-2013 11:02 AM

My constraint is that I cannot impute the missing values. Let's suppose the missing is random and the variable is numeric. I actually can get the code run, if I have the specification I mentioned in the inital post, but I am not sure whether the results are reliable.

Solution

09-06-2013
03:59 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

09-06-2013 03:59 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

09-19-2013 10:52 AM

Thanks, kessler. You are right. That is how SAS does it behind the scene.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

09-06-2013 06:20 PM

If you want to capture the fact that the average of missing x2 values might not be zero, you could try fitting your model this way :

**proc model data=yx1x2;**

**parameters b0 b1 b2 bz;**

**z = missing(x2);**

**if z then x2=0;**

**y = b0 + x1*b1 + x2*b2 + z*bz;**

**fit y;**

**run;**

Parameter **b0** will account for the overall intercept (you may remove it later if it is not significant) and **bz** will account for the average effect of missing x2 values.

PG

PG

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

09-19-2013 10:55 AM

Thanks, PG. That is a good way to work around the issue.