主要是上置信区间与下置信区间的线，您要自行创立，
我提供一个较简单做法，您参考一下。
【就非参文献而言，这种置信区间是很烂的，有学问的要执行wild bootstrap下的confidence interval】

webuse motorcycle,clear
lpoly accel time, degree(3) kernel(gau) generate(x s) se(se) level(90)
g ub=s+1.645*se
g lb=s-1.645*se
kernreg1 accel time, b(2.4) k(4) np(100) gen(m2p4 g2p4)
sort x g2p4
twoway (line m2p4  g2p4 ) (line s x)   (line ub x)  (line lb x)

each figure presents the estimated Engel curve for four models: thick black lines indicate 3SLS estimates of the symmetry-restricted exact model; thick gray lines indicate 3SLS estimates of the symmetry-restricted quadratic (having a second- instead of fifth-order polynomial iny) exact model; thin black lines indicate approximate SUR estimates of the model; and medium black lines indicate exact 3SLS estimates of the asymmetric model (italics denote shorthand). Estimates are computed at each percentile of log expenditures in the data, and estimated 90 percent confidence intervals (computed via the delta method) are displayed with small crosses for the exact 3SLS symmetry-restricted model at each fifth percentile of the expenditure distribution.

这个图有什么stata实现的方法，还是必须用R软件？

唉，惭愧，都忘得差不多了，帮不上忙

您附上的那个图，一看就知道是用R绘制的。
我对3SLS不是很熟，按文中所述，3SLS的一些模型变体，我个人不是很了解，
在Stata方面，您应当参考 reg3指令，查查看吧！【reg3里面可以估SUR】

至于R 我想，您google一下，3SLS in R，应该有不少东西。

祝 顺利 ＆ 开心

g ub=s+1.645*se
g lb=s-1.645*se
加减1.645的依据是什么？