First, let us fix some parameters.
Construct a four dimensional hierarchical Gumbel copula with parameters (theta_0,theta_1,theta_2) = (2,3,4)
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
Need a high-speed mirror for your open-source project?
Contact our mirror admin team at info@clientvps.com.
This archive is provided as a free public service to the community.
Proudly supported by infrastructure from VPSPulse , RxServers , BuyNumber , UnitVPS , OffshoreName and secure payment technology by ArionPay.