Errors in running code in vignettes:
when running code in ‘gstat.Rnw’
  ...

> library(sp)

> data(meuse)

> class(meuse)
[1] "data.frame"

> names(meuse)
 [1] "x"       "y"       "cadmium" "copper"  "lead"    "zinc"    "elev"   
 [8] "dist"    "om"      "ffreq"   "soil"    "lime"    "landuse" "dist.m" 

> coordinates(meuse) = ~x + y

> class(meuse)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"

> summary(meuse)
Object of class SpatialPointsDataFrame
Coordinates:
     min    max
x 178605 181390
y 329714 333611
Is projected: NA 
proj4string : [NA]
Number of points: 155
Data attributes:
    cadmium           copper            lead            zinc       
 Min.   : 0.200   Min.   : 14.00   Min.   : 37.0   Min.   : 113.0  
 1st Qu.: 0.800   1st Qu.: 23.00   1st Qu.: 72.5   1st Qu.: 198.0  
 Median : 2.100   Median : 31.00   Median :123.0   Median : 326.0  
 Mean   : 3.246   Mean   : 40.32   Mean   :153.4   Mean   : 469.7  
 3rd Qu.: 3.850   3rd Qu.: 49.50   3rd Qu.:207.0   3rd Qu.: 674.5  
 Max.   :18.100   Max.   :128.00   Max.   :654.0   Max.   :1839.0  
                                                                   
      elev             dist               om         ffreq  soil   lime   
 Min.   : 5.180   Min.   :0.00000   Min.   : 1.000   1:84   1:97   0:111  
 1st Qu.: 7.546   1st Qu.:0.07569   1st Qu.: 5.300   2:48   2:46   1: 44  
 Median : 8.180   Median :0.21184   Median : 6.900   3:23   3:12          
 Mean   : 8.165   Mean   :0.24002   Mean   : 7.478                        
 3rd Qu.: 8.955   3rd Qu.:0.36407   3rd Qu.: 9.000                        
 Max.   :10.520   Max.   :0.88039   Max.   :17.000                        
                                    NA's   :2                             
    landuse       dist.m      
 W      :50   Min.   :  10.0  
 Ah     :39   1st Qu.:  80.0  
 Am     :22   Median : 270.0  
 Fw     :10   Mean   : 290.3  
 Ab     : 8   3rd Qu.: 450.0  
 (Other):25   Max.   :1000.0  
 NA's   : 1                   

> coordinates(meuse)[1:5, ]
       x      y
1 181072 333611
2 181025 333558
3 181165 333537
4 181298 333484
5 181307 333330

> bubble(meuse, "zinc", col = c("#00ff0088", "#00ff0088"), 
+     main = "zinc concentrations (ppm)")

> print(bubble(meuse, "zinc", col = c("#00ff0088", "#00ff0088"), 
+     main = "zinc concentrations (ppm)"))

> data(meuse.grid)

> summary(meuse.grid)
       x                y              part.a           part.b      
 Min.   :178460   Min.   :329620   Min.   :0.0000   Min.   :0.0000  
 1st Qu.:179420   1st Qu.:330460   1st Qu.:0.0000   1st Qu.:0.0000  
 Median :179980   Median :331220   Median :0.0000   Median :1.0000  
 Mean   :179985   Mean   :331348   Mean   :0.3986   Mean   :0.6014  
 3rd Qu.:180580   3rd Qu.:332140   3rd Qu.:1.0000   3rd Qu.:1.0000  
 Max.   :181540   Max.   :333740   Max.   :1.0000   Max.   :1.0000  
      dist        soil     ffreq   
 Min.   :0.0000   1:1665   1: 779  
 1st Qu.:0.1193   2:1084   2:1335  
 Median :0.2715   3: 354   3: 989  
 Mean   :0.2971                    
 3rd Qu.:0.4402                    
 Max.   :0.9926                    

> class(meuse.grid)
[1] "data.frame"

> coordinates(meuse.grid) = ~x + y

