Skip to content
Snippets Groups Projects
Unverified Commit e75338b9 authored by Reuben Pengelly's avatar Reuben Pengelly Committed by GitHub
Browse files

Update README.md

parent 7e8118eb
No related branches found
No related tags found
No related merge requests found
...@@ -122,24 +122,18 @@ Running LDMAP ...@@ -122,24 +122,18 @@ Running LDMAP
------------- -------------
#### ####
LDMAP is run using the following syntax:
Interaction with the core LDMAP binary, `ldmapper1` is most user friendly through the shell wrapper script `ldmap`. Running `ldmap` brings up the following menu: `./ldmapper1 [input.tped] [intermediate.txt] [job] [out.ldmap] [out.log] [MAF] [HWE]`
L D M A P where
Construction of linkage disequilibrium maps from diplotype data...
VERSION 2.0 November 2014
`[intermediate.txt]` is an intermediate file that will be created;\
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> `[job]` is an existing file containing run parameters (included, does not generally need edits);\
Enter number to select function `[out.ldmap]` will be created and contain the final LD map;\
1 Construct LD map in segments from intermediate file `[out.log]` will be created and contain detailed run information;\
2 quit `[MAF]` is the desired minor allele frequency cutoff (suggest 0.05);\
`[HWE]` is the desired Hardy-Weinberg equilibrium deviation p-value cutoff (Suggest 0.001).
####
The wrapper script `ldmap` works by requesting parameters from the user, and writing a small executable job files (`.cp`), which then runs `ldmapper1` with these desired parameters when executed. If automation of processes is desired, the required commands may be entered into a `.cp` file directly. An example `.cp` file may look like:
`./ldmapper1 input.tped int job example.map example.log 0.05 0.001`
where `int` is a temporary file used in LD map generation, `job` is a file containing run parameters (eponymous file included with source code), `example.map` will contain the final LD map, and `example.log` will contain detailed run information, `0.05` is the desired AF cutoff, and `0.001` is the desired HWE cutoff.
Output files Output files
------------ ------------
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment