首页 > 其他分享 >R语言广义加性混合模型(GAMM)分析长沙气象因子、空气污染、PM2.5浓度、显著性检验、逐日变化可视化|附代码数据

R语言广义加性混合模型(GAMM)分析长沙气象因子、空气污染、PM2.5浓度、显著性检验、逐日变化可视化|附代码数据

时间:2024-07-12 18:53:57浏览次数:10  
标签:显著性 模型 颗粒物 PM2.5 逐日 GAMM 浓度

全文链接:https://tecdat.cn/?p=32981

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

气候变化和空气污染对现代社会产生了越来越大的影响。在这种背景下,研究气象和空气污染之间的关系以及其对PM2.5浓度的影响变得非常重要。为了更好地理解和解释这些关系,广义加性混合模型(GAMM)成为一种强大的工具。

长沙作为湖南省的省会城市,其气象条件和空气质量一直备受关注。通过分析长沙地区的气象数据、空气污染指标和PM2.5浓度,可以更加全面地了解该地区的空气质量状况,并揭示气象因素对其变化的影响。

本研究旨在使用R语言中的广义加性混合模型(GAMM)方法,帮助客户来探索长沙地区气象因素与空气污染之间的关系,并进一步研究它们对PM2.5浓度变化的影响。我们将采集长沙地区的逐日气象数据、空气污染指标以及PM2.5浓度数据,应用GAMM模型进行分析。

在分析过程中,我们将考虑多种气象因素,如温度、湿度、风速等,并结合空气污染指标,如PM10、SO2、NO2等,来建立相应的GAMM模型。通过显著性检验,我们将评估各个因素对PM2.5浓度的影响程度,并进一步进行逐日变化的可视化分析,以呈现其动态变化规律。

查看数据

image.png

image.png

读取数据

   

airquality=read.csv("长沙气象站逐日气象+空气污染 数据.csv")  
  
Region=read.csv("长沙市年各个站点的PM2.5 每日浓度.csv",skip=1)

image.png

image.png

Table 2. Significance tests of PM2.5levels for different season, region, daytime and day of week. 对不同季节、地区、白天和周日的PM2.5水平进行显著性检验。

这段内容提到了对不同季节、地区、白天和周日的PM2.5水平进行显著性检验。显著性检验是一种统计方法,用于确定两个或多个样本之间是否存在显著差异。

在这个情境中,研究人员可能希望了解不同季节、地区、白天和周日的PM2.5水平是否存在显著差异。PM2.5是指空气中直径小于或等于2.5微米的颗粒物,它对人体健康和环境质量有重要影响。

为了进行显著性检验,研究人员可能会收集来自不同季节、地区、白天和周日的PM2.5水平数据,并使用统计方法来比较这些数据。常用的显著性检验方法包括t检验、方差分析和卡方检验等。

通过进行显著性检验,研究人员可以确定不同季节、地区、白天和周日的PM2.5水平是否存在显著差异,从而为环境保护和公共健康管理提供重要的科学依据。

Season

   
 
kruskal.test(PM2.5.ug.m3. ~ Season, data = airquality)

image.png

Region

image.png

Daytime

image.png

Season

   
pairwise.t.tes=="summ

image.png

image.png

image.png

image.png

image.png

image.png

image.png

Region

image.png

Day of week

   
Weekends1=seq(3,nrow(airquality),by=7)  
  
Weekends2=seq(4,nrow(airquality),by=7)

