Multiple axis labels for scale_discrete

I think CRAN - Package ggh4x might help you, if you consider the yes, no a different facet. You could do something like facet_nested(~marital+year)and leave the labels just by year.