Post by xnr94 on Jun 5, 2020 17:12:35 GMT
Hello All.
When i'm running GCTA-GREML analysis with one variable (weight) I've got this results:
system("./gcta64 --grm SS10 --pheno SS10.weight --reml --autosome --autosome-num 29 --out REML --thread-num 10")
-V(G) 115.85 +- 83
-V(e) 4661 +- 183
-Vp 4777.64 +- 176.56
V(G)/Vp = 0.024 +- 0.017
But when I run the bivariate I got the next results (ADDING survival as 0 and 1) I get the next:
system("./gcta64 --reml-bivar 1 2 --grm SS10 --pheno pheno.txt --out BiVar")
-V(G)_tr1 144.04 +- 90
-V(e) 4661 +- 183
-Vp_tr1 4869.86 +- 182.47
V(G)/Vp_tr1 = 0.0249 +- 0.018
Does this models consider each variables in a joint way? I mean, if this genetic variance estimate is only for trait 1, then shouldn't they be equal ? Why does that change? And why would phenotypic variance change too in the same dataset?
Thanks for helping me out.
When i'm running GCTA-GREML analysis with one variable (weight) I've got this results:
system("./gcta64 --grm SS10 --pheno SS10.weight --reml --autosome --autosome-num 29 --out REML --thread-num 10")
-V(G) 115.85 +- 83
-V(e) 4661 +- 183
-Vp 4777.64 +- 176.56
V(G)/Vp = 0.024 +- 0.017
But when I run the bivariate I got the next results (ADDING survival as 0 and 1) I get the next:
system("./gcta64 --reml-bivar 1 2 --grm SS10 --pheno pheno.txt --out BiVar")
-V(G)_tr1 144.04 +- 90
-V(e) 4661 +- 183
-Vp_tr1 4869.86 +- 182.47
V(G)/Vp_tr1 = 0.0249 +- 0.018
Does this models consider each variables in a joint way? I mean, if this genetic variance estimate is only for trait 1, then shouldn't they be equal ? Why does that change? And why would phenotypic variance change too in the same dataset?
Thanks for helping me out.