Computes the maximum likelihood estimates of the symmetric NNTS parameters
Usage
nntsestimationSymmetric(M = 0, data, maxit = 500)
Arguments
M
Number of components in the NNTS
data
Vector of angles in radians
maxit
Maximum number of iterations in the optimization algorithm
Value
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
Note
For the maximization of the loglikelihood function the function constrOptim from the package stats is used
Author(s)
Juan Jose Fernandez-Duran y Maria Mercedes Gregorio-Dominguez
References
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
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(CircNNTSR)
Attaching package: 'CircNNTSR'
The following object is masked from 'package:grDevices':
trans3d
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/CircNNTSR/nntsestimationSymmetric.Rd_%03d_medium.png", width=480, height=480)
> ### Name: nntsestimationSymmetric
> ### Title: NNTS Symmetric Coefficient estimation
> ### Aliases: nntsestimationSymmetric
>
> ### ** Examples
>
> 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)
>
>
>
>
>
> dev.off()
null device
1
>