tests/testthat/nestedFiles/calibfit.R

# The 'isoscape' example without the additional calibration fits
if(IsoriX::getOption_IsoriX("example_maxtime") > 10) {
  GNIPDataDEagg <- prepsources(data = GNIPDataDE)
  
  ## We fit the models for Germany:
  GermanFit <- isofit(data = GNIPDataDEagg,
                      mean_model_fix = list(elev = TRUE, lat_abs = TRUE))
  
  ## We fit the calibration model:
  CalibAlien <- calibfit(data = CalibDataAlien, isofit = GermanFit)
  
  ## We display minimal information:
  CalibAlien
  
  ## We display more information:
  summary(CalibAlien)
  
  ## We plot the calibration function:
  plot(CalibAlien)
}

Try the spaMM package in your browser

Any scripts or data that you put into this service are public.

spaMM documentation built on Aug. 30, 2023, 1:07 a.m.