copy files
DESCRIPTION
Copy files. Go to Faculty \marleen\Boulder2012\ Multivariate Copy all files to your own directory Go to Faculty \kees\Boulder2012\ Multivariate Copy all files to your own directory. Introduction to Multivariate Genetic Analysis (1) Marleen de Moor, Kees-Jan Kan & Nick Martin. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/1.jpg)
March 7, 2012 M. de Moor, Twin Workshop Boulder 1
Copy files
• Go to Faculty\marleen\Boulder2012\Multivariate• Copy all files to your own directory
• Go to Faculty\kees\Boulder2012\Multivariate• Copy all files to your own directory
![Page 2: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/2.jpg)
Introduction to Multivariate
Genetic Analysis (1)
Marleen de Moor, Kees-Jan Kan & Nick Martin
March 7, 2012 2M. de Moor, Twin Workshop Boulder
![Page 3: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/3.jpg)
March 7, 2012 M. de Moor, Twin Workshop Boulder 3
Outline
• 11.00-12.30– Lecture Bivariate Cholesky Decomposition– Practical Bivariate analysis of IQ and attention problems
• 12.30-13.30 LUNCH• 13.30-15.00
– Lecture Multivariate Cholesky Decomposition– PCA versus Cholesky– Practical Tri- and Four-variate analysis of IQ, educational
attainment and attention problems
![Page 4: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/4.jpg)
March 7, 2012 M. de Moor, Twin Workshop Boulder 4
Outline
• 11.00-12.30– Lecture Bivariate Cholesky Decomposition– Practical Bivariate analysis of IQ and attention problems
• 12.30-13.30 LUNCH• 13.30-15.00
– Lecture Multivariate Cholesky Decomposition– PCA versus Cholesky– Practical Tri- and Four-variate analysis of IQ, educational
attainment and attention problems
![Page 5: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/5.jpg)
Aim / Rationale multivariate models
Aim:To examine the source of factors that make traits correlate or co-vary
Rationale:• Traits may be correlated due to shared genetic factors (A or D) or shared environmental factors (C or E)• Can use information on multiple traits from twin pairs to partition covariation into genetic and environmental components
March 7, 2012 M. de Moor, Twin Workshop Boulder 5
![Page 6: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/6.jpg)
Example
• Interested in relationship between ADHD and IQ
• How can we explain the association– Additive genetic factors (rG)
– Common environment (rC)
– Unique environment (rE)
March 7, 2012 M. de Moor, Twin Workshop Boulder 6
ADHD
E1 E2
A1 A2
a11
IQ
rG
C1 C2
rC
1 1 1 1
1 1
rE
c11 a22 c22
e11 e22
![Page 7: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/7.jpg)
March 7, 2012 M. de Moor, Twin Workshop Boulder 7
![Page 8: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/8.jpg)
March 7, 2012 M. de Moor, Twin Workshop Boulder 8
Bivariate Cholesky
Multivariate Cholesky
![Page 9: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/9.jpg)
Sources of information
• Two traits measured in twin pairs• Interested in:
– Cross-trait covariance within individuals = phenotypic covariance
– Cross-trait covariance between twins = cross-trait cross-twin covariance
– MZ:DZ ratio of cross-trait covariance between twins
March 7, 2012 M. de Moor, Twin Workshop Boulder 9
![Page 10: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/10.jpg)
Observed Covariance Matrix: 4x4
Phenotype 1 Phenotype 2 Phenotype 1 Phenotype 2
Phenotype 1 Variance P1
Phenotype 2 Covariance P1-P2
Variance P2
Phenotype 1 Within-traitP1
Cross-trait Variance P1
Phenotype 2 Cross-trait Within-traitP2
Covariance P1-P2
Variance P2
Twin 1
Twin
1Tw
i n 2
Twin 2
![Page 11: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/11.jpg)
Observed Covariance Matrix: 4x4
Phenotype 1 Phenotype 2 Phenotype 1 Phenotype 2
Phenotype 1 Variance P1
Phenotype 2 Covariance P1-P2
Variance P2
Phenotype 1 Within-traitP1
Cross-trait Variance P1
Phenotype 2 Cross-trait Within-traitP2
Covariance P1-P2
Variance P2
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covariance
![Page 12: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/12.jpg)
Observed Covariance Matrix: 4x4
Phenotype 1 Phenotype 2 Phenotype 1 Phenotype 2
Phenotype 1 Variance P1
Phenotype 2 Covariance P1-P2
Variance P2
Phenotype 1 Within-traitP1
Cross-trait Variance P1
Phenotype 2 Cross-trait Within-traitP2
Covariance P1-P2
Variance P2
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covarianceCross-twin covariance
![Page 13: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/13.jpg)
Cholesky decomposition
March 7, 2012 M. de Moor, Twin Workshop Boulder 13
Twin 1Phenotype 1
A1
E1
a11 a21
e11 e21
1
1
C1
c11
c21
A2
E2
a22
e22
1
1
Twin 1Phenotype 2
C2
c22
11
Twin 2Phenotype 1
A1 A2
E1 E2
a11 a21 a22
e11 e21 e22
1
1
1
1
Twin 2Phenotype 2
C1 C2
c11
c21
c22
11
1/.5
1
1/.5
1
![Page 14: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/14.jpg)
Now let’s do the path tracing!
March 7, 2012 M. de Moor, Twin Workshop Boulder 14
Twin 1Phenotype 1
A1 A2
E1 E2
a11 a21 a22
e11 e21 e22
1
1
1
1
Twin 1Phenotype 2
C1 C2
c11
c21
c22
11
Twin 2Phenotype 1
A1 A2
E1 E2
a11 a21 a22
e11 e21 e22
1
1
1
1
Twin 2Phenotype 2
C1 C2
c11
c21
c22
11
1/.5 1/.5
11
![Page 15: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/15.jpg)
Within-Twin Covariances (A)
A1 A2
a11 a21 a22
11
Phenotype 1 Phenotype 2
Phenotype 1 a112+c11
2+e112
Phenotype 2a11a21+c11c21+e11e21 a22
2+a212+c22
2+c212+
e222+e21
2
Twin 1
Twin
1
Twin 1Phenotype 1
Twin 1Phenotype 2
![Page 16: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/16.jpg)
Within-Twin Covariances (A)
A1 A2
a11 a21 a22
11
Phenotype 1 Phenotype 2
Phenotype 1 a112+c11
2+e112
Phenotype 2a11a21+c11c21+e11e21 a22
2+a212+c22
2+c212+
e222+e21
2
Twin 1
Twin
1
Twin 1Phenotype 1
Twin 1Phenotype 2
![Page 17: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/17.jpg)
Within-Twin Covariances (A)
A1 A2
a11 a21 a22
11
Phenotype 1 Phenotype 2
Phenotype 1 a112+c11
2+e112
Phenotype 2a11a21+c11c21+e11e21 a22
2+a212+c22
2+c212+
e222+e21
2
Twin 1
Twin
1
Twin 1Phenotype 1
Twin 1Phenotype 2
![Page 18: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/18.jpg)
Within-Twin Covariances (A)
A1 A2
a11 a21 a22
11
Phenotype 1 Phenotype 2
Phenotype 1 a112+c11
2+e112
Phenotype 2a11a21+c11c21+e11e21 a22
2+a212+c22
2+c212+
e222+e21
2
Twin 1
Twin
1
Twin 1Phenotype 1
Twin 1Phenotype 2
![Page 19: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/19.jpg)
Within-Twin Covariances (C)
C1 C2
c11 c21 c22
11
Phenotype 1 Phenotype 2
Phenotype 1 a112+c11
2+e112
Phenotype 2a11a21+c11c21+e11e21 a22
2+a212+c22
2+c212+
e222+e21
2
Twin 1
Twin
1
Twin 1Phenotype 1
Twin 1Phenotype 2
![Page 20: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/20.jpg)
Within-Twin Covariances (E)
Phenotype 1 Phenotype 2
Phenotype 1 a112+c11
2+e112
Phenotype 2a11a21+c11c21+e11e21 a22
2+a212+c22
2+c212+
e222+e21
2
Twin 1
Twin
1
Twin 1Phenotype 1
E1 E2
e11 e21 e22
11
Twin 1Phenotype 2
![Page 21: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/21.jpg)
Cross-Twin Covariances (A)
Twin 1Phenotype 1
A1 A2
a11 a21 a22
11
Twin 1Phenotype 2
Twin 2Phenotype 1
A1 A2
a11 a21 a22
11
Twin 2Phenotype 2
1/0.5 1/0.5
Phenotype 1 Phenotype 2
Phenotype 1
Phenotype 2 +c11c21 +c222+c21
2
Twin 1
Twin
2
![Page 22: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/22.jpg)
Cross-Twin Covariances (A)
Twin 1Phenotype 1
A1 A2
a11 a21 a22
11
Twin 1Phenotype 2
Twin 2Phenotype 1
A1 A2
a11 a21 a22
11
Twin 2Phenotype 2
1/0.5 1/0.5
Phenotype 1 Phenotype 2
Phenotype 1 1/0.5a112+
Phenotype 2 +c11c21 +c222+c21
2
Twin 1
Twin
2
![Page 23: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/23.jpg)
Cross-Twin Covariances (A)
Twin 1Phenotype 1
A1 A2
a11 a21 a22
11
Twin 1Phenotype 2
Twin 2Phenotype 1
A1 A2
a11 a21 a22
11
Twin 2Phenotype 2
1/0.5 1/0.5
Phenotype 1 Phenotype 2
Phenotype 1 1/0.5a112+c11
2
Phenotype 2 1/0.5a11a21+c11c21 +c222+c21
2
Twin 1
Twin
2
![Page 24: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/24.jpg)
Cross-Twin Covariances (A)
Twin 1Phenotype 1
A1 A2
a11 a21 a22
11
Twin 1Phenotype 2
Twin 2Phenotype 1
A1 A2
a11 a21 a22
11
Twin 2Phenotype 2
1/0.5 1/0.5
Phenotype 1 Phenotype 2
Phenotype 1 1/0.5a112+c11
2
Phenotype 2 1/0.5a11a21+c11c21 1/0.5a222+1/0.5a21
2+c222+c21
2
Twin 1
Twin
2
![Page 25: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/25.jpg)
Cross-Twin Covariances (C)
Twin 1Phenotype 1
Twin 1Phenotype 2
C1 C2
c11
c21
c22
11
Twin 2Phenotype 1
Twin 2Phenotype 2
C1 C2
c11
c21
c22
11
1 1
Phenotype 1 Phenotype 2
Phenotype 1 1/0.5a112+c11
2
Phenotype 2 1/0.5a11a21+c11c21 1/0.5a222+1/0.5a21
2+c222+c21
2
Twin 1
Twin
2
![Page 26: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/26.jpg)
Predicted Model
Phenotype 1 Phenotype 2 Phenotype 1 Phenotype 2
Phenotype 1 a112+c11
2+e112
Phenotype 2a11a21+c11c21+e
11e21
a222+a21
2+c222+
c212+e22
2+e212
Phenotype 1 1/.5a112+c11
2 a112+c11
2+e112
Phenotype 21/.5a11a21+
c11c21
1/.5a222+1/.5
a212+c22
2+c212
a11a21+c11c21+e
11e21
a222+a21
2+c222+
c212+e22
2+e21
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covarianceCross-twin covariance
![Page 27: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/27.jpg)
Predicted Model
Phenotype 1 Phenotype 2 Phenotype 1 Phenotype 2
Phenotype 1Variance
P1
Phenotype 2Covariance
P1-P2Variance
P2
Phenotype 1Within-trait
P1Cross-trait Variance
P1
Phenotype 2Cross-trait Within-trait
P2Covariance
P1-P2Variance
P2
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covarianceCross-twin covariance
![Page 28: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/28.jpg)
Phenotype 1 Phenotype 2 Phenotype 1 Phenotype 2
Phenotype 1Variance
P1
Phenotype 2Covariance
P1-P2Variance
P2
Phenotype 1Within-trait
P1Cross-trait Variance
P1
Phenotype 2Cross-trait Within-trait
P2Covariance
P1-P2Variance
P2
Predicted Model
Variance of P1 and P2the same across twinsand zygosity groups
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covarianceCross-twin covariance
![Page 29: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/29.jpg)
Phenotype 1 Phenotype 2 Phenotype 1 Phenotype 2
Phenotype 1Variance
P1
Phenotype 2Covariance
P1-P2Variance
P2
Phenotype 1Within-trait
P1Cross-trait Variance
P1
Phenotype 2Cross-trait Within-trait
P2Covariance
P1-P2Variance
P2
Predicted Model
Covariance of P1 and P2the same across twinsand zygosity groups
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covarianceCross-twin covariance
![Page 30: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/30.jpg)
Phenotype 1 Phenotype 2 Phenotype 1 Phenotype 2
Phenotype 1Variance
P1
Phenotype 2Covariance
P1-P2Variance
P2
Phenotype 1Within-trait
P1Cross-trait Variance
P1
Phenotype 2Cross-trait Within-trait
P2Covariance
P1-P2Variance
P2
Predicted Model
Cross-twin covariancewithin each trait differsby zygosity
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covarianceCross-twin covariance
![Page 31: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/31.jpg)
Phenotype 1 Phenotype 2 Phenotype 1 Phenotype 2
Phenotype 1Variance
P1
Phenotype 2Covariance
P1-P2Variance
P2
Phenotype 1Within-trait
P1Cross-trait Variance
P1
Phenotype 2Cross-trait Within-trait
P2Covariance
P1-P2Variance
P2
Predicted Model
Cross-twin cross-traitcovariance differs byzygosity
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covarianceCross-twin covariance
![Page 32: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/32.jpg)
Example covariance matrix MZ
ADHD IQ ADHD IQ
ADHD 1
IQ -0.26 1
ADHD 0.64 -0.21 1
IQ -0.25 0.70 -0.31 1
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covarianceCross-twin covariance
![Page 33: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/33.jpg)
Example covariance matrix DZ
ADHD IQ ADHD IQ
ADHD 1
IQ -0.31 1
ADHD 0.20 -0.12 1
IQ -0.12 0.53 -0.27 1
Twin 1
Twin
1Tw
i n 2
Twin 2
Within-twin covariance
Within-twin covarianceCross-twin covariance
![Page 34: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/34.jpg)
Kuntsi et al. study
March 7, 2012 M. de Moor, Twin Workshop Boulder 34
![Page 35: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/35.jpg)
Summary
• Within-twin cross-trait covariance (phenotypic covariance) implies common aetiological influences
• Cross-twin cross-trait covariances >0 implies common aetiological influences are familial
• Whether familial influences are genetic or common environmental is shown by MZ:DZ ratio of cross-twin cross-trait covariances
![Page 36: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/36.jpg)
Specification in OpenMx?
March 7, 2012 M. de Moor, Twin Workshop Boulder 36
Twin 1Phenotype 1
A1
A2
E1
E2
a11 a21 a22
e11 e21 e22
1
1
1
1
Twin 1Phenotype 2
C1 C2
c11
c21
c22
11
Twin 2Phenotype 1
A1
A2
E1
E2
a11 a21 a22
e11 e21 e22
1
1
1
1
Twin 2Phenotype 2
C1 C2
c11
c21
c22
11
1/.5 1/.5
11 OpenMx script….
![Page 37: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/37.jpg)
Within-Twin Covariance (A)
P1 P2
A1 A2
a11 a21 a22
11 Path Tracing:
Lower 2 x 2 matrix:
P1
P2
a1 a2
a11 a21 a22
A a* aT
A = a%*%t(a)
![Page 38: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/38.jpg)
Within-Twin Covariance (A)
A a* aT
A = a%*%t(a)
Vars <- c("FSIQ","AttProb") nv <- length(Vars) aLabs <- c("a11“, "a21", "a22")
pathA <- mxMatrix(name = "a", type = "Lower", nrow = nv, ncol = nv, labels = aLabs)covA <- mxAlgebra(name = "A", expression = a %*% t(a))
OpenMx
![Page 39: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/39.jpg)
Within-Twin Covariance (A+C+E)
Using matrix addition, the total within-twin covariancefor the phenotypes is defined as:
A = a %*% t(a)
C = c %*% t(c)
E = e %*% t(e)
V
V
![Page 40: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/40.jpg)
OpenMx Matrices & Algebra
Vars <- c("FSIQ","AttProb") nv <- length(Vars)
aLabs <- c("a11“, "a21", "a22")cLabs <- c("c11", "c21", "c22")eLabs <- c("e11", "e21", "e22")
# Matrices a, c, and e to store a, c, and e Path CoefficientspathA <- mxMatrix(name = "a", type = "Lower", nrow = nv, ncol = nv, labels = aLabs)pathC <- mxMatrix(name = "c", type = "Lower", nrow = nv, ncol = nv, labels = cLabs)pathE <- mxMatrix(name = "e", type = "Lower", nrow = nv, ncol = nv, labels = eLabs)
# Matrices generated to hold A, C, and E computed Variance ComponentscovA <- mxAlgebra(name = "A", expression = a %*% t(a))covC <- mxAlgebra(name = "C", expression = c %*% t(c))covE <- mxAlgebra(name = "E", expression = e %*% t(e))
# Algebra to compute total variances and standard deviations (diagonal only)covPh <- mxAlgebra(name = "V", expression = A+C+E)matI <- mxMatrix(name= "I", type="Iden", nrow = nv, ncol = nv)invSD <- mxAlgebra(name ="iSD", expression = solve(sqrt(I*V)))
OpenMx
![Page 41: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/41.jpg)
MZ Cross-Twin Covariance (A)
P1 P2
A1 A2
a11 a21 a22
11
P1 P2
A1 A2
a11 a21 a22
11
1 1
Twin 1 Twin 2
Cross-twin within-trait:P1-P1 = 1*a11
2
P2-P2 = 1*a222+1*a21
2
Cross-twin cross-trait:P1-P2 = 1*a11a21
P2-P1 = 1*a21a11
)(%*% ataA
![Page 42: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/42.jpg)
DZ Cross-Twin Covariance (A)
P1 P2
A1 A2
a11 a21 a22
11
P1 P2
A1 A2
a11 a21 a22
11
0.5 0.5
Twin 1 Twin 2
Cross-twin within-trait:P1-P1 = 0.5a11
2
P2-P2 = 0.5a222+0.5a21
2
Cross-twin cross-trait:P1-P2 = 0.5a11a21
P2-P1 = 0.5a21a11
))(%*%%(%5.0%%5.0 ataxAx
![Page 43: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/43.jpg)
MZ/DZ Cross-Twin Covariance (C)Twin 1 Twin 2
P1 P2
C1 C2
c11 c21 c22
11
P1 P2
C1 C2
c11 c21 c22
11
1 1
)(%*% ctcC Cross-twin within-trait:P1-P1 = 1*c11
2
P2-P2 = 1*c222+1*c21
2
Cross-twin cross-trait:P1-P2 = 1*c11c21
P2-P1 = 1*c21c11
![Page 44: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/44.jpg)
Covariance Model for Twin Pairs
# Algebra for expected variance/covariance matrix in MZexpCovMZ <- mxAlgebra(name = "expCovMZ",
expression = rbind (cbind(A+C+E, A+C), cbind(A+C, A+C+E) ) )
# Algebra for expected variance/covariance matrix in DZexpCovDZ <- mxAlgebra(name = "expCovDZ",
expression = rbind (cbind(A+C+E, 0.5%x%A+C), cbind(0.5%x%A+C, A+C+E) ) )
OpenMx
![Page 45: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/45.jpg)
Unstandardized vs standardized solution
Twin 1Phenotype 1
A1
E1
a11 a21
e11 e21
1
1
C1
c11
c21
A2
E2
a22
e22
1
1
Twin 1Phenotype 2
C2
c22
11
Twin 2Phenotype 1
A1
A2
E1
E2
a11 a21 a22
e11 e21 e22
1
1
1
1
Twin 2Phenotype 2
C1 C2
c11
c21
c22
11
1/.5
1
1/.5
1
Twin 1Phenotype 1
A1
E1
a11
e11
1
1
C1
c11
A2
E2
a22
e22
1
1
Twin 1Phenotype 2
C2
c22
11
Twin 2Phenotype 1
A1
A2
E1
E2
a11 a22
e11 e22
1
1
1
1
Twin 2Phenotype 2
C1 C2
c11 c22
11
1/.5
1
1/.5
1
![Page 46: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/46.jpg)
Genetic correlation
• It is calculated by dividing the genetic covarianceby the square root of the product of the geneticvariances of the two variables
![Page 47: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/47.jpg)
1/.5 1/.5
A1 A2
a11
P11 P21
a22a21
A1 A2
a11
P12 P22
a22a21
Twin 1 Twin 2
rg a21a11
a112 *(a21
2 a222 )
Genetic correlation
Genetic covariance
Sqrt of the product of the genetic variances
![Page 48: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/48.jpg)
1/.5 1/.5
A1 A2
a11
P11 P21
a22
A1 A2
a11
P12 P22
a22
Twin 1 Twin 2
gr
Standardized Solution = Correlated Factors Solution
gr
![Page 49: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/49.jpg)
Genetic correlation – matrix algebra
corA <- mxAlgebra(name ="rA", expression = solve(sqrt(I*A))%*%A%*%solve(sqrt(I*A)))
OpenMx
![Page 50: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/50.jpg)
Contribution to phenotypic correlation
A1 A2
a11
P11 P21
a22
Twin 1
rg
If the rg = 1, the two sets of genes overlap completely
If however a11 and a22 are near to zero, genes do not contribute much to the phenotypic correlation
The contribution to the phenotypic correlation is a function of both heritabilities and the rg
![Page 51: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/51.jpg)
Contribution to phenotypic correlation
A1 A2
aP12
P2
rg
1 1
aP22
Proportion of rP1,P2 due to additive genetic factors:
2,1
22
21 **
PP
PgP
r
araP1
A1 A2
a11
P2
rg
1 1
a22
P1
a21112111211121
1121
eeccaa
aa
![Page 52: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/52.jpg)
Contribution to phenotypic correlation
ACEcovMatrices <- c("A","C","E","V","A/V","C/V","E/V")ACEcovLabels <-("covComp_A","covComp_C","covComp_E","Var","stCovComp_A","stCovComp_C","stCovComp_E")formatOutputMatrices(CholACEFit,ACEcovMatrices,ACEcovLabels,Vars,4)
OpenMx
Proportion of the phenotypic correlation due to geneticeffects
Proportion of the phenotypic correlation due to unshared environmental effects
Proportion of the phenotypic correlation due to shared environmental effects
![Page 53: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/53.jpg)
Summary / Interpretation
• Genetic correlation (rg) = the correlation between two latent genetic factors– High genetic correlation = large overlap in genetic
effects on the two phenotypes
• Contribution of genes to phenotypic correlation = The proportion of the phenotypic correlation explained by the overlapping genetic factors– This is a function of the rg and the heritabilities of the
two traits
![Page 54: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/54.jpg)
March 7, 2012 M. de Moor, Twin Workshop Boulder 54
Outline
• 11.00-12.30– Lecture Bivariate Cholesky Decomposition– Practical Bivariate analysis of IQ and attention problems
• 12.30-13.30 LUNCH• 13.30-15.00
– Lecture Multivariate Cholesky Decomposition– Practical Tri- and Four-variate analysis of IQ, educational
attainment and attention problems
![Page 55: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/55.jpg)
Practical
• Replicate findings from Kuntsi et al.
• 126 MZ and 126 DZ twin pairs from Netherlands Twin Register
• Age 12
• FSIQ• Attention Problems (AP) [mother-report]
March 7, 2012 M. de Moor, Twin Workshop Boulder 55
![Page 56: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/56.jpg)
Practical – exercise 1
• Script CholeskyBivariate.R• Dataset Cholesky.dat
• Run script up to saturated model
March 7, 2012 M. de Moor, Twin Workshop Boulder 56
![Page 57: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/57.jpg)
Practical – exercise 1
March 7, 2012 M. de Moor, Twin Workshop Boulder 57
MZ FSIQ1 AP1 FSIQ2 AP2
FSIQ1 1
AP1 1
FSIQ2 1
AP2 1
DZ FSIQ1 AP1 FSIQ2 AP2
FSIQ1 1
AP1 1
FSIQ2 1
AP2 1
• Fill in the table with correlations:
![Page 58: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/58.jpg)
Practical – exercise 1 - Questions
• Are correlations similar to those reported by Kuntsi et al.?
• What is the phenotypic correlation between FSIQ and AP?
• What are the MZ and DZ cross-twin cross-trait correlations?
• What are your expectations for the common aetiological influences?– Are they familial?– If yes, are they genetic or shared environmental?
March 7, 2012 M. de Moor, Twin Workshop Boulder 58
![Page 59: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/59.jpg)
Practical – exercise 2
• Run Bivariate ACE model in the script• Look whether you understand the output. If not,
ask us!• Adapt the first submodel such that you drop all C• Compare fit of AE model with ACE model
Script: CholeskyBivariate.R
March 7, 2012 M. de Moor, Twin Workshop Boulder 59
![Page 60: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/60.jpg)
Practical – exercise 2
• Fill in the table with fit statistics:
• Question:– Is C significant?
Script: CholeskyBivariate.R
March 7, 2012 M. de Moor, Twin Workshop Boulder 60
-2LL df chi2 ∆df P-value
ACE model
- - -
AE model
![Page 61: Copy files](https://reader031.vdocument.in/reader031/viewer/2022012919/56815259550346895dc08d25/html5/thumbnails/61.jpg)
Practical – exercise 3
• Now try to fill in the estimates for all paths in the path model (grey boxes):
March 7, 2012
M. de Moor, Twin Workshop Boulder
61
FSIQ
E1 E2
A1 A2
AP
1 1
1 1