R/victory_int.R

In CorReg: Linear Regression Based on Linear Structure Between Variables

victory_int<-function(x=x,y=y,col=c("black","red","green"),alpha=0.2,what=which.min){
   quimin=apply(y,1,what)
   lim <- par("usr")
   minloc=quimin[1]
   xdeb=x[1]
   for(i in 2:length(x)){
      if(quimin[i]!=minloc & i!=length(x)){
         couleur_loc=c(t(col2rgb(col[quimin[i-1]]))/255)
         couleur_loc=rgb(red=couleur_loc[1],green=couleur_loc[2],blue=couleur_loc[3],alpha)
         rect(xdeb, lim[3]-1, (x[i-1]+x[i])/2, lim[4]+1, border = couleur_loc, col = couleur_loc)
         xdeb=(x[i-1]+x[i])/2
         minloc=quimin[i]
      }else if(quimin[i]!=minloc & i==length(x)){#changement et on est au dernier, on ferme le precedent et on ajoute le rectange final
         couleur_loc=c(t(col2rgb(col[quimin[i-1]]))/255)
         couleur_loc=rgb(red=couleur_loc[1],green=couleur_loc[2],blue=couleur_loc[3],alpha)
         rect(xdeb, lim[3]-1, (x[i-1]+x[i])/2, lim[4]+1, border = couleur_loc, col = couleur_loc)
         xdeb=(x[i-1]+x[i])/2
         minloc=quimin[i]
         couleur_loc=c(t(col2rgb(col[quimin[i]]))/255)
         couleur_loc=rgb(red=couleur_loc[1],green=couleur_loc[2],blue=couleur_loc[3],alpha)
         rect(xdeb, lim[3]-1, x[i], lim[4]+1, border = couleur_loc, col = couleur_loc)
      }else if(i==length(x)){#fini donc on trace
         couleur_loc=c(t(col2rgb(col[quimin[i]]))/255)
         couleur_loc=rgb(red=couleur_loc[1],green=couleur_loc[2],blue=couleur_loc[3],alpha)
         rect(xdeb, lim[3]-1, x[i], lim[4]+1, border = couleur_loc, col = couleur_loc)
      }
   }
}