step-by-step creation of models with accumulation of variables

Can you please tell me if step linear models of the independent variable "ols_step_both_p()" (R) are possible with the accumulation of predictors in the amount of 58, 220 and 299, naturally with the Bonferoni corrections, which are also accumulated? can all intermediate results be analyzed in the article or not? or is it allowed to publish only the result of one stage (any)? The main problem is that with the gradual addition of variables to the model, the reliable components of the models change due to the negative interaction between the variables and the increasing Bonferroni correction factor to p-value. Please help me write code to find a linear relationship model by choosing from 297 predictors (and if possible + 2 more categorical predictors)