现在做断点回归需要控制时间固定效应和行业固定效应， 不控制固定效应的stata代码为
rd Y z, deg(4) kernel (triangle) mbw (100) cluster(id).  Y是结果变量，z是处理变量，求问如果想控制上述两种固定效应，stata代码应该怎么写。
谢谢解答。