一、实验目的：
（1）理解一元线性回归和多元线性回归的数学原理，能够利用sklearn中相关库解决现实世界中的各类回归问题；
（2）掌握利用matplotlib对一元线性回归模型进行可视化的方法，并分析模型的优劣；
（3）掌握利用statsmodels进行线性回归并进行模型评价的方法；
（4）掌握利用回归模型对现实问题进行分析和解释的能力。
二、实验仪器及实验环境
（1）硬件：PC机；
（2）软件：Anaconda Jupyter Notebook，Spyder，Pandas
三、实验内容：
（一）基于伊春市木材剩余物数据利用回归模型预测木材剩余物
伊春林区位于黑龙江省东北部。全区有森林面积218.9732万公顷，木材蓄积量为2.324602亿m3。森林覆盖率为62.5%，是我国主要的木材工业基地之一。1999年伊春林区木材采伐量为532万m3。按此速度44年之后，1999年的蓄积量将被采伐一空。所以目前亟待调整木材采伐规划与方式，保护森林生态环境。为缓解森林资源危机，并解决部分职工就业问题，除了做好木材的深加工外，还要充分利用木材剩余物生产林业产品，如纸浆、纸袋、纸板等。因此预测林区的年木材剩余物是安排木材剩余物加工生产的一个关键环节。下面，利用一元线性回归模型预测林区每年的木材剩余物。显然引起木材剩余物变化的关键因素是年木材采伐量。
给出伊春林区16个林业局1999年木材剩余物和年木材采伐量数据见“木材剩余物.csv”。
1.读取伊春市木材剩余物数据集，并显示输入
OSError: Initializing from file failed
解决方法：http://t.csdnimg.cn/afyEh
import numpy as np
import pandas as pdfilepath=r"D:\木材剩余物.csv"
df=pd.read_csv(filepath,sep=",",index_col=0,engine='python')#engine='python'
columns=["剩余物","采伐量"]
df.columns=columns
df.head()
#因为\在python中时转义的意思，这里所表示的路径在编码时无法被正确识别，加上r可强制不转义。
#因为\在python中时转义的意思，这里所表示的路径在编码时无法被正确识别，加上r可强制不转义。
index_col=0告诉Pandas使用第一列作为索引。
2、利用分割X和Y数据集。
X=df[["采伐量"]]#这里如果不加两层中括号没有表格
Y=df["剩余物"]
X.head()
3、对数据进行可视化显示。
import matplotlib.pyplot as plt#调用模块进行可视化显示
plt.scatter(X,Y)
4、利用sklearn中的线性回归模型建立回归模型，对模型进行训练，输出模型参数。
from sklearn.linear_model import LinearRegression
regr=LinearRegression()
regr.fit(X,Y)
regr.coef_,regr.intercept_#regr.coef代表y=ax+b中的a,权值，而regr.intercept代表截距，就是b
5、假设乌伊岭林业局2000年计划采伐木材20万m3，求木材剩余物的点预测值。
通过计算，置信度为0.95的2000年平均木材剩余物E(y2000)的置信区间是
从而得出预测结果，2000年若采伐木材20万m3，产生木材剩余物的点估计值是7.3231万m3。平均木材剩余物产出量的置信区间估计是在 [5.8736, 8.7726] 万m3之间。从而为恰当安排2000年木材剩余物的加工生产提供依据。
6、利用statsmodels实现线性回归并对模型进行评估。
import statsmodels.api as sm
X2=sm.add_constant(X)
est=sm.OLS(Y,X2).fit()
print(est.summary())
代码解读：
sm.add_constant()函数用于增加截距项，也就是增加一个常数项。
sm.OLS()函数使用OLS（普通最小二乘法）建立线性回归模型est。
est.summary()用于输出模型评估结果。
