BookmarkSubscribeRSS Feed
🔒 This topic is solved and locked. Need further help from the community? Please sign in and ask a new question.
CatM
Obsidian | Level 7

Hello, 

 

I am trying to do a multilinear regression of income data as the outcome with disability as the main predictor and controlling for confounding of age education and gender. the income data is from two different years where i am trying to see if there have been any improvements or socioeconomic gains for people with disabilities. In order to look at this from a more relative angle, it was suggested that i convert the income data into centiles (100 by 1)  and use linear regression. However, i am running into issues because the income data includes negative and 0 income values which we do not wish to take out but analyze all together. These are being grouped into one centile which makes the data skewed when running the regression analysis. Is there any ideas on how this data should be analyzed instead. there is a weight for the data however this is not helping meet the assumption of normality. 

 

proc rank data=incdis groups=100 out=cent_incdis ;
var atinc42;
ranks atinc42_centile;
by year;
run;

proc surveyreg data=cent_incdis;
class ecsex99_n (ref='0') education (ref='0');
model atinc42_centile = disabs26_n ecsex99_n education ecage26 / solution;
weight icswt26;
run;

 

2000.PNG2011.PNG

1 ACCEPTED SOLUTION

Accepted Solutions
Reeza
Super User
Doh! You have a lot of tied data...look at the options for handling tied data, it’s under the PROC RANK statement documentation I believe.

View solution in original post

4 REPLIES 4
Reeza
Super User

The assumption is that residuals are normally distributed, not necessarily the variables. 

 

I think centiles is great idea, not sure how it can be skewed when by definition, 1% is in each group. Do you mean the residuals become skewed?

Reeza
Super User
Doh! You have a lot of tied data...look at the options for handling tied data, it’s under the PROC RANK statement documentation I believe.
CatM
Obsidian | Level 7

Ok great thank you! i just checked the residuals and they look fine. i'll look at tied data documents

Reeza
Super User
I’m not sure that I’d be that concerned anyways, the idea behind centiles is to see if the groups improvement is occurring at same rate as general population, so if that group has A larger number that’s reflective of reality. When this happens though youre more likely to see something like Simpsons Paradox occur due to an imbalance so be careful if you do subgroup analysis.

Ready to join fellow brilliant minds for the SAS Hackathon?

Build your skills. Make connections. Enjoy creative freedom. Maybe change the world. Registration is now open through August 30th. Visit the SAS Hackathon homepage.

Register today!
What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 4 replies
  • 820 views
  • 2 likes
  • 2 in conversation