Build regression model from a set of candidate predictor variables by entering predictors based on p values, in a stepwise manner until there is no variable left to enter any more.

ols_step_forward_p(model, ...)

# S3 method for default
ols_step_forward_p(model, penter = 0.3, progress = FALSE, details = FALSE, ...)

# S3 method for ols_step_forward_p
plot(x, model = NA, print_plot = TRUE, ...)

Arguments

model

An object of class lm; the model should include all candidate predictor variables.

...

Other arguments.

penter

p value; variables with p value less than penter will enter into the model

progress

Logical; if TRUE, will display variable selection progress.

details

Logical; if TRUE, will print the regression result at each step.

x

An object of class ols_step_forward_p.

print_plot

logical; if TRUE, prints the plot else returns a plot object.

Value

ols_step_forward_p returns an object of class "ols_step_forward_p". An object of class "ols_step_forward_p" is a list containing the following components:

model

final model; an object of class lm

steps

number of steps

predictors

variables added to the model

rsquare

coefficient of determination

aic

akaike information criteria

sbc

bayesian information criteria

sbic

sawa's bayesian information criteria

adjr

adjusted r-square

rmse

root mean square error

mallows_cp

mallow's Cp

indvar

predictors

Deprecated Function

ols_step_forward() has been deprecated. Instead use ols_step_forward_p().

References

Chatterjee, Samprit and Hadi, Ali. Regression Analysis by Example. 5th ed. N.p.: John Wiley & Sons, 2012. Print.

Kutner, MH, Nachtscheim CJ, Neter J and Li W., 2004, Applied Linear Statistical Models (5th edition). Chicago, IL., McGraw Hill/Irwin.

See also

Examples

# stepwise forward regression model <- lm(y ~ ., data = surgical) ols_step_forward_p(model)
#> #> Selection Summary #> ------------------------------------------------------------------------------ #> Variable Adj. #> Step Entered R-Square R-Square C(p) AIC RMSE #> ------------------------------------------------------------------------------ #> 1 liver_test 0.4545 0.4440 62.5119 771.8753 296.2992 #> 2 alc_heavy 0.5667 0.5498 41.3681 761.4394 266.6484 #> 3 enzyme_test 0.6590 0.6385 24.3379 750.5089 238.9145 #> 4 pindex 0.7501 0.7297 7.5373 735.7146 206.5835 #> 5 bcs 0.7809 0.7581 3.1925 730.6204 195.4544 #> ------------------------------------------------------------------------------
# stepwise forward regression plot model <- lm(y ~ ., data = surgical) k <- ols_step_forward_p(model) plot(k)
# final model k$model
#> #> Call: #> lm(formula = paste(response, "~", paste(preds, collapse = " + ")), #> data = l) #> #> Coefficients: #> (Intercept) liver_test alc_heavy enzyme_test pindex bcs #> -1178.330 58.064 317.848 9.748 8.924 59.864 #>