Uses the "almost classical lambda" (=matrix where the mean of each row is equal to the classical lambda) to create a robust lambda by using M estimation.
determine_robust_lambda(
almost_classical_lambda,
fastoption = TRUE,
fastoption2 = FALSE
)M-estimator of location of the parameter, by minimizing sum of rho()
matrix where the mean of each row is equal to the classical lambda
Uses nlm() instead of optim(). This is faster.
experimental parameter: can speed nlm() up (10%), but loses accuracy. May benefit from finetuning.