|
|
A large quad-variate data set was simulated assuming an underlying Cholesky model, with 5000 observations per trait and high polygenicity (20,000 SNPs per genetic factor). The data can be downloaded [here](https://gitlab.gwdg.de/beate.stpourcain/grmsem_external):
|
|
|
|
|
|
```{r eval = FALSE}
|
|
|
[1] "Defined genetic factor loadings lambda_A:"
|
|
|
> print(lambda_A)
|
|
|
[,1] [,2] [,3] [,4]
|
|
|
[1,] 0.5477226 0.0000000 0.0000000 0.0000000
|
|
|
[2,] 0.4472136 0.6324555 0.0000000 0.0000000
|
|
|
[3,] 0.4472136 0.4472136 0.4472136 0.0000000
|
|
|
[4,] 0.3162278 0.4472136 0.0000000 0.6324555
|
|
|
|
|
|
[1] "Defined genetic variance A:"
|
|
|
> print(h2)
|
|
|
[,1] [,2] [,3] [,4]
|
|
|
[1,] 0.3000000 0.2449490 0.2449490 0.1732051
|
|
|
[2,] 0.2449490 0.6000000 0.4828427 0.4242641
|
|
|
[3,] 0.2449490 0.4828427 0.6000000 0.3414214
|
|
|
[4,] 0.1732051 0.4242641 0.3414214 0.7000000
|
|
|
|
|
|
[1] "Defined residual factor loadings lambda_E:"
|
|
|
> print(lambda_E)
|
|
|
[,1] [,2] [,3] [,4]
|
|
|
[1,] 0.8366600 0.0000000 0.0000000 0.0000000
|
|
|
[2,] 0.3162278 0.5477226 0.0000000 0.0000000
|
|
|
[3,] 0.4472136 0.0000000 0.4472136 0.0000000
|
|
|
[4,] 0.0000000 0.0000000 0.0000000 0.5477226
|
|
|
|
|
|
[1] "Defined residual variance E:"
|
|
|
> print(e2)
|
|
|
[,1] [,2] [,3] [,4]
|
|
|
[1,] 0.7000000 0.2645751 0.3741657 0.0
|
|
|
[2,] 0.2645751 0.4000000 0.1414214 0.0
|
|
|
[3,] 0.3741657 0.1414214 0.4000000 0.0
|
|
|
[4,] 0.0000000 0.0000000 0.0000000 0.3
|
|
|
|
|
|
[1] "Defined phenotypic variance V:"
|
|
|
> print(sigma_y)
|
|
|
[,1] [,2] [,3] [,4]
|
|
|
[1,] 1.0000000 0.5095241 0.6191147 0.1732051
|
|
|
[2,] 0.5095241 1.0000000 0.6242641 0.4242641
|
|
|
[3,] 0.6191147 0.6242641 1.0000000 0.3414214
|
|
|
[4,] 0.1732051 0.4242641 0.3414214 1.0000000
|
|
|
|
|
|
[1] "Simulated mean mu:"
|
|
|
[1] -8.120861e-18 -2.446934e-17 7.609713e-18 -1.669021e-17
|
|
|
|
|
|
[1] "Simulated phenotypic variance V:"
|
|
|
> print(var(ph))
|
|
|
Y1 Y2 Y3 Y4
|
|
|
Y1 1.0000000 0.5061277 0.6278417 0.1792348
|
|
|
Y2 0.5061277 1.0000000 0.6167281 0.4352479
|
|
|
Y3 0.6278417 0.6167281 1.0000000 0.3331679
|
|
|
Y4 0.1792348 0.4352479 0.3331679 1.0000000
|
|
|
|
|
|
``` |
|
|
\ No newline at end of file |