R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据

全文链接:http://tecdat.cn/?p=30914

原文出处:拓端数据部落公众号

我们正和一位朋友讨论如何在R软件中用GLM模型处理全国的气候数据。本文获取了全国的2021年全国的气候数据。

采样时间:2021年1月1号~2021年12月31号
采样地点:全国各地。

本次调查搜集了2021年全国不同地区的风向、降雨量、风速、风速变化、最大风速、最大降雨量、闪电概率等数据。并对不同变量之间的相关性进行了调查,对国家数据预测的错误率进行了GLM模型拟合。

读取数据

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第1张图片

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第2张图片

library(car)
library(MuMIn)
head(data)

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第3张图片

读取因变量

numberFaults=data$numblts

head(data1)     

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第4张图片

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第5张图片

# 相关分析

调查的出的各指标数据用R软件进行处理并且用箱图进行对比显示。

# 部分指标的箱线图

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第6张图片

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第7张图片

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第8张图片

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第9张图片

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第10张图片

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第11张图片

查看各变量之间的相关系数

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第12张图片

有显著的相关关系。从变量相关关系图和矩阵,可以看到temperatureMin和temperatureMax,windChillMin,windChillMax,以及gustSpeed和windSpeed之间,rainfall和changeInRainfall,以及lightningRisk和lightningCategory之间都有教明显的线性相关关系。yearDay和windChill之间也有一定的相关关系。

glm 线性回归模型

summary(glm.po)

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第13张图片

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第14张图片

检验是否存在多重共线性问题

kappa(cor(data[,c(1:15,17:20)]), exact=T)
## [1] 3.020456e+18

判断多重共线性变量

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第15张图片

进一步模型优化

step(glm.po2)

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第16张图片

summary(glm.step)

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第17张图片

vif

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第18张图片

从模型中变量的VIF值来看,大多数变量之间不存在较强的多重共线性关系。 从残差拟合图来看,大部分样本拟合值分布在0周围,说明拟合结果较理想。981,2331和524号样本可能为异常点。从正态分布qq图来看,大部分点分布在图中直线附近。说明样本点服从正态分布。同样,拟合值的标准残差也分布在红线周围,说明拟合效果较好。同样,大部分样本的cook’ distance距离在正常范围内,392,624,622号样本的cook’ distance较大,可能会对模型产生较大的影响。

全子集回归来选出最优的模型

全子集回归,即基于全模型获得可能的模型子集,并根据AIC值等对子集排序以从中获取最优子集。

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第19张图片

重新拟合模型

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第20张图片

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第21张图片

优化模型

avg(ms1, subset = delta < 10,fit=T,rank = "AIC")

残差图

plot(pre-numberFaults)

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第22张图片

计算R-squre值,查看模型拟合情况

Rsquare=ssr/sst

R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第23张图片

从逻辑回归结果来看,模型中部分自变量对因变量的影响较为明显,达到了0.01的显著性水平,具有一定的理论意义。然而从AIC的值来看,达到了4024.881,数值较大,说明模型的拟合度较差,有进一步改进的空间。 一般认为计算条件数kappa(X),k<100,说明共线性程度小,如果1001000,存在严重的多重共线性。 从结果来看,kappa值远远大于1000,因此判断该模型存在严重的共线性问题,即线性回归模型中的解释变量之间由于存在精确相关关系或高度相关关系而使模型估计失真或难以估计准确。changintemp,changinrainfall,windspeedmin,lightningcategory ,rainfallMax之间存在较高的相关关系会对模型的拟合失真并且难以用以估计和预测,因此,删掉这些变量后重新对模型进行拟合。 删除部分共线性程度高的变量后可以看到模型的AIC降低了,因此,模型的拟合程度提高了。 大部分样本的cook’ distance距离在正常范围内,392,624,622号样本的cook’ distance较大,可能会对模型产生较大的影响。


R语言广义线性模型(GLM)、全子集回归模型选择、检验分析全国风向气候数据_第24张图片

最受欢迎的见解

1.[](http://tecdat.cn/r%e8%af%ad%e...)R语言多元Logistic逻辑回归 应用案例

2.[](http://tecdat.cn/r%e8%af%ad%e...)面板平滑转移回归(PSTR)分析案例实现

3.[](http://tecdat.cn/r%e8%af%ad%e...)matlab中的偏最小二乘回归(PLSR)和主成分回归(PCR)

4.[](http://tecdat.cn/r%e8%af%ad%e...)R语言泊松Poisson回归模型分析案例

5.[](http://tecdat.cn/r%e8%af%ad%e...)R语言回归中的Hosmer-Lemeshow拟合优度检验

6.[](http://tecdat.cn/r%e8%af%ad%e...)r语言中对LASSO回归,Ridge岭回归和Elastic Net模型实现

7.[](http://tecdat.cn/r-%e8%af%ad%...)在R语言中实现Logistic逻辑回归

8.[](http://tecdat.cn/r%e8%af%ad%e...)python用线性回归预测股票价格

9.[](http://tecdat.cn/r%e8%af%ad%e...)R语言如何在生存分析与Cox回归中计算IDI,NRI指标

你可能感兴趣的:(数据挖掘深度学习人工智能算法)