nntsestimationSymmetric | R Documentation |
Computes the maximum likelihood estimates of the symmetric NNTS parameters
nntsestimationSymmetric(M = 0, data, maxit = 500)
M |
Number of components in the NNTS |
data |
Vector of angles in radians |
maxit |
Maximum number of iterations in the optimization algorithm |
coef |
Vector of length M+1. The first M components are the squared moduli of the c parameters, and the last number is the mean of symmetry |
loglik |
Optimum log-likelihood value |
AIC |
Value of Akaike's Information Criterion |
BIC |
Value of Bayesian Information Criterion |
convergence |
An integer code: zero indicates successful convergence; error codes are the following: one indicates that the iteration limit maxit has been reached, and 10 indicates degeneracy of the Nelder-Mead simplex |
For the maximization of the loglikelihood function the function constrOptim from the package stats is used
Juan Jose Fernandez-Duran y Maria Mercedes Gregorio-Dominguez
Fernandez-Duran, J.J., Gregorio-Dominguez, M.M. (2009) Symmetric Circular Distributions Based on Nonnegative Trigonometric Sums. Working Paper, DE-C09.12, Department of Statistics, ITAM, Mexico
b<-c(runif(10,3*pi/2,2*pi-0.00000001),runif(10,pi/2,pi-0.00000001))
estS<-nntsestimationSymmetric(2,b)
nntsplotSymmetric(estS$coef,2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.