R/TV.R

Defines functions TV

TV<-function(z,S){
  if(length(z)==length(which(z==0))){ return(0) }
  kik=which(z!=0)
  return(t(z[kik])%*%ginv(S[kik,kik])%*%z[kik])
}
roqe/T2GA documentation built on March 9, 2021, 9:38 p.m.