Skip to contents

Estimates concentration parameter, angular standard deviation, and confidence limit.

Usage

fisher_statistics(x, w = NULL, conf.level = 0.95, na.rm = TRUE)

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.

conf.level

numeric. Level of confidence.

na.rm

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

Value

list, with

"k"

estimated concentration parameter \(\kappa\) for the von Mises-Fisher distribution

"csd"

estimated angular standard deviation enclosing 63% of the orientation data. Angle is in degrees if x is a spherical object, and radians if otherwise.

"alpha"

Confidence limit for given conf.level. Angle is in degrees if x is a spherical object, and radians if otherwise.

Examples

set.seed(20250411)
x <- rvmf(100, mu = Line(120, 50), k = 5)
fisher_statistics(x)
#> $k
#> [1] 4.416173
#> 
#> $csd
#> [1] 38.54446
#> 
#> $csd_2s
#> [1] 66.62006
#> 
#> $alpha
#> [1] 7.640117
#>