R/residuals.R

In jslefche/piecewiseSEM: Piecewise Structural Equation Modeling

#' Residual values from fit models
#'
#' @param object a \code{\link{psem}} object
#' @param ... additional arguments to residuals
#' 
#' @return a \code{data.frame} of residuals for endogenous variables as columns
#' 
#' @export
#' 
residuals.psem <- function(object, ...) {
  
  object <- removeData(object, formulas = 1)
  
  #get residuals of relationships
  ret <- lapply(object, residuals, ...)
  
  #get column names
  formulaList <- listFormula(object)
  
  resp <- sapply(formulaList, function(x) all.vars_notrans(x)[1])
  
  names(ret) <- paste(resp, "residuals", sep="_")
  
  do.call(cbind, ret)
  
}