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", "Ray", or "Plane", where the rows are the observations and the columns are the coordinates.

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.667941 3.969685
#> 
#> $a95
#> [1] 8.571166 7.026677
#> 
#> $beta
#> [1] 0.6720788
#>