数据分享|R语言生态学种群空间点格局分析:聚类泊松点过程对植物、蚂蚁巢穴分布数据可视化...

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

点模式分析(点格局分析)是一组用于分析空间点数据的技术。在生态学中,这种类型的分析可能在客户的几个情境下出现,但对数据生成方式做出了特定的假设,因此让我们首先看看哪些生态数据可能与点模式分析相关或不相关点击文末“阅读原文”获取完整代码数据)。

相关视频

哪些数据适用于点模式分析?

点模式分析的最重要假设是点的数量和位置都需要是随机的。此外,我们需要知道客户的采样区域(即所谓的窗口)。以下是适用于点模式分析的示例:

  • 森林样地中树木的位置

  • 草地区域中蚂蚁巢穴的分布

以下是不适合进行点模式分析的示例:

  • 在较大的样地中事先定义好的子样方形成的规则网格中的群落组成

  • 记录了给定区域内固定数量的鸟巢的位置

以下示例可能适用于点模式分析,也可能不适用:

  • 动物移动的无线电跟踪数据(请参阅针对此特定类型数据的众多技术)

  • 每年记录森林样地中树木位置,形成复制的点模式

在 R 中进行点模式分析

我将使用一个数据集来展示分析。该数据集包含沿海沙丘系统中植物和蚂蚁巢穴的位置查看文末了解数据免费获取方式

从现有数据创建点模式

#加载数据集
dat <- read.table("daf.csv",sep=";",head=TRUE)
#将坐标转换为米
dat$X <- dat$X/100
dat$Y <- dat$Y/100#创建点模式
(all_pp <- ppp(

3da8fb97a9498e013c19e68a3c355839.png

plot(

11d3ba68fb37294b961c16927bb3c800.png

对点模式对象进行数据操作

#一个首次操作是给每个点添加信息,即所谓的标记
#在这个示例中,我们可以添加植物的物种名称
marks(all_pp) <-#第二个操作可以是删除任何重复点
#我们可以使用简单的行索引来做到这一点
all_pp <- unique#然后添加坐标单位
summary(all_pp)

1285db628c6bd73a1bcf2d6b2dc95b73.png

#我们可以使用标记对点模式进行子集化
ant_pp <- subset(all_pp,marks=="Tetramorium_caespitum")
#在这种情况下,我们不再需要标记
ant_pp <- unmark(ant_pp)

标记的概念非常重要,所以我会多花几句话来介绍。标记可以是与点模式长度相同的数字或因子向量,它们是为每个点收集的额外信息。在本示例中,这是记录的植物和蚂蚁的物种名称,但也可以是树木高度或鸟巢中的蛋数。

基于窗口的第二个有趣的操作集是根据特定窗口对点模式进行子集化:

w <- hexon(entre=c(5,5))
plot(antp[w])

b9572a957cb70cc8c1b2dd2cc8dff180.png

#基于物种名称拆分
split_pp <- si(all_pp)
class(spitpp)
## [1] "splitppp" "ppplist"  "solist"   "list"
as.matrix(lapplysplit_p,npis),ncol=1)

22546b0e9c0564c66786e8a1d10ed894.png

#也可以使用:by(all_p,maksall_p,npoints)来拆分#基于窗口进行拆分
spl_ant_pp <- splitn_p,f=w)
summary(plt_at_pp)

1f83ea0f3bc59f8d64ed2b9a25fb9143.png

d637d982d0aca1bae2bdb0c45e1274c3.png

点格局的探索性分析

这是任何点格局分析中非常重要的一步,这一步可以帮助你:(i) 探索密度,以及 (ii) 观察点格局是否偏离随机期望。

den_all <- desit(slit_pp)
plot(dns_ll)

432c9127cf947bb4591bb42980cf1515.png

首先要找出点格局是否由一个强度函数生成,如果是,则点格局是均匀的;如果点格局是由多个强度函数生成的,则点格局是不均匀的。这是点格局分析中的一个重要的第一步,因为大多数函数和模型默认假设点格局是均匀的。我将展示两种推断点格局均匀性的方法:(i) 模拟和 (ii) 方块计数。


点击标题查阅往期内容

2b45058ab909779775833006c0db2ab0.jpeg

生态学建模:增强回归树(BRT)预测短鳍鳗生存分布和影响因素

outside_default.png

左右滑动查看更多

outside_default.png

01

outside_default.png

02

outside_default.png

03

outside_default.png

04

outside_default.png

第一种方法是根据观察到的点格局的平均强度,在空间上完全随机模拟点格局。如果观察到的点格局和模拟的点格局的密度估计相似,那么我们就有证据表明点格局是均匀的。

#将观察到的密度与基于强度的随机模拟进行比较#在图中为观测数据选择一个随机位置
pos <- sale(1:16,1)
#模拟15个CSR点
simp <- roisp(lmbda= intesityant_pp,win  Windw(nt_p),nm=15)
#将模拟集合中的第pos个位置替换为观察数据集
tmp <- sp[pos]]
sim[pos]] <- ant_pp
sip16]] <- tmp
naes(imp)[16] <- "Siulio 16"
#计算密度估计
densp <- dsiymp)
#绘图,你能分辨出哪一个是观测数据集吗?
pr(mfow=c(,4)mar=c(0,0,4,4))
plt(as.isto(densp),zlim=rage(unlst(lapl(despran))))

