[R-sig-Geo] Huge kriging weights when using kriging with external drift in gstat
Tobias Vetter
vetter at pik-potsdam.de
Tue Jul 2 16:16:31 CEST 2013
Thank you.
So as long as all my distances are smaller than the range I will not
have this problem?
Tobias
> Tobias, try using another variogram model: the linear model with sill is
> not positive definite in more than one dimension (and yes: a warning
> about this would have been very nice!)
>
> On 07/02/2013 03:37 PM, Tobias Vetter wrote:
>> Dear all,
>> I am interpolating automatically daily air temperature values using
>> kriging with external drift (altitude) and the library gstat. From time
>> to time (quite random) I get unrealistic predictions like -400°C with
>> very high kriging weights like 2000.
>> Looking at the data everything seems to be ok. Then I changed the
>> variogram ranges a little bit and the problem disappears.
>> How can this behavior be explained?
>>
>> I have attached a small example script for demonstration and some data.
>>
>> Thanks you
>> Tobias
>>
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
More information about the R-sig-Geo
mailing list