xnr94
New Member
Posts: 4
|
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.
|
|
ax
New Member
Posts: 27
|
Post by ax on Jun 22, 2020 6:04:42 GMT
Hi xnr94,
Is that because your input two phenotypes are no 100% overlapped, so the software only take the intersect of them?
Cheers, Angli
|
|
|
Post by prerce on Feb 17, 2022 0:18:42 GMT
buy ivermectin uk <a href="https://stromectolreef.com/#">https://stromectolreef.com/</a>
|
|
|
Post by prerce on Feb 17, 2022 0:20:45 GMT
buy ivermectin uk <a href="https://stromectolreef.com/#">https://stromectolreef.com/</a>
|
|