dnestlog {BMAmevt}  R Documentation 
Pairwise Beta (PB) and Nested Asymmetric Logistic (NL) distributions
Description
Likelihood function (spectral density) and random generator in the Pairwise Beta and NL models.
Usage
dnestlog(
x = rbind(c(0.1, 0.3, 0.6), c(0.3, 0.3, 0.4)),
par = c(0.5, 0.5, 0.2, 0.3),
log = FALSE,
vectorial = TRUE
)
dpairbeta(
x,
par = c(1, rep(2, choose(4, 2) + 1)),
log = FALSE,
vectorial = TRUE
)
rnestlog(
n = 5,
par = c(0.2, 0.3, 0.4, 0.5),
threshold = 1000,
return.points = FALSE
)
rpairbeta(n = 1, dimData = 3, par = c(1, rep(1, 3)))
Arguments
x 
An angular data set (may be reduced to a single point).
A 
par 
The parameter for the Pairwise Beta or the Nested Logistic density.

log 
Logical. Should the density be returned on the log scale ? 
vectorial 
Logical. Should a vector or a single value be returned ? 
n 
The number of points on the simplex to be generated. 
threshold 
The radial threshold
to hold. 
return.points 
logical: should the censored vectorial dataset corresponding to the angular one be returned ? 
dimData 
the dimension of the sample space, which is 
Details
Applies to angular data sets. The density is given with respect to the Lebesgue measure on R^{p1}
, where p
is the number of columns in x
(or the length of x
, if the latter is a single point).
Value
The value returned by the likelihood function is imposed (see
e.g. posteriorMCMC
.
In contrast, the random variable have unconstrained output format.

dpairbeta
returns the likelihood as a single number ifvectorial ==FALSE
, or as a vector of sizenrow(x)
containing the likelihood of each angular data point. Iflog == TRUE
, the loglikelihood is returned instead.rpairbeta
returns a matrix withn
rows anddimData
columns. 
dnestlog
returns the likelihood as a single number ifvectorial ==FALSE
, or as a vector of sizenrow(x)
containing the likelihood of each angular data point. Iflog == TRUE
, the loglikelihood is returned instead.rnestlog
returns a matrix withn
rows anddimData
columns ifreturn.points==FALSE
(the default). Otherwise, a list is returned, with two elements:
Angles
: The angular data set 
Points
: The full trivariate data set abovethreshold
(i.e.Angles
multiplied by the radial components)