outside_default.png

如果你能在模拟数据中找到真实数据集,那么就证据表明一种不均匀的过程生成了数据。在这种情况下,你需要使用特殊的函数。第二种方法是将窗口划分成方块,并统计每个方块中的点的数量。通过卡方检验,我们可以推断点格局是均匀的(p > 0.05)还是不均匀的(p < 0.05)。

quar.tst(ant_pp)

outside_default.png

输出结果告诉我们,点模式的零假设是完全随机生成的空间过程,被拒绝了。我们有一些证据表明,点模式是不均匀或非平稳的。如果点模式遵循完全随机性(CSR),则这个计数值(K)和考虑的距离(r)之间存在已知的关系。在R中,对蚂蚁巢穴点模式执行这个操作的代码是:

ee_iso <- eneope(antpun=Kst, nar=liscrectio="rder"), gloal=TRUE

outside_default.png

plot(eeso)

outside_default.png

在这里,我围绕模拟的随机点过程的期望(Ktheo(r))从中推导出一个包络。然而,正如我们之前所看到的,我们有一些证据表明蚂蚁巢穴的点模式是不均匀的,因此我们应该通过使用Ripley的K函数的修改版本来考虑这一点:

eeihm < elop(ant_pp, fun=Khom,glba=TUE)

outside_default.png

plot(ee_iom)

outside_default.png

这次观察到的曲线在较大距离上低于预期曲线,说明巢穴中的分散程度超过了考虑点模式非均匀性时的预期CSR。

从这个探索性分析中我们得出以下结论:

  • 蚂蚁巢穴显示出不均匀性模式

  • 有一些证据表明在较大距离上,蚂蚁巢穴的间距比从估计的强度函数中预期的要大。

现在我们可以进入下一步,对我们的点模式进行建模。

构建点模式模型

在我们的蚂蚁巢穴数据示例中,我们可能对巢穴密度是否取决于特定植物物种的密度感兴趣。

#拟合一个仅包含截距的泊松点过程模型
m0 <- ppm(tp ~ 1)
m0

outside_default.png 这是最简单的拟合模型,该模型基本上告诉我们强度(蚁巢的密度)在观察区域内始终为e^(-0.36)=0.70。现在我们可以使用坐标作为预测变量来建立一个模型:

m1 <- ppm(tpp~ plnm(x,y,2))
m1

outside_default.png

plt(m,se=FS,ho="iae")

outside_default.png

该模型拟合了以下关系:log(λ)=x+y+x^2+y^2+x:y,基本上是每个坐标轴的二次关系加上一个交互项。图示了模型预测的强度(λ)与观察到的蚁巢的关系。这个模型在一些区域拟合效果较差。对于每个模型,重要的下一步是模型验证,可以绘制许多重要的模型诊断图:

diospm(m1whih = "ooh")

outside_default.png

outside_default.png

默认情况下,生成四个图,我只请求绘制平滑残差图,以便我们可以确定模型在观察到的点格局上拟合得不好的区域。该模型在某些区域的拟合效果较差。也可以使用拟合强度在Kinhom函数中查看观察到的点格局是否比模型拟合预期的更或更少聚集:

