Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2012

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

Search the Archive

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
>



  • Prev by Date: Re: Options Question ver 8
  • Next by Date: Re: filling
  • Previous by thread: Re: "Bug" in SurvivalDistribution?
  • Next by thread: Options Question ver 8