¡Hola!

#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 parametros de otro año y tipo de manejo
#Gracias