prob: prob

In flr/kobe: Tools for Scientific Fisheries Management Advice

prob

Description

Calculates the probability of an obervations occurring in a 2D cell using HPDregionplot Given a sample calculates the bivariate region of highest marginal posterior density for two variables, using kde2d from MASS to calculate a bivariate density.

Usage

probFn2(x, y, prob = c(0.5, 0.75, 0.95), n = 21,
  h = c(bandwidth.nrd(x), bandwidth.nrd(y)), lims = NULL,
  na.rm = FALSE)

Arguments

x	a vector
y	a vector
prob	probability levels
n	number of points at which to evaluate the density grid
h	bandwidth of 2D kernel smoother (previous default value was c(1,1), which worked poorly with some plots with very small scales; if not specified, defaults to values in kde2d)
lims	limits, specified as (x.lower,x.upper,y.lower,y.upper) (passed to kde2d)
na.rm	logical; if true, any NA and NaN's are removed from x before calculations
...	any other argument

Value

a data.frame with three variables x, y coordinates of the grid points, vectors of length n. level contours corresponding to prob

Examples

## Not run: 
   y=rnorm(20)
   x  =rnorm(20)
   prob(x,y)
## End(Not run)

flr/kobe documentation built on Sept. 10, 2023, 9:46 a.m.