How can I get the percentile estimate in the gamlss model with categorical variable?

Question

there is a subject I want to ask you about centiles prediction.

Model1<-gamlss (y~cs(x,df=5)+il, sigma.formula=~cs(x,df=3)+il,nu.formula=~cs(x,df=1),tau.formula=~c