[R-sig-Geo] fitting variogram in regrssion kriging

[Tomislav Hengl]

Helen,

Most likely you need to subset the missing values by using e.g. 
"precip[!is.na(precip$annaul_2001),]", but it is difficult to see if you 
maybe have a typo with variable names.

Try using "traceback()" and check the amount of missing values by using 
e.g. "summary(is.na(precip$annaul_2001))" and/or 
"summary(complete.cases(precip))".

HTH,

T. (Tom) Hengl
Researcher @ ISRIC - World Soil Information
Url: http://www.wageningenur.nl/en/Persons/dr.-T-Tom-Hengl.htm
Network: http://profiles.google.com/tom.hengl
Publications: http://scholar.google.com/citations?user=2oYU7S8AAAAJ


On 25-9-2014 7:57, Helen Chen wrote:
> Hi,
>
> I am following the code provide from this link:
> http://spatial-analyst.net/wiki/index.php?title=Regression-kriging_guide)
> with my own data.
>
> but I got stuck in fitting variogram, the error message is as follow:
>
>> vgm_precip_r <- fit.variogram(variogram(annaul_2001~mask_10m.asc,
> precip), model=null.vgm)
> Error in na.fail.default(list(annaul_2001 = c(945.134, 648.716, 559.562,  :
>    missing values in object
>
> But I have checked that there is no missing value in my rainfall data. Any
> help would be very much appreciated.
>
> Helen
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>