求教各位大佬，我做了一个回归，让x和年份进行交互，想计算出核心解释变量在每一年产生的影响，想把结果进行图示，用margins然后marginsplot，但margins一直算不出来，所以想请问各位大佬有什么好的方法。
代码如下：
reghdfe lnpergdp c.lnseg##i.year $control, a(city_code year) vce(cl city_code)
margins, dydx(lnseg) at(year=(2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019))
希望的效果见附件，横坐标就是时间