数学建模——线性回归模型-LMLPHP

目录

1.线性回归模型的具体步骤和要点: 

1.收集数据:

2.探索性数据分析:

3.选择模型:

4.拟合模型:

5.评估模型:

1.R平方(R-squared):

2.调整R平方(Adjusted R-squared):

3.残差分析:

4.方差膨胀因子(VIF):

6.解释结果:

7.预测与应用:

8.检验假设:

2.线性回归模型公式分析包括以下几个方面:

​编辑

3.模型代码实现

1.代码_python

2.图形


1.线性回归模型的具体步骤和要点: 

1.收集数据

首先,需要收集与研究问题相关的数据。这些数据应包括一个或多个自变量(特征)和一个因变量(目标)。

2.探索性数据分析

在建立模型之前,通常会对数据进行探索性分析,包括可视化和描述性统计分析,以了解数据的分布、相关性和异常值等情况。

3.选择模型

根据问题的特点选择合适的线性回归模型。如果只有一个自变量,可以使用简单线性回归模型;如果有多个自变量,可以使用多元线性回归模型。

4.拟合模型

利用最小二乘法或其他拟合方法来估计模型的参数。最小二乘法是一种常用的方法,它通过最小化观测值与模型预测值之间的残差平方和来确定参数。

5.评估模型

评估模型的好坏以及对数据的拟合程度。常用的评估指标包括R平方、调整R平方、均方误差等。

1.R平方(R-squared)

R平方是一个衡量模型拟合优度的指标,表示因变量的变异中能被自变量解释的比例。R平方越接近1,说明模型对数据的拟合越好。

2.调整R平方(Adjusted R-squared)

调整R平方考虑了自变量的数量和样本量,相比于R平方更可靠。

3.残差分析

分析残差是否呈现出随机分布,检查是否满足模型假设。

4.方差膨胀因子(VIF)

用于检测自变量之间的多重共线性问题。

6.解释结果

分析模型的参数估计,理解自变量与因变量之间的关系。通过检查参数的符号和大小,可以了解自变量对因变量的影响方向和程度。

7.预测与应用

利用拟合好的模型进行预测或者应用。可以使用模型对新的数据进行预测,也可以利用模型进行决策支持或政策制定等。

8.检验假设

在应用模型时,需要检验模型的假设是否成立,例如线性关系、常数方差、独立误差等。如果假设不成立,可能需要对模型进行修正或者选择其他的建模方法。

2.线性回归模型公式分析包括以下几个方面:

数学建模——线性回归模型-LMLPHP

3.模型代码实现

具体的需要根据具体数据磨合 

1.代码_python

import numpy as np
import statsmodels.api as sm
import matplotlib.pyplot as plt

# 准备数据
np.random.seed(0)
X = np.random.rand(100, 2)  # 两个自变量
y = 2 * X[:,0] + 3 * X[:,1] + np.random.randn(100)  # 因变量

# 添加常数项
X = sm.add_constant(X)

# 拟合线性回归模型
model = sm.OLS(y, X).fit()

# 绘制残差图
plt.figure(figsize=(12, 6))
plt.subplot(1, 2, 1)
plt.scatter(model.fittedvalues, model.resid)
plt.xlabel('Fitted values')
plt.ylabel('Residuals')
plt.title('Residuals vs Fitted')

# 绘制预测值与观测值的散点图
plt.subplot(1, 2, 2)
plt.scatter(model.fittedvalues, y)
plt.xlabel('Fitted values')
plt.ylabel('Observed values')
plt.title('Observed vs Fitted')

# 添加拟合直线
plt.plot(model.fittedvalues, model.fittedvalues, color='red')

plt.tight_layout()
plt.show()

# 绘制参数估计的置信区间
plt.figure(figsize=(8, 6))
model_params = model.params
conf_int = model.conf_int()
plt.errorbar(model_params.index, model_params, yerr=model_params - conf_int[:, 0], fmt='o')
plt.axhline(0, color='red', linestyle='--')
plt.xlabel('Parameters')
plt.ylabel('Estimate')
plt.title('Parameter Estimates with Confidence Intervals')
plt.xticks(rotation=45)
plt.show()

2.图形

数学建模——线性回归模型-LMLPHP

05-15 07:18