Re: "Bug" in SurvivalDistribution?
- To: mathgroup at smc.vnet.net
- Subject: [mg127063] Re: "Bug" in SurvivalDistribution?
- From: Bob Hanlon <hanlonr357 at gmail.com>
- Date: Wed, 27 Jun 2012 04:12:54 -0400 (EDT)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- References: <201206260848.EAA06907@smc.vnet.net>
I see the same behavior on my Mac. Even with RandomReal replaced by
RandomVariate, the same problem occurs.
Bob Hanlon
On Tue, Jun 26, 2012 at 4:48 AM, Mark Coleman <markspcoleman at gmail.com> wrote:
> Greetings,
>
> I've bumped into an unusual condition using the built-in SurvivalDistribution function in Mathematica v8.0.4. I'm hoping another list member could run this code on their system and let me know if they get the same behavior.
>
> I've been doing some basic exploratory survival analysis and attempted to generate a SurvivalDistribution object for n=1900 event times. On my system (Dell T3400 with 8 Gb Ram running 64-bit Windows 7 Enterprise edition, SP1), the following code works for n<1900:
>
> n = 1899;
> data0 = RandomReal[LogNormalDistribution[1, 5], n];
> survDist = SurvivalDistribution[data0]
> Plot[SurvivalFunction[survDist, z], {z, 1, 100}]
>
> but when n>= 1900, my call to SurvivalDistribution does not evaluate properly.
>
> I've duiplicated this on a much larger Mathematica v8 Linux server. I've also reported it to WRI.
>
> Any idea on what's going on would be most welcomed. My final datasets contain much larger samples, n>100000.
>
> Thanks,
>
> Mark
>
--
Bob Hanlon
- References:
- "Bug" in SurvivalDistribution?
- From: Mark Coleman <markspcoleman@gmail.com>
- "Bug" in SurvivalDistribution?