Main Content


Weibull negative log-likelihood


nlogL = wbllike(params,data)
[logL,AVAR] = wbllike(params,data)
[...] = wbllike(params,data,censoring)
[...] = wbllike(params,data,censoring,freq)


nlogL = wbllike(params,data) returns the Weibull log-likelihood. params(1) is the scale parameter, A, and params(2) is the shape parameter, B.

[logL,AVAR] = wbllike(params,data) also returns AVAR, which is the asymptotic variance-covariance matrix of the parameter estimates if the values in params are the maximum likelihood estimates. AVAR is the inverse of Fisher's information matrix. The diagonal elements of AVAR are the asymptotic variances of their respective parameters.

[...] = wbllike(params,data,censoring) accepts a Boolean vector, censoring, of the same size as data, which is 1 for observations that are right-censored and 0 for observations that are observed exactly.

[...] = wbllike(params,data,censoring,freq) accepts a frequency vector, freq, of the same size as data. freq typically contains integer frequencies for the corresponding elements in data, but can contain any nonnegative values. Pass in [] for censoring to use its default value.

The Weibull negative log-likelihood for uncensored data is


where f is the Weibull pdf.

wbllike is a utility function for maximum likelihood estimation.


This example continues the example from wblfit.

r = wblrnd(0.5,0.8,100,1);
[logL, AVAR] = wbllike(wblfit(r),r)
logL =
  0.0048  0.0014
  0.0014  0.0040


[1] Patel, J. K., C. H. Kapadia, and D. B. Owen. Handbook of Statistical Distributions. New York: Marcel Dekker, 1976.

Version History

Introduced before R2006a