R/mon_predict_pred2.R

In CorReg: Linear Regression Based on Linear Structure Between Variables

# ' reconstitution du predictif
mon_predict_pred2<-function(mod_expl=mod_expl,mod_inj=mod_inj,newdata=newdata,B=B){
  coef=as.matrix(c(mod_expl$coefficients,mod_inj$coefficients))
  Atilde=mod_expl$coefficients
  quiI1etoile=which(duplicated(c(names(Atilde),names(newdata))))-length(Atilde)
  A2=mod_inj$coefficients
  quiI2etoile=which(duplicated(c(names(A2),names(newdata))))-length(A2)
  Z=B[-1,-1]
  Z=Z-diag(diag(Z))
  Z[Z!=0]=1
  quiI2=which(colSums(Z)!=0)
  quiI2etoile=which(duplicated(c(names(A2),names(newdata[,quiI2]))))-length(A2)
  if(names(mod_expl$coefficients)[1]=="(Intercept)"){
    X=cbind(rep(1,times=nrow(newdata)),as.matrix(newdata[,quiI1etoile]),as.matrix(extract_Xtilde(X=newdata,B=B)[,quiI2etoile]))
  }else{
    X=cbind(as.matrix(newdata[,quiI1etoile]),extract_Xtilde(X=newdata,B=B))
  }
  return(list(hatY=as.matrix(X)%*%as.matrix(coef),hatA=as.matrix(coef)))
}