[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