concordance rate between SNP observed and SNP not observed
0
0
Entering edit mode
4.4 years ago

Hi, i am just beginning with imputation and i try to understand the concordance rate between SNP observed and SNP not observed. can you tell me if the SNP observed are referred to reference panel ?? and what's the difference between SNP observed and not observed phasing and SNP observed and not observed not phasing ?? please if you can help me!

SNP • 1.8k views
0
Entering edit mode

Please provide some context... such as:

• which program you're using?
• what data have you got?
0
Entering edit mode

Hi, i will use Beagle sofware to impute missing data contained in a genotyping chip concerning a Malagasy population (sample) using a reference panel of 1000 genome. the imputation will concern chromosome 20. and i will repeat this operation with an european, asian and mexican population.The purpose of this study is to better understand their history, biological adaptation processes, and the prevalence of certain genetic risk factors.

i guess when i do the imputation, I need to measure the concordance rate between my imputed SNPs and the SNPs observed in the reference panel! But i don't know how to explain this part! can you bring me some HELP please?!

i want to install Beagle beagle V 4.1 on Ubuntu. If you have some advices about installation and use of Beagle I'm all ears!!

0
Entering edit mode

Okay, let's just get the program installed, first. Can you try to follow this? - https://francoismichonneau.net/2014/05/how-to-install-beagle-on-ubuntu/

0
Entering edit mode

just to know before trying to install Beagle, first I have to install java on Ubuntu is not it? because Beagle is implemented in java is not it ?? i will try to follow this and, I come back to you thank you for helping me Sir! I am a real beginner.

0
Entering edit mode

Yes, but that is done via the following command:

sudo apt-get install build-essential autoconf automake \
libtool subversion pkg-config openjdk-7-jdk


openjdk-7-jdk is a collection of JAVA libraries for run-time and development.

0
Entering edit mode

Ah Okay! thank you, so i will try it. And I will write to you tomorrow if you do not mind !

0
Entering edit mode

Sure thing!

0
Entering edit mode

Hi! the version of my Ubuntu is 16.04 and the site you advised me to follow is for Ubuntu 14.04, will that be fine ?? thank you.

0
Entering edit mode

I believe it should be okay. I also use Ubuntu 16.04

0
Entering edit mode

I type the firt line in my terminal

sudo apt-get install build-essential autoconf automake \
libtool subversion pkg-config openjdk-7-jdk


and the result is this:

 E: Unable to find package automake libtool
E: Can not find package pkg-donfig
E: The package "openjdk-7-jdk" has no version that can be installed


is there a solution? thank you

0
Entering edit mode

maybe i should install all this package separately

0
Entering edit mode

Doing that, the most updated version of java jdk will be installed

0
Entering edit mode

i followed this website: https://github.com/beagle-dev/beagle-lib/wiki/LinuxInstallInstructions and and he installed me libBeagle version 3.0 i want to install the version 4.1 of beagle beacause is the most accurate version of phasing and imputation!

0
Entering edit mode

0
Entering edit mode

i try to download this jar file "beagle.27Jan18.7e1.jar" ? when I click on it I have something incomprehensible. but i found some thing like that on google,is that it? Sorry for my questions!

#!/bin/bash
if [ ! -f beagle.27Jan18.7e1.jar ]; then
echo
wget http://faculty.washington.edu/browning/beagle/beagle.27Jan18.7e1.jar
fi
if [ ! -f bref.27Jan18.7e1.jar ]; then
echo
wget http://faculty.washington.edu/browning/beagle/bref.27Jan18.7e1.jar
fi

0
Entering edit mode

I'm assuming you're not conversant with java or jar files. Please Google how to execute a jar file and follow any of the Stack Overflow links.

0
Entering edit mode

You then execute it with java -jar beagle.27Jan18.7e1.jar

java -jar beagle.27Jan18.7e1.jar
beagle.27Jan18.7e1.jar (version 4.1)
Copyright (C) 2014-2015 Brian L. Browning
Command line syntax: java -jar beagle.27Jan18.7e1.jar [arguments]

data input/output parameters ...
gt=<VCF file: use GT field>                        (optional)
gl=<VCF file: use GL/PL field>                     (optional)
gtgl=<VCF file: use GT (preferred) or GL/PL field> (optional)
ref=<VCF file with phased genotypes>               (optional)
out=<output file prefix>                           (required)
excludesamples=<file with 1 sample ID per line>    (optional)
excludemarkers=<file with 1 marker ID per line>    (optional)
map=<PLINK map file with cM units>                 (optional)
chrom=<[chrom] or [chrom]:[start]-[end]>           (optional)
maxlr=<max GL/PL likelihood ratio>                 (default=5000)

general parameters ...
lowmem=<use low-memory algorithm (true/false)>     (default=false)
window=<markers per window>                        (default=50000)
overlap=<overlap between windows>                  (default=3000)
seed=<random seed>                                 (default=-99999)

phasing and imputation parameters ...
niterations=<number of phasing iterations>         (default=5)
impute=<impute ungenotyped markers (true/false)>   (default=true)
gprobs=<print GP field for imputed markers>        (default=false)
ne=<effective population size>                     (default=1000000)
err=<allele miscall rate>                          (default=0.0001)
cluster=<max cM in a marker cluster>               (default=0.005)

IBD parameters ...
ibd=<perform IBD detection (true/false)>           (default=false)
ibdlod=<min LOD score of reported IBD segments>    (default=3.0)
ibdcm=<min cM length of reported IBD segments>     (default=1.5)
ibdscale=<model scale factor for Refined IBD>      (default: data-dependent)
ibdtrim=<markers at each segment end>              (default=40)

0
Entering edit mode

sorry for the way of writing my messages!! I'm working on Ubuntu when i click on th link "beagle.27Jan18.7e1.jar" found in "https://faculty.washington.edu/browning/beagle/beagle.html#introduction" i have this :

---- binary junk deleted --------

0
Entering edit mode

What about right-click, and then 'Save As' or 'Save Target As?

Your browser may be trying to open the JAR file as plain text

0
Entering edit mode

thank you i did right -click and then "save target as" on my desk and then i did on my terminal "Java -jar beagle.27Jan18.7e1.jar" but it tell me "unable to access jarfile beagle.27Jan18.7e1.jar" !!

0
Entering edit mode

Because the directory may not be in your \$PATH. Make sure you cd to Desktop if you save the file there.

And then try: java -jar ./beagle.27Jan18.7e1.jar

0
Entering edit mode

Hi , I had problems with my ubuntu it is blocked but finally I was able to install beagle in another computer. that's why I was absent. thanks for your previous help.

Now, how to impute on Beagle 4.1 my missing data using a referecne panel. is there a script or some commands to do this operation ?? I can not find a very precise document on that. Can you help me please ?!

0
Entering edit mode

Hello, I recommend that you open a new thread specifically about running Beagle. This thread's issue has already been solved, and it would be good to bring in other users to the discussion.

0
Entering edit mode

ok, thank you i did it

0
Entering edit mode

Thank you. I will check again tomorrow to see if there are any responses to the new thread.

0
Entering edit mode
#!/bin/bash


0
Entering edit mode

Please use the formatting bar (especially the code option) to present your post better. I've done it for you this time.

0
Entering edit mode

0
Entering edit mode
 #!/bin/bash
if [ ! -f beagle.27Jan18.7e1.jar ]; then
echo
wget http://faculty.washington.edu/browning/beagle/beagle.27Jan18.7e1.jar
fi
if [ ! -f bref.27Jan18.7e1.jar ]; then
echo
wget http://faculty.washington.edu/browning/beagle/bref.27Jan18.7e1.jar
fi