Schuster_update: Computes the Newton-Raphson update
Description
Computes both gradient and hessian, and then solves the system of equations
Usage
Schuster_update(n, marginal_pi, kappa, v)
Value
the vector of updates
Arguments
- n
matrix of observed counts
- marginal_pi
expected proportions for each category
- kappa
current value of kappa coefficient
- v
symmetry matrix