OpenMx Structural Equation Modeling
cffalk2
Joined: 10/06/2025
mxComputePenaltySearch
Happy to see this implemented (thank you!). Questions about how it works. I've seen the vignette here: https://cran.r-project.org/web/packages/OpenMx/vignettes/regularization.html
But, that includes use of mxPenaltySearch. The one above (but not mxPenaltySearch) has an option to change ebicGamma, which is desired. I don't see documentation or examples for mxComputePenaltySearch, but would guess it could form part of a custom compute plan.
- Read more about mxComputePenaltySearch
- 1 comment
- Log in or register to post comments
cffalk2
Joined: 10/06/2025
mxComputePenaltySearch
Happy to see this implemented (thank you!). Questions about how it works. I've seen the vignette here: https://cran.r-project.org/web/packages/OpenMx/vignettes/regularization.html
But, that includes use of mxPenaltySearch. The one above (but not mxPenaltySearch) has an option to change ebicGamma, which is desired. I don't see documentation or examples for mxComputePenaltySearch, but would guess it could form part of a custom compute plan.
- Read more about mxComputePenaltySearch
- Log in or register to post comments
cffalk2
Joined: 10/06/2025
mxComputePenaltySearch
Happy to see this implemented (thank you!). Questions about how it works. I've seen the vignette here: https://cran.r-project.org/web/packages/OpenMx/vignettes/regularization.html
But, that includes use of mxPenaltySearch. The one above (but not mxPenaltySearch) has an option to change ebicGamma, which is desired. I don't see documentation or examples for mxComputePenaltySearch, but would guess it could form part of a custom compute plan.
- Read more about mxComputePenaltySearch
- Log in or register to post comments
Manuel Rein
Joined: 04/01/2024
Vivian
Joined: 04/30/2025
IvanVoronin
Joined: 08/18/2013
Mixed effects model - matrix specification
Hello,
I am looking for advice with regard to the implementation of mixed effects model in OpenMx.
I am starting with a simple model where one dependent variable Y is predicted by four correlated variables X1-X4. The participants in the sample are clustered in pairs (they are twins).
- Read more about Mixed effects model - matrix specification
- 3 comments
- Log in or register to post comments
Ilaria
Joined: 07/22/2024
Problems with Constrain expected Means and Variances to be equal across Twin Order
Hi all,
I am trying to fit several bivariate [between my exposure at age 16 and my outcomes at age 21 (accounting for age and sex)] and trivariate [my exposure at age 16 and my outcomes at ages 21 and 26 (also accounting for age and sex)] ACE models. However, I’ve started by fitting univariate saturated and submodels, followed by ACE and submodels, to eventually fit my bivariate models. The saturated models run fine across all, but when I try to fit the constrained models (by twin order and zygosity), I get the following error message:
didenursahin
Joined: 03/20/2023
outdated-solved
outdated-solved
- Read more about outdated-solved
- Log in or register to post comments
Tugce Yildiz
Joined: 01/15/2024
NA CIs for Bivariate Model
Hi, I am trying to run a bivariate model with two variables and 2 covariates. However, I get NA !!! for BivA, BivC and BivE. Can you help me solve this issue?
Here is the script:
### SELECT VARIABLES HERE
data_file <- "Cerebellum_volumes.csv"
variable_input1 <- "lsas_total"
variable_input2 <- "Vermis_VII"
cov1 <- "Total_Cerebel_Vol"
cov2 <- "age"
### Import data
dat <- read.csv(file=data_file, header=TRUE, sep = ";")
- Read more about NA CIs for Bivariate Model
- 1 comment
- Log in or register to post comments
lior abramson
Joined: 07/21/2017
checking total h2 at time 2 in a longitudinal choelsky model with paths script
Hello,
I am conducting a longitudinal Cholesky model with one variable in two time points.
I use the script written below with paths.
How can I modify this script to find for the total h2 at time 2? that is, I want to find the heritability explained by both the genetic effect from time 1 and the genetic effect from time 2, and the confidence interval of this estimate.
Thank you very much for your help