pmvlog(q, dep, d = 2, mar = c(0, 1, 0))
rmvlog(n, dep, d = 2, mar = c(0, 1, 0))
d
or a matrix with d
columns, in which case the distribution is evaluated across
the rows.pmvlog
gives the distribution function and rmvlog
generates random deviates.d
dimensional multivariate logistic distribution
function (Gumbel, 1960) with parameter $\code{dep} = r$ is
$$G(z) = \exp\left[-(y_1^{1/r}+\ldots+y_d^{1/r})^r\right]$$
where $0 < r \leq 1$ and
$$y_i = {1+s_i(z_i-a_i)/b_i}^{-1/s_i}$$
for $1+s_i(z_i-a_i)/b_i > 0$ and
$i = 1,\ldots,d$.
Different parameters on each margin are not implemented, so
$\code{mar} = (a_i,b_i,s_i)$
for every $i$.
If $s_i = 0$ then $y_i$ is defined by continuity.
This is a special case of the multivariate asymmetric logistic
model.
The univariate marginal distributions are generalized extreme value.rbvlog
, rmvalog
,
rgev
pmvlog(matrix(rep(0:4,5), ncol=5), .7, d = 5)
pmvlog(rep(4,5), .7, d = 5)
rmvlog(10, .7, d = 5)
Run the code above in your browser using DataLab