Problem creating saturated output from script: MultivariateTwinAnalysis_MatrixRaw.R
Posted on
trinewa
Joined: 11/30/2009
Forums
I am running the script: MultivariateTwinAnalysis_MatrixRaw.R from the TC24 homepage, on my own data.
I don't seem to be able to create any saturated output from the following commands, although the multiTwinSatFit is not empty:
parameterSpecifications(multiTwinSatFit)
expectedMeansCovariances(multiTwinSatFit)
tableFitStatistics(multiTwinSatFit)
I get the following error messages:
Error: could not find function "parameterSpecifications"
Error: could not find function "expectedMeansCovariances"
Error: could not find function "tableFitStatistics"
I am just beginning to learn to run twin analyyses through OpenMx - very happy if anyone could tell me what I am doing wrong.
I think you need to source
I think you need to source the Genetic Epi Helper Functions: http://www.vipbg.vcu.edu/~vipbg/Tc24/GenEpiHelperFunctions.R. The link is at the bottom right of the TC24 homepage. But I admit there's no instructions that you need to include this file.
Log in or register to post comments
In reply to I think you need to source by mspiegel
Thank you, that helped! One
Thank you, that helped!
One more question concerning this script:
The parts of the script that concern fitting the multivariate saturated model, and fitting multivariate model with equal means & variances across twin order and zygosity, fail when running the script on my data.
In the first I get the following error message:
Running multiTwinSat
Error in FUN(X[[3L]], ...) :
Could not find function with name : and 2 arguments,
and in the second the error message is tha object "uniTwinSatFit" cannot be found.
I have edited the script by setting correct starting values. Are there any other adaptations that need to be made in order to run the script on my own data?
Log in or register to post comments
In reply to Thank you, that helped! One by trinewa
I think upgrading to the
I think upgrading to the latest version of OpenMx will make the first error go away. The ":" operator for matrix selection was implemented in the past month.
Log in or register to post comments
In reply to I think upgrading to the by mspiegel
Thanks once again for useful
Thanks once again for useful help! Ugrading OpenMx fixed the first problem. However, I discovered an error in the example script, in that it calls for a non-existing object uniTwinSatFit that needs to be changed to multiTwinSatFit.
Log in or register to post comments