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>