Behavioral Genetics Models

Hi
Is there any way in OpenMx or R to get a graphical representation of a multivariate common pathways heterogeneity model, with standardized path coefficients? Running the omxGraphviz produces the error message that the model is not a RAM type model.

Is there any way to get the 'adjusted BIC' statistic output in OpenMx?

Hi ,

I am running ACE and nested submodels on 100 twin pairs ( 50 MZ and 50 DZ twin pairs) to estimate A,C and E for 10,000 genes. Do I need multiple hypothesis correction on these estimates?

Thanks in advance.

Sabha

Hi,

I'm wondering how to amend an ACE script to allow for different means across zygosity? (means modelling of my phenotype showed poor fits when trying to constrain MZ/DZ means).

Can someone point me in the right direction?

Paul

For anyone interested in running multivariate genetic models with sex as a moderating variable, here is a script
based on the correlated factors approach to ensure that the order of the variables does not affect the ability of the model to account for the DZOS data.
(Ref: Neale et al., Multivariate genetic analysis of sex-lim and G x E interaction, Twin Research & Human Genetics, 2006 ).

The author is dr.Frühling Rijsdijk, who has approved of the sharing of the script.

Hi all,

I am trying to fit an ACE model on a latent variabel using path specification. I started out with a twin model in which each individual has one common factor explaining the variance on 7 observed variables. These factors are obviously correlated within an MZ or DZ pair. I am now removing this correlation between twins and adding a ACE decomposition of the factor variance. Basically i am fitting a common pathway model.

I am not sure if this topic has been covered elsewhere in this forum. How to incorporate age as continuous var into the twin univariate model and calculating the significance of age along with components of A,C ,E and testing the deterioration in model fit after each component is dropped from the full model.

It would be helpful if there any example scripts available to follow.

Thanks in advance.

Sabha

Hi,

I need help with how to estimate the A, C,E for multiple variables and store them in a table format.
I have a example data with two values for each twin pair for 5 genes ( A,B,C,D,E).

I started out in following way from the Twin Univariate model script but it's not working.
Data<- read.table("ExampleData.txt",sep="\t",header=T) # Example Data attached.
Varsi<-c("GeneA_","GeneB_","GeneC_","GeneD_","GeneE_")
for(i in 1:5){
nv <- 1
selVars<- paste(Varsi[i],c(rep(1,nv),rep(2,nv)),sep="")
mzDatai<-subset(Data,Type=="MZ",selVars)
dzDatai<-subset(Data,Type=="DZ",selVars)

I have the following problem.

i estimated a single common factor model in twin pairs:

correlated factors within twinpairs, correlated residuals between measurements done at equal timepoints within twinpairs.

now i want to model the factor means by intoducing a definition variable.

i got the single factor twin model working, which essentially is a two factor model one factor for each twin. Now i need to introduce a definition variable that influences the factor mean. I have found some MX scripts for this but i need it to work in openMX.

Can anyone help me with the function or syntex needed to calculate the 95% CI for A,C,E estimates in the existing
UnivariateTwinAnalysis_MatrixRaw.R script.

Thanks in advance.