De un modelo lineal mixto (lme), me interesa obtener los intervalos de confianza (IC) de todos los niveles de los factores fijos. Con la función `intervals()` consigo el intervalo del intercepto, como referente, y sólo de algunos niveles, pero no todos ¿Saben cómo puedo obtener los IC de todos los niveles?

El modelo es el siguiente:

    Ca <- lme(fixed = Ca ~ ano+manejo, data = mc1, random = ~ 1 | parcela/subparcela, control=list(maxIter=1000, niterEM=1000))

    intervals(Ca, which = "fixed")
    Approximate 95% confidence intervals
    
     Fixed effects:
                    lower      est.     upper
    (Intercept)  6113.837  6960.000 7806.1632
    anosiete    -2067.245 -1503.333 -939.4221
    manejoO     -2112.795 -1154.167 -195.5381

faltarían los parámetros de otro año y tipo de manejo