e <- envepe(an_ppKinom,unargs  list(lda=1),obl=TUE)

outside_default.png

plt(em)

outside_default.png

在这里,我们可以看到观察到的点格局比模型预期更加聚集。一种解决方法是使用聚类泊松点过程模型:

m2 <- kp(ntp ~ lyo(x,,2))
plotm,wht="satistc",pase=ASE)

outside_default.png

虚线绿线显示了基于模型中的预测变量的预期K值,实线黑线将拟合的聚类过程添加到预测变量中(在这种情况下是一个Thomas过程。从配适模型中模拟点格局也很容易,我们将在这里使用它来查看观察到的点格局和模拟的点格局之间是否存在明显差异:

#一个随机位置
pos <- sample(1,1:16)
#从该模型中模拟15个点格局
sims <- (2,nim = 15)

outside_default.png

#将观察到的点格局放在随机位置
mp < sims[[pos]]
sims[s<- tmulation 16"
#计算密度估计
den <- nsity(sims)
#绘图,你能分辨出哪一个是观察到的数据集吗?
par(row=c(4,),mar=c(0,0,4,4))
plot(aslistf(dnsp)zliangunlis(lpply(nsange))))

outside_default.png

我无法真正辨认出观察到的格局,因此该模型相当好。

点格局的预测变量也可以是像素图像(或"im")对象,我们将使用一个植物物种Senecio_inaequidens的密度作为预测变量:

m3 <- kpm(an_p ~ Seco_jcobaa,dta=desall)
#让我们来看一下预期的K值
plt(m3,wat="aisti",pus=FALSE)

这个模型显然比之前的模型好。可以使用effectfun函数绘制协变量的效果,例如:

#查看Senecio jacobaea的效应
plot(effetfun(m3Seecio_jacaea"))

outside_default.png

这个模型显然比之前的模型好。可以使用effectfun函数绘制协变量的效果,例如:

#查看Senecio jacobaea的效应
plot(efcfun(m3"Seeciojacobea"))

outside_default.png

可以绘制更酷炫的地图:

#绘制一个有趣的透视图
pp <- predctm3)
M<-persp(s_al$,lin=pox=FALSE,visble=TREapron = TRUE,et55,ph=,exnd=6,mi=" density")
perspPins(ant_pp,Zdn_lSeneco_cobaea,M=M,pc=20)

outside_default.png

图中的高度表示Senecio jacobaea苗的密度,颜色表示蚁巢的拟合强度,点表示实际观察到的蚁巢。

数据获取

在公众号后台回复“蚂蚁”,可免费获取完整数据。

outside_default.png

本文中分析的数据分享到会员群,扫描下面二维码即可加群!

outside_default.png


outside_default.png

点击文末“阅读原文”

获取全文完整代码数据资料。

本文选自《R语言生态学种群空间点格局分析:聚类泊松点过程对植物、蚂蚁巢穴分布数据可视化》。

outside_default.png

outside_default.png

点击标题查阅往期内容

Python决策树、随机森林、朴素贝叶斯、KNN(K-最近邻居)分类分析银行拉新活动挖掘潜在贷款客户

R语言逻辑回归(Logistic Regression)、回归决策树、随机森林信用卡违约分析信贷数据集

R语言基于Bagging分类的逻辑回归(Logistic Regression)、决策树、森林分析心脏病患者

R语言样条曲线、决策树、Adaboost、梯度提升(GBM)算法进行回归、分类和动态可视化

R语言用主成分PCA、 逻辑回归、决策树、随机森林分析心脏病数据并高维可视化

matlab使用分位数随机森林(QRF)回归树检测异常值

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

R语言中使用线性模型、回归决策树自动组合特征因子水平

R语言中自编基尼系数的CART回归决策树的实现

Python对商店数据进行lstm和xgboost销售量时间序列建模预测分析

R语言基于树的方法:决策树,随机森林,Bagging,增强树

R语言实现偏最小二乘回归法 partial least squares (PLS)回归

R语言多项式回归拟合非线性关系

R语言逻辑回归(Logistic回归)模型分类预测病人冠心病风险

R语言用局部加权回归(Lowess)对logistic逻辑回归诊断和残差分析

R语言混合效应逻辑回归(mixed effects logistic)模型分析肺癌数据

outside_default.png

outside_default.png

outside_default.png

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/139893.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【LeetCode】——双指针(快慢指针)/多指针

个人主页 代码仓库 C语言专栏 初阶数据结构专栏 Linux专栏 前言 大家好&#xff01;这是新开的LeetCode刷题专栏&#xff0c;这个专栏不只是随便的拿一些我练过的题讲解&#xff0c;而是总结我在刷题中的一些方法适用于一大类的题&#xff0c;是给大家提供这一大类题的解题…

爬虫使用Selenium生成Cookie

在爬虫的世界中&#xff0c;有时候我们需要模拟登录来获取特定网站的数据&#xff0c;而使用Selenium登录并生成Cookie是一种常见且有效的方法。本文将为你介绍如何使用Selenium进行登录&#xff0c;并生成Cookie以便后续的爬取操作。让我们一起探索吧&#xff01; 一、Seleni…

深度学习自学笔记一:神经网络和深度学习

神经网络是一种模拟人脑神经元之间相互连接的计算模型&#xff0c;它由多个节点&#xff08;或称为神经元&#xff09;组成&#xff0c;并通过调整节点之间的连接权重来学习和处理数据。深度学习则是指利用深层次的神经网络进行学习和建模的机器学习方法。 假设有一个数据集&a…

企业架构LNMP学习笔记53

PHP扩展安装&#xff1a; server01和server03上安装redis扩展&#xff1a; 解压编译安装&#xff1a; shell > tar xvf redis-4.3.0.tgz shell > cd redis-4.3.0 shell > phpize shell > ./configure && make && make install 配置文件php.ini&…

HTML5day02综合案例2

案例展示 <!DOCTYPE html> <html lang"en"> <head><meta charset"UTF-8"><meta name"viewport" content"widthdevice-width, initial-scale1.0"><title>注册信息</title> </head> &l…

记一次MySQL安装过程中遇到的问题

由于太久没用MySQL&#xff0c;今天在重装MySQL时遇到一个问题&#xff0c;被卡了近2个小时。。。。。。 由于我本人原先安装过MySQL&#xff0c;所以在重装的时候必须要先卸载原先的MySQL。 下面先给出正确的卸载流程&#xff08;作者就是在卸载的时候操作失误导致安装过程被…

使用python处理MNIST数据集

文章目录 一. MNIST数据集1.1 什么是MNIST数据集1.2MNIST数据集文件格式1.3使用python访问MNIST数据集文件内容 附录程序源码 一. MNIST数据集 1.1 什么是MNIST数据集 MNIST数据集是入门机器学习/识别模式的最经典数据集之一。最早于1998年Yan Lecun在论文:[Gradient-based l…

MyBatis基础之SqlSession

SqlSession 线程安全问题 当你翻看 SqlSession 的源码时&#xff0c;你会发现它只是一个接口。我们通过 MyBatis 操作数据库&#xff0c;实际上就是通过 SqlSession 获取一个 JDBC 链接&#xff0c;然后操作数据库。 SqlSession 接口有 3 个实现类&#xff1a; #实现类1Defa…

【Verilog语法】比较不同计数器的运算方式,其中有一个数是延迟打一拍的效果,目的是使得两个计数器的结果相同。

比较不同计数器的运算方式&#xff0c;其中有一个数是延迟打一拍的效果&#xff0c;目的是使得两个计数器的结果相同。 1&#xff0c;第一种2&#xff0c;第二种3&#xff0c;第三种 第三种方案&#xff0c;完成实现。 1&#xff0c;第一种 &#xff08;1&#xff09;RTL modu…

索引(含B树、B+树)

1、索引&#xff08;index&#xff09; 索引是在数据库表的字段上添加的&#xff0c;是为了提高查询效率存在的一种机制。 一张表的一个字段可以添加一个索引&#xff0c;当然&#xff0c;多个字段联合起来也可以添加索引。 索引相当于一本书的目录&#xff0c;是为了缩小扫描…

数据仓库整理

数仓 olap vs oltp OLTP主要用于支持日常的业务操作&#xff0c;如银行交易、电子商务等&#xff0c;强调数据的准确性、实时性和并发性。OLAP主要用于支持复杂的数据分析&#xff0c;如数据仓库、决策支持等&#xff0c;强调数据的维度、聚合和可视化。 将OLTP数据库的数据…

Intellij idea 2023 年下载、安装教程、亲测可用

文章目录 1 下载与安装IDEA2 常用设置设置 Java JDK 版本自动导入包、移除包IDEA 自动生成 author 注释签名java.io.File 类无法自动提示导入&#xff1f;高亮显示与选中字符串相同的内容IDEA 配置 MavenIDEA 连接 Mysql 数据库 3 参考文章 1 下载与安装IDEA 首先先到官网下载…

硬件知识积累 网口接口 百兆,千兆,万兆 接口介绍与定义 (RJ45 --简单介绍)

1. 百兆网口 1.1百兆网的定义 百兆网的意思是是100Mb/S&#xff0c;中文叫做100兆位/秒。 1.2百兆网口的常用连接器 1.1.1 一般百兆网口的连接器一般是RJ45 下面是 实物图&#xff0c; 原理图&#xff0c;封装图。 1.3 百兆网口连接线的介绍 1.3.1 百兆需要使用的线的定义 百…

基于Vue+ELement搭建登陆注册页面实现后端交互

&#x1f389;&#x1f389;欢迎来到我的CSDN主页&#xff01;&#x1f389;&#x1f389; &#x1f3c5;我是Java方文山&#xff0c;一个在CSDN分享笔记的博主。&#x1f4da;&#x1f4da; &#x1f31f;推荐给大家我的专栏《ELement》。&#x1f3af;&#x1f3af; &#x1…

代码随想录算法训练营 动态规划part17

一、回文子串 647. 回文子串 - 力扣&#xff08;LeetCode&#xff09; class Solution {public int countSubstrings(String s) {boolean[][] dp new boolean[s.length()][s.length()];int ans 0;for (int j 0; j < s.length(); j) {for (int i 0; i < j; i) {if …

stm32无人机-飞行力学原理

惯性导航&#xff0c;是一种无源导航&#xff0c;不需要向外部辐射或接收信号源&#xff0c;就能自主进行确定自己在什么地方的一种导航方法。 惯性导航主要由惯性器件计算实现&#xff0c;惯性器件包括陀螺仪和加速度计。一般来说&#xff0c;惯性器件与导航物体固连&#xf…

EMQX Enterprise 5.2 发布:Flow 设计器,Amazon Kinesis,Azure Event Hubs

EMQX Enterprise 5.2.0 版本现已正式发布&#xff01; 新版本带来了一系列重磅更新&#xff0c;最令人瞩目的是可拖拽的可视化 Flow 设计器&#xff0c;它可以帮助企业快速创建、测试和部署数据集成。同时&#xff0c;我们新增了对 Amazon Kinesis 和 Azure Event Hubs 的支持…

Flowable主要子流程介绍

1. 内嵌子流程 &#xff08;1&#xff09;说明 内嵌子流程又叫嵌入式子流程&#xff0c;它是一个可以包含其它活动、分支、事件&#xff0c;等的活动。我们通常意义上说的子流程通常就是指的内嵌子流程&#xff0c;它表现为将一个流程&#xff08;子流程&#xff09;定…

【再识C进阶3(上)】详细地认识字符串函数、进行模拟字符串函数以及拓展内容

小编在写这篇博客时&#xff0c;经过了九一八&#xff0c;回想起了祖国曾经的伤疤&#xff0c;勿忘国耻&#xff0c;振兴中华&#xff01;加油&#xff0c;逐梦少年&#xff01; 前言 &#x1f493;作者简介&#xff1a; 加油&#xff0c;旭杏&#xff0c;目前大二&#xff0c;…

Linux 链表示例 LIST_INIT LIST_INSERT_HEAD

list(3) — Linux manual page 用Visual Studio 2022创建CMake项目 * CmakeLists.txt # CMakeList.txt : Top-level CMake project file, do global configuration # and include sub-projects here. # cmake_minimum_required (VERSION 3.12)project ("llist")# I…