Skip to contents

Elliptical concentration and confidence cone estimation

Source: R/vector_stats.R
bingham_statistics.Rd

Elliptical concentration and confidence cone estimation

Usage

bingham_statistics(x, w = NULL, na.rm = TRUE)

Source

Borradaile, G. (2003). Spherical-Orientation Data. In: Statistics of Earth Science Data. Springer, Berlin, Heidelberg. https://doi.org/10.1007/978-3-662-05223-5_10

Arguments

x

object of class "Vec3", "Line", or "Plane".

w

numeric. Optional weights for each observation.

na.rm

logical. Whether NA values should be removed before the computation proceeds.

Value

list

k

two-column vector containing the estimates for the minimum (\(k_\text{min}\)) and maximum concentration (\(k_\text{max}\)).

a95

two-column vector containing the estimates for the minimum and maximum 95% confidence cone.

beta

The shape factor of the distribution given by the ratio \(\frac{k_\text{min}}{k_\text{max}}\).

See also

inertia_tensor.spherical()

Examples

set.seed(1234)
x <- rfb(100, mu = Line(120, 50), k = 15, A = diag(c(-5, 0, 5)))

stereoplot()
stereo_point(x)


bingham_statistics(x)
#> $k
#> [1] 2.591337 4.034218
#> 
#> $a95
#> [1] 8.696932 6.970250
#> 
#> $beta
#> [1] 0.6423395
#>