[WP34S] Chi-Square Functions Broken in 2935
#1

Sorry to harp on the old stats distributions again, but this is a fundamentally serious bug, not discourse on improving accuracy, etc.

In short, the CDF and PDF return zero, irrespective of the arguments entered. The inverse CDF returns a roughly close result, quickly--presumably the unrefined initial estimate.

This occurs in both the real calc and the emulator and is new. It was working fine in 2932. The regularized incomplete gamma function still works fine, FWIW, so, at least in the CDF, the problem doesn't lie there.

Thought you'd like to know ASAP about a glitch so basic and important!

Cheers--

Les

#2

I can see the problem here.

The cause eludes me for the moment. The chi^2 code hasn't changed between those two versions which means something more subtle is at fault.


- Pauli

#3

Quote:
The chi^2 code hasn't changed between those two versions which means something more subtle is at fault.

Uh oh...

#4

I've changed some of the inner workings in the command dispatcher (there was a nasty bug lurking) and it is possible that it has broken things. I'll check.

#5

I can't see a problem with my Windows emulator here:

(Tracing is still not perfect as is the display of <= but that's not the point here.)

#6

Put 0 in register 00 and it will work! There was a typing error in chi2.wp34s, confounding x[<=]0? with x[<=]? 00. You can see it in the trace after the BSRF 001 command (replace the space by <=).

Will be fixed in the next build.

Edited: 3 May 2012, 4:42 a.m.

#7

My error :-( Nice find too.

I was suspecting something uninitialised being the root cause here. At home I wasn't reproducing the problem....


- Pauli

#8

Man that is subtle!

Sure enough, Reg 00 in both my calcs AND the emulator were nonzero due to some detritus left by user programs of my own I was running.

I verily thought I was losing my mind. Glad the fix was so easy, if tough to find.

Les

P.S. I am delighted to see work on the upper-tail distributions in progress. When computing inverses for cumulative probabilities close to 1, it will be much more accurate to work with an upper tail probability close to 0. One can already do this for the normal, t, and F. The first two are symmetrical and it boils down to working with the a probability less than .5 and flipping the sign later if necessary. The F distributions tails are related by swapping the numerator and denominator dfs and taking the reciprocal of the quantile. The challenging one with the existing calculator configuration is chi-squared, but this is simply a special case of the right-sided incomplete regularized gamma, which in the upper tail is very rapidly and accurately computed by its continued fraction. Indeed, I vote that the calculator offer the right-sided incomplete gamma explicitly so that it can be calculated accurately deep into the upper-tail if desired. I have no expertise with the other less common distributions so can't comment. My 2 euros.

Edited: 3 May 2012, 1:22 p.m.



Possibly Related Threads…
Thread Author Replies Views Last Post
  Broken keys on HP Prime? Tarcisi C 0 1,433 12-10-2013, 05:59 PM
Last Post: Tarcisi C
  [HP-Prime CAS] "Warning, ^ (Command) Is ambiguous on non square matrices"?? CompSystems 1 2,141 12-07-2013, 07:15 PM
Last Post: CompSystems
  [41CL] New Extra Functions version Monte Dalrymple 0 1,085 11-08-2013, 04:32 PM
Last Post: Monte Dalrymple
  HP Prime: in need of help with defining functions Alberto Candel 14 4,205 10-27-2013, 10:48 AM
Last Post: Alberto Candel
  Is my HP Prime broken ? Michael de Estrada 10 3,011 10-13-2013, 11:16 AM
Last Post: Mike Powell
  HP Prime spreadsheet functions SanS 0 1,877 10-04-2013, 04:23 AM
Last Post: SanS
  Stats functions on the HP34S Nicholas van Stigt 5 1,952 09-24-2013, 02:45 AM
Last Post: Nick_S
  Trig Functions Howard Owen 11 3,412 09-16-2013, 02:53 PM
Last Post: Fred Lusk
  50g piecewise functions Kurt Fankhauser 6 2,296 09-15-2013, 08:01 PM
Last Post: Kurt Fankhauser
  A brand new calculator benchmark: "middle square method seed test" Pier Aiello 25 6,771 09-13-2013, 01:58 PM
Last Post: Pier Aiello

Forum Jump: