I would try PROC GLIMMIX. It appears that you have an ordinal multinomial response, one continuous variable, one classification variable (country), and one variable that is going to cause problems (year). I don't think year has the same meaning in each country, and in some sense, it is a repeated measurement on the country. Something like this (untested) might be a good start: proc glimmix data=yourdata; class country year; model happiness= income country income*country/dist=mult link=clogit; random year/subject=country residual type=un; run; Options for odds ratios could be added. Example 40.13 in the GLIMMIX documentation may help as well. Steve Denham
... View more