Entering edit mode
4.7 years ago
ssadegi42
▴
20
I want use Beagle to impute my data; How I can imputed low density (5k) to 50K by beagle software?
I want use Beagle to impute my data; How I can imputed low density (5k) to 50K by beagle software?
thanks a lot @ dr. Fatima and Dr.Kevin I use GATK software to convert un-phased.vcf to phased.vcf for ref input
also I could find a simple way to convert unphased to phased using Beagle
thanks
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Please confirm that you have first tried to do this on your own (?). There is no way to provide a quick / simple answer for you. Have you looked at the Beagle manual, for example, and other online tutorials?
I create vcf file by plink for reference (50k) file and test (5k) file by bfile; ./plink --file ref --cow --recode vcf --out ref ./plink --file test --cow --recode vcf --out test and i want run a code to impute missing genotype with MY reference panel (50k) java -Xmx20g -jar beagle5.jar gt=test.vcf ref=ref.vcf out=out_myfile
but ref.vcf file not found???????????????
I haven't used any of these tools, but if ref.vcf is not found then your first command is not working.
Which version of plink are you using?
https://www.cog-genomics.org/plink2/input
https://www.cog-genomics.org/plink/2.0/output
I think you can also use --recode-beagle instead of --recode to convert your file to something recognizable by beagle:
hello dear Fatima
I use beagle 5 version and this version uses --recode vcf to run the input file of beagle software, recode-beagle use to beagle 4 version. i run .map and .ped files by plink and achieve the output file of VCF for ref and test data. vcf file format is:
Are beagle and ref.vcf in the same folder? If not, are you using the correct path to the vcf file?
https://faculty.washington.edu/browning/beagle/beagle_5.1_08Nov19.pdf https://faculty.washington.edu/browning/beagle/run.beagle.21Mar20.d65.example
I run beagle software but I encounter this error. please help me to solve
Command line:
I think that the problem is that your reference file must have phased genotypes.
Yes, as @kevin mentioned the error is related to phase.
(I deleted my previous comment about path to the file since it sounded irrelevant now that you're getting another error rather than ref.vcf not found. )
You can also google these error to find how to solve them.
thanks a lot, dear Kevin for your guide; when I use the Plink to create VCF file for "ref" data, the phased allele separator was ‘/’
since I'm a beginner
I use R software and replace ‘|’ instead of ‘/’ in Ref.VCF and Beagle was executed. But, I don't know that this way is true or false?
also this comment "No genetic map is specified: using 1 cM = 1 Mb" still observed; and I don't know whether this comment is a problem or not?
also I convert the imputed Beagle output (by this way) by plink to create Ped and Map file, but the most of imputed genotypes are hemozygote.
please help me to create ref.VCF using ‘|’ separator.
That is bad practice, and immediately puts error in your data. Can you please take a look at the difference between phased and un-phased genotypes? For Beagle, you will require a phased reference dataset.
The 1000 Genomes Phase III data is phased, by the way. You can download it here: Produce PCA bi-plot for 1000 Genomes Phase III - Version 2
Dr. Kevin thanks for your help and i hope to know you meaning. I studied phased and un-phased for target and REF dataset in beagle manual. I want use of my Ref file (50k) not 1000genome project as the reference dataset. so i want convert the un-phased my Ref.Vcf to phased Ref.vcf to impute target file. can you help me and introduce a software to convert un-phased to phased file? I hope that I have been able to convey the concept
I don't know if using tools like these ones can help or puts error in your data:
https://github.com/broadinstitute/gatk-docs/blob/master/gatk3-methods-and-algorithms/Purpose_and_operation_of_Read-backed_Phasing.md
https://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html
http://manpages.ubuntu.com/manpages/bionic/man1/vcftools.1.html
Dear Fatima;
Thank you for taking the time to help me, I will try to use these softwares.