OpenMx General Help
parallel computation in linux cluster
I want to conduct several 10,000-simulations for an OpenMx model in Linux cluster. I have over 20 parameters in the model and for each parameter, I need 1000-bootstrap confidence interval to test the coverage. I've add lines in below screenshot, and there are 16 cores available for OpenMx. Even though, the simulation may take about one month at my best guess. Could someone kindly advise some other ways to speed up my simulations? Thanks in advance.
- Read more about parallel computation in linux cluster
- 4 comments
- Log in or register to post comments
Create a function to generate mxMatrix and mxAlgebra
I am writing an R function which can create and run an OpenMx model. I can do it in a simple way. Partial codes are shown below. However, I have over 20 definition variables and for each I need corresponding mxAlgebra. I may want to write a function/a loop to generate mxMatrix and mxAlgebra. I've tried a couple of methods, yet none of them worked. The mainly problem is that in the expression of mxAlgebra, I need call the name of corresponding mxMatrix. Thank you in advance!
Performance issues many algebras
I am currently implementing a new model specification language for longitudinal panel models. As "backend", I use OpenMx. Everything works but it is painfully slow.
- Read more about Performance issues many algebras
- 21 comments
- Log in or register to post comments
Single Common Factor Model with Some Residual Errors Correlated
Power estimation for the detection of rG and rE
Specifically, the analysis employed a trivariate Cholesky (AE providing the best fit, with all significant rGs, and one significant rE). The CIs have been calculated for all estimates. The sample is on the small side: 200 same-sex pairs (half MZ, half DZ) and may have been underpowered for the other smaller rEs but I’m not sure what size effect I had enough power to detect.
How to set the value of a computed parameter to that computed value?
Background and use-case: OpenMx allows us to use labels to set the value of a parameter (matrix cell). Here's a RAM model where I insert the value of covXY into the row 2, column1 of the symmetric path (S) matrix:
m1$S
$labels
x y
x "x_with_x" "covXY[1,1]"
y "covXY[1,1]" "y_with_y"
$values
x y
x 1.198632 0.000000
y 0.000000 2.018483
Rmpi
Error in { : task 18 failed - "job.num is at least 2."
Calls: %dopar% ->
Execution halted
--------------------------------------------------------------------------
mpirun has exited due to process rank 0 with PID 1077 on
node compute-0-11.local exiting improperly. There are two reasons this could occur:
1. this process did not call "init" before exiting, but others in
the job did. This can cause a job to hang indefinitely while it waits
- Read more about Rmpi
- 6 comments
- Log in or register to post comments
Standard Errors versus Confidence Intervals
If I treat the ordinal data as if it were quantitative, I have no trouble fitting a common factor model.
- Read more about Standard Errors versus Confidence Intervals
- 6 comments
- Log in or register to post comments
How to use "mxFitFunctionAlgebra" instead of "mxAlgebraObjective"
###############
models_1_v4 <- mxModel("Models_1",
Model_1_t1_v4,
Model_1_t2_v4,
mxAlgebra(
Model_1_t1_v4.objective+
Model_1_t2_v4.objective,
name="multi"),
# mxAlgebraObjective("multi"),
mxFitFunctionAlgebra(algebra = "multi")
)
m1_fit_v4 <- mxRun(models_1_v4)
##################
Output from mxFactorScores - Order of rows?
The rows are in the order of the _sorted_ data.
I'm not sure what this means. The data I have is in order by date and it's clear that the estimated factor scores for "ML" are no longer in the same order. There are two correlated factors each with 2 indicators (and nothing else). How can I get the factor scores into the same order as the rest of the data?
- Read more about Output from mxFactorScores - Order of rows?
- 5 comments
- Log in or register to post comments
Pagination
- Previous page
- Page 7
- Next page