Comments (4)
Hi, you have to phase, MSMC does not work on unphased data.
The command line you input is correct in principle, but if the VCFs are not phased, the output will be hard to interpret.
The other point here is that you are using way too many samples. You should first try one phased diploid from each population, then two from each at most. Eight individuals is the absolute maximum that I have ever run, but it takes excessive memory and requires even better phasing. Please try 2 phased diploids first!
Stephan
On 16 Nov 2016, at 04:15, dzaccook [email protected] wrote:
Hi,
I have two populations and each include many individual, I do not know the relationship of each individual in corresponding population, I had called SNP using GATK, but I can not phase the vcf. Can I generate multihetsep files to fed to MSMC (just determine the population size) like this:
./generate_multihetsep.py
--mask=22_chr1_mask.bed.gz (population 1)
--mask=23_chr1_mask.bed.gz (population 1)
--mask=24_chr1_mask.bed.gz (population 1)
--mask=25_chr1_mask.bed.gz (population 1)
--mask=1_chr1_mask.bed.gz (population 2)
--mask=2_chr1_mask.bed.gz (population 2)
--mask=3_chr1_mask.bed.gz (population 2)
--mask=4_chr1_mask.bed.gz (population 2)
--mask=chr1.mask.bed.gz
22_chr1.vcf.gz
23_chr1.vcf.gz
24_chr1.vcf.gz
25_chr1.vcf.gz
1_chr1.vcf.gz
2_chr1.vcf.gz
3_chr1.vcf.gz
4_chr1.vcf.gz >multihetsep.txt,
or any idea?
Thank you!—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub #25, or mute the thread https://github.com/notifications/unsubscribe-auth/AAbQmqcZlXXGhT5XbF4iLyfrENYA79pZks5q-nVNgaJpZM4KzTQf.
from msmc.
Thank you very much!
from msmc.
Hi,
@stschiff , I have only one sample and have not reference haplotypes.
So, I do not use shapeit to phase the vcf.
How can deal with it?
Maybe, can I use the software beagle deal with it?
from msmc.
If you have just one sample you don't need to phase.
from msmc.
Related Issues (20)
- AssertError bootstrap HOT 2
- another Segmentation Fault HOT 1
- error running msmc HOT 2
- running msmc without final result HOT 1
- install error HOT 1
- MSMC running with the divergent individuals HOT 1
- MSMC error HOT 12
- decode time units HOT 3
- Effect of introgression and selfing for estimation of population split HOT 1
- Abnormal results of MSMC, and need help ~ HOT 2
- trying to replicate Fig2 of Schiffels & Durbin (2014) HOT 6
- How to install MSMC HOT 1
- powell exceeding maximum iterations HOT 6
- precompiled binaries not working on macOS Catalina HOT 1
- "The sequence "chr1" was not found" HOT 1
- MSMC2- install error HOT 2
- Error in installing software
- Unexpected 'g' when converting from type char[] to type ulong HOT 7
- AssertionError while using generate_multihetsep.py HOT 1
- using the phased-vcf files as input for MSMC2 HOT 1
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 msmc.