Use of constants in formulae supplied as arguments to mxConstraint
> onelocus<-mxModel("onelocus", mxAlgebraObjective("NegativeLogLikelihood"),
mxMatrix("Full", nrow=1, ncol=1, free=TRUE, values=c(.3333), name="P"), # P, freq of allele 1
mxMatrix("Full", nrow=1, ncol=1, free=TRUE, values=c(.3333), name="Q"), # Q, freq of allele 2
mxMatrix("Full", nrow=1, ncol=1, free=TRUE, values=c(.3333), name="R"), # R, freq of allele 3
mxModel, remove = T
example code
library(OpenMx)
thisOMxModel <- mxModel("testModel", mxMatrix(name="Afac", type="Full",nrow=8,ncol=5))
testit <- mxModel(thisOMxModel, mxMatrix(name="Afac", type="Full",nrow=8,ncol=5), remove=T)
# Error in as.character.default(
# no method for coercing this S4 class to a vector
- Read more about mxModel, remove = T
- 2 comments
- Log in or register to post comments
Mx parser for multiple groups
Michael-Neales-Macbook-Pro-2% ~/mx/openmx/OpenMx/trunk/inst/tools/mxAlgebraParser.py < ~/mx/openmxscripts/testalgebraparser.txt > ~/mx/openmxscripts/testalgebraparser.out
Traceback (most recent call last):
File "/Users/neale/mx/openmx/OpenMx/trunk/inst/tools/mxAlgebraParser.py", line 149, in
print(parser.parse(input))
File "build/bdist.macosx-10.3-i386/egg/ply/yacc.py", line 265, in parse
- Read more about Mx parser for multiple groups
- 7 comments
- Log in or register to post comments
Can't seem to get Mx1Parser to work, svn 956
Traceback (most recent call last):
File "./inst/tools/mxAlgebraParser.py", line 45, in
import ply.lex as lex
ImportError: No module named ply.lex
- Read more about Can't seem to get Mx1Parser to work, svn 956
- 3 comments
- Log in or register to post comments
Something funky with sudo make build install OS X after svn up to 956
OpenMx doesn't work with R 2.10.0
- Read more about OpenMx doesn't work with R 2.10.0
- 14 comments
- Log in or register to post comments
crash: when manifests not in data?
I think it is because the dataset in the model doesn't match the manifests list (salthouse vs demoOneFactor), but shouldn't crash....
require(OpenMx)
salthouse <- read.table("http://www.subjectpool.com/data/2009/andrea/salthouse.txt", header=TRUE, sep="\t", as.is=c(TRUE))
manifests <- names(salthouse)
latents <- c("G")
factorModel <- mxModel("One Factor", type="RAM",
manifestVars = manifests,
latentVars = latents,
mxPath(from=latents, to=manifests),
mxPath(from=manifests, arrows=2),
- Read more about crash: when manifests not in data?
- 2 comments
- Log in or register to post comments
replace trac subversion browser
- Read more about replace trac subversion browser
- 2 comments
- Log in or register to post comments
Version 0.2.2-951 is released.
- Read more about Version 0.2.2-951 is released.
- 4 comments
- Log in or register to post comments