Nothing
################### the main gel function ###################
########## The functions is to avoid having to build model objects
gel4 <- function (g, x=NULL, theta0=NULL,lambda0=NULL, getVcov=FALSE,
gelType = c("EL","ET","EEL","HD", "REEL","ETEL","ETHD"),
vcov = c("MDS","iid","HAC"), grad=NULL,
vcovOptions=list(), centeredVcov = TRUE,
cstLHS=NULL, cstRHS=NULL, lamSlv=NULL,
rhoFct=NULL, initTheta=c("gmm", "theta0"),
data = parent.frame(),
coefSlv=c("optim","nlminb","constrOptim"),
lControl=list(), tControl=list())
{
Call <- match.call()
vcov <- match.arg(vcov)
gelType <- match.arg(gelType)
coefSlv <- match.arg(coefSlv)
initTheta <- match.arg(initTheta)
if (initTheta=="theta0" & is.null(theta0))
stop("theta0 is required when initTheta='theta0'")
model <- gelModel(g, x, gelType, rhoFct, theta0, grad, vcov, vcovOptions,
centeredVcov, data=data)
if (initTheta == "theta0")
{
initTheta <- "modelTheta0"
names(theta0) = model@parNames
} else {
theta0 <- NULL
}
if (!is.null(cstLHS))
{
model <- restModel(model, cstLHS, cstRHS)
spec <- modelDims(model)
if (!is.null(theta0))
theta0 <- theta0[(names(theta0) %in% spec$parNames)]
}
fit <- modelFit(model=model, initTheta=initTheta, theta0=theta0,
lambda0=lambda0, vcov=getVcov, coefSlv=coefSlv,
lamSlv=lamSlv, tControl=tControl, lControl=lControl)
fit@call <- Call
fit
}
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.