Unit-specific linear time trends Goodman-Bacon decomposition?

[Robert Garcia]

Hi

Is there any way to control for unit-specific linear time trends when performing a Goodman-Bacon decomposition in Stata?

My DiD model is

xtreg outcome treatment i.year i.unit##c.year, fe cluster(unit)

I tried to run

bacondecomp outcome treatment i.unit##c.year, cluster(unit)

but got the error message

factor-variable and time-series operators not allowed
r(101);

Is there any way to create the interaction term for unit-specific linear time trends besides

i.unit##c.year

?

I am using Stata 15.1

All input is appreciated!