Selección de modelo para modelos con efectos mixtos basados en AIC.

Me gustaría ajustar un modelo de efectos mixtos y obtener una lista de los diez principales modelos candidatos clasificados por AIC. Para modelos de efectos fijos, utilizo la función glmulti:glmulti para ajustar el modelo y la función glmulti:weightable para obtener una lista de modelos candidatos clasificados por AIC, lo que da como resultado esto:

modelo     aicc      pesos

1          resultado ~ 1 + ist + hwt 188.3154 0.6044616289

2 resultado ~ 1 + ist + hwt + troadn 190.3552 0.2179837514

3                resultado ~ 1 + hwt 191.7862 0.1065837756

4       resultado ~ 1 + hwt + troadn 193.5704 0.0436764670

5                resultado ~ 1 + ist 195.8793 0.0137685009

6       resultado ~ 1 + ist + troadn 196.3339 0.0109688767

7             resultado ~ 1 + troadn 200.0687 0.0016949662

8                      resultado ~ 1 201.4210 0.0008620334

¿Existe un paquete equivalente para completar estos mismos pasos con un modelo de efectos mixtos? Usé la función lme4::glmer para ajustar un modelo de efectos mixtos, pero esto no produce una lista de los principales modelos en función del puntaje AIC. ¿Existe una función o paquete alternativo, o no se utiliza este proceso con modelos de efectos mixtos?

Me gustaría obtener una lista de los 10 mejores modelos candidatos para:

(finalidad ~ sexo + temporada + año + (1 | nombre_obs) + (1 | pájaro), datos = detectar, familia = binomial)