Controlling False Discovery Rate (FDR) in Regression with Correlated IVs

I am looking for a procedure/code to control FDR in a multiple regression analysis, including 1 DV and 28 IVs (some of them interaction terms). Some IVs are correlated. Where can I find the appropriate FDR procedure/code for this problem? Thanks. Salar