[Date Index]
[Thread Index]
[Author Index]
Re: weibull plot on weibull scaled paper
*To*: mathgroup at smc.vnet.net
*Subject*: [mg116722] Re: weibull plot on weibull scaled paper
*From*: DrMajorBob <btreat1 at austin.rr.com>
*Date*: Fri, 25 Feb 2011 06:34:42 -0500 (EST)
You're probably right.
Bobby
On Thu, 24 Feb 2011 05:21:14 -0600, Darren Glosemeyer
<darreng at wolfram.com> wrote:
> On 2/23/2011 4:24 AM, Bill Rowe wrote:
>> On 2/22/11 at 1:12 PM, btreat1 at austin.rr.com (DrMajorBob) wrote:
>>
>>> This gave the FindRoot:jsing error on the first try just now... and
>>> also the third try.
>>> dist = WeibullDistribution[7, 200];
>>> data = RandomVariate[dist, 300];
>>> ProbabilityScalePlot[data, "Weibull"]
>>> QuantilePlot[data, dist]
>>> So the error is not unusual.
>> Hmmm... If I copy and paste just
>>
>> dist = WeibullDistribution[7, 200];
>> data = RandomVariate[dist, 300];
>> ProbabilityScalePlot[data, "Weibull"]
>>
>> into a single cell then execute it, I see the error you
>> reported. Alternatively, if I type the above by hand into a
>> single cell I still see the error you reported consistently. Or
>> I can copy each line one by one pasting each into a cell then
>> execute each line before pasting the next and still get the error.
>>
>> But if I enter each line by hand into individual cells, I don't
>> see the error. So, there appears to be a bug somewhere. That is,
>> the various ways I've described above of doing the computation
>> should all behave identically.
>>
>>
>
> ProbabilityScalePlot and its friends call the internal parameter
> estimation code for FindDistributionParameters/EstimatedDistribution to
> get the estimates. The reason for the message is that the equations in
> the estimation code for WeibullDistribution were less numerically stable
> than they could have been and so convergence fails more often than it
> should when parameter values are somewhat large. This has been fixed for
> the next release.
>
> I think the fact that the example worked without error when you typed
> the inputs in separate cells is a red herring. The example will converge
> for some data sets and not for others. I suspect the random values you
> got in the case that worked just happened to give a convergent result
> while the random values you'd gotten in the other attempts did not.
>
> Darren Glosemeyer
> Wolfram Research
>
--
DrMajorBob at yahoo.com
Prev by Date:
**Re: weibull plot on weibull scaled paper**
Next by Date:
**Exporting video of simulation**
Previous by thread:
**Re: weibull plot on weibull scaled paper**
Next by thread:
**Re: weibull plot on weibull scaled paper**
| |