How to calculate effect size f-squared in PLS-SEM using R programming

Dear R experts,

I am conducting a research using PLS-SEM approach using plspm package. However, I have no idea how to calculate effect size f-squared using R script in this case. Can anyone share some knowledge or any technical paper regarding this issue? I knew smart-pls software can calculate the f2 values, however, I do not know how to run the smart-pls. Preferable R script if possible. If unavailable, any idea how to calculate manually one by one for each path relationships? Thanks for sharing in advance. Good day.