£ÒÈ÷˺Ͽ - µ»ö°ìÍ÷
- Åê¹Æ¼Ô¡§ ¤ß¤å
- ¥«¥Æ¥´¥ê¡§ ¤Ê¤·
- Í¥ÀèÅÙ¡§ ÉáÄÌ
- ¾õÂÖ¡§ ´°Î»
- Æü»þ¡§ 2009ǯ06·î30Æü 09»þ05ʬ30ÉÃ
- ¶õ´ÖÅý·× ¢
- ¥é¥¤¥Ö¥é¥ê¡Úgstat¡Û¤ò»È¤¦
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))
£ÒÈ÷˺Ͽ /¾õÂÖ¶õ´Ö¥â¥Ç¥ê¥ó¥°/donlp2/¤½¤Î¾¤Î¥á¥â