R/axialnntsdensity.R

In CircNNTSRaxial: Axial Data using NNTS Models

axialnntsdensity<-function (data, cpars = 1/sqrt(pi), M = 0)
{
  data <- as.matrix(data)
  n <- nrow(data)
  R <- 1
  if (R != length(M))
    return("Error: M must correspond to a univariate variable")
  if (length(cpars) != (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 norm of componentes does not satisfy condition")
  if (sum(M) == 0)
    return(t(rep(1/(pi), n)))
  gridcomb <- 0:M
  statisticsmatrix <- matrix(0, nrow = M + 1, ncol = n)
  statisticsmatrix <- exp((0+1i) * gridcomb %*% t(2*data))
  aux <- t(cpars) %*% statisticsmatrix
  res <- aux * Conj(aux)
  return(t(Re(res)))
}