R/getExpectedCov.R

In brandmaier/semper: Structural Equation Modeling, power, effect size, and reliability

#
# obtain expected covariance matrix from a sempower model
#
#'
#' @export

getExpectedCovariance <- function(model) {


  if (inherits(model,"semper")) {  
   omx <- getOpenMxRepresentation(model)
  } else if (inherits(model,"MxRAMModel")) {
    omx <- model
  } else {
   stop("Unknown model type!")
  }

 # data <- simulateData(omx,10)
  # orsome fake data
  data <- data.frame(matrix(rnorm(100*length(omx$manifestVars)),
                            nrow = 100,ncol=length(omx$manifestVars)))
  names(data) <- omx$manifestVars
  
  
  
  omx <- mxModel(omx, mxData(observed = data,type="raw"))
  omx <- omxSetParameters(omx, labels=names(omxGetParameters(omx)), free=FALSE)
  run <- mxRun(omx,silent = TRUE)
  sm <- summary(run)
  dataMat <- attr(run$fitfunction$result,"expCov")
  return(dataMat)
}