R/compare_zero.R

In CorReg: Linear Regression Based on Linear Structure Between Variables

# ' compare 0 values in two vectors
# ' @param trueA first vector
# ' @param Aalgo second vector
# ' @param taux boolean. Computes the ratio of each statistic or not.
compare_zero<-function(trueA=trueA,Aalgo=Aalgo,taux=FALSE){
  quivrai0=which(trueA==0)
  if(length(quivrai0)>0){
    nbbon0=length(which(Aalgo[quivrai0]==0))
    nbbon1=length(which(Aalgo[-quivrai0]!=0))
    nbfaux0=length(which(Aalgo[-quivrai0]==0))
    nb0mank=length(quivrai0)-nbbon0
  }else{
    nbbon0=0
    nbbon1=length(which(Aalgo!=0))
    nbfaux0=length(which(Aalgo==0))
    nb0mank=0
  }
  if(taux){
    tauxbon0=nbbon0/length(quivrai0)
    tauxfaux0=nbfaux0/length(which(Aalgo==0))
    return(list(true0=nbbon0,true1=nbbon1,false0=nbfaux0,false1=nb0mank,ratio_true0=tauxbon0,ratio_false0=tauxfaux0))
  }else{
    return(list(true0=nbbon0,true1=nbbon1,false0=nbfaux0,false1=nb0mank))    
  }
}