トップ 差分 一覧 ソース 検索 ヘルプ RSS ログイン

BugTrack-R備忘録/39

R備忘録 /状態空間モデリング/donlp2/その他のメモ

R備忘録 - 記事一覧

空間統計 - ライブラリ【gstat】を使う

  • 投稿者: みゅ
  • カテゴリ: なし
  • 優先度: 普通
  • 状態: 完了
  • 日時: 2009年06月30日 09時05分30秒

内容

  • 空間統計
  • ライブラリ【gstat】を使う

メウス川のデータ

  • meuse
library(gstat)
library(lattice)
# - メウス川
data(meuse)
data(meuse.grid)
levelplot(log(zinc)~x+y, data=meuse, cuts=99, col.regions=rainbow(100,start=0.5))

空間の相関(空間従属性)を調べる

  • 空間統計 で述べた様な、二次定常性と等方性を仮定し、空間の相関を調べる.
x <- variogram(object=log(zinc)~1, locations=~x+y, data=meuse, cutoff=3000, cloud=TRUE)
plot(x)
x <- variogram(object=log(zinc)~1, locations=~x+y, data=meuse, cutoff=3000)
model1 <- vgm(psill=0.5, model="Sph", range=600, nugget=0.1)
plot(x, model=model1)
  • フィットさせる
model2 <- fit.variogram(object=x, model=vgm(psill=0.5, model="Sph", range=600, nugget=0.1))
plot(x, model=model2)
model2
  model      psill    range
1   Nug 0.06268121   0.0000
2   Sph 0.55768569 884.9763

クリギング

g <- gstat(id="logzn", formula=log(zinc)~1, locations=~x+y, data=meuse, model=model2)
xo <- predict.gstat(g, meuse.grid[,1:2])
levelplot(logzn.pred~x+y, data=xo, cuts=99, col.regions=rainbow(100,start=0.5))

NEXT