Copyright © 2007-2024 The OpenMx Project
Wiki home page [1]Ideas and example functions that extend OpenMx, encapsulate tedious work, and make scripts easier to write or more compact.
You will probably define helper functions, especially for summarising the output of model you use frequently.
If you have questions not answers, then add those here: That's how a wiki works. Please add material here as you learn...
source [2]
readLowerTriangle <- function(file, nrows, fill=TRUE) { xvector <- scan(file) X <- matrix(NA, nrows, nrows) i <- 1 for(row in 1:nrows) { for(col in 1:nrows) { if(col>row) next X[row,col] <- xvector[i] i <- i + 1 if (fill) X[col,row] <- X[row,col] } } return(X) }
See also read.moments() in http://cran.r-project.org/web/packages/sem/sem.pdf
require(sem) # install.packages("sem", dep=T)
read.moments(file = "", diag = TRUE,
names = as.character(paste("X", 1:n, sep = "")))