Re: "Bug" in SurvivalDistribution?
- To: mathgroup at smc.vnet.net
- Subject: [mg127044] Re: "Bug" in SurvivalDistribution?
- From: Mark Coleman <markspcoleman at gmail.com>
- Date: Wed, 27 Jun 2012 04:06:19 -0400 (EDT)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- References: <201206260848.EAA06907@smc.vnet.net> <CAEtRDSeXh9Atxybz2dPJ2tzkSuRwv-OdspLK6kRQfTL10V+rqQ@mail.gmail.com>
Hi Bob,
I heard back from a contact at WRI. It does seem to be a bug. He sent me a
work around for the issue:
Unprotect[DataDistribution];
DataDistribution /:
HoldPattern[
DistributionParameterQ[
DataDistribution["Survival", {_, _}, _, _]]] := True
I ran this code on my Windows machine and it works correctly, even for very
large sample sizes.
Best,
Mark
On Tue, Jun 26, 2012 at 9:41 AM, Bob Hanlon <hanlonr357 at gmail.com> wrote:
> 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
> 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
> >
>
>
>
> --
> Bob Hanlon
>
- References:
- "Bug" in SurvivalDistribution?
- From: Mark Coleman <markspcoleman@gmail.com>
- "Bug" in SurvivalDistribution?