Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.soft-sys.math.maple > #260

Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs!

Show all headers | View raw

jdm <james.d.mclaughlin@gmail.com> writes:

> with(Statistics):
> Phiinverse := (x) -> Quantile(Normal(0, 1), x) :
>
> evalf(Phiinverse(1-1/2^24));
>
> and the result is 5.691523310-2.147372008*I
>
> Now, I can get a real-valued solution to this by using the equivalent
> functionality from stats instead of Statistics.
>
> with(stats):
> PhiinvClassic:= (x) -> statevalf[icdf, normald[0, 1]](x):
>
> evalf(PhiinvClassic(1-1/2^24));
>
> However, this is not a perfect solution either:
>
> evalf(PhiinvClassic(1-1/2^24.5));
>
> Float(undefined)
>
> This may seem like an odd thing to use Maple for, but it is relevant
> to some cryptanalytic research I'm doing, and I'd be interested to
> know if this is a known bug and if anyone knows how to work around it.
>
> (Anyone who knows whether Mathematica also has this problem, please
> post!)

Try the following,

with(Statistics):
Phiinverse := (x) -> Quantile(Normal(0, 1), x):

f1 := Phiinverse(x);
                                                       1/2
                     f1 := RootOf(-erf(_Z) - 1 + 2 x) 2

f2 := subsindets(f1
                 , 'specfunc(anything,RootOf)'
                 , rootof -> unapply('fsolve'(op(rootof), _Z), x)
                );
                                                            1/2
               f2 := (x -> fsolve(-erf(_Z) - 1 + 2 x, _Z)) 2


evalf(f2(1-1/2^24));
                                  5.294704084


-- 
Joe Riel

Back to comp.soft-sys.math.maple | Previous | Next — Previous in thread | Next in thread | Find similar

Thread

Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! jdm <james.d.mclaughlin@gmail.com> - 2011-12-05 09:09 -0800
  Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! Joe Riel <joer@san.rr.com> - 2011-12-05 10:18 -0800
    Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! jdm <james.d.mclaughlin@gmail.com> - 2011-12-05 12:08 -0800
      Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! Joe Riel <joer@san.rr.com> - 2011-12-05 12:51 -0800
        Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! jdm <james.d.mclaughlin@gmail.com> - 2011-12-05 13:28 -0800
          Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! jdm <james.d.mclaughlin@gmail.com> - 2011-12-05 13:46 -0800
            Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! jdm <james.d.mclaughlin@gmail.com> - 2011-12-05 14:43 -0800
              Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! Joe Riel <joer@san.rr.com> - 2011-12-05 15:32 -0800
                Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! jdm <james.d.mclaughlin@gmail.com> - 2011-12-05 16:04 -0800
                Re: Inverse CDF of standard normal is (according to Maple 15) complex valued for some positive inputs! Herman Rubin <hrubin@skew.stat.purdue.edu> - 2011-12-07 02:06 +0000

csiph-web