¥È¥Ã¥× º¹Ê¬ °ìÍ÷ ¥½¡¼¥¹ ¸¡º÷ ¥Ø¥ë¥× RSS ¥í¥°¥¤¥ó

BugTrack-£ÒÈ÷˺Ͽ/39

£ÒÈ÷˺Ͽ /¾õÂÖ¶õ´Ö¥â¥Ç¥ê¥ó¥°/donlp2/¤½¤Î¾¤Î¥á¥â

£ÒÈ÷˺Ͽ - µ­»ö°ìÍ÷

¶õ´ÖÅý·× ­¢ - ¥é¥¤¥Ö¥é¥ê¡Ú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