Dears,
I am writing sice having question about your R package "copula".
I am using this package to fit my data with some of the 22 copulas mentioned
in Nelsen 1999.
However, when trying random number generator I have got some difficulties.
For example when I have familie number A13 where cupola is equal to
C(u,v) = exp(1-((1-log(u))^theta+(1-log(v))^theta-1)^(1/theta))
devuv = (exp(1 - (-1 + (1 - log(u))^theta + (1 -
log(v))^theta)^theta^(-1))*(1 - log(u))^(-1 + theta)*(-1 + theta + (-1 + (1
- log(u))^theta + (1 - log(v))^theta)^theta^(-1))*(-1 + (1 - log(u))^theta +
(1 - log(v))^theta)^(-2 + theta^(-1))*(1 - log(v))^(-1 + theta))/(u*v)
devu = (exp(1 - (-1 + (1 - log(u))^theta + (1 -
log(v))^theta)^theta^(-1))*(1 - log(u))^(-1 + theta)*(-1 + (1 -
log(u))^theta + (1 - log(v))^theta)^(-1 + theta^(-1)))/u
phi = (1-log(t))^theta-1
devphi = - ((theta*(1 - log(t))^(-1 + theta))/t)
K(t) = t - ((1-log(t))^theta-1) / (-((theta*(1 - log(t))^(-1 + theta))/t))
I am always meeting some difficulies to random generate data. Or to make a
program code to count inverse function of my phi which should help in the
algorithm to generate data.
It would be grateful if you could help me.
Thanks a lot.
--
Mgr. Rastislav Matúš
PhD. Student
Department of Land and Water Resources Management
Faculty of Civil Engineering
Slovak University of Technology
Blok C, 12. posch.,
Radlinského 11,
813 68 Bratislava 15,
EUROPE
BUREAU: Tel.: +421 2 59274 622, fax: +421 2 52923575
PERSONAL: +421 908 623450
E-mail: matus@svf.stuba.sk
www.kvhk.sk
[[alternative HTML version deleted]]
