rancourtCR: Rancourt CR Method

In heike/rotations: Tools for working with rotational data

Description

Find the radius of a 100α% confidence region for the projected mean Rancourt et al. (2000)

Usage

  rancourtCR(Qs, a)

  ## S3 method for class 'Q4'
 rancourtCR(Qs, a)

  ## S3 method for class 'SO3'
 rancourtCR(Rs, a)

Arguments

Rs,Qs	A n-by-p matrix where each row corresponds to a random rotation in matrix (p=9) or quaternion form (p=4)
a	The alpha level desired, e.g. 0.05 or 0.10

Details

This works in the same way as done in Bingham et al. (2009) which assumes rotational symmetry and is therefore conservative.

Value

radius of the confidence region centered at the projected mean

References

Bingham M, Nordman D and Vardeman S (2009). "Modeling and Inference for Measured Crystal Orientations and a Tractable Class of Symmetric Distributions for Rotations in three Dimensions." _Journal of the American Statistical Association_, *104*(488), pp. 1385-1397.

Rancourt D, Rivest L and Asselin J (2000). "Using orientation statistics to investigate variations in human kinematics." _Journal of the Royal Statistical Society: Series C (Applied Statistics)_, *49*(1), pp. 81-94.

Examples

Qs<-ruars(20,rcayley,kappa=100,space='Q4')
region(Qs,method='prentice',alpha=0.9)

heike/rotations documentation built on May 17, 2019, 3:24 p.m.