Como eu poderia extrair coeficientes (b0 e b1) com seus respectivos erros padrão para cada unidade experimental (plot) em um modelo linear misto como este:

Melhor ajuste para um modelo linear

com esse mesmo conjunto de dados (df) e para o modelo ajustado (fitL1): como obter um quadro de dados como este ...

   plot    b0      b0_se   b1    b1_se 
    1    2898.69   53.85   -7.5  4.3

   ...    ...       ...     ...   ...