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