eryn
New Member
Posts: 1
|
Post by eryn on Dec 21, 2016 15:30:31 GMT
Hello,
I'm sorry if this has been asked somewhere, and if I just haven't found the answer.
I have been partitioning variance where I made a separate grm for each chromosome, and then asked how much of the variance each explained.
./gcta64 --mgrm multi_grm.txt --pheno phen --reml --reml-maxit 500 --covar covar_sex.txt --out test
Where, multi_grm.txt is chr1_grm chr2_grm...
Depending on the trait this does seem to work, but sometimes it only works if I drop the smallest chromosomes from the analysis. Then I have an answer for some chromosomes, but not all of them.
What I would really like to do is the equivalent (I think) of a mlma-loco, where I use all of the chromosomes except one to get a grm while I estimate how much of the variance is explained by the focal chromosome. I'm not as keen on the GWAS that pops out of the mlma-loco, I just want a chromosome specific heritability estimate.
Thanks so much for any help!
Eryn
|
|
|
Post by Jian Yang on Dec 22, 2016 4:38:43 GMT
If your purpose is to estimate h2SNP for a particular chromosome, you can fit a model y = c + g_minus_c + e where c is the genetic effect for the focal chromosome and g_minus_c is the genetic effect for all chromosomes except the focal chromosome. For example, the focal chr is chr1, then the multi_grm.txt should be
chr1_grm all_but_chr1_grm
all_but_chr1_grm can be computed by merging the GRMs of the other 21 chromosomes.
|
|
mswat
New Member
Posts: 1
|
Post by mswat on Oct 8, 2018 14:29:18 GMT
I also interested in doing this LOCO approach for partitioning variance across chromosomes. For each chromosome, however, the variance explained is very small (e.g. 0.001-0.01) and does not come close to adding up to the total heritability of the trait. I realize this could be due to covariance among chromosomes so I have filtered for LD 0.8 before calculating the "all_but_chr*_grm"'s. Are there other explanations for why the variance explained by each chromosome is so low? Would relatedness among individuals cause this?
|
|