Dear Statalisters,

Currently I am working on a panel dataset where my probit model predicts either 1 or 0 if the probability is higher than a certain percentage. These values correspond to dates which I would like to graph as a shaded area underneath the timeseries data per panel. Unfortunately I can't find an option to do this I know this options exists in the normal twoway graph using the 'rarea' option but this doesn't seem to work in a panel data setting. Any help is much appreciated!

Code:
xtprobit Boom lRent lGDP UNE Long_interest Short_interest lM3 CPI
predict Boomhat, pr
gen Boomhat1 = 1 if Boomhat>0.35
xtline NHP Boomhat1 if country ~=3, byopts(title(Nominal House Price Indicator)l1title("") b1title("") note("Source: OECD Nominal House Price Indicator")  rescale)  i(Country) t(date) tlabel(#3) ylabel(#2) ylabel(, angle(horizontal))