close
Warning:
Can't synchronize with repository "(default)" (/var/svn/tolp does not appear to be a Subversion repository.). Look in the Trac log for more information.
- Timestamp:
-
Dec 21, 2010, 1:36:51 PM (15 years ago)
- Author:
-
Víctor de Buen Remiro
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v11
|
v12
|
|
80 | 80 | [[LatexEquation( \frac{\partial L\left(\beta\right)}{\partial\beta_{i}\partial_{j}}=-\underset{k}{\sum}w_{k}\frac{\left(1+e^{-x_{k}^{t}\beta}\right)e^{-x_{k}^{t}\beta}x_{ki}x_{kj}-\left(e^{-x_{k}^{t}\beta}\right)^{2}x_{ki}x_{kj}}{\left(1+e^{-x_{k}^{t}\beta}\right)}=-\underset{k}{\sum}x_{ki}x_{kj}w_{k}\pi_{i}\left(1-\pi_{i}\right) )]] |
81 | 81 | |
82 | | Desde el punto de vista del cálculo numérico en aritmética discreta hay que tener en cuenta que |
| 82 | From the standpoint of arithmetic discrete numerical calculation must take into account that[[BR]] |
83 | 83 | |
84 | 84 | [[LatexEquation( e^{710}=\infty )]] [[BR]] [[BR]] |
85 | 85 | |
86 | | por lo que hay que tratar de forma cuidadosa con las expresiones que contienen la exponencial. |
87 | | En este caso concreto se hará uso de las siguientes igualdades asintóticas |
| 86 | For this reason we must carefully try to contain the exponential expressions. |
| 87 | In this case it will use the following asymptotic equalities[[BR]] |
88 | 88 | |
89 | 89 | [[LatexEquation( \ln\left(1+e^{-z}\right)\;\overset{z\rightarrow-\infty}{\longrightarrow}\;-z )]] [[BR]] [[BR]] |
90 | | |
91 | 90 | [[LatexEquation( \frac{e^{-z}}{1+e^{-z}}\;\overset{z\rightarrow-\infty}{\longrightarrow}\;1 )]] [[BR]] [[BR]] |
92 | 91 | |
… |
… |
|
112 | 111 | [[LatexEquation( \frac{\partial L\left(\beta\right)}{\partial\beta_{i}\partial_{j}}=-\underset{k}{\sum}w_{k}\phi\left(x_{k}\beta\right)\left(y_{k}\frac{z\Phi\left(x_{k}\beta\right)+\phi\left(x_{k}\beta\right)}{\Phi^{2}\left(x_{k}\beta\right)}+\left(1-y_{k}\right)\frac{-z \Phi\left(-x_{k}\beta\right)+\phi\left(x_{k}\beta\right)}{\Phi\left(-x_{k}\beta\right)^{2}}\right)x_{ik}x_{jk} )]] |
113 | 112 | |
| 113 | To avoid numerical problems will use the following equality |
114 | 114 | |
| 115 | [[LatexEquation( \ln\left(\Phi\left(z\right)\right)=\ln\left(1-erf\left(\frac{-z}{\sqrt{2}}\right)\right)-\ln2 )]] [[BR]] [[BR]] |
| 116 | |
| 117 | The function logarithm of complemetary error function |
| 118 | |
| 119 | [[LatexEquation( \ln\left(1-erf\left(u\right)\right) )]] |
| 120 | |
| 121 | is implemented as |
| 122 | [http://www.gnu.org/software/gsl/manual/html_node/Log-Complementary-Error-Function.html gsl_sf_log_erfc] |
| 123 | that is available in TOL. |
| 124 | |
| 125 | The Hazard function [[BR]] [[BR]] |
| 126 | |
| 127 | [[LatexEquation( h\left(z\right)=\frac{\phi\left(z\right)}{\Phi\left(z\right)} )]] |
| 128 | |
| 129 | decreases rapidly as [[LatexEquation(z)]] approaches [[LatexEquation(-\infty)]] and asymptotes to |
| 130 | [[LatexEquation(h\left(z\right) \sim z )]] as [[LatexEquation(z)]] approaches [[LatexEquation(+\infty.)]] |
| 131 | |
| 132 | Hazard function is implemented as [http://www.gnu.org/software/gsl/manual/html_node/Probability-functions.html gsl_sf_hazard ] |
| 133 | that is also available in TOL. |