OpenMx General Help

fyi,

got the following error message after mxRun():

Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
BLAS/LAPACK routine 'DPOTRF' gave error code -4

to save y'all some time, turns out that i was running an old script that had an MxAlgebra object with the expression ".5 * VA" where VA was a matrix. to fix it, i just had to change this to ".5 %x% VA."

suggest you trap the DPOTRF error and return a more intelligible error message.

greg

the default likelihood function for observed data of type=cov or cor multiplies
log(det(predictedCov)) + trace(observedCov %*% inverse(predictedCov))
by (numObs - 1). my understanding is that both the derivations from the multivariate normal pdf and the wishart pdf lead to this being multiplied by numObs and not (numObs - 1). (note that math/stat treatment of the wishart talks of "degrees of freedom" but the actual algebra defines this quantity as numObs).

i stand to be corrected here.
greg

Hi

I am new to OpenMx and thought that I would try things out by working through some of the examples.

Does anyone know where I can get the dataset myTwinData.txt that is used in some of them?

Thankyou

Karin

anyone explain this?

> model2 <- mxModel(model1, mxAlgebra(.5*CTf, name="CTDZf"))
> model2 <- mxRun(model2)
Running SameMZDZf
Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
BLAS/LAPACK routine 'DPOTRF' gave error code -4

> model2 <- mxModel(model1, mxAlgebra(.5 %x% CTf, name="CTDZf"))
> model2 <- mxRun(model2)
Running SameMZDZf
Warning message:
In model 'SameMZDZf' NPSOL returned a non-zero status code 1. blah blah ...

yet
> mat <- matrix(rnorm(16), 4)
> .5 * mat == .5 %x% mat
[,1] [,2] [,3] [,4]
[1,] TRUE TRUE TRUE TRUE

Does anybody know how can we download the source codes of OpenMx. Looks like no links to the source codes from the OpenMx Website.

Hello,