Description Usage Arguments Details Value Examples
Bivariate asymmetric Gumbel and Galambos copulas
1 2 3 4 5 6 7 8 9 | pasymgum(u,v,cpar) # C(u,v)
pcondasymgum21(v,u,cpar) # C_{2|1}(v|u)
pcondasymgum12(u,v,cpar) # C_{1|2}(u|v))
dasymgum(u,v,cpar) # c(u,v)
asymgum.cpar2tau(cpar) # Kendall's tau
asymgum.cpar2rhoS(cpar) # Spearman's rho
Basymgum(ww,cpar,mxiter=30,eps=1.e-7,iprint=F) # B(w)=A(w,1-w)
# similar to above with 'gum' replaced by 'gal' or 'gumMO' except
AasymgumMO(x,y,cpar)
|
u |
value in interval 0,1; could be a vector |
v |
value in interval 0,1; could be a vector |
cpar |
copula parameter vector, for asymgum, the two parameters are each in (0,1); for asymgal, the two parameters are each negative; for asymgumMO, the 3 parameters are delta>1, pi1 in (0,1] and pi2 in (0,1] |
ww |
vector of values in (0,1) |
mxiter |
maximum number of iterations |
eps |
tolerance for convergence |
iprint |
print flag for intermediate results |
x |
vector of positive values |
y |
vector of positive values |
The distribution in asymgum is also known as bilogistic in the multivariate extreme value literature – it is the bivariate Gumbel when the two parameters are equal and the Gumbel parameter is the reciprocal. The distribution in asymgal is also known as negative bilogistic in the multivariate extreme value literature – it is the bivariate Galambos when the two parameters are equal and the Galambos parameter is the negative reciprocal. The distribution in asymgumMO becomes bivariate Gumbel when the second and third parameters are both 1.
cdf for pasymgum
conditional cdf for pcondasymgum21 and pcondasymgum12
pdf for dasymgum
Kendall's tau for asymgum.cpar2tau
Spearman's rho for asymgum.cpar2rhoS
list with $Bfn (function values), $Bder (first derivatives), $Bder2 (second derivatives) for Basymgum
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | ze=.3; eta=.2
cpar=c(ze,eta)
heps=1.e-5
u=.3; v=.8
chkcopderiv(u,v,cpar,bcdf=pasymgum,pcond=pcondasymgum21,bpdf=dasymgum,
str="asymgum",eps=heps)
tau=asymgum.cpar2tau(cpar)
rho=asymgum.cpar2rhoS(cpar)
cat(cpar,tau,rho,"\n")
# special case of Gumbel
cpar=c(.5,.5)
cdf=pasymgum(.4,seq(.1,.9,.2),cpar)
cdf2=pgum(.4,seq(.1,.9,.2),1/cpar[1])
print(cbind(cdf,cdf2))
tau=asymgum.cpar2tau(cpar)
rho=asymgum.cpar2rhoS(cpar)
cat(tau,rho,"\n")
tau=gum.cpar2tau(1/cpar[1])
rho=gum.cpar2rhoS(1/cpar[1])
cat(tau,rho,"\n")
# asymmetric Galambos
cpar=c(-.5,-.5)
cdf=pasymgal(.4,seq(.1,.9,.2),cpar)
cdf2=pgal(.4,seq(.1,.9,.2),-1/cpar[1])
print(cbind(cdf,cdf2))
#and gumMO to add
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.