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



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;

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;




Accepted Solutions
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

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?

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.
Obsidian | Level 7

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

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.



Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.

If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website. 

Register now!

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.

Get the $99 certification deal.jpg



Back in the Classroom!

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

View all other training opportunities.

Discussion stats
  • 4 replies
  • 2 in conversation