//¤³¤Î¥Ú¡¼¥¸¤¬É½¼¨¤µ¤ì¤ëÊý¤Ï¡¢URL¤«¤é¡Öaction=SOURCE&¡×¤òºï½ü¤·¤Æ¤ß¤Æ¤¯¤À¤µ¤¤
[[£ÒÈ÷˺Ͽ - µ­»ö°ìÍ÷]]

!!!¶õ´ÖÅý·× ­¢ - ¥é¥¤¥Ö¥é¥ê¡Úgstat¡Û¤ò»È¤¦
*Åê¹Æ¼Ô¡§ ¤ß¤å
*¥«¥Æ¥´¥ê¡§ ¤Ê¤·
*Í¥ÀèÅÙ¡§ ÉáÄÌ
*¾õÂÖ¡§ ´°Î»
*Æü»þ¡§ 2009ǯ06·î30Æü 09»þ05ʬ30ÉÃ
//{{bugstate}}
!!ÆâÍÆ
*¶õ´ÖÅý·× ­¢
*¥é¥¤¥Ö¥é¥ê¡Ú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))
{{ref_image ss_2_1.png}}
!!¶õ´Ö¤ÎÁê´Ø¡Ê¶õ´Ö½¾Â°À­¡Ë¤òÄ´¤Ù¤ë
*[¶õ´ÖÅý·× ­¡|http://phoenixx.sakura.ne.jp/R/wiki.cgi?page=BugTrack%2D%A3%D2%C8%F7%CB%BA%CF%BF%2F38]¤Ç½Ò¤Ù¤¿Íͤʡ¢Æó¼¡Äê¾ïÀ­¤ÈÅùÊýÀ­¤ò²¾Äꤷ¡¢¶õ´Ö¤ÎÁê´Ø¤òÄ´¤Ù¤ë¡¥
 x <- variogram(object=log(zinc)~1, locations=~x+y, data=meuse, cutoff=3000, cloud=TRUE)
 plot(x)
{{ref_image ss_2_2.png}}
 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)
{{ref_image ss_2_3.png}}
*¥Õ¥£¥Ã¥È¤µ¤»¤ë
 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
{{ref_image ss_2_4.png}}
!!¥¯¥ê¥®¥ó¥°
 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))
{{ref_image ss_2_5.png}}
!!NEXT
*[[BugTrack-£ÒÈ÷˺Ͽ/40]] ¡§ ¶õ´ÖÅý·× ­£ - ¡Úgstat¡ÛÉáÊ×·¿¥¯¥ê¥®¥ó¥°