Logistic Log-Likelihood Functions and their Partials

Logistic 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{{{T}_{i}}-\mu }{\sigma }}}}{\sigma {{(1+{{e}^{\tfrac{{{T}_{i}}-\mu }{\sigma }}})}^{2}}} \right)-\underset{i=1}{\mathop{\overset{S}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime }}\ln (1+{{e}^{\tfrac{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{T_{{{L}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}}-\frac{1}{1+{{e}^{\tfrac{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 logistic shape parameter (unknown a priori, the first of two parameters to be found)

• [math]\displaystyle{ \eta }[/math] is the logistic 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 group

• [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{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{{{e}^{\tfrac{{{T}_{i}}-\mu }{\sigma }}}}{1+{{e}^{\tfrac{{{T}_{i}}-\mu }{\sigma }}}}+\frac{1}{\sigma }\underset{i=1}{\mathop{\overset{S}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime }}\frac{{{e}^{\tfrac{T_{i}^{^{\prime }}-\mu }{\sigma }}}}{1+{{e}^{\tfrac{T_{i}^{^{\prime }}-\mu }{\sigma }}}} \\ & & -\frac{\underset{i=1}{\mathop{\overset{{{F}_{I}}}{\mathop{\mathop{}_{}^{}}}\,}}\,N_{i}^{^{\prime \prime }}}{\sigma }+\frac{1}{\sigma }\underset{i=1}{\mathop{\overset{{{F}_{I}}}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime \prime }}\left( \frac{{{e}^{\tfrac{T_{{{L}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}}{1+{{e}^{\tfrac{T_{{{L}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}}+\frac{{{e}^{\tfrac{T_{{{R}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}}{1+{{e}^{\tfrac{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{{{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{{{T}_{i}}-\mu }{\sigma }{{e}^{\tfrac{{{T}_{i}}-\mu }{\sigma }}}}{1+{{e}^{\tfrac{{{T}_{i}}-\mu }{\sigma }}}} \\ & & +\frac{1}{\sigma }\underset{i=1}{\mathop{\overset{S}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{^{\prime }}\frac{\tfrac{T_{i}^{^{\prime }}-\mu }{\sigma }{{e}^{\tfrac{T_{i}^{^{\prime }}-\mu }{\sigma }}}}{1+{{e}^{\tfrac{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{T_{{{L}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }{{e}^{\tfrac{T_{{{L}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}}{1+{{e}^{\tfrac{T_{{{L}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}}+\frac{\tfrac{T_{{{R}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }{{e}^{\tfrac{T_{{{R}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}}{1+{{e}^{\tfrac{T_{{{R}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}} \\ & & -\frac{\tfrac{T_{{{R}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }{{e}^{\tfrac{T_{{{R}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}-\tfrac{T_{{{L}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }{{e}^{\tfrac{T_{{{L}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}}{{{e}^{\tfrac{T_{{{R}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}-{{e}^{\tfrac{T_{{{L}_{i}}}^{^{\prime \prime }}-\mu }{\sigma }}}}) \end{align} }[/math]