R/extract_Xtilde.R

In CorReg: Linear Regression Based on Linear Structure Between Variables

# ' extrait Xtilde
extract_Xtilde<-function(X=X,B=B){
  Z=B[-1,-1]
  Z=as.matrix(Z)
  Z=Z-diag(diag(Z))
  n=nrow(X)
  Z[Z!=0]=1
  quiI2=which(colSums(Z)!=0)
  quiI1=(1:ncol(Z))[-quiI2]
  X=as.matrix(X)
  return(X[,quiI2]-cbind(rep(1,times=n),X[,quiI1])%*%B[c(1,quiI1+1),quiI2])
}