R/lav_representation.R

In lavaan: Latent Variable Analysis

# user visible function to add 'matrix' entries in the parameter table
lavMatrixRepresentation <- function(partable, representation = "LISREL",
                                    add.attributes = FALSE,
                                    as.data.frame. = TRUE) {

    # check parameter table
    partable <- lav_partable_complete(partable)

    # get model matrices
    if(representation == "LISREL") {
        REP <- lav_lisrel(partable, target = NULL, extra = add.attributes)
    } else if(representation == "RAM") {
        REP <- lav_ram(partable, target = NULL, extra = add.attributes)
    } else {
        stop("lavaan ERROR: representation must either \"LISREL\" or \"RAM\".")
    }

    partable$mat <- REP$mat
    partable$row <- REP$row
    partable$col <- REP$col

    if(as.data.frame.) {
        partable <- as.data.frame(partable, stringsAsFactors=FALSE)
        class(partable) <- c("lavaan.data.frame", "data.frame")
    }

    if(add.attributes) {
        if(representation == "LISREL") {
            attr(partable, "ov.dummy.names.nox") <- attr(REP, "ov.dummy.names.nox")
            attr(partable, "ov.dummy.names.x")   <- attr(REP, "ov.dummy.names.x")
        } else if(representation == "RAM") {
            attr(partable, "ov.idx") <- attr(REP, "ov.idx")
        }
        attr(partable, "mmNames")  <- attr(REP, "mmNames")
        attr(partable, "mmNumber") <- attr(REP, "mmNumber")
        attr(partable, "mmRows")   <- attr(REP, "mmRows")
        attr(partable, "mmCols")   <- attr(REP, "mmCols")
        attr(partable, "mmDimNames")  <- attr(REP, "mmDimNames")
        attr(partable, "mmSymmetric") <- attr(REP, "mmSymmetric")
    }

    partable
}