[R-sig-Geo] error from universal kriging with krige()

[Edzer Pebesma]

Is it a recent R/gstat version? If yes, could you provide me with a 
reproducable example?
--
Edzer

Sven Burbeck wrote:
> Hello,
>
> when I try to apply universal kriging by krige() I produce the following
> error:
>   
>> myVg_2 <- variogram(log10(M_BL) ~  GEO_UNIT - 1, data = myBlt_utm)
>> myVgm <- vgm(psill = 0.05,model = "Sph",range = 25000,nugget = 0.1)
>> myVgmf <- fit.variogram(object = myVg_2,model = myVgm,fit.sills = TRUE,
>>     
> fit.ranges = TRUE,fit.method = 7)
>   
>> myBlt_grid_UK <- krige(log10(M_BL) ~  GEO_UNIT - 1, locations =
>>     
> myBlt_utm, newdata = myPointGrid_utm, model = myVgmf)
> Error in gstat.formula.predict(d$formula, newdata, na.action = na.action,  :
>   NROW(locs) != NROW(X): this should not occur
> In addition: Warning messages:
> 1: 'newdata' had 2701 rows but variable(s) found have 911 rows
> 2: 'newdata' had 2701 rows but variable(s) found have 911 rows
>
>
> I know there was a posting about this error. There it was a coordinate
> problem but still I don't know how to solve the problem.
>
> Can anybody give me a hint?
>
> Thanks in advance,
> Sven
>
> _______________________________________________
> 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/