CalibDataBat: Calibration dataset for bat species

Description

This dataset contains deuterium delta values of fur keratin from sedentary bat species captured between 2005 and 2009 from Popa-Lisseanu et al. (2012). The data can be used as an example to fit a calibration model using the function calibfit.

Arguments

Format

A dataframe with 178 observations on 7 variables:

[, 1]	siteID	(`Factor`)	Identification of the sampling site
[, 2]	long	(`numeric`)	Longitude coordinate [decimal degrees]
[, 3]	lat	(`numeric`)	Latitude coordinate [decimal degrees]
[, 4]	elev	(`numeric`)	Elevation asl [m]
[, 5]	animalID	(`Factor`)	Identification of the sampled animal
[, 6]	species	(`Factor`)	Species name
[, 7]	tissue.value	(`numeric`)	Deuterium delta value of the tissue

Details

Users who wish to use their own dataset for calibration should create a dataframe of similar structure than this one (only the column 'species' can be dropped). The columns should possess the same names as the ones described above. If the elevation is unknown at the sampling sites, elevation information can be extracted from a high resolution elevation raster using the function extract. In this dataset, we retrieved elevations from the Global Multi-resolution Terrain Elevation Data 2010. Note that the original study used a different source of elevation data.

References

Popa-Lisseanu, A. G., Soergel, K., Luckner, A., Wassenaar, L. I., Ibanez, C., Kramer-Schadt, S., Ciechanowski, M., Goerfoel, T., Niermann, I., Beuneux, G., Myslajek, R. W., Juste, J., Fonderflick, J., Kelm, D., Voigt, C. C. (2012). A triple isotope approach to predict the breeding origins of European bats. PLoS ONE 7(1):e30388.

Examples

Run this code

# NOT RUN {
head(CalibDataBat)
str(CalibDataBat)

## The following example require to have downloaded
## a large elevation raster with the function getelev()
## and will therefore not run unless you type:
## example(CalibDataBat, run.dontrun=TRUE)

# }
# NOT RUN {
if(require(raster)){
    ## We delete the elevation data
    CalibDataBat$elev <- NULL

    ## We reconstruct the elevation data using an elevation raster
    ## (see ?getelev for details on how to get the tif file)
    elevationrasterbig <- raster("gmted2010_30mn.tif")
    CalibDataBat$elev <- extract(
        elevationrasterbig,
        cbind(CalibDataBat$long, CalibDataBat$lat))
    head(CalibDataBat)
}

# }
# NOT RUN {

# }