title: Extracting coefficients and their standard error for each unit in an lme model fit
body: <p>How could I extract coefficients (b0 and b1) with their respectively standard errors for each experimental unit (plot )in a linear mixed model such as this one: </p>
  
  <p><a href="https://stackoverflow.com/questions/23812039/better-fits-for-a-linear-model">Better fits for a linear model</a></p>
  
  <p>with this same dataset(df), and for the fitted model (fitL1): how could I get a data frame as this one...</p>
  
  <pre><code>   plot    b0      b0_se   b1    b1_se 
      1    2898.69   53.85   -7.5  4.3
  
     ...    ...       ...     ...   ...
  </code></pre>