Re: "Bug" in SurvivalDistribution?
- To: mathgroup at smc.vnet.net
- Subject: [mg127048] Re: "Bug" in SurvivalDistribution?
- From: Darren Glosemeyer <darreng at wolfram.com>
- Date: Wed, 27 Jun 2012 04:07:42 -0400 (EDT)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- References: <201206260848.EAA06907@smc.vnet.net>
On 6/26/2012 3:48 AM, Mark Coleman wrote:
> 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
>
This appears to be fixed in the development version. As a workaround for
now, you can evaluate the following (it seems the object being created
is not being flagged as a valid distribution):
Unprotect[DataDistribution];
DataDistribution /:
HoldPattern[
DistributionParameterQ[
DataDistribution["Survival", {_, _}, _, _]]] := True
After that,
n = 1900;
data0 = RandomReal[LogNormalDistribution[1, 5], n];
survDist = SurvivalDistribution[data0];
Plot[SurvivalFunction[survDist, x], {x, 1, 100}]
will work just fine.
Darren Glosemeyer
Wolfram Research
- References:
- "Bug" in SurvivalDistribution?
- From: Mark Coleman <markspcoleman@gmail.com>
- "Bug" in SurvivalDistribution?