GWAS Analysis in R: A Working Guide
How to run a genome-wide association study in R end to end, from VCF to QC to association testing to polygenic scores, and what you can and cannot do with a single genome.
2 posts
How to run a genome-wide association study in R end to end, from VCF to QC to association testing to polygenic scores, and what you can and cannot do with a single genome.
A post hoc |log2FC| > 1 filter is the wrong tool. Test against a fold-change threshold directly with DESeq2's lfcThreshold or edgeR's glmTreat, shrink effect sizes before ranking, and pick the threshold from your own replicate structure.