[R-sig-Geo] Robust estimation with INLA

Lucy Mburu mburulou at yahoo.com
Wed Sep 2 01:08:28 CEST 2015


Hi all,

I am trying to fit a spatial model with INLA, but I need to incorporate 
the weights of influential points. I have estimated the weights after 
running a robust regression model and stored the list of weights as a 
named vector. Is there a possibility for including this vector into the 
call for slm.inla?

Using the columbus example:

mod1<-rlm(CRIME~INC+HOVAL,data=columbus,psi="psi.bisquare")
xweight<-mod1$w

How can I now incorporate "xweight" into the example on page 13 of the 
INLABMA documentation 
(https://cran.r-project.org/web/packages/INLABMA/INLABMA.pdf)?


Thanks for your ideas.

Lucy Mburu.



More information about the R-sig-Geo mailing list