R/axialnntsloglik.R

In CircNNTSRaxial: Axial Data using NNTS Models

axialnntsloglik<-function (data, cpars = 1/sqrt(pi), M = 0)
{
  if (M == 0)
    return(-length(data) * log(pi))
  size <- length(cpars)
  if (size != M + 1)
    return("Length of cpars must be equal to M+1")
  if (abs(sum(Mod(cpars)^2) - 1/(pi)) > 1e-10)
    return("sum of the squared norms of componentes greater than condition")
  y <- sum(log(axialnntsdensity(data, cpars, M)))
  return(y)
}