Nothing
marg.likelihood <- function(penden.env,pen.likelihood,temp=FALSE) {
#help <- eigen(get("DDD.sum",penden.env),symmetric=TRUE)
help <- get("eigen.pen.mat",penden.env)
index <- get("index.eigen.pen.mat",penden.env)
diag.help2 <- Matrix(diag(help$values[index]),sparse=FALSE,doDiag=TRUE)
Utilde <- Matrix(get("Utilde.eigen.pen.mat",penden.env))
evalues <- help$values[index]
t.Utilde <- t(Utilde)
k1 <- 0.5*sum(log(get("lambda",penden.env)*evalues))
k2 <- pen.likelihood
#browser()
if(!temp) eneu <- eigen(t(Utilde)%*%-get("Derv2.pen",penden.env)%*%Utilde)$values
else eneu <- eigen(t(Utilde)%*%-get("Derv2.pen.temp",penden.env)%*%Utilde)$values
if(is.complex(eneu)) eneu <- Re(eneu)
if(any(eneu<=0)) print("eneu <=0")
k3 <- -0.5*sum(log(eneu))
if(!temp) assign("marg.log.like",k1+k2+k3,penden.env)
else assign("marg.log.like.temp",k1+k2+k3,penden.env)
}
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.