wilcox.test( airquality$PM2.5.ug.m3.[airquality$Weekends==0] , airquali

image.png

Fig 3. Day to day variations of PM2.5in different seasons, 不同季节PM2.5的逐日变化

不同季节PM2.5的逐日变化指的是在不同季节中,每天空气中PM2.5浓度的变化情况。PM2.5是指空气中直径小于或等于2.5微米的颗粒物,它们可以悬浮在空气中并被人们吸入。PM2.5来源于工业排放、交通尾气、燃煤、木材燃烧等。

在不同季节中,PM2.5的逐日变化可能会有所不同。一般来说,春季和夏季的PM2.5浓度相对较低,而秋季和冬季的PM2.5浓度相对较高。

在春季,天气温暖,风力较大,降雨较多,这些因素有助于减少空气中的PM2.5浓度。此时,植物生长活跃,能够吸收大量的二氧化碳和颗粒物,从而净化空气。

夏季的PM2.5浓度一般较低,主要是因为高温和强烈的阳光可以加速颗粒物的分解和扩散。此外,夏季也是降雨量较多的季节,雨水可以洗刷空气中的颗粒物。

秋季的PM2.5浓度开始上升,主要是因为天气变冷,大气稳定,风力较小,这有利于颗粒物的积累。此外,秋季还是农作物收获季节,农田燃烧残留物会释放出大量的颗粒物。

冬季是PM2.5浓度最高的季节。在冬季,天气寒冷,大气层稳定,风力较小,这导致颗粒物在空气中滞留时间较长。此外,冬季是取暖季节,燃煤和木材燃烧释放的颗粒物也会增加。

总的来说,不同季节PM2.5的逐日变化受到多种因素的影响,包括气温、风力、降雨量、植物生长和人类活动等。了解这些变化有助于我们采取相应的措施来减少PM2.5的污染。

   
points(which(airquality$Season=="spring" & airquality$Weekends==1),airq

image.png

image.png

image.png

image.png

不同季节PM2.5的逐日变化 Diurnal variations of PM2.5in different months

不同季节PM2.5的逐日变化是指在不同季节中,每天PM2.5浓度的变化情况。PM2.5是指空气中直径小于或等于2.5微米的颗粒物,它们对人体健康有害。这些颗粒物主要来自于工业排放、交通尾气、燃烧过程和自然来源等。

在不同季节中,PM2.5的逐日变化可能受到多种因素的影响。例如,冬季由于采暖需求增加,燃煤和燃气的使用量增加,导致PM2.5浓度较高。而夏季由于气温升高,光化学反应增多,加上交通尾气和工业排放的贡献,PM2.5浓度也可能较高。

此外,气象条件也会对PM2.5的逐日变化产生影响。例如,风速、风向、降雨等天气因素都会影响PM2.5的扩散和沉降,从而影响其浓度。

因此,不同季节中PM2.5的逐日变化是一个动态的过程,受到多种因素的综合影响。通过对这种变化的研究和监测,可以更好地了解和应对空气质量问题。

   
grep( date,airquality$日期)),]$PM2.5.ug.m3.  
plot(Mar,type="p",col=j,ylab="95% CI")

3-5月份

image.png

6-8月份

image.png

9-11 月份

1111.png

12-2 月份

image.png

Correlation coefficient and COD versus distance between the stations 相关系数和COD随站间距离的变化

相关系数是用来衡量两个变量之间的线性关系强度和方向的统计指标。它的取值范围在-1到1之间,其中1表示完全正相关,-1表示完全负相关,0表示无相关性。

相关系数和站间距离的变化之间可能存在一定的关系。如果两个变量之间的相关系数较高(接近1或-1),则说明它们之间存在较强的线性关系。在这种情况下,站间距离的变化可能较小,因为变量之间的关系比较稳定。相反,如果两个变量之间的相关系数较低(接近0),则说明它们之间的线性关系较弱或不存在。在这种情况下,站间距离的变化可能较大,因为变量之间的关系较为不稳定。

   
abline(lm(cordata~distdata),col="red")

image.png

image.png

GAMM模型分析温度  湿度,风速 能见度 气压对PM2.5的影响

GAMM(Generalized Additive Mixed Models)是一种灵活的统计模型,结合了广义可加模型(GAM,Generalized Additive Models)和混合模型(Mixed Models)的优点。它可以用于解决包含非线性关系和随机效应的数据建模问题。

GAMM模型主要用于处理长期观测数据和重复测量数据,其中数据可能受到时间、空间或其他相关因素的影响。GAMM模型具有以下特点:

  1. 广义可加模型(GAM)的优点:GAM模型可以通过非线性平滑函数来建模解释变量与响应变量之间的复杂关系,适用于非线性关系的建模。
  2. 混合模型(Mixed Models)的优点:GAMM模型可以处理数据中的随机效应,比如个体间的差异或层级结构中的变化。
  3. 可解释性:GAMM模型通过可加函数模型来描述数据,可以直观地解释模型中的每个平滑项的效应。
  4.  
   
