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]) }
Add the following code to your website.
REMOVE THIS Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.