Nothing
## -----------------------------------------------------------------------------
library(geex)
library(dplyr)
myee <- function(data){
Y1 <- data$Y1
function(theta){
c(Y1 - theta[1],
(Y1 - theta[1])^2 - theta[2])
}
}
## -----------------------------------------------------------------------------
mybasis <- new("m_estimation_basis",
.estFUN = myee,
.data = geexex)
## -----------------------------------------------------------------------------
slotNames(mybasis)
## -----------------------------------------------------------------------------
mybasis@.psiFUN_list[1:2]
## ---- eval=FALSE--------------------------------------------------------------
# m <- nrow(geexex)
# lapply(split(geexex, f = 1:m), function(O_i){
# myee(O_i)
# })
## -----------------------------------------------------------------------------
mybasis@.GFUN
## -----------------------------------------------------------------------------
rootSolve::multiroot(
f = mybasis@.GFUN,
start = c(0, 0))
## ----setup--------------------------------------------------------------------
mycontrol <- new('geex_control', .root = setup_root_control(start = c(1, 1)))
mybasis@.control <- mycontrol
roots <- mybasis %>%
estimate_GFUN_roots()
roots
## -----------------------------------------------------------------------------
-numDeriv::jacobian(func = mybasis@.psiFUN_list[[1]], x = roots$root)
## -----------------------------------------------------------------------------
mats <- mybasis %>%
estimate_sandwich_matrices(.theta = roots$root)
# Compare to the numDeriv computation above
grab_bread_list(mats)[[1]]
## -----------------------------------------------------------------------------
mats %>%
{compute_sigma(A = grab_bread(.), B = grab_meat(.))}
## -----------------------------------------------------------------------------
m_estimate(
estFUN = myee,
data = geexex,
root_control = setup_root_control(start = c(0, 0))
)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.