This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Rehh package input files

Dear respected members, How to make "map.inp" files for running rehh? I have .ped and .map files from plink. I have got the haplotype information from fastPHASE (hapguess_switch.out).

Thank you in advance.

r rehh map.inp
• 2,001 views • link
written by kas • kas

Did you read the manual?

It is a simple text file, example from the manual:

# show first 6 lines of file
cat(readLines("map.inp", n = 6), sep = "\n")
> F0100190 1 113642 T A
> F0100220 1 244699 C G
> F0100250 1 369419 G C
> F0100270 1 447278 A T
> F0100280 1 487654 T A
> F0100290 1 524507 C G
• 0 views • link
written by zx8754 zx8754

Thank you for your response. I've read the manual. Here totally 5 columns are there. Last 2 for ancestral and derived alleles. But I don't have that information. How to get those ancestral and derived alleles?

If I load my map file in R, It's showing like this 

Error in read.table(map_file, row.names = 1, colClasses = "character") :
  duplicate 'row.names' are not allowed.

This is my map file format.

rs41609588  chr1    135098
rs108982244 chr1    267940
rs43703977  chr1    393248
rs110887329 chr1    471078
rs29015852  chr1    516404
rs110936909 chr1    571340
rs108981857 chr1    845494
• 0 views • link
updated by zx8754 zx8754 • written by kas • kas

0 answers

No answers yet.

Log in to answer this question.

More posts like this