gamm(PM2.5.ug.m3.~s(平均温度.摄

image.png

image.png

image.png

   
summary(b$gam) # gam style summary of fitted model

image.png

simple checking plots 模型检验

image.png

用AR1残差项拟合模型

image.png

lme拟合的细节

image.png

image.png

拟合模型的Gam风格总结

image.png

模型检验

image.png


QQ截图20230122145443.png

最受欢迎的见解

1.R语言多元Logistic逻辑回归 应用案例

2.面板平滑转移回归(PSTR)分析案例实现

3.matlab中的偏最小二乘回归(PLSR)和主成分回归(PCR)

4.R语言泊松Poisson回归模型分析案例

5.R语言混合效应逻辑回归Logistic模型分析肺癌

6.r语言中对LASSO回归,Ridge岭回归和Elastic Net模型实现

7.R语言逻辑回归、Naive Bayes贝叶斯、决策树、随机森林算法预测心脏病

8.python用线性回归预测股票价格

9.R语言用逻辑回归、决策树和随机森林对信贷数据集进行分类预测

标签:显著性,模型,颗粒物,PM2.5,逐日,GAMM,浓度
From: https://www.cnblogs.com/tecdat/p/18299213

相关文章

  • 【Mathematical Model】基于Python的相关性/显著性分析&成图
        很久之前编写的代码了,当时是用来分析遥感波段组合对于某地物反演的相关性分析。今天正好整理数据时一块分享出来。原创作者:RS迷途小书童博客地址:https://blog.csdn.net/m0_56729804?type=blog1相关性的概念        “相关性”是统计学中的一个基本......
  • Python统计实战:一题搞定一元线性回归的回归系数、显著性及预测值计算
    为了解决特定问题而进行的学习是提高效率的最佳途径。这种方法能够使我们专注于最相关的知识和技能,从而更快地掌握解决问题所需的能力。(以下练习题来源于《统计学—基于Python》。联系获取完整数据和Python源代码文件。)练习题随机抽取10家航空公司,对其最近一年的航班准点率......
  • 显著性目标检测(弱监督):Mutual Information Regularization for Weakly-supervised
    文章地址:MutualInformationRegularizationforWeakly-SupervisedRGB-DSalientObjectDetection|IEEEJournals&Magazine|IEEEXplore摘要:1.引入一个信息上界和一个互信息最小化正则项,鼓励每个模态的解纠缠表示用于SOD;2.运用非对称的特征提取器;3.引入多模态的......
  • 探索气象数据的多维度三维可视化:PM2.5、风速与高度分析
    探索气象数据的多维度可视化:PM2.5、风速与高度分析摘要在现代气象学中,数据可视化是理解复杂气象模式和趋势的关键工具。本文将介绍一种先进的数据可视化技术,它能够将PM2.5浓度、风速和高度等多维度数据以直观和动态的方式展现出来。一、视角与复位气象数据可视化首先需......
  • 【python】中国五大城市PM2.5 数据分析(含数据集及完整代码)
    目录问题描述(数据集、项目代码已资源绑定)前瞻 数据预处理数据可视化北京上海广州成都沈阳总览对比图北京部分代码问题描述(数据集、项目代码已资源绑定,见文章开头) 数据集前瞻         随着工业化和城市化的迅速发展,全球多个城市都面临着空气污......
  • R语言空气污染数据的地理空间可视化和分析:颗粒物2.5(PM2.5)和空气质量指数(AQI)|附代码数
    原文链接:http://tecdat.cn/?p=23800最近我们被客户要求撰写关于空气污染数据的研究报告,包括一些图形和统计输出。由于空气污染对公众健康的不利影响,人们一直非常关注。世界各国的环境部门都通过各种方法(例如地面观测网络)来监测和评估空气污染问题介绍全球的地面站及时测量了许......
  • 基于单片机环境监测温湿度PM2.5系统设计
    **单片机设计介绍,基于单片机环境监测温湿度PM2.5系统设计文章目录一概要二、功能设计设计思路三、软件设计原理图五、程序六、文章目录一概要  基于单片机环境监测温湿度PM2.5系统是一个集成了传感器技术、单片机控制以及数据处理与显示功能的综合性系统......
  • 基于51单片机的空气质量检测【温湿度,PM2.5,上位机】(仿真)
    设计一个简单、实用的基于单片机的空气质量检测仪。仪器特点:①pm2.5检测功能:当空气质量高于设置最小值低于最大值时,亮黄灯提示轻度污染,当高于设定值最大值时,亮红灯通过蜂鸣器警告,同时启动风扇进行排风;②温度检测功能:当环境温度低于设定最小值时,亮红灯蜂鸣器警告,高于最大设定......
  • [Unity/Gamma校正]通俗解释Gamma空间
    经常接触纹理制作或是shader编写应该多多少少都听说过Gamma校正Gamma校正本身很好理解,\(Gamma\=\Linear^{2.2}\)通过\(pow2.2\)这条曲线,将线性的颜色值映射到非线性 耳熟能详的案例是中灰值 人眼对灰度的感知不是线性的,如果线性地显示深度会觉得中灰偏亮这是因为人......
  • 计算降水和ENSO指数的相关系数或者回归系数,并做显著性检验
    '''Description:计算降水和ENSO指数的相关系数或者回归系数,并做显著性检验-----------------------------------------Time:2024/02/1910:42:04Author:ForxdVersion:1.0'''#%%importxarrayasxrimportnumpyasnpi......