Nothing
marg.likelihood <- function(penden.env,pen.likelihood,temp=FALSE) {
help <- eigen(get("DDD.sum",penden.env))
index <- which(help$values>1e-08)
diag.help2 <- diag(help$values[index])
if(is.null(get("ind.val",penden.env))) {
Utilde <- get("Utilde.eigen.pen.mat",penden.env)
evalues <- help$values[index]
t.Utilde <- t(Utilde)
}
if(!is.null(get("ind.val",penden.env))) {
ind.val<-get("ind.val",penden.env)
pen.mat <- get("DDD.sum",penden.env)[-ind.val,-ind.val]
help2 <- eigen(pen.mat)
index <- which(help2$values>1e-16)
Utilde <- help2$vectors[,index]
evalues <- help2$values[index]
t.Utilde <- t(Utilde)
}
if(!temp) k1 <- 0.5*sum(log(get("lambda",penden.env)[1]*evalues))
if(temp) k1 <- 0.5*sum(log(get("lambda.temp",penden.env)[1]*evalues))
k2 <- pen.likelihood
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
ind<-which(eneu<0&eneu>-1e-8)
eneu<-eneu[-ind]
if(is.complex(eneu)) eneu <- Re(eneu)
if(any(eneu<=0)) print("eneu <=0")
k3 <- -0.5*sum(log(eneu))
#print(c(k1,k2,k3))
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.