[R-sig-Geo] Kriging with External Drift, multivariate experience?

[Edzer Pebesma]

Tobin Cara wrote:
> Hello,
>
> Thank you all for your previous help. I began using R 2 weeks ago, and I am getting somewhere finally. I have been able to run universal kriging with a Digital Elevation Model trend. Now, does anyone have experience with kriging more than one trend variable?
>
> I assume you have to have the other variables defined at the same locations as the original grid and then at the new grid. However, can more than one new grid be defined?
>
> For example,
>
> elev_ked <- krige(meanPrec~Z+LAM, meanPrec, newdata=elev, fitted_vario)
>
> In this example, I want to krig the mean Precipitation relative to the elevation Z and a Limited Area Model (LAM). There is only one way to define the newdata that I know and in this case I am interpolating back to the Digital Elevation Model 'elev'. If I want to interpolate back to the LAM grid, I presume I change to newdata=LAM?
>
> However, is there an approach that uses the entire newdata grid of elev and a LAM grid for interpolation? 
yes; make them into a single object, e.g. by

elev$LAM = LAM[[1]]

which only make sense if elev and LAM have the same topology (grid
structure).
> Do they have to be at the same resolution?
>   
Yes.
> Thank you so much!
>
> Cara
>
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>   

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de/
http://www.springer.com/978-0-387-78170-9 e.pebesma at wwu.de