Dear Stata users,

I have generated restricted cubic splines (4 knots) for A, and then fit the median regression model. Now I would like to calculate 95% CI for median regression line, but do not know how to do this on Stata. I would really appreciate if anyone may help me with this.

Codes:
PHP Code:

mkspline A_ 
v nknots(4cubic

lqreg B A_1 A_2 A_3
quantiles(50nodots seed(123)

lqregpred adjs, for(A_1 A_2 A_3plotvs(A
Thank you so much.
Best wishes.