相关软件（R代码）

# 使用MASS包中的Pima数据
# install.packages("MASS")
library(MASS)
data(Pima.tr)
# 建立线性模型
model <- lm(glu ~ bmi + age, data = Pima.tr)
summary(model)
# 模型诊断图
par(mfrow = c(2, 2))
plot(model)
# 提取系数和置信区间
confint(model, level = 0.95)
# 预测新数据
new_data <- data.frame(bmi = c(25, 30, 35), age = c(30, 45, 60))
predict(model, newdata = new_data, interval = "confidence")
# 检查多重共线性
# install.packages("car")
library(car)
vif(model)

2.3.2 极大似然估计法