d <- read.csv(‘https://stats.dip.jp/01_ds/data/winequality-red.csv’) datatable(d)

fit.full <-lm(fixed.acidity~volatile.acidity+citric.acid+residual.sugar+chlorides+free.sulfur.dioxide+total.sulfur.dioxide+density+pH+sulphates+alcohol,data = d)

library(sjPlot) tab_model(fit.full, show.stat = T, show.aic = T)

fit.manual <-lm(fixed.acidity~volatile.acidity+citric.acid+residual.sugar+chlorides+free.sulfur.dioxide+total.sulfur.dioxide+density+pH+sulphates+alcohol,data = d)

tab_model(fit.manual, show.stat = T, show.aic = T)

library(MASS) fit.aic <- stepAIC(lm(fixed.acidity~(volatile.acidity+citric.acid+residual.sugar+chlorides+free.sulfur.dioxide+total.sulfur.dioxide+density+pH+sulphates+alcohol)^2,data = d), trace =0)

tab_model(fit.aic, show.stat = T, show.aic = T)

tab_model(fit.full, fit.manual, fit.aic, show.stat = T, show.aic = T, dv.labels = c(‘quality除きフルモデル’, ‘手動選択モデル’, ‘ステップワイズ変数選択モデル’))