The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.

EstPhiNestedGumbel

library(VecDep)

First, let us fix some parameters.

q = 4
dim = c(2,2)
n = 1000

Construct a four dimensional hierarchical Gumbel copula with parameters (theta_0,theta_1,theta_2) = (2,3,4)

hac = gethac(dim,c(2,3,4),type = 1)

Draw a sample of size n from the copula

sample = suppressWarnings(HAC::rHAC(n,hac))
#> [1] "X1 <-> 1"
#> [1] "X2 <-> 2"
#> [1] "X3 <-> 3"
#> [1] "X4 <-> 4"

Estimate the mutual information between two random vectors of size two via pseudo maximum likelihood estimation

estphi(sample,dim,list("hac",type = 1,M = 10000),function(t){t * log(t)})
#> Warning in hac2nacopula(hac): NAs introduced by coercion
#> [1] "X1 <-> 1"
#> [1] "X2 <-> 2"
#> [1] "X3 <-> 3"
#> [1] "X4 <-> 4"
#> [1] 0.4854667

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.