MathGroup Archive 2001

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

Search the Archive

Re: weibull distribution


In a message dated 2001/10/31 3:46:48 AM, chandrika at ittc.ukans.edu writes:

>I am trying to fit data to weibul distribution. Is there a way to do
>using mathematica?
>

Needs["Statistics`ContinuousDistributions`"];

data = RandomArray[WeibullDistribution[5., 5.], 100];

logPDF[x_] := Evaluate[PowerExpand[Log[PDF[WeibullDistribution[a, b], x]]]];

logProd = Plus @@ (logPDF[#]& /@ data);

eqnMLE = ((D[logProd, #] == 0)& /@ {a, b});

soln = FindRoot[eqnMLE, {a, 1.}, {b, Mean[data]}]

{a -> 4.99725689108074, 
  b -> 5.033979933866775}


Bob Hanlon
Chantilly, VA  USA


  • Prev by Date: Re: graphics on the same line
  • Next by Date: Re: graphics on the same line
  • Previous by thread: Re: weibull distribution
  • Next by thread: Re: weibull distribution