Services & ResourcesWolfram ForumsMathGroup Archive

[Date Index] [Thread Index] [Author Index]

"Bug" in SurvivalDistribution?

• To: mathgroup at smc.vnet.net
• Subject: [mg127036] "Bug" in SurvivalDistribution?
• From: Mark Coleman <markspcoleman at gmail.com>
• Date: Tue, 26 Jun 2012 04:48:00 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com

Greetings,

I've bumped into an unusual condition using the built-in SurvivalDistributi=
on 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

• Follow-Ups:
  • Re: "Bug" in SurvivalDistribution?
    • From: Darren Glosemeyer <darreng@wolfram.com>
  • Re: "Bug" in SurvivalDistribution?
    • From: Bob Hanlon <hanlonr357@gmail.com>