如下面所示数据，对应每一个code的公司观察值，虚拟变量D打分需要做如下修改：
自从第一次被观察到D=1后（最早是2012年，不确定哪一年被观察到D=1，以300294为例，2013年第一次被观察到D=1；600969则是2014年第一次被观察到D=1）），后续年度D都赋值为1，而不是目前这种偶尔带0的情况（600969在2014被观察到D=1后，2015一直到2020都应该D=1，但是目前有些年度D=0）。
请问如何编程实现上述目标？
谢谢大神！