find. hyper ()function gamlss

I am trying to create a model to estimate the centiles of the angle variable by the age variable (age). However, I can't run the following model because I don't know how to generate p [1] and I don't know what p [1 or 2 or 3] means in the example Can anybody help me?

Tage<-(Homens$idade)^(m0$power)
mOD <- gamlss(angulo~cs(Tage,df=op$par[1]), 
             sigma.formula=~cs(Tage,df=op$par[2]),
             nu.formula=~cs(Tage,df=op$par[3]),c.spar=(-1.5,2.5), family=BCCGo,data=Homens)