R/count_2triangle.R

In ericstrobl/CCI: Cyclic Causal Inference

count_2triangle <- function(pag,a,b,c){
  
  indD = which( (pag[a,]==1 | pag[a,]==2) & (pag[b,]==1 | pag[b,]==3) &
                  (pag[,b]==1 | pag[,b] ==3) )
  
  indD = setdiff(indD,c);
  
  indDf=c();
  
  for(iD in seq_len(length(indD))){
    
    if ( length(setdiff(which(pag[,c]!=0),b)>1) ){
      
      if (isPossible_undirected_path_fast(pag,indD[iD],c,b)){
        indDf = c(indDf, indD[iD]);
      }
      
      
    }
  }
  
  return(indDf)
  
}