The estimation of functional β(t)^ by bΣ is done in step (2), and separated from the estimation of dependence structure ρ in step (3). This allows for the flexibility of including any type of dependence structures designed by the user (for instance, using particular restrictions) that are typically not included in the usual packages (like nlme).