ca
for fixmahal dependent on
the number of points and variables of the current fixed point cluster
(FPC).
This is experimental and only thought for use in fixmahal.cmahal(n, p, nmin, cmin, nc1, c1 = cmin, q = 1)ca is computed. For smaller FPC sizes, ca is set to
the value for nmin.ca.ca=c1.cmahal.
Value for ca for FPC size equal to nc1.ca as function of the FPC size. Should presumably always be 1.n, giving the values for ca
for all FPC sizes smaller or equal to n.ca should
decrease with increasing FPC size and increase with increasing
p in fixmahal. This is to prevent too small
meaningless FPCs while maintaining the significant larger
ones. cmahal with q=1 computes ca in such a way
that as long as ca>cmin, the decrease in n is as steep
as possible in order to maintain the validity of the convergence
theorem in Hennig and Christlieb (2002).fixmahalplot(1:100,cmahal(100,3,nmin=5,cmin=qchisq(0.99,3),nc1=90),
xlab="FPC size", ylab="cmahal")Run the code above in your browser using DataLab