How to donate letters to significant differences in R?

Hi Shahram,
Check out the cld() (compact letter display) function in the emmeans package. There's more information about that in this past question.

Good luck!