|
|
GCTA `GREML` estimates were fitted as follows:
|
|
|
|
|
|
```{r eval = FALSE}
|
|
|
gcta64 --reml-bivar 1 2 --grm-gz large.gcta --pheno large.gcta.phe
|
|
|
--reml-bivar-lrt-rg 0 --thread-num 4 --out Y1_Y2 --reml-maxit 1000 >
|
|
|
large.gcta.Y1Y2.log`
|
|
|
gcta64 --reml-bivar 1 2 --grm-gz large.gcta --pheno large.gcta.phe --reml-bivar-lrt-rg 0 --thread-num 4 --out Y1_Y2 --reml-maxit 1000 > large.gcta.Y1Y2.log`
|
|
|
```
|
|
|
|
|
|
which yielded the estimates:
|
... | ... | |