R/temp.resid.R

In multicon: Multivariate Constructs

temp.resid <-
function(template, y.set, nomiss=.8) {
  x.set <- matrix(template, ncol=length(template), nrow=nrow(y.set), byrow=T)
  coefs <- Profile.reg(x.set=x.set, y.set=y.set, center="none", nomiss=nomiss)
  res <- data.frame(y.set - (coefs[,1] + coefs[,2]*x.set))
  names(res) <- paste(colnames(y.set), sep="", ".res")
  return(res)
}