genden.data: Create density surface from animal location data

Description Usage Arguments Details Value Note Author(s) See Also Examples

Description

If user possesses location data (from aerial photography, telemetry, etc.) of animal distribution within the study area, this information can be used to generate realizations of this population for subsequent simulation.

Usage

1
genden.data(reg, xylocations, nint.x = 100, nint.y = 100, gam.param = 1)

Arguments

reg

a region onto which the density is mapped

xylocations

Data frame consisting of only x and y columns. These are the animal locations. Note that range of x- and y-values must fall within the boundaries specified by the argument reg

nint.x

number of x-intervals in the density matrix.

nint.y

number of y-intervals in the density matrix.

gam.param

Gamma parameter in gam smooth; lower this value <1 to cause a more 'patchy' fit to your data; a value >1.5 will cause a higher degree of smoothness to the surface

Details

For users wishing to be able to more realistically simulate the population of interest to them, this function has been created to allow density surfaces that more closely approximate the dispersion of animals within a study area. If location data are available, a data frame consisting of the x- and y-coordinates can be used as data for a simple fitting exercise.

This function is intended as a substitute for add.hotspot and set.stripe when data about animal locations are available.

Value

An object of class 'density.population´ which has the following elements:

n.interval.x

number of x-intervals

n.interval.y

number of y-intervals

matrix

the density matrix.

Note

The method by which the density surface is a simple smooth of the location data, using only the x- and y- coordinates as predictors. This is done using gams (general additive models) from the mgcv package. Consult that package for more details regarding fitting smooths.

Author(s)

Eric Rexstad, RUWPA, Univ. of St. Andrews.

See Also

add.hotspot, set.stripe, run.FShake3D

Examples

1
2
3
4
5
6
	library(wisp)
	data(beach); data(seal.pop)
#			Borrow animal locations from the 'seal.pop' dataset, and fit a gam to those locations	
	gam.seal <- genden.data(beach, xylocations=data.frame(x=seal.pop$posx, y=seal.pop$posy), 
						  gam.param=0.7)
	plot.density.population(gam.seal, method='image')

DistanceDevelopment/WiSP documentation built on Sept. 18, 2020, 2:55 p.m.