FAQFAQ   SearchSearch   MemberlistMemberlist   UsergroupsUsergroups 
 ProfileProfile   PreferencesPreferences   Log in to check your private messagesLog in to check your private messages   Log inLog in 
Forum index » Science and Technology » Math » num-analysis
Is there a program for numerically integrating a PDF to get CDF
Post new topic   Reply to topic Page 1 of 1 [3 Posts] View previous topic :: View next topic
Author Message
glare22@gmail.com
science forum beginner


Joined: 12 Jul 2006
Posts: 7

PostPosted: Wed Jul 12, 2006 8:38 pm    Post subject: Re: Is there a program for numerically integrating a PDF to get CDF Reply with quote

Thanks a lot.

I have tried to use dqk61 to estimate the value in a small range of
0.0125 where the peak exists, I get a very large number. I will try
dqagi and dqags.


Peter Spellucci wrote:
Quote:
In article <1152695303.927180.168540@m79g2000cwm.googlegroups.com>,
glare22@gmail.com writes:
I want to get the CDF from the integration of a given PDF. Is there a
program I can use for this purpose? Currently I use qk15i from the
quadpack. But the result is not good. For the PDF with a sharp peak, it
gives totally wrong results. Can anyone help me? Thanks.


Liang


you are using the wrong code:
1. quadpack/qk15i.f ( plus dependencies )
gams: H2a3a2,H2a4a2
for: estimate 1D integral on (semi)infinite interval using a 15 point quadrature rule,non-automatic
prec: single
Score: 100%

first of all, switch to an adaptive automatic integrator. then you also must be
able to estimate the witdth of possible peaks in order to set the tolerances \so
tight such that the adaptive mechanism is forced the refine the grid
correspondingly. this is not too hard to do. with prior knowledge of
the integrand you can outwit any numerical integrator: it "sees" the integrand
only at the grid which is constructed automatically using some rules.
In your case, where you want to construct a table of a indefinite integral,
you can proceed stepwise, using a semiinfinte interval first with a strong
tolerance and then using finite intervals using other adaptive code from
quadpack.
dqagi and dqags would be a choice
hth
peter
Back to top
Peter Spellucci
science forum Guru


Joined: 29 Apr 2005
Posts: 702

PostPosted: Wed Jul 12, 2006 10:27 am    Post subject: Re: Is there a program for numerically integrating a PDF to get CDF Reply with quote

In article <1152695303.927180.168540@m79g2000cwm.googlegroups.com>,
glare22@gmail.com writes:
Quote:
I want to get the CDF from the integration of a given PDF. Is there a
program I can use for this purpose? Currently I use qk15i from the
quadpack. But the result is not good. For the PDF with a sharp peak, it
gives totally wrong results. Can anyone help me? Thanks.


Liang


you are using the wrong code:
1. quadpack/qk15i.f ( plus dependencies )
gams: H2a3a2,H2a4a2
for: estimate 1D integral on (semi)infinite interval using a 15 point quadrature rule,non-automatic
prec: single
Score: 100%

first of all, switch to an adaptive automatic integrator. then you also must be
able to estimate the witdth of possible peaks in order to set the tolerances \so
tight such that the adaptive mechanism is forced the refine the grid
correspondingly. this is not too hard to do. with prior knowledge of
the integrand you can outwit any numerical integrator: it "sees" the integrand
only at the grid which is constructed automatically using some rules.
In your case, where you want to construct a table of a indefinite integral,
you can proceed stepwise, using a semiinfinte interval first with a strong
tolerance and then using finite intervals using other adaptive code from
quadpack.
dqagi and dqags would be a choice
hth
peter
Back to top
glare22@gmail.com
science forum beginner


Joined: 12 Jul 2006
Posts: 7

PostPosted: Wed Jul 12, 2006 9:08 am    Post subject: Is there a program for numerically integrating a PDF to get CDF Reply with quote

I want to get the CDF from the integration of a given PDF. Is there a
program I can use for this purpose? Currently I use qk15i from the
quadpack. But the result is not good. For the PDF with a sharp peak, it
gives totally wrong results. Can anyone help me? Thanks.


Liang
Back to top
Google

Back to top
Display posts from previous:   
Post new topic   Reply to topic Page 1 of 1 [3 Posts] View previous topic :: View next topic
The time now is Thu Dec 14, 2017 4:56 pm | All times are GMT
Forum index » Science and Technology » Math » num-analysis
Jump to:  

Similar Topics
Topic Author Forum Replies Last Post
No new posts How to solve linear program with matrix variable ? Fan num-analysis 4 Thu Jul 13, 2006 5:55 am
No new posts Is there a program for numerically integrating a PDF to g... glare22@gmail.com num-analysis 5 Wed Jul 12, 2006 9:02 am
No new posts DONT let the CHILD RAPING EVIL AMERICAN GOVT PSYCHOPATHS ... McGinnn Math 2 Sun Jul 09, 2006 10:28 pm
No new posts Integrating t^(-3/2) * e^(1/t) Konrad Viltersten Math 7 Wed Jul 05, 2006 2:27 pm
No new posts eCheminfo Drug Discovery Community of Practice Meeting, B... The Ferryman Chem 1 Tue Jun 27, 2006 2:17 pm

Copyright © 2004-2005 DeniX Solutions SRL
Other DeniX Solutions sites: Electronics forum |  Medicine forum |  Unix/Linux blog |  Unix/Linux documentation |  Unix/Linux forums  |  send newsletters
 


Powered by phpBB © 2001, 2005 phpBB Group
[ Time: 0.0193s ][ Queries: 20 (0.0036s) ][ GZIP on - Debug on ]