## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  eval = FALSE,
  fig.width=8, 
  fig.height=5,
  warning = FALSE,
  message = FALSE)

## ----Load packages------------------------------------------------------------
#  
#  library(intSDM)
#  library(INLA)
#  

## ----startWorkflow------------------------------------------------------------
#  
#  Rich <- startWorkflow(Species = c("Lolium perenne L.","Rubus caesius L.",
#                                    "Rosa spinosissima L.",
#                                    "Poa trivialis L.",
#                                    "Galium verum L.",
#                                    "Tanacetum vulgare L.",
#                                    "Viola tricolor L.",
#                                    "Epilobium L."),
#          Projection =  '+proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=km +no_defs',
#                        Save = FALSE, Richness = TRUE,
#                        saveOptions = list(projectName = 'Richness'))
#  

## ----addArea------------------------------------------------------------------
#  
#  Ned <- giscoR::gisco_get_countries(country = 'Netherlands', resolution = 60)
#  Ned <- st_cast(st_as_sf(Ned), 'POLYGON')
#  Ned <- Ned[which.max(st_area(Ned)),]
#  Ned <- rmapshaper::ms_simplify(Ned, keep = 0.5)
#  Rich$addArea(Ned)
#  

## ----addGBIF------------------------------------------------------------------
#  
#  Rich$addGBIF(datasetName = 'DVD',
#               datasetKey = '740df67d-5663-41a2-9d12-33ec33876c47',
#               datasetType = 'PA', generateAbsences = TRUE)
#  
#  Rich$addGBIF(datasetName = 'iNat',
#               datasetKey = '50c9509d-22c7-4a22-a47d-8c48425ef4a7')
#  

## ----Mesh---------------------------------------------------------------------
#  
#  Rich$addMesh(max.edge = c(5, 10))
#  Rich$plot(Mesh = TRUE)
#  

## ----priors-------------------------------------------------------------------
#  
#  Rich$specifySpatial(prior.range = c(0.2,0.1),
#                      prior.sigma = c(2, 0.1))
#  
#  Rich$biasFields('iNat', prior.range = c(0.1, 0.1),
#                  prior.sigma = c(0.2, 0.1))
#  
#  Rich$specifyPriors(priorIntercept = list(prior = 'pc.prec', param = c(0.02, 0.01)))
#  

## ----modelFormula-------------------------------------------------------------
#  
#  Rich$addCovariates(worldClim = c('tavg'), res = 10, Function = scale)
#  Rich$modelFormula(covariateFormula = ~ tavg + I(tavg^2))
#  

## ----specRich-----------------------------------------------------------------
#  
#  Rich$modelOptions(ISDM = list(Offset = 'sampleSizeValue'),
#                    Richness = list(predictionIntercept = 'DVD'))
#  

## ----workflow-----------------------------------------------------------------
#  
#  Rich$workflowOutput('Maps')
#  
#  RichModel <- sdmWorkflow(Rich, inlaOptions = list(verbose = TRUE))
#  

## ----Rich---------------------------------------------------------------------
#  
#  ggplot() + gg(RichModel$Richness$Richness, aes(col = q0.025))
#  ggplot() + gg(RichModel$Richness$Richness, aes(col = q0.5))
#  ggplot() + gg(RichModel$Richness$Richness, aes(col = q0.975))
#  

## ----prob---------------------------------------------------------------------
#  
#  ggplot() + gg(RichModel$Richness$Probabilities$Galium_verum_L., aes(col = mean))
#  ggplot() + gg(RichModel$Richness$Probabilities$Tanacetum_vulgare_L., aes(col = mean))
#  
#