Behavioral Genetics Models
Why mxTryHardOrdinal()dosen't report start values?
When i use
SatFit <- mxTryHardOrdinal(SatModel, intervals=TRUE,extraTries = 15,finetuneGradient=FALSE,bestInitsOutput=TRUE)
,the result is Solution found! Final fit=22505.957 (started at 22819.562) (16 attempt(s): 16 valid, 0 errors)
. I don't know why it dosen't report start values from the best fit.saturated univariate SEM model of a binary variable
I have some problems of fitting saturated univariate SEM model of a binary variable "whtr_g". I would appreciate it if anyone could help me.
1. I always get warning message when fitting saturated SEM model(see warning message.png), whether using SLSQP, NPSOL or CSOLNP optimizer. What’ more, the result of three optimizers are different. But I don’t know if there are errors in my code. Please help me to check it.
Moderation model produces weird results
- Read more about Moderation model produces weird results
- 36 comments
- Log in or register to post comments
How to calculate bivariate h2(or c2, e2) in bivaraite SEM model with two binary traits?
I have encountered several problems in doing bivariate SEM models, which is an important part of my graduation project.The main problems are as follows:
G*E interaction with a moderator
Also, is it appropriate to see the unstandardized variance components by SES moderation, rather than the standardized variance components?
- Read more about G*E interaction with a moderator
- 4 comments
- Log in or register to post comments
Nonlinear transformation and estimates
Thank you in advance!
Julia
- Read more about Nonlinear transformation and estimates
- 3 comments
- Log in or register to post comments
Question about Confidence Interval for Standardized path coefficients (multivariate model)
Hi everyone, I am new here and learning OpenMX. I have a problem of getting confidence interval for standardized path coefficients. Could anyone please provide me the OpenMx language about this? I also provide my scripts below and hope someone can let me know what the problem is (the script about 'ciACE <- mxCI("VC[1:5, 16:30]")' ) or how to read the output about standardized CI.
# Load Data
nl <- read.csv(file="C:/R/TwinPairACE.csv", header=TRUE, sep=",")
Fit tests for ACE models
I am running an ACE model for an ordered categorical outcome using script adapted from [Maes](http://ibg.colorado.edu/cdrom2016/maes/UnivariateAnalysis/one/oneACEo.R). I seem to be able to generate most indicators of model fit, including AIC, BIC, -2LL and p-values associated with likelihood-ratio tests for nested models. However, I cannot seem to generate an associated chi-square value. Is there a simple way to generate chi-square value for these tests? Thanks,
Chance
- Read more about Fit tests for ACE models
- 2 comments
- Log in or register to post comments
Pagination
- Previous page
- Page 11
- Next page