Running ‘RFsimulate.R’ [12s/12s]
  Running ‘krige.R’
  Running ‘lgcp.R’ [10s/10s]
  Running ‘lgm.R’ [17s/18s]
  Running ‘lgmRaster.R’ [199s/121s]
  Running ‘likfitLgm.R’ [12s/12s]
  Running ‘matern.R’
  Running ‘maternGmrfPrec.R’ [11s/11s]
  Running ‘profLlgm.R’ [15s/12s]
  Running ‘simLgcp.R’ [17s/23s]
  Running ‘stackRasterList.R’
Running the tests in ‘tests/lgcp.R’ failed.
Complete output:
  > options("rgdal_show_exportToProj4_warnings"="none") 
  > if(Sys.info()['sysname'] =='Linux' &
  +   requireNamespace("INLA", quietly=TRUE)) {   
  +   INLA::inla.setOption(inla.call = 
  +       system.file(paste(
  +           "bin/linux/",          
  +           ifelse(
  +             .Machine$sizeof.pointer == 4, 
  +             "32", "64"),
  +           'bit/inla.static', sep=''),
  +         package="INLA")) 
  + }
  > 
  > 
  > havePackages = c(
  +     'INLA' = requireNamespace('INLA', quietly=TRUE)
  + )
  > 
  > print(havePackages)
  INLA 
  TRUE 
  > 
  > # as in example
  > require('geostatsp')
  Loading required package: geostatsp
  Loading required package: Matrix
  Loading required package: raster
  Loading required package: sp
  > myPoints = SpatialPoints(cbind(rbeta(100,2,2), rbeta(100,3,4)))
  > myPoints@bbox = cbind(c(0,0), c(1,1))
  > 
  > mycov = raster(matrix(rbinom(100, 1, 0.5), 10, 10), 0, 1, 0, 1)
  > names(mycov)="x1"
  > 
  > if(all(havePackages)) {
  +   
  + res = lgcp(
  +     formula=~factor(x1),
  +     data=myPoints, 
  +     grid=20, 
  +     covariates=mycov,
  + 		priorCI=list(sd=c(u=0.1, alpha = 0.01), range=c(0.4, 0.41))
  + )
  + 
  + knitr::kable(res$parameters$summary, digits=3)
  + 
  + plot(res$raster[["predict.exp"]])
  + plot(myPoints,add=TRUE,col="#0000FF30",cex=0.5)
  + 
  + 	
  + # intercept only
  + 
  + res = lgcp(
  +     data=myPoints, 
  +     grid=20, 
  +     covariates=mycov,
  + 		formula=~1,
  + 		priorCI=list(sd=c(0.9, 1.1), range=c(0.4, 0.41))
  + )
  + 
  + knitr::kable(res$parameters$summary, digits=3)
  + 
  + plot(res$raster[["predict.exp"]])
  + plot(myPoints,add=TRUE,col="#0000FF30",cex=0.5)
  + 
  + if(interactive()  | Sys.info()['user'] =='patrick') {
  + 
  + # some missing values
  + 
  + temp = values(mycov)
  + temp[1:4] = NA
  + values(mycov) = temp
  + 
  + res = lgcp(data=myPoints, grid=20, covariates=mycov,
  + 		formula=~factor(x1),
  + 		priorCI=list(sd=c(0.9, 1.1), range=c(0.4, 0.41))
  + )
  + 
  + plot(res$raster[["predict.exp"]])
  + plot(myPoints,add=TRUE,col="#0000FF30",cex=0.5)
  + 
  + data('murder')
  + data('torontoPop')
  + myCov = list(
  +     pop=torontoPdens,
  +     inc = log(torontoIncome)
  + )
  + 
  + formula = ~ inc + offset(pop, log=TRUE)
  + 
  + resL=lgcp(formula, data=murder, 
  +     grid=squareRaster(murder, 30),
  +     covariates=myCov,
  +     border=torontoBorder)
  + 
  + resO = lgcp( ~ inc + pop, 
  +     data=murder, 
  +     grid=squareRaster(murder, 30),
  +     covariates=list(inc=myCov$inc, pop=log(myCov$pop)),
  +     border=torontoBorder)
  + 
  + rbind(resL$param$summary, resO$param$summary)
  + 
  + 
  + 
  + }	
  + }
  Error in plot.window(...) : need finite 'xlim' values
  Calls: plot -> plot -> plot.default -> localWindow -> plot.window
  In addition: Warning messages:
  1: In kable_pipe(x = character(0), caption = NULL, escape = TRUE) :
    The table should have a header (column names)
  2: In min(x) : no non-missing arguments to min; returning Inf
  3: In max(x) : no non-missing arguments to max; returning -Inf
  4: In min(x) : no non-missing arguments to min; returning Inf
  5: In max(x) : no non-missing arguments to max; returning -Inf
  Execution halted