连老师，您好：

我有一个cross-sectional firm level data, 做了OLS regression with three fixed effects (country, industry, year)，即加入这三种变量的dummies。

模型是：Y = a*PD^2 + b*PD + c*X+e
Y是firm performance measured by ROA (return on assets),
PD是key variable (product diversification),
X是control variables和three fixed effects (country, industry, year) 的dummies (country_*, year_* and industry_*)

regression的结果是beta系数a负显著，b正显著，control variables的beta系数也都显著。a负b正说明PD和Y有倒U型的关系。这也符合当今主流的product diversification-performance relationship literature的观点。即适中的产品多样化的对公司绩效有最大正作用，很小或者过度的产品多样化对公司绩效有负作用。

现在我想把PD和Y的倒U型的关系画在直角坐标系上，x轴是PD，y轴是Y。这曲线图最好是已经控制了control variables对Y的影响的情况下，PD和Y的曲线关系。

我尝试用以下stata code, 但画出来的两种曲线图都不理想（产生的图里面同时带有散点和曲线，而我只需要曲线）。也不知道下面的code正不正确。
reg Y PD PD^2 X country_* year_* industry_*, r
predict Y_hat
twoway (scatter Y_hat PD) (qfit Y_hat PD)
graph save graph_1, replace
twoway (scatter Y PD) (qfit Y PD)
graph save graph_2, replace


请问如何用stata画出这种曲线图(curvilinear relationship)？

请您抽空指导一下，谢谢。

学生：Rooney