R語言計算moran I

install.packages("

maptools
")#畫地圖的包
install.packages("
spdep
")#空間統計，moran'i
install.packages("
tripack")
install.packages(
"rann")
library(
"maptools")
library(
"spdep")
library(
"tripack")
library(
"rann")
#讀空間資料
rdata <- readshapepoly("
d:/data/cairo.shp
")#讀取shp面資料
names(rdata)#
顯示資料字段名字
head(rdata@data)#
顯示屬性表
#按照字段兩幅畫圖
spplot(rdata[c("
tfr96_03
","tfr86_03
") ],
main = "
spatial distribute of tfr96_03
",#圖名
xlab = "
x coords
",#橫座標名字
ylab = "
y coords
",#縱座標名字
cut = 30#分段)
#按照字段一幅畫圖
spplot(rdata["
tfr96_03"],
main = "
spatial distribute of tfr96_03
",#圖名
xlab = "
x coords
",#橫座標名字
ylab = "
y coords
",#縱座標名字
cut = 30#分段)
#按邊鄰接角鄰接生成鄰居（方式一）
queen_nb <- poly2nb(rdata, queen = true)#
有8個rook_nb <- poly2nb(rdata, queen = false)#
有4個#
獲取中心點座標編號
coords <-coordinates(rdata)
#地圖資料轉資料框-->生成id
ids <-row.names(as.data.frame(rdata))
#設定畫引數
oopar <- par(mfrow = c(1,2),
mar = c(3,3,1,1)+0.1)
#畫邊界
plot(rdata, border = "
grey
", main = "
queen-style")
#add修改上乙個圖
plot(queen_nb, coords, col = "
dodgerblue",
add = true, pch = 19, cex = 0.5)
#生成鄰接關係（方式二）knearneigh()定義k-near
k4_nb 
<- knn2nb(knearneigh(coords, k = 4), row.names =ids) 
#識別鄰接關係是否對稱
is.symmetric.nb(k4_nb, verbose = false, force =true)
#補全鄰接關係 
n.comp.nb(k4_nb)$nc k4_w <-nb2listw(k4_nb)
#計算moran'i 
moran.test(rdata$tfr96_03, listw =k4_w) 
#計算moran'i (蒙特卡洛方法)
moran.mc(rdata$tfr96_03, listw = k4_w, nsim = 999)

GeoDa計算全域性Moran I

1 匯入包含資料的.shp檔案 2 建立權重矩陣，點選weight manger，再點選create，weights file id variable 其中包含的數值要唯一 或者 add id variable.選擇計算的距離方式 3 接下來space univariate moran s i 4 ...

R語言平行計算程式設計

r語言並行程式設計例項，將資料集rcuters 21578的50個文件賦值100,000次，得到包含500萬個文件的資料庫，查詢正規表示式。tm軟體包是r語言中為文字挖掘提供處理的package，提供資料匯入 語料庫處理 預處理 元資料管理 建立term document矩陣等功能。平行計算 lib...

R語言樣本量計算

連續型變數 統計功效與樣本量間資料視覺化 一般的研究設計中，臨床結局變數可分為連續性或二分型別，不同型別的設計往往對應不同的結局變數。樣本量計算方法也有所不同。二分類結局是指結局為二分類且二者互斥，如實施某種 後，患者是否 況。服藥後病情是否改善。給定組1 n 28 某事件發生概率為30 組2 n ...