how to extract random effects from a gamlssZadj object in R?

Question

I fitted the following model using the function gamlssZadj from the R package gamlss.inf:

fit <- gamlssZadj(y = score,
    mu.formula =~ bs(age) + gender + ran