[R-br] Como fazer predição para pontos não amostrados?

Jacqueline Cantú jacqueline.cantu em hotmail.com
Segunda Julho 28 10:20:52 BRT 2014


Olá pessoal..Tenho outra dúvida..Construí o semivariograma cruzado, ajustei um modelo e quero fazer a predição para pontos que não foram amostrados para isso construí uma malha fina...Minha dúvida é: como fazer a predição espacial para os pontos não amostrados?Lembrando que são 2 variáveis, sendo que para predizer os valores de uma variável será utilizado os valores da outra variável  (cokrigagem)
Segue rotina..  require (geoR)library(sp)require(gstat)require(splancs)require(classInt)
#lendo os dadosh <- read.table("phca.txt ",head=T)names(h)class(h)coordinates(h) = ~x+y
g = gstat(NULL, "ph", ph~1, h)g = gstat(g, "ca", ca~1, h)plot(variogram(g))vm = variogram(g)variogram(g)
##Ajuste do modelo esfericosph <- fit.lmc(vm,g, vgm(0.3, "Sph", 200, 0.2))sph
#modelo ajustado plot(vm,sph)
##Predição para pontos não amostrados#construção de uma malha finaborda <- read.table("borda.txt") borda plot(borda)polygon(borda)apply(borda,2,range) gr1 <- expand.grid(x=seq(239182.8, 241027.9, by=5),y=seq(7236744, 7238593, by=5))gi <- polygrid(gr1, bor=borda)plot(gi)
Att, Jacqueline







 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140728/6cac5cf0/attachment.html>


Mais detalhes sobre a lista de discussão R-br