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. |