OpenMx plotting issue:

[andrew0302]

Working on an SEM model with cross-nested random intercepts using OpenMx. The model runs in OpenMx, however, semPaths() throws the following error when I try to display it.

I'm a bit of a noob, so forgive me if this is the wrong venue to report issues like this. Any potential solution?

[emstruong]

I am also having troubles, but this seems to be a different way that semPaths doesn't work with OpenMx.

library(OpenMx)
library(semPlot)

# Instructions from https://brandmaier.github.io/semtree/articles/getting-started.html
# Simulated Data

set.seed(23)
N <- 1000
M <- 5
icept <- rnorm(N, 10, sd = 4)
slope <- rnorm(N, 3, sd = 1.2)
p1 <- sample(c(0, 1), size = N, replace = TRUE)
loadings <- 0:4
x <-
  (slope + p1 * 5) %*% t(loadings) + 
  matrix(rep(icept, each = M), byrow = TRUE, ncol = M) + 
  rnorm(N * M, sd = .08)
growth.data <- data.frame(x, p1)
names(growth.data) <- c(paste0("X", 1:M), "P1")

manifests <- names(growth.data)[1:5]
growthCurveModel <- mxModel("Linear Growth Curve Model Path Specification",
                            type="RAM",
                            manifestVars=manifests,
                            latentVars=c("intercept","slope"),
                            mxData(growth.data, type="raw"),
                            # residual variances
                            mxPath(
                              from=manifests,
                              arrows=2,
                              free=TRUE,
                              values = c(.1, .1, .1, .1, .1),
                              labels=c("residual","residual","residual","residual","residual")
                            ),
                            # latent variances and covariance
                            mxPath(
                              from=c("intercept","slope"),
                              arrows=2,
                              connect="unique.pairs",
                              free=TRUE,
                              values=c(2, 0, 1),
                              labels=c("vari", "cov", "vars")
                            ),
                            # intercept loadings
                            mxPath(
                              from="intercept",
                              to=manifests,
                              arrows=1,
                              free=FALSE,
                              values=c(1, 1, 1, 1, 1)
                            ),
                            # slope loadings
                            mxPath(
                              from="slope",
                              to=manifests,
                              arrows=1,
                              free=FALSE,
                              values=c(0, 1, 2, 3, 4)
                            ),
                            # manifest means
                            mxPath(
                              from="one",
                              to=manifests,
                              arrows=1,
                              free=FALSE,
                              values=c(0, 0, 0, 0, 0)
                            ),
                            # latent means
                            mxPath(
                              from="one",
                              to=c("intercept", "slope"),
                              arrows=1,
                              free=TRUE,
                              values=c(1, 1),
                              labels=c("meani", "means")
                            )
) # close model

GrowthCurveRun <- mxRun(growthCurveModel)
#> Running Linear Growth Curve Model Path Specification with 6 parameters

semPaths(GrowthCurveRun)
#> Error in MeanStd[v.idx[v]] <- stdM$Std.Value[stdM$col == v.idx[v]]: replacement has length zero

^{Created on 2023-05-30 with reprex v2.0.2}