Nothing
#predict fucntions for family(exp)
merlin_p_exp_h <- function(gml,t=NULL)
{
if (length(t)==0) t = merlin_util_timevar(gml)
return(exp(merlin_util_xzb(gml,t)))
}
merlin_p_exp_ch <- function(gml,t=NULL)
{
if (length(t)==0) t = merlin_util_timevar(gml)
if (gml$NI[[gml$modelind]]) {
N = gml$Nobs[gml$modtouse]
gq = merlin_get_GK()
qp = sweep((t %*% (gq[[1]]/2)),1,(t/2),"+")
if (gml$Nlevels & !gml$fixedonly) {
result = matrix(0,nrow=N,ncol=gml$Ndim[gml$Nlevels])
}
else result = matrix(0,nrow=N,ncol=1)
for (q in 1:15) {
result = result + sweep(merlin_p_exp_h(gml,qp[,q]),1,gq[[2]][q] * t/2,"*")
}
return(result)
}
else return(sweep(exp(merlin_util_xzb(gml,t)),1,t,"*"))
}
merlin_p_exp_logch <- function(gml,t=NULL)
{
if (length(t)==0) t = merlin_util_timevar(gml)
return(log(merlin_p_exp_ch(gml,t)))
}
merlin_p_exp_s <- function(gml,t=NULL)
{
if (length(t)==0) t = merlin_util_timevar(gml)
return(exp(-merlin_p_exp_ch(gml,t)))
}
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.