[R-br] Como fazer cokrigagem
Jacqueline Cantú
jacqueline.cantu em hotmail.com
Quarta Julho 2 14:49:22 BRT 2014
Muito obrigada Luis Paulo Braga.
Att, Jacqueline
Date: Wed, 2 Jul 2014 13:29:18 -0300
From: lpbraga em geologia.ufrj.br
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br] Como fazer cokrigagem
Co-Krigagem
(Modelo Linear de Co-Regionalização)
A
Co-Krigagem é recomendada quando uma variável secundária correlacionada
espacialmente à variável primária contribui para a melhoria das predições. Isso
ocorre particularmente para o caso em que o modelo para cada variável e para a
variável cruzada é um Modelo Linear de Corregionalização.
data(meuse)
g <-
gstat(id = "ln.zinc", formula = log(zinc)~1, locations = ~x+y, data =
meuse)
g <-
gstat(g, id = "ln.lead", formula = log(lead)~1, locations = ~x+y,
data = meuse)
#carregou os dados de zinco e chumbo
vm<-variogram(g)
#obteve os semivariogramas
experimentais univariados e cruzados
vm.fit <-
fit.lmc(vm,g, vgm(1, "Sph", 800, 1))
#ajustou
plot(vm,vm.fit)
ck<-predict(vm.fit,meuse.grid)
spplot(ck["ln.zinc.pred"],col.regions=bpy.colors(20),
key.space=list(x=0.1,y=.95,corner=c(0,1)),
main = "co kriging predictions(ln(Zn))", xlab="Xloc",ylab=
"Yloc ")
Em 2 de julho de 2014 10:52, Jacqueline Cantú <jacqueline.cantu em hotmail.com> escreveu:
Olá pessoal,
Alguém saberia me ajudar em como fazer
cokrigagem..
Deste já agradeço a atenção.
Grata, Jacqueline
_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140702/015671dd/attachment.html>
Mais detalhes sobre a lista de discussão R-br