plotpopn.bvn: Plot Population with Elliptical Home Ranges

Description Usage Arguments Details Value See Also Examples

Description

An secr popn object with attributes 's2xy' and 'theta' is plotted (over an existing plot if add = TRUE).

Usage

1
plotpopn.bvn(popn, radius = 2.45, add = FALSE, region = NULL, ...)

Arguments

popn

object of class 'popn'

radius

radius corresponding to chosen p values (see examples)

add

logical; if FALSE a new plot is started by calling plot(popn)

region

polygon(s) to which ellipses should be truncated (optional)

...

other arguments passed to polygon

Details

Coordinates of each ellipse are generated and plotted.

region may be an sp object that inherits from SpatialPolygons or a 2-column matrix of boundary coordinates.

Value

None. The function is used for its side effect of adding to an existing plot.

See Also

simpopn.bvn

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
tempgrid <- make.grid()
circular.pop <- simpopn.bvn(s2xy = c(625,625), core = tempgrid,
    buffer = 100, D = 1)
elliptical.pop <- simpopn.bvn(s2xy = c(625/3,625*3), core = tempgrid,
    buffer = 100, D = 1)
par(mfrow = c(1,2), par(xpd = TRUE))
plot(circular.pop)   ## this step determines the axes, frame etc.
plotpopn.bvn(circular.pop, add = TRUE)
plotpopn.bvn(elliptical.pop)

## demo aligned ellipses
par(mfrow=c(2,2))
for (i in 1:4) {
    aligned.pop <- simpopn.bvn(s2xy = c(625/3,625*3), theta=-1,
        core = tempgrid, buffer = 100, D = 1)
    plotpopn.bvn(aligned.pop, col='lightblue')
}

## demo variable size (up to 9-fold)
par(mfrow=c(1,1))
pop <- simpopn.bvn(s2xy = c(625,625), 
    core = tempgrid, buffer = 100, D = 1)
attr(pop, 's2xy') <- attr(pop, 's2xy') * (0.5 + runif(nrow(pop)*2))^2
plot(pop, cex=0)
plotpopn.bvn(pop, add = TRUE, col='grey')

MurrayEfford/secrBVN documentation built on May 26, 2019, 4:38 p.m.