[WP34s] Regularized incomplete Beta function
05-08-2014, 11:29 PM
Post: #41
 Paul Dale Senior Member Posts: 1,725 Joined: Dec 2013
RE: [WP34s] Regularized incomplete Beta function
(05-08-2014 10:58 PM)Dieter Wrote:  The user may enter dof values that lead to problems and strange results. At the moment I am looking at the Chi² quantile. For Chi² = n the CDF should be near 0,5 (slightly larger). Try this with n=2000 and you get 0,5042... (correct). Then try n=3000 and the 34s returns zero (!).

This is a bug, not a gradual loss of accuracy

The series expansion for the incomplete gamma function (as opposed to the continued fraction variation) isn't converging in 500 iterations for x around the degrees of freedom. For the n = 2000 case it is taking 440 iterations to converge. For n = 3000, it wants 532 iterations. For n = 250000, 4591 iterations are required.

Perhaps I got the threshold for switching between the continued fraction and the series expansion wrong.

This might also explain the gradual loss of digits you are seeing at time. That would match the behaviour of the continued fraction path which attempts to return something based on the best guess. The series code just returns zero when its iteration count expires.

Pauli
 « Next Oldest | Next Newest »

 Messages In This Thread [WP34s] Regularized incomplete Beta function - Dieter - 05-01-2014, 07:17 PM RE: [WP34s] Regularized incomplete Beta function - Thomas Klemm - 05-01-2014, 07:48 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-01-2014, 08:00 PM RE: [WP34s] Regularized incomplete Beta function - Thomas Klemm - 05-01-2014, 08:28 PM RE: [WP34s] Regularized incomplete Beta function - walter b - 05-01-2014, 09:08 PM RE: [WP34s] Regularized incomplete Beta function - walter b - 05-01-2014, 09:23 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-01-2014, 09:26 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-01-2014, 09:55 PM RE: [WP34s] Regularized incomplete Beta function - walter b - 05-01-2014, 10:14 PM RE: [WP34s] Regularized incomplete Beta function - Thomas Klemm - 05-01-2014, 11:13 PM RE: [WP34s] Regularized incomplete Beta function - Thomas Klemm - 05-01-2014, 11:01 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-02-2014, 08:48 PM RE: [WP34s] Regularized incomplete Beta function - walter b - 05-02-2014, 09:25 PM RE: [WP34s] Regularized incomplete Beta function - Thomas Klemm - 05-02-2014, 10:22 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-02-2014, 10:34 PM RE: [WP34s] Regularized incomplete Beta function - walter b - 05-02-2014, 10:46 PM RE: [WP34s] Regularized incomplete Beta function - Thomas Klemm - 05-02-2014, 11:29 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-02-2014, 09:03 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-04-2014, 06:29 AM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-04-2014, 05:26 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-05-2014, 02:20 AM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-05-2014, 09:13 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-05-2014, 09:57 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-05-2014, 08:26 AM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-11-2014, 01:50 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-11-2014, 09:51 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-12-2014, 06:18 AM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-12-2014, 06:58 AM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-12-2014, 07:47 PM RE: [WP34s] Regularized incomplete Beta function - walter b - 05-12-2014, 08:14 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-12-2014, 09:59 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-22-2014, 08:51 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-25-2014, 09:43 PM RE: [WP34s] Regularized incomplete Beta function - Manolo Sobrino - 05-03-2014, 01:38 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-05-2014, 11:02 AM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-05-2014, 08:30 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-05-2014, 09:49 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-05-2014, 10:04 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-05-2014, 10:16 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-05-2014, 11:09 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-06-2014, 12:54 AM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-06-2014, 11:01 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-07-2014, 12:12 AM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-08-2014, 06:26 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-08-2014, 10:28 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-08-2014, 10:58 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-08-2014 11:29 PM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-09-2014, 09:06 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-10-2014, 12:03 AM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-11-2014, 01:23 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-11-2014, 10:33 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-05-2014, 11:15 AM RE: [WP34s] Regularized incomplete Beta function - Dieter - 05-05-2014, 08:14 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-05-2014, 09:40 PM RE: [WP34s] Regularized incomplete Beta function - Paul Dale - 05-12-2014, 06:05 AM

User(s) browsing this thread: 1 Guest(s)