视频网站如何做微信营销,网站型与商城型有什么区别吗,video.js wordpress,办公室装修设计图片一、实验目的#xff1a;
#xff08;1#xff09;理解一元线性回归和多元线性回归的数学原理#xff0c;能够利用sklearn中相关库解决现实世界中的各类回归问题#xff1b;
#xff08;2#xff09;掌握利用matplotlib对一元线性回归模型进行可视化的方法#xff0c…一、实验目的
1理解一元线性回归和多元线性回归的数学原理能够利用sklearn中相关库解决现实世界中的各类回归问题
2掌握利用matplotlib对一元线性回归模型进行可视化的方法并分析模型的优劣
3掌握利用statsmodels进行线性回归并进行模型评价的方法
4掌握利用回归模型对现实问题进行分析和解释的能力。
二、实验仪器及实验环境
1硬件PC机
2软件Anaconda Jupyter NotebookSpyderPandas
三、实验内容 一基于伊春市木材剩余物数据利用回归模型预测木材剩余物
伊春林区位于黑龙江省东北部。全区有森林面积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 pd
filepathrD:\木材剩余物.csv
dfpd.read_csv(filepath,sep,,index_col0,enginepython)#enginepython
columns[剩余物,采伐量]
df.columnscolumns
df.head() #因为\在python中时转义的意思这里所表示的路径在编码时无法被正确识别加上r可强制不转义。 #因为\在python中时转义的意思这里所表示的路径在编码时无法被正确识别加上r可强制不转义。 index_col0告诉Pandas使用第一列作为索引。 2、利用分割X和Y数据集。
Xdf[[采伐量]]#这里如果不加两层中括号没有表格
Ydf[剩余物]
X.head() 3、对数据进行可视化显示。
import matplotlib.pyplot as plt#调用模块进行可视化显示
plt.scatter(X,Y) 4、利用sklearn中的线性回归模型建立回归模型对模型进行训练输出模型参数。
from sklearn.linear_model import LinearRegression
regrLinearRegression()
regr.fit(X,Y)
regr.coef_,regr.intercept_#regr.coef代表yaxb中的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
X2sm.add_constant(X)
estsm.OLS(Y,X2).fit()
print(est.summary()) 代码解读 sm.add_constant()函数用于增加截距项也就是增加一个常数项。 sm.OLS()函数使用OLS普通最小二乘法建立线性回归模型est。 est.summary()用于输出模型评估结果。
