How to find the scoring vector and the Hessian array in R
to apply the Newton Raphson
method in the code below:
rm(list=ls())
cat("4")
#Função para simular variáveis aleatórias de uma modelo de regressão exponencial.
simula.exponencial <- function(formula, beta) {
X <- model.matrix(formula)
lambda <- exp(-X %*% beta)
y <- rexp(nrow(X), lambda)
return(data.frame(y = y, X))
}
set.seed(123)
n=10
cov <- seq(0, 5, length=n)
dados1 <- simula.exponencial(~cov, c(2,0.5))
dados1
## Função escore
escore <- function(par, formula, dados){
mf <- model.frame(formula, dados)
X <- model.matrix(formula, data=mf)
esco <- ?????????
return(drop(esco))
}
## Hessiana
hessiano <- function(par, formula, dados){
X <- model.matrix(formula, data=dados)
mat <- matrix(0, nrow(X), nrow(X))
diag(mat) <- ???????
H <- ??????????
return(H)
}
Thank you in advance!