Newton-Raphson method: Logistic regression (R)

Question

Have an exercise to write three function: the score function S<- function(theta, y, X){...} and the Fisher information matrix I <- function(theta,