Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 07-26-2013 11:10 AM
(1673 views)

Dear all,

I am evaluating a double integral using the quad function but i got some errors. The integral is f(Z).g(Q).T(Z,Q) where Z is N(0,1) and Q follows a scaled chi distribution .The program is attached and here is a part of the error:" Convergence could not be attained over the subinterval ( .M , .P)". Could you help me to figure out the problem?

Thanks

Aya

17 REPLIES 17

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Be sure to test and debug each function (integrand) before you try to integrate them. Graphing the functions are a good way to ensure that they are implemented correctly.

An integrand will return a scalar quantity when it receives a scalar input. It looks like the INNER function is returning a vector of values:

Q = 1;

y = inner(0);

print y; /* should be a scalar */

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

I am sorry that I am asking too many questions concerning this code.But it is really very important in my research. I have been trying to correct the errors concerning the singularity of a matrix, so i used ginv function to avoid these errors. However, I got other errors related to

Convergence could not be attained over the subinterval

( 0 , .P)

Convergence could not be attained over the subinterval

( .M , .P)

Could you please help me to point out the reasons behind these errors?

Waiting for your reply,

Aya

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

The plots are attached. Each plot is for a certain QQ value (1,3,5) and different ZED values(-5,-4,-3,-2,-1,0,1,2,3,4,5). I have made many trials to have this integration done, but I am always faced with errors. Sometimes related to the singularity of a matrix, so i tried the " ginv" function, so I got errors related to the convergence of the integral over the specified interval. I really can't figure out what the problem is. I tried running the syntax for other integration limits; ZEDinterval (-1,1) and Qinterval (0,1), It gave me an answer without errors " 1.659E-59". The correct answer over the whole interval (-inf,inf) and(0,inf) is supposed to be 339.4!!

Aya

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

I think the convergence problem may be due to scaling. The y axis ranges from near 0 to 10^15. I would suggest integrating the log of the function, and applying the chain rule to get back to the integrated value of the basis function.

Steve Denham

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

While the OP doesn't say f(x) > 0, the plots that were attached certainly were defined only on f(x) > 0. So, here's my thinking. We want Integral f(x) dx. We find Integral log(f(x)) dx, by taking the log of the response variable, and integrating. Then I was thinking because Integral (log x) dx was x log x - x, we could just substitute. But after putting pencil to paper, that's probably not going to work, as the the original function won't be unaffected by the integration.

Steve Denham

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

It's always hard to understand complicated code that someone else wrote, but at the end of the INNER module you have

normal=(1/sqrt(2*pi))*exp(-0.5*(ZED**2));

b=(v*avrl)cc*normal;

return(b);

To me, it looks like NORMAL is simply the standard normal PDF. The quantity c=(v*avl) is a constant that depends on Q. Therefore the

integral of INNER over the domain (-Infinity, Infinity) is c, because the integral of the PDF is unity. Therefore it is not necessary to call the QUAD function for INNER.

It looks like OUTER might be handled similarly. If OUTER is some multiple k of the chi-square distribution, then the integral on [0, infinity) is simply k.

Incidentally, if you are interested in other limits of integration, you can use the CDF function to compute integrals of a PDF.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Thank you so much for your efforts. i really appreciate it. My integration is the product of three quantities

1-standard normal pdf "normal"

2- a scaled chi distribution with degrees of freedom "df"

3- The quantity (v*avrl) which is a function of both variables QQ and ZED where ZED is the variable following the standard normal distribution and QQ is the variable following a scaled chi distribution.

The inner module is used to compute the quantity "v*avrl". After computing the quantitty "v*avrl" which depends on both ZED and QQ, the inner integrand is computed by multipying "v*avrl" by the normal density.

The result of the inner integration should be a function of QQ which is called "w" in our syntax. This w is multiplied by the scaled chi distribution: **(1**/cf)*(**1**/sqrt(df))*(**2****(**1**-(df/**2**)))*(QQ**(df-**1**))*(exp((-QQ****2**)/**2**))/gamma(df/**2**);

Thus, i think I have to call QUAD for both the inner and outer integrals since "v*avrl" depends on both QQ and ZED, see this part in the syntax:

QQstar1=(QQ*z[i, ])+(ZED/(sqrt(mm)));

QQstar=QQstar1-shift;

kone=QQstar+QQ*fun1;

ktwo=QQstar+QQ*fun2;

pij=probnorm(kone)-probnorm(ktwo);

Thus "v*avrl"is not a constant with respect to the inner integral. Am I correct??

Aya

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Yep. I missed that v*avrl depends on ZED. Too bad.

I don't suppose you know anything analytically about the function f(z)=v(z)*avrl(z), such as values at +/- infinity and derivative info? If you do, you might be able to integrate by parts using the fact that phi(x) has a known antiderivative, the normal CDF.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

No, I don't know the antiderivative information. What I know is that the v*avrl is the average run length of a certain control chart. I am computing the average run length using the markov chain technique. The martrix "r" is the transitional probability matrix and I have t transiet states.

Aya

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

I recommend the following: factor out c = f(0)=v(0)*avrl(0) from the integral. Then you have the integral of

f(z)*phi(z) where phi is the standard normal PDF. By scaling the integral, you might get convergence.

You can call QUAD on this bounded integrand, and then multiply the result of QUAD by c=f(0) in the optimization code.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Sorry for repeating my question, but I don't know what do you mean by scaling the integral? I tried to understand the procedure you explained but I couldn't get the point.

Thanks,

Aya Emara

Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.

**If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website. **

Multiple Linear Regression in SAS

Learn how to run multiple linear regression models with and without interactions, presented by SAS user Alex Chaplin.

Find more tutorials on the SAS Users YouTube channel.