How to handle with negative variance components when genomic analysis is interrupted (reacta, gcta)
Hi all!
One question regarding regional heritability mapping analysis of 50K sheep genotypes using reacta and gcta software.
The sample size is very low and as a result I get errors such as: "the matrix V becomes negative-definite because of one of the variance component is negative" .
The suggested solution is to run the analysis without the --reml-no-constrain option where I get the error: "analysis stopped because more than half of the variance components are constrained. The result would be unreliable"
Any suggestions to perform the analysis successfully even with low number of individuals?
Many thanks in advance! Valentina
• 1,469 views
•
link
written
by
valentinatsar •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Blog: The Complete Guide to Metatranscriptome Sequencing
written by Novogene •
-
Tutorial: Single-cell RNA-seq: Preprocessing: Preprocessing the data of scRNA-seq with omicverse-1
written by Julia Ma •
-
Problems using ERCC spike-ins for normalization in DESeq2
written by manuel.fernandez •
-
Tutorial: GEO data mining (III) - Gene pathway enrichment analysis with complete code sharing
written by Novogene •
-
Forum: Sharing frequently asked questions about microbial genome sequencing II
written by Novogene •
-
Extremely low pve in a GWAS study?
written by maya123z •
-
EBSeq explained to a biologist
written by Rik Verdonck •
-
Adjusting for population substructure in heterozygosity using regression with PCs
written by Vivek
-
Tutorial: InsideDNA: Benchmarking seven most popular genome assemblers
written by shilparaopradeep •
-
Forum: Help Me Gain A Little More Insight On What I'M Seeing Here Please
written by proxify •
