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? Do they have to be at the same resolution?
Thank you so much!
Cara
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
