Open xiangzhu opened 6 years ago
I wonder if we could use a warnings() (or even a stop()) when people try to feed summary data from MHC region to the program.
warnings()
stop()
The main reason is that LD pattern of MHC is complicated, and may be hard to estimated reliably from 1000 Genomes data.
Related notes: https://crerecombinase.github.io/PolygenicRSS/New_WholeGenome.html#mhc
I wonder if we could use a
warnings()(or even astop()) when people try to feed summary data from MHC region to the program.The main reason is that LD pattern of MHC is complicated, and may be hard to estimated reliably from 1000 Genomes data.
Related notes: https://crerecombinase.github.io/PolygenicRSS/New_WholeGenome.html#mhc