How to obtain confidence intervals for predicted value in using natural spline in geeglm

Question

Below is what I have done and the error message I am getting each time

gee_fit <- geeglm(rslt~ ns(day,3), data = epidem,id =patid, family = binomial, corst