data("RedWine")
full_formula <- as.formula(quality ~ fixed.acidity + volatile.acidity + citric.acid
+ residual.sugar + chlorides + free.sulfur.dioxide +
total.sulfur.dioxide + density + pH + sulphates + alcohol)
fullmodel <- polr(formula = full_formula,data=RedWine, method = "probit")
select_model <- update(fullmodel, formula. = ". ~ . - fixed.acidity -
citric.acid - residual.sugar - density")
surr_rsq_select <- surr_rsq(select_model, fullmodel, data = RedWine, avg.num = 30)
# surr_rsq_ci(surr_rsq_select, alpha = 0.05, B = 1000) # Not run, it takes time.
Run the code above in your browser using DataLab