gen Flnincome1=F.lnincome1
xtset Stkcd y
xi:qregpd Flnimcome1 lnPatent $control, i.industry i.y, q(0.1) identifier(Stkcd) fix(y)
报错：invalid 'q'
想请问大家，stata中做面板数据的分位数回归，如何把行业和年度固定效应加进去？如何在行业层面聚类？