myee <- function(data){
function(theta){
c(data$Y1 - theta[1],
(data$Y1 - theta[1])^2 - theta[2])
}
}
results <- m_estimate(
estFUN = myee,
data = geexex,
root_control = setup_root_control(start = c(1,1)))
head(grab_meat_list(results@sandwich_components))
Run the code above in your browser using DataLab