白話空間統計十二 R語言對點資料分析的實現(1)

2021-07-04 18:43:05 字數 2812 閱讀 9680

前段時間,有人批評我寫白話空間統計的部落格是在寫軟文給

arcgis

平台和esri

打廣告,話說這個實在是太恭維我了。如果讀到空間統計,而且還能讀懂的人,不可能沒有聽說過

arcgis

軟體吧,這種情況到底是先有雞還是先有蛋,自然一目了然了。

雖然蝦神在

esri

中國幹了好多年了,給公司打打廣告也是理所當然的事情。但是寫部落格的時候確實還真沒有這個想法,不過既然有人批評了,那麼蝦神我先挑明一下這系列白話空間統計的整個寫作背景:

首先肯定是對空間統計進行乙個科普,畢竟這個東西比較冷門,能找到的資料不多,所以我也萌發了寫個部落格來科普一下的想法,而且也是對自己知識的乙個整理的過程。

其次,我盡量的把所有演算法都描述了一下,如果編碼能力比較強的,我覺得按照演算法都可以自己編寫分析的過程的,也就不必要用到任何的軟體。

第三,製作示例,作圖和出示分析結果的時候,為了省事,確實用的是

arcgis

軟體,先不說這是蝦神我的份內之事,再說「師夷長技以制夷」嘛,沒必要上綱上線。

第四,整個寫作的動力,都是自己熱血沸騰,心潮澎湃的結果,如果這是公司給的任務,我大週末的蹲家裡幹點啥不行?再說用業餘時間寫作能給蝦神我算加班麼?

俗話說:知錯能改,善莫大焉。既然有人批評了,我就有則改之,無則加勉,今天內容就不用

arcgis

了,我們用

r語言來實現一下我前面所說的那些方法和內容。

今天先來說說用

r語言來實現這兩天說的幾個點資料的分析演算法,至於莫蘭指數什麼的,以後有機會也會慢慢得講講如何用

r語言來實現。

當然,為了省事,我就不用介紹

r語言的基本使用和安裝了,網路上成堆(ps

r語言的包非常多,目前公開的已經超過了

6000

個,你能想到的分析,它肯定有,你沒想到的,它自然也有。

今天我們介紹乙個用於點分析的包:

aspace

這個包是加拿大的多倫多大學的randybui, ron n. buliung, tarmo k. remmel在

2012

年公布的,專用於對點資料的探索和分析以及視覺化的專用包。

這個包主要可以用於對空間資料的中心模式相關的計算(包括了標準距離,標準差橢圓,中心點、中位數等功能)。而且他最大的特點是他支援對

esri

的shapefile

檔案的讀寫(這個真心不是打廣告),也就是說,他可以把計算出來的結果寫成

shape file。

今天我們先對他的基礎應用進行乙個簡單的了解,如下:

首先對這個包要進行安裝,安裝的方法就直接通過

r語言提供的安裝指令碼實現就可以了:

我們直接用上次做的長江太湖流域傷寒的病例資料,資料描述如下:

將它匯出成

shape file

資料,然後在

r語言裡面來使用。

首先載入

aspace

包,命令如下:

>library(aspace)

載入完成之後,會提示一堆東西,表示你目前要使用這個包需要的各種環境,如下:

然後因為我們要讀取

shapefile

進行計算,最簡單的

shapefile

讀取包就是

r語言裡面的

maptools

包,當然,我也要按照和載入這個包,如下:

安裝完了之後,載入到

r環境中去:

讀取shapefile

檔案為點資料集,然後顯示出來看看:

然後我們來計算一下平均中心:

> xy

>mean_centre(id=1,filename="e:/output.txt",

+ weighted=false,weights=null,points=xy)

結果如下:

平均中心就計算出來了,至於

r語言裡面的各個語句是什麼意思,大家有興趣的可以查閱

r語言的語法,當然,這個

aspace

包功能還很多,今天我們就簡單的做了個

hello world

,下次我們繼續說它的其他功能。

R語言空間面板資料統計分析

install.packages splm library splm library spdep library sp library matrix library sf library spdata install.packages ecdat data produc,package ecdat ...

通過R語言對資料進行聚類分析

聚類分析是將資料根據相似特徵分類到不同的類或簇,用以對資料進行區分。1.聚類 k means 通過資料中每個數值間的距離進行聚類,是一種縱向聚類,將樣本資料據成型別,用以對資料的特徵進行區分。1 讀取資料 data read.csv cludata data complete.cases data ...

R語言對爬取的天氣資料進行整理

1.批量輸出每個省的每個月份出現頻次最多的天氣情況 2.批量輸出每個省每個月的氣溫狀況 包括對類似5 10 資料的處理 天氣情況彙總 rm list ls all true gc library readxl library lubridate library data.table library ...