> class(meuse.grid)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"

> gridded(meuse.grid) = TRUE

> class(meuse.grid)
[1] "SpatialPixelsDataFrame"
attr(,"package")
[1] "sp"

> image(meuse.grid["dist"])

> title("distance to river (red = 0)")

> library(gstat)

> zinc.idw = idw(zinc ~ 1, meuse, meuse.grid)
[inverse distance weighted interpolation]

> class(zinc.idw)
[1] "SpatialPixelsDataFrame"
attr(,"package")
[1] "sp"

> spplot(zinc.idw["var1.pred"], main = "zinc inverse distance weighted interpolations")

> print(spplot(zinc.idw["var1.pred"], main = "zinc inverse distance weighted interpolations"))

> plot(log(zinc) ~ sqrt(dist), meuse)

> abline(lm(log(zinc) ~ sqrt(dist), meuse))

> lzn.vgm = variogram(log(zinc) ~ 1, meuse)

> lzn.vgm
    np       dist     gamma dir.hor dir.ver   id
1   57   79.29244 0.1234479       0       0 var1
2  299  163.97367 0.2162185       0       0 var1
3  419  267.36483 0.3027859       0       0 var1
4  457  372.73542 0.4121448       0       0 var1
5  547  478.47670 0.4634128       0       0 var1
6  533  585.34058 0.5646933       0       0 var1
7  574  693.14526 0.5689683       0       0 var1
8  564  796.18365 0.6186769       0       0 var1
9  589  903.14650 0.6471479       0       0 var1
10 543 1011.29177 0.6915705       0       0 var1
11 500 1117.86235 0.7033984       0       0 var1
12 477 1221.32810 0.6038770       0       0 var1
13 452 1329.16407 0.6517158       0       0 var1
14 457 1437.25620 0.5665318       0       0 var1
15 415 1543.20248 0.5748227       0       0 var1

> lzn.fit = fit.variogram(lzn.vgm, model = vgm(1, "Sph", 
+     900, 1))

> lzn.fit
  model      psill    range
1   Nug 0.05066243   0.0000
2   Sph 0.59060780 897.0209

> plot(lzn.vgm, lzn.fit)

> print(plot(lzn.vgm, lzn.fit))

> lznr.vgm = variogram(log(zinc) ~ sqrt(dist), meuse)

> lznr.fit = fit.variogram(lznr.vgm, model = vgm(1, 
+     "Exp", 300, 1))

> lznr.fit
  model      psill    range
1   Nug 0.05712231   0.0000
2   Exp 0.17641559 340.3201

> plot(lznr.vgm, lznr.fit)

> print(plot(lznr.vgm, lznr.fit))

> lzn.kriged = krige(log(zinc) ~ 1, meuse, meuse.grid, 
+     model = lzn.fit)
[using ordinary kriging]

> spplot(lzn.kriged["var1.pred"])

> print(spplot(lzn.kriged["var1.pred"]))

> lzn.condsim = krige(log(zinc) ~ 1, meuse, meuse.grid, 
+     model = lzn.fit, nmax = 30, nsim = 4)
drawing 4 GLS realisations of beta...
[using conditional Gaussian simulation]

> spplot(lzn.condsim, main = "four conditional simulations")

> print(spplot(lzn.condsim, main = "four conditional simulations"))

> lzn.condsim2 = krige(log(zinc) ~ sqrt(dist), meuse, 
+     meuse.grid, model = lznr.fit, nmax = 30, nsim = 4)
drawing 4 GLS realisations of beta...
[using conditional Gaussian simulation]

> spplot(lzn.condsim2, main = "four UK conditional simulations")

> print(spplot(lzn.condsim2, main = "four UK conditional simulations"))

> lzn.dir = variogram(log(zinc) ~ 1, meuse, alpha = c(0, 
+     45, 90, 135))

  When sourcing ‘gstat.R’:
Error: replacement has 30 rows, data has 15
Execution halted