R/Choleski_factors.R

Defines functions Choleski_factors

Choleski_factors <-
function(E,rho)
{
  m=nrow(E);
  n=ncol(E);
  if (m>=n)
    U=chol(t(E)%*%E+rho*diag(n)) else 
      U=chol(diag(m)+1/rho*(E%*%t(E)))
  
  L=t(U);
  return(list(L,U))
}
umich-biostatistics/ECLasso documentation built on Jan. 29, 2020, 9:12 p.m.