ヘッセ行列
Using this score function and Hessian matrix, the partial likelihood can be maximized using the Newton-Raphson algorithm. The inverse of the Hessian matrix, evaluated at the estimate of "β", can be used as an approximate variance-covariance matrix for the estimate, and used to produce approximate standard errors for the regression coefficients.