R/marg.likelihood.r

In penDvine: Flexible Pair-Copula Estimation in D-Vines using Bivariate Penalized Splines

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)
}