R/axialnntsrandominitial.R

In CircNNTSRaxial: Axial Data using NNTS Models

axialnntsrandominitial<-function (M = 1)
{
  if (M < 0)
    return("M must be nonnegative")
  res <- complex(M + 1)
  aux <- rnorm(2 * (M + 1))
  aux <- sqrt(1/(pi)) * (aux/sqrt(sum(aux^2)))
  res[1] <- Mod(aux[1] + (0+1i) * aux[M + 2])
  if (M > 0) {
    for (k in 2:(M + 1)) {
      res[k] <- aux[k] + (0+1i) * aux[k + M + 1]
    }
  }
  else res <- 1/(pi)
  return(res)
}