IsoriX (version 0.5)

calibdata: Calibration dataset

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.

Usage

data(calibdata)

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.

See Also

calibfit to fit a calibration model

Examples

# NOT RUN {
data(calibdata)
head(calibdata)
str(calibdata)

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

# }
# NOT RUN {
if(require(raster)){
    ### DELETE AND RECREATE ELEVATION DATA
    calibdata$elev <- NULL  ## we delete them

    ## 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")
    calibdata$elev <- extract(
        elevationrasterbig,
        cbind(calibdata$long, calibdata$lat))
    head(calibdata)
}

# }
# NOT RUN {
# }