I'm examining whether the rates of increase in income within a certain profession are significantly different in different parts of the income distribution, to see if the income gap is significantly stretching or closing.
THE QUANTREG MODEL
I have performed a quantile regression in SPSS (I'm new to coding, only have vey basic knowledge of Python so I need your help). The dependent variable is indexed income, the independent variables are time (quarter in this dataset), demographic groups, segments of the profession. I have also added the interaction terms of each dummy with the time variable.
So (at least the way I see it), this model allows comparison of changes in income on three levels:
- How does belonging to a certain demographic group or segment impact income (e.g., compared with data entry jobs: data analytics adds 100€, data science adds 200€)
- How are the effects of each different category/dummy changing over time (e.g., compared with data entry jobs, the positive effect of being a data scientist has increased 10% and now adds 220€)
- How do these changing effects differ between different parts of the income distribution (e.g., the coefficient of time*data_scientist is much larger in the 90%Q than in the 10%Q, indicating that the higher-earning data scientists have seen a bigger increase in income over time than the lower-earning data scientists)
QUESTION
So I've got my output of this quantile regression on SPSS, a huge table with all the coefficients and their significance and confidence intervals.
Now I want to find out whether the differences between the 90%Q and the 10%Q are statistically significant, in order to make statements about whether the income gap in this profession has significantly increased or decreased. I thought to do this on Python instead of SPSS, I've searched how to cut the data into quantiles, and how to perform a quantile regression. But how should one get on with testing the significance of difference between 90%Q and 10%Q?