Comments (4)
Hi Wang Min,
In practice, we often find that the additive genetic variance is much more important than dominant/recessive effects. This is typically the case for complex phenotypes in human populations, and is often the case in other organisms/populations. Incorporating dominance effects into the LMM is challenging (at least if we want the computation to be fast), and is generally not used much in human genetics, therefore we have not put much effort into this direction, but we will keep your request open.
For an LMM method that does estimate for dominance effects, I recommend QTLRel, although note that this package is about 10x slower than GEMMA.
Thank you for your interest in GEMMA,
Peter
from gemma.
Closed due to inactivity. Please reopen if someone wants to work on it.
from gemma.
I'm faced with a similar challenge. Would it be possible to calculate kinship as normally, and then recode as @kumine describes above? Or does this break something downstream? Probably a naive question!
from gemma.
Hmm, Peter is of course correct, but once you have a locus or SNP of interest, it should be trivial to compute the dominance deviation on a per-SNP basis WITHOUT any correction for kinship. This estimate will often be "good enough" for most users. One problem with dominance deviations is that they depend greatly on scale. Log2 expression data will generally show a dominance deviation that is simply due to scaling.
from gemma.
Related Issues (20)
- Parsing input file '' failed in function AnalyzeBimbamGXE
- Enforce failed for number of analyzed individuals equals 0. in src/param.cpp at line 2073 in ProcessCvtPhen HOT 2
- long haplotype as population structure
- basic usage doesn't work
- Gemma run-time for Estimate Relatedness Matrix calculation
- very low number of analyzed SNPs/var
- The input files includes 26 linkage groups but the output files only includes 22 linkage group
- Categorical covariates
- Pve is .999
- gemma-wrapper:279:in `<main>': Expected GEMMA -a genotype file switch (RuntimeError)
- Got a nan in the p-value results HOT 1
- gsl: newton.c:88: ERROR: function value is not finite 3.08%
- GSL ERROR: function value is not finite in brent.c at line 58 errno 9
- WARNING: Brent did not converge
- Why the diagonal of relatedness matrix is not equal to 1 or a constant value? HOT 1
- The problem of filtering snp parameters
- ERROR: Enforce failed for Problem reading FAM file (phenotypes do not match geno file) in src/gemma_io.cpp at line 600 in ReadFile_fam GEMMA 0.98.3 (2020-11-28) by Xiang Zhou and team (C) 2012-2020
- Identical P-Values at end of run HOT 1
- BIMBAM format (preferred)
- ##ERROR: FAILED: Parsing input file 'mikk_small_hapVariants_updated_noheaders.csv' failed in function ReadFile_geno in src/gemma_io.cpp at line 744, for GEMMA 0.98.5 (2021-08-25) by Xiang Zhou, Pjotr Prins and team (C) 2012-2021
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from gemma.