Learn R Programming

CircNNTSRaxial (version 0.1.0)

axialnntssimulation: NNTS axial density simulation function

Description

Simulation for the density function for the NNTS axial model

Usage

axialnntssimulation(nsim=1, cpars = 1/pi, M = 0)

Value

A list with two components:

simulations

The function generates a vector with nsim simulated random values from the NNTS axial density function

conteo

Number of uniform random numbers used to obtain the nsim simulations by the acceptance-rejection simulation method

Arguments

nsim

Number of simulations

cpars

Parameters of the model. A vector of complex numbers of dimension M+1. The sum of the squared moduli of the c parameters must be equal to 1/pi

M

Number of components in the NNTS axial model

Author

Juan Jose Fernandez-Duran and Maria Mercedes Gregorio-Dominguez

Examples

Run this code
set.seed(1234567)
ccoef<-axialnntsrandominitial(M=4)
data<-axialnntssimulation(100,ccoef,M=4)$simulations
hist(data,breaks=seq(0,pi,pi/7),xlab="Orientations (radians)",freq=FALSE,
ylab="",main="",ylim=c(0,1))
axialnntsplot(ccoef,M=4,add=TRUE)

Run the code above in your browser using DataLab