R/compare_A_Z.R

In CorReg: Linear Regression Based on Linear Structure Between Variables

# ' description de la repartition des A vis a vis de la structure
compare_A_Z <- function(A = A, Z = Z)
{
  if (length(A) == ncol(Z) + 1) {
    A = A[-1]
  }
  quiA = which(A != 0)
  quip2 = which(colSums(Z) != 0)
  p2 = length(quip2)
  nbA = length(quiA)
  nbAp2 = sum(duplicated(c(quiA, quip2)))
  quip1 = which(rowSums(Z) != 0)
  p1 = length(quip1)
  nbAp1 = sum(duplicated(c(quiA, quip1)))
  p3 = ncol(Z) - p1 - p2
  nbAp3 = length(which(A != 0)) - nbAp1 - nbAp2
  
  return(list(p2 = p2, p1 = p1, p3 = p3, nbAp2 = nbAp2, nbAp1 = nbAp1,
              nbAp3 = nbAp3))
}