Learn R Programming

Directional (version 7.0)

Density of the Wood bimodal distribution on the sphere: Density of the Wood bimodal distribution on the sphere

Description

Density of the Wood bimodal distribution on the sphere.

Usage

dwood(y, param, logden = FALSE)

Value

A vector with the (log) density values of y.

Arguments

y

A matrix containing two columns. The first one is the latitude and the second is the longitude, both expressed in degrees.

param

A vector with the 5 parameters, in the order they are returned by the wood.mle function. That is, \((\gamma, \delta, \alpha, \beta, \kappa)\).

logden

If you the logarithm of the density values set this to TRUE.

Author

Michail Tsagris.

R implementation and documentation: Michail Tsagris mtsagris@uoc.gr.

Details

The density of the spherical Wood distribution is computed.

References

Wood A.T.A. (1982). A bimodal distribution on the sphere. Journal of the Royal Statistical Society, Series C, 31(1): 52--58.

See Also

dkent, desag, wood.mle

Examples

Run this code
x <- rvmf(100, rnorm(3), 15)
x <- euclid.inv(x)
mod <- wood.mle(x)
d <- dwood(x, mod$info[, 1])

Run the code above in your browser using DataLab