The Loglogistic Log-Likelihood Functions and their Partials

The Loglogistic Log-Likelihood Functions and their Partials

This log-likelihood function is composed of three summation portions:

[math]\displaystyle{ \begin{align} \ln (L)= & \Lambda =\underset{i=1}{\mathop{\overset{{{F}_{e}}}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,{{N}_{i}}\ln \left( \frac{{{e}^{\tfrac{\ln ({{T}_{i}})-\mu }{\sigma }}}}{\sigma t{{(1+{{e}^{\tfrac{\ln ({{T}_{i}})-\mu }{\sigma }}})}^{2}}} \right) \\ & -\underset{i=1}{\mathop{\overset{S}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime }}\ln (1+{{e}^{\tfrac{\ln (T_{i}^{^{\prime }})-\mu }{\sigma }}}) \\ & +\underset{i=1}{\mathop{\overset{{{F}_{I}}}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime \prime }}\ln \left( \frac{1}{1+{{e}^{\tfrac{\ln (T_{{{L}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}}-\frac{1}{1+{{e}^{\tfrac{\ln (T_{{{R}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}} \right) \end{align} }[/math]

where:

• [math]\displaystyle{ {{F}_{e}} }[/math] is the number of groups of times-to-failure data points

• [math]\displaystyle{ {{N}_{i}} }[/math] is the number of times-to-failure in the [math]\displaystyle{ {{i}^{th}} }[/math] time-to-failure data group

• [math]\displaystyle{ \mu }[/math] is the loglogistic shape parameter (unknown a priori, the first of two parameters to be found)

• [math]\displaystyle{ \sigma }[/math] is the loglogistic scale parameter (unknown a priori, the second of two parameters to be found)

• [math]\displaystyle{ {{T}_{i}} }[/math] is the time of the [math]\displaystyle{ {{i}^{th}} }[/math] group of time-to-failure data

• [math]\displaystyle{ S }[/math] is the number of groups of suspension data points

• [math]\displaystyle{ N_{i}^{\prime } }[/math] is the number of suspensions in [math]\displaystyle{ {{i}^{th}} }[/math] group of suspension data points

• [math]\displaystyle{ T_{i}^{\prime } }[/math] is the time of the [math]\displaystyle{ {{i}^{th}} }[/math] suspension data group

• [math]\displaystyle{ FI }[/math] is the number of interval failure data groups,

• [math]\displaystyle{ N_{i}^{\prime \prime } }[/math] is the number of intervals in [math]\displaystyle{ {{i}^{th}} }[/math] group of data intervals

• [math]\displaystyle{ T_{Li}^{\prime \prime } }[/math] is the beginning of the [math]\displaystyle{ {{i}^{th}} }[/math] interval

• and [math]\displaystyle{ T_{Ri}^{\prime \prime } }[/math] is the ending of the [math]\displaystyle{ {{i}^{th}} }[/math] interval

For the purposes of MLE, left censored data will be considered to be intervals with [math]\displaystyle{ T_{Li}^{\prime \prime }=0. }[/math]

The solution of the maximum log-likelihood function is found by solving for ([math]\displaystyle{ \widehat{\mu },\widehat{\sigma }) }[/math] so that [math]\displaystyle{ \tfrac{\partial \Lambda }{\partial \mu }=0,\tfrac{\partial \Lambda }{\partial \sigma }=0. }[/math]

[math]\displaystyle{ \begin{align} \frac{\partial \Lambda }{\partial \mu }= & -\frac{\underset{i=1}{\mathop{\overset{{{F}_{e}}}{\mathop{\mathop{}_{}^{}}}\,}}\,{{N}_{i}}}{\sigma }+\frac{2}{\sigma }\underset{i=1}{\mathop{\overset{{{F}_{e}}}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,{{N}_{i}}\frac{{{e}^{\tfrac{\ln ({{T}_{i}})-\mu }{\sigma }}}}{1+{{e}^{\tfrac{\ln ({{T}_{i}})-\mu }{\sigma }}}} \\ & +\frac{1}{\sigma }\underset{i=1}{\mathop{\overset{S}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime }}\frac{{{e}^{\tfrac{\ln (T_{i}^{^{\prime }})-\mu }{\sigma }}}}{1+{{e}^{\tfrac{\ln (T_{i}^{^{\prime }})-\mu }{\sigma }}}}-\frac{{{F}_{I}}}{\sigma } \\ & +\frac{1}{\sigma }\underset{i=1}{\mathop{\overset{{{F}_{I}}}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime \prime }}\left( \frac{{{e}^{\tfrac{\ln (T_{{{L}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}}{1+{{e}^{\tfrac{\ln (T_{{{L}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}}+\frac{{{e}^{\tfrac{\ln (T_{{{R}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}}{1+{{e}^{\tfrac{\ln (T_{{{R}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}} \right) \end{align} }[/math]

[math]\displaystyle{ \begin{align} \frac{\partial \Lambda }{\partial \sigma }= & -\underset{i=1}{\mathop{\overset{{{F}_{e}}}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,{{N}_{i}}\frac{\ln ({{T}_{i}})-\mu }{{{\sigma }^{2}}}-\frac{1}{\sigma }\underset{i=1}{\mathop{\overset{{{F}_{e}}}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,{{N}_{i}}+\frac{2}{\sigma }\underset{i=1}{\mathop{\overset{{{F}_{e}}}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,{{N}_{i}}\frac{\tfrac{\ln ({{T}_{i}})-\mu }{\sigma }{{e}^{\tfrac{\ln ({{T}_{i}})-\mu }{\sigma }}}}{1+{{e}^{\tfrac{\ln ({{T}_{i}})-\mu }{\sigma }}}} \\ & +\frac{1}{\sigma }\underset{i=1}{\mathop{\overset{S}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime }}\frac{\tfrac{\ln (T_{i}^{^{\prime }})-\mu }{\sigma }{{e}^{\tfrac{\ln (T_{i}^{^{\prime }})-\mu }{\sigma }}}}{1+{{e}^{\tfrac{\ln (T_{i}^{^{\prime }})-\mu }{\sigma }}}} \\ & \frac{1}{\sigma }\underset{i=1}{\mathop{\overset{{{F}_{I}}}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime \prime }}(\frac{\tfrac{\ln (T_{{{L}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }{{e}^{\tfrac{\ln (T_{{{L}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}}{1+{{e}^{\tfrac{\ln (T_{{{L}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}}+\frac{\tfrac{T_{{{R}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }{{e}^{\tfrac{\ln (T_{{{R}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}}{1+{{e}^{\tfrac{\ln (T_{{{R}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}} \\ & -\frac{\tfrac{\ln (T_{{{R}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }{{e}^{\tfrac{\ln (T_{{{R}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}-\tfrac{\ln (T_{{{L}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }{{e}^{\tfrac{\ln (T_{{{L}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}}{{{e}^{\tfrac{\ln (T_{{{R}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}-{{e}^{\tfrac{\ln (T_{{{L}_{i}}}^{^{\prime \prime }})-\mu }{\sigma }}}}) \end{align} }[/math]