Dear all,
I´m geneticist. Hence, sorry for the dummy questions.
I have to generate a new variable using a set of parameters.
The new variable is generated as follows:
((z^y)/(gamma(y,z)))* integral of x^y-1 * exp^-zx dx from c to infinity,
where
z, y, c are values previously calculated
gamma is the Gamma distribution with parameters y and z.
How do I implement such formula using Stata?
I have already checked "functions" commands, but still feel unsure
regarding the correct implementation.
I´ll be glad for any help. Thanks in advance.
Tiago
--
Tiago V. Pereira
Heart Institute - InCor
Laboratory of Cardiovascular Genetics
Federal University of São Paulo
University of São Paulo Medical School
Av. Dr. Eneas de Carvalho Aguiar, 44;
Cerqueira César - CEP 05403-000
Sao Paulo, SP Brazil.
Tel./fax: +55 11 3069 5068
email: tiago.pereira@incor.usp.br
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/