空间分析与R语言系列之聚集性探测中的Besag&Newell方法
2016...
上次为大家讲解了一下空间聚集性探测中应用较广的空间扫描统计量,其实常用的空间聚集性探测方法还有本次要讲的Besag&Newell方法。
Besag&Newell方法主要用于罕见病的研究,它会首先定义一个固定的病例数k,然后以发生病例的每个区域为中心,按距离远近收集邻近区域的病例,计算累积病例数不少于k所需的区域数。所需区域数越少的位置,越有可能存在聚集的位置。下面通过一个实例来说明在R语言中如何利用Besag&Newell方法实现空间聚集性探测。
数据依然采用maptools工具包自带的sids数据,包括美国北卡罗来纳州1974年100个地区的经纬度坐标、婴儿猝死综合症病例数据和出生人口数据(Cressie and Chan 1989, and Cressie 1993)。所用R包为RCluster,要求的数据格式包括婴儿猝死综合症病例、预期病例数(各地区风险人口数*卡罗莱纳州病例总数/卡罗莱纳州总风险人口数)、风险人口数和坐标信息。
数据获取代码:
###所用数据是maptools工具包自带的sids数据
nc_file
关注 空间分析
微信扫一扫关注公众号