MLINK is one of the options of LINKAGE (LCP). MLINK allows you to find two point lodscores for your markers. Before you run LCP, make sure you have the following files in your directory:
- A Recoded *.ped File
- The datain.dat File
These step-by-step instructions will guide you through the program. On the left side of the table are instructions and on the right is what should appear on your screen. Typed commands are in bold text.
1. At the user prompt, type lcp. Press (RETURN) user@humgen% lcp 2. A screen entitled LINKAGE CONTROL PROGRAM will appear. On this screen, you must specify the input files for the program. Use the arrows on your keyboard to move up and down the list. Use the (DELETE) key when editing information in LCP. DO NOT USE THE BACKSPACE KEY. If you accidentally use the wrong key, you will end up on another screen. Just hit backspace unitl you return to the right screen. Change the PEDIGREE file name to the file name of your recoded data. The default PARAMETER file name is usually correct. Make sure it matches the name of your datain.dat file (if you named it something other than datain.dat). In this example, r-pedigree1.ped is used.
After you are finished, press (CTRL) N
COMMAND file name [pedin] :
LOG file name [final.out] :
STREAM file name [stream.out] :
PEDIGREE file name [pedin.dat] :
PARAMETER file name [datain.dat] :
Secondary PEDIGREE file name [] :
Secondary PARAMETER file name [] :pedin
final.out
stream.out
r-pedigree1.ped
datain.dat
3. At the Pedigree Options screen select General Pedigrees. Press (CTRL) N
General pedigrees :
Three-generation pedigrees :
Experimental cross pedigrees :<- 4. At the General Pedigree Analysis Options screen, choose MLINK. Press (CTRL) N
LODSCORE :
ILINK :
LINKMAP :
MLINK :
<-5. At the MLINK Test Options screen, choose the appropriate option. In this example, a Lod score table will be generated. Press (CTRL) N
Specific evaluation :
Lod score table :
Multiple pairwise Lod table :
<-6. At the MLINK Sex Difference Options screen, the only option is No sex difference. Press (CTRL) N
No sex difference : <- 7. Next, you will arrive at the Command Screen. Here, you must specify loci pairs for lodscoring. Select Locus order and type in a pair of loci separated by a space. You can only enter one pair at a time by typing a pair, then pressing (CTRL) N. You will not advance to another screen. Instead, if the program has accepted your specifications, the cursor will jump to he first number in your locus pair. In this example, the disease locus (1) will be paired with markers 1 (2) and 2 (3). After you have entered the last locus pair, press (CTRL) Z and the user prompt should appear.
Locus order [] :
Recombination Fractions [.0] :
Recombination varied [1] :
Other recomb. [.01 .05 .1 .2 .3 .4] :
Locus order [] :
Recombination Fractions [.0] :
Recombination varied [1] :
Other recomb. [.01 .05 .1 .2 .3 .4] :1 2
.0
1
.01 .05 .1 .2 .3 .4
1 3
.0
1
.01 .05 .1 .2 .3 .4user@humgen%
The next step is to run the pedin program. The pedin program will use all the specifications you entered in the LCP program to calculate two-point lodscores for your data. At the user prompt, type the command pedin. If this does not work, type ./pedin. Press (RETURN). If the run is successful, the screen will fill up with a series of data sets that look something like this:---------------------------
THETAS 0.400
---------------------------
PEDIGREE | LN LIKE | LOG 10 LIKE
------------------------
2 -60.736561 -26.377497 LOD= -0.108518
TOTALS -60.736561 -26.377497
-2 LN(LIKE) = 1.1214e+02 LOD SCORE = -0.108518
Now you can view lodscores with the Linkage Report Program. This program interprets the output of the pedin program. It allows you to see the two-point lodscores in a simple format.
1. At the user prompt, type lrp. Press (RETURN) user@humgen% lrp 2. A screen entitled LINKAGE REPORT PROGRAM will appear. On this screen, you can create a title for the report. Use the arrows on your keyboard to move up and down the list. Use the (DELETE) key when editing information in LRP. DO NOT USE THE BACKSPACE KEY. If you accidentally use the wrong key, you will end up on another screen. Just hit backspace unitl you return to the right screen. Do not change the default setting for STREAM file name. Type in the title you wish to give the report (optional). After you are finished, press (CTRL) N
STREAM file name [stream.out] :
REPORT title [] :stream.out
Pedigree One3. At the Report Options screen select General pedigree reports. Press (CTRL) N
General pedigree reports :
Three-generation pedigree reports :
Experimental cross pedigree reports :
Full stream file reports :<- 4. At the General Pedigree Report Options screen, choose MLINK. Press (CTRL) N
Two-point lodscore report (LODSCORE) :
Multi-point order report (ILINK) :
Location score report (LINKMAP) :
Lod table report (MLINK) :
<-5. Select a format. In this example, the table format will be used. Press (CTRL) N
Table format :
Full format :<- 5. At the Lod Table Report (MLINK) Options screen, you can specify what will be shown in the data table. In this example, the default settings will be used. Press (CTRL) N
Display LOG 10 results [Yes] :
Display LOG e results [No] :
Include Pedigrees [No] :Yes
No
No6. At the Report Output Options screen, you may chose to just view the data or save it to a file. Since it is possible to return to this screen after viewing the data, select Output report to the screen. Press (CTRL) N
Output report to the screen :
Output report to a file :<- 7. The two-point lodscore report should now appear on your screen. To save this report to a file, press (CTRL) R. This will take you back to the Report Output Options screen (In step #6). Select Output report to a file and press (CTRL) N. Create a report file name (the default is report.txt) and leave all the other default settings unchanged. Press (CTRL) N. You will return to the Report Output Options screen again. Press (CTRL) Z to exit the program.
L O D T A B L E R E P O R T
File: stream.out Pedigree One Screen: 1 of 1
_____________________________________________________
Order
-----
1=2
1=30.0
-----
-2.99
-5.120.01
-----
-2.30
-3.830.05
-----
-1.24
-2.200.1
-----
-0.68
-1.400.2
-----
-0.19
-0.660.3
-----
-0.02
-0.300.4
-----
0.02
-0.11= = Test Interval
a = LOD Scores
b = LOG 10 LikelihoodsYou should now have a *.txt file in your directory. This file can be edited in emacs. To print the file, at the user prompt type the command lpr followed by the name of the *.